# Exploring the Mechanisms of Recovery in Anorexia Nervosa through a Translational Approach: From Original Ecological Measurements in Human to Brain Tissue Analyses in Mice.

Philibert Duriez, Ida a K Nilsson, Ophelia Le Thuc, David Alexandre, Nicolas Chartrel, Carole Rovere, Christophe Chauveau, Philip Gorwood, Virginie Tolle, Odile Viltart

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

## To cite this version:

Philibert Duriez, Ida a K Nilsson, Ophelia Le Thuc, David Alexandre, Nicolas Chartrel, et al. Exploring the Mechanisms of Recovery in Anorexia Nervosa through a Translational Approach: From Original Ecological Measurements in Human to Brain Tissue Analyses in Mice. Nutrients, 2021, 13 (8), pp.2786. 10.3390/nu13082786. [hal-03667160](https://hal-univ-lille.fr/hal-03667160)

HAL Id: hal-03667160

https://hal.univ-lille.fr/hal-03667160

Submitted on 13 May 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Exploring the Mechanisms of Recovery in Anorexia Nervosa through a Translational Approach: From Original Ecological Measurements in Human to Brain Tissue Analyses in Mice

Philibert Duriez <fnref n="1" />,<fnref n="2" /> ID, Ida A. K. Nilsson <fnref n="3" /> ID, Ophelia Le Thuc <fnref n="4" /> ID, David Alexandre <fnref n="5" /> ID, Nicolas Chartrel <fnref n="5" />, Carole Rovere <fnref n="4" />, Christophe Chauveau <fnref n="6" />,<fnref n="7" />, Philip Gorwood <fnref n="1" />,<fnref n="2" /> ID, Virginie Tolle <fnref n="1" /> and Odile Viltart <fnref n="1" />,<fnref n="7" />,* ID

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

check for
updates

**Citation:** Duriez, P.; Nilsson, I.A.K.; Le Thuc, O.; Alexandre, D.; Chartrel, N.; Rovere, C.; Chauveau, C.; Gorwood, P.; Tolle, V.; Viltart, O. Exploring the Mechanisms of Recovery in Anorexia Nervosa through a Translational Approach: From Original Ecological Measurements in Human to Brain Tissue Analyses in Mice. *Nutrients* **2021**, *13*, 2786. https://doi.org/10.3390/nu13082786

Academic Editors:
Odysseas Androutsos,
George Mastorakos and
Evangelia Charmandari

Received: 13 July 2021
Accepted: 9 August 2021
Published: 13 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Anorexia nervosa (AN) is a severe eating disorder where caloric restriction, excessive physical activity and metabolic alterations lead to life-threatening situations. Despite weight restoration after treatment, a significant part of patients experience relapses. In this translational study, we combined clinical and preclinical approaches. We describe preliminary data about the effect of weight gain on the symptomatology of patients suffering from acute AN (n = 225) and partially recovered (n = 41). We measured more precisely physical activity with continuous cardiac monitoring in a sub-group (n = 68). Using a mouse model, we investigated whether a long-term food restriction followed by nutritional recovery associated or not with physical activity may differentially impact peripheral and central homeostatic regulation. We assessed the plasma concentration of acyl ghrelin, desacyl ghrelin and leptin and the mRNA expression of hypothalamic neuropeptides and their receptors. Our data show an effect of undernutrition history on the level of physical activity in AN. The preclinical model supports an important role of physical activity in the recovery process and points out the leptin system as one factor that can drive a reliable restoration of metabolic variables through the hypothalamic regulation of neuropeptides involved in feeding behavior.

**Keywords:** physical activity; heart rate; anorexia nervosa; body mass index; hypothalamus; neuropeptides; ghrelin; leptin; animal model

## 1. Introduction

Anorexia nervosa (AN) is a complex psychiatric disorder emerging most commonly in adolescence and young adulthood whose etiology has remained poorly understood up to now, even if metabolic and endocrine aspects are now considered as key pathophysiological determinants [1–5]. The severe weight loss due to food restriction is usually combined with excessive physical activity [6–8]. The identification of factors influencing the onset, duration and quality of the recovery process is a research priority for eating disorders [9]. One out of two patients with AN relapses within a year following inpatient treatment

<footnote n="1">¹ Institute of Psychiatry and Neuroscience of Paris (IPNP), University of Paris, INSERM UMR-S 1266, F-75014 Paris, France; p.duriez@ghu-paris.fr (P.D.); p.gorwood@ghu-paris.fr (P.G.); virginie.tolle@inserm.fr (V.T.)</footnote>

<footnote n="2">² GHU Paris Psychiatry and Neurosciences, Hospital Sainte-Anne, F-75014 Paris, France</footnote>

<footnote n="3">³ Department of Molecular Medicine & Surgery, Karolinska Institutet, Centre for Eating Disorders Innovation (CEDI), Medical University, Karolinska Institutet, S-17176 Stockholm, Sweden; ida@ki.se</footnote>

<footnote n="4">⁴ CNRS UMR 7275, Institute of Molecular and Cellular Pharmacology (IPMC), University of Nice-Sophia Antipolis, F-06560 Valbonne, France; ophelia.lethuc@helmholtz-muenchen.de (O.L.T.); rovere@ipmc.cnrs.fr (C.R.)</footnote>

<footnote n="5">5 INSERM U1239, Laboratory of Neuronal and Neuroendocrine Differentiation and Communication, University of Normandie, UNIROUEN, F-76821 Mont-Saint-Aignan, France; david.alexandre@univ-rouen.fr (D.A.); nicolas chartrel@univ-rouen.fr (N.C.)</footnote>

<footnote n="6">6 Marrow Adiposity and Bone Laboratory (MABLab), University of Littoral Côté d'Opale, CHRU Lille, F-62327 Boulogne sur Mer, France; christophe.chauveau@univ-littoral.fr</footnote>

<footnote n="7">7 Faculty of Sciences and Technologies, University of Lille, F-59650 Villeneuve d'Ascq, France</footnote>

<footnote>* Correspondence: odile.viltart@univ-lille.fr; Tel.: +33-6-76-88-05-06</footnote>
and approximately 20% of patients experience either recurrent patterns of remission and relapses or have a chronic outcome [10–12]. Weight gain results in an initial improvement in AN symptomatology and constitutes a major first step in the medical treatment of AN. Weight gain stabilization, within a range of global symptomatology improvement, is a therapeutic priority, which has to be maintained for months. Partial recovery is probably a frequent evolution after an acute relapse, with persistent qualitative food restriction, problematic physical exercise and/or body image alteration [13].

Previous studies aimed at determining biological markers of the remission processes. Leptin and ghrelin are key candidates well known for their role in regulating feeding behavior, energy metabolism and motivational and rewarding processes associated with feeding and running [4,5]. Leptin is an adipokine, mainly expressed in adipocytes [14], for which plasma concentrations increase in proportion to body fat stores [15]. This hormone is involved in weight regulation and various biological functions such as immune function, reproduction, bone formation, anxiety and physical activity [5,16–18]. Ghrelin is a pleiotropic hormone, mainly released by the X/A-like endocrine cells in the oxyntic glands of the gastric fundus and initially described as a growth hormone secretagogue [19,20]. Similar to leptin, ghrelin is known to participate both in metabolic regulations such as food initiation, glucose homeostasis or gastric emptying and in complex behaviors such as anxiety, the reward process or compulsive physical activity [21–25]. Two ghrelin isoforms have been described: acylated ghrelin (AG), which has an n-octanoylated serine in position three and desacyl ghrelin (DAG, [26,27]). Both forms have distinct biological functions [21,28].

At the clinical level, basal plasma leptin concentrations have recently been consid-
ered as a reliable marker of stable recovery with a better predictivity than body mass
index [29]. Furthermore, in a pilot study of three patients with restrictive AN, a two-week
treatment with metreleptin (recombinant human leptin) generated positive outcomes such
as reduction in physical activity, of weight phobia and of negative emotional feelings [30].
In patients with AN, plasma variations of AG and DAG during the weight restoration
period might be considered as reliable biomarkers of a successful recovery [31]. Indeed, we
recently showed that DAG plasma concentration after weight restoration is associated with
instable recovery and an increase in physical activity after discharge [32]. However, an
accurate measurement of physical activity in AN remains difficult and data are still lacking,
limiting translational benefits from preclinical models such as the activity-based rodent
model of AN [33].

Both ghrelin and leptin target brain structures involved in the homeostatic control of food intake. More specifically, the hormones receptors are expressed on neurones located in the hypothalamus, namely the orexigenic neurones co-expressing agouti-related peptide (AgRP) and neuropeptide Y (NPY) and the anorexigenic neurones expressing proopiomelanocortin (POMC). These neurones constitute critical nodes of a circuit that sense key metabolic cues to regulate energy metabolism and feeding behavior [34–38]. Recently, both AgRP and POMC neurones were shown to also regulate circuitries controlling non-feeding behavior, such as the reward pathway, anxiety, stereotypic behaviors, compulsive exercise and delay discounting all related to AN [34,39–45]. These neurones are precisely located in the ventromedial part of the hypothalamic arcuate nucleus (ARC) and are considered as primary targets of convergent peripheral hormones [46,47]. Recently, Florent et al. (2019), using magnetic resonance imaging, clearly demonstrated alterations in the hypothalamic response to a meal in patients with AN compared to control and subjects with constitutional thinness [48].

The current imaging techniques are not precise enough to capture the changes oc-
curring at subcellular levels. For this purpose, using preclinical models is a prerequisite
to better understand mechanisms involved in the remission or relapse of AN. Several
animal models are now available to assess the various aspects of AN [32,49–52]. Using the
"activity based-anorexia" (ABA) model, Verhagen et al. (2011) demonstrate that central
leptin injections in the lateral ventricle or local injections of leptin into the ventral tegmental
area suppress the running wheel activity [53]. In particular, using a mouse model that reca-
pitulates metabolic signatures of AN, we recently showed that plasma DAG concentrations were positively correlated with diurnal physical activity during a two-week nutritional recovery [32]. These data let us hypothesize that chronic food restriction, associated with a moderate physical activity, might differentially impact the expression of hypothalamic neuropeptides targeted by ghrelin and leptin.

In this translational study, we combined a clinical and preclinical approach to find potential explanations to remission or relapse in partially recovered patients with AN. We first describe initial original data about physical activity in both acutely ill and partially recovered patients with AN in order to highlight the candidate mechanisms of the recovery process. Secondly, using the FRW (food restricted with wheel) model, we investigate whether short (15 days) or long-term (65 days) food restriction associated or not with physical activity, may differentially impact the mRNA expression of hypothalamic neuropeptides and their receptors. As we will also measure the plasma concentration of AG, DAG and leptin, we will present the hypothalamic mRNA expression of the ghrelin receptor and leptin receptor.

## 2. Materials and Methods

### 2.1. *Experiment 1: Clinical Assessment*

All female patients with a current diagnosis of AN (n = 266) were enrolled in the present study after referring to a specialized center for eating disorders, the CMME (Clinique des Maladies Mentales et de l'Encéphale, Paris, France). All patients signed a written consent for their participation in research protocol, according to the Helsinki Declaration. The diagnosis of AN was established by a clinical interview conducted by a trained clinician using a validated instrument (MINI-S), following DSM-5 diagnostic criteria [54]. According to DSM-5, patients who previously met all criteria for AN but for a sustained period were not low in body weight were diagnosed as partial recovery AN (prAN) [54].

Clinical and biological data were obtained after a two-day multidisciplinary assessment of the patient. We gathered age, current BMI, adult minimal BMI and illness duration. Minimal BMI reflects the severity of weight loss and undernutrition state during the course of AN according to DSM-5 [54]. We also collected classical nutrition markers: albumin (normal range: 35–52 g/L) and ferritin (normal range: 15–150 µg/L) blood levels.

Eating disorder symptoms were assessed using the Eating Disorder Inventory-2 [55]. The EDI2 consists of 91 items measuring 11 psychopathological dimensions. High scores represent greater levels of psychopathology. Depression and anxiety scores were measured with the self-report Hospital and Anxiety Depressive Scale (HADS) [56]. A high score (max 21) indicates a high level of anxiety or depression, respectively. Compulsivity and obsessionality were monitored with the Yale Brown Scale [57]. A high score (max 20) shows a high level of obsessionality and/or compulsivity.

We recently developed the evaluation of cardiac activity through a 72 h monitoring, using electrocardiogram Holter monitoring (ePatch®). On the 68 patients evaluated, 55 were patients with acute AN (aAN) and 13 were patients with prAN. The Holter monitor was fitted on day 1 and removed three days later. Maximum heart rate (HRmax) was calculated for each patient with the formula: 220-age. We calculated proportion of time when the heart rate was in four ranges defined as: bradycardia (below 60 beats/min), physiological (between 60 beat/min and 50% of HRmax), moderate activity (between 50% and 60% of HRmax), high activity (higher than 60% of HRmax).

### 2.2. *Experiment 2: Preclinical Approach*

#### 2.2.1. Animals

Adult C57BL/6J female mice (7 wk old, average initial body weight 18.8 ± 0.15 g; Charles River Laboratories, L'Arbresle, France) were housed two per cage in all the procedures. Mice are social animals; thus, in an ethological perspective and for the welfare of the animals, we decided to maintain our mice two per cage to avoid stress isolation and hypothermia induced by our long-term caloric restriction protocol. They were kept
in a pathogen-free barrier facility maintained at 21.5 °C with a 12:12 h dark-light cycle (lights on at 07:30 a.m.). During one week of acclimation, mice were weighed every day to become used to handling and had free access to water and standard chow diet (4.30% fat, 22.30% protein and 51.20% carbohydrate; Special Diet Service RM3; Dietex, Essex, UK). Mice were randomized into three groups: The experimental group "food restriction and wheel" (group FRW), mice were placed in a cage equipped with a free running wheel (INTELLIBIO, Seichamps, France) and exposed first to a 30% quantitative food restriction for three days then followed by 50% quantitative food restriction for all the duration of the protocol. This restriction was calculated every day from the total food consumed by each mouse in the group of mice fed ad libitum (group AL) the previous day, by weighing the whole pellets in the feeder. Food (one pellet per mouse) was distributed directly into the cage every day at 6:30 p.m. This restriction was also applied to another experimental group: A pair-fed group or group FR (food-restricted). Body weights were monitored. Control of body weight was conducted to ensure that two mice in the same cage had similar feeding behavior and activity. If evident differences were noted between the two mice, the data obtained from this cage were excluded from analysis. The daily locomotor activity of FRW mice was assessed in their home cage equipped with a wheel (diameter: 230 mm; width: 50 mm; 1 revolution = 0.72 m) and linked to a computer system (ActiviWheel Software; INTELLIBIO, Seichamps, France). No mice died in any protocol. All experiments were carried out in accordance with the European Communities Council Directive (86/609/EEC) and approved by the Regional Ethics Committee of Nord-Pas de Calais of Lille, France (protocol no. CEEA 392012).

### 2.2.2. Protocols

Three protocols were performed. At the end of each protocol, the mice were sacrificed. The short-term protocol (ST) consisted of food restriction for two weeks (FR, *n* = 6; FRW, *n* = 6; AL, *n* = 6; 3 cages per group) while in the long-term protocol (LT), the food restriction was maintained for ten weeks (FR, *n* = 8; FRW, *n* = 8; AL, *n* = 8; 4 cages per group). Finally, in another set of animals, a two-week nutritional recovery phase was performed, following the long-term protocol (LT-Rec; FR, *n* = 6; FRW, *n* = 6; AL, *n* = 6; 3 cages per group). For this purpose, FR and FRW had ad libitum access to standard diet. Free access to the running wheel was maintained for FRW mice. Body weight was measured every two days in all protocols while food intake was measured daily. Body weight gain was calculated to evaluate the percentage of loss between D0 and D65 (100-((BW at D65/BW at D0) × 100)) and the percentage of gain after two weeks of recovery (100-((BW at D86/BW at D65) × 100)). The locomotor activity was assessed daily from the running wheel linked to a computer system that measured interval counts (10 min) per mean wheel revolution (ActiviWheel Software; INTELLIBIO, Seichamps, France). Data were extracted with an excel macro (Microsoft Office Standard, 2016) to obtain cumulative activity or day/night activity.

### 2.2.3. Euthanasia and Tissue Collection

At the end of the protocols, mice were sacrificed in the morning between 8:00 a.m. and 11:00 a.m. They were deeply anesthetized with an overdose of ketamine (100 mg/kg) and xylazine mix (20 mg/kg). Blood was collected through cardiac puncture with a 1 mL syringe and kept at 4 °C for ~2 h in neutral tubes containing p-hydroxy-mercuribenzoic acid (PHMB 0.4 mM final), a serine protease inhibitor, until centrifugation (8000 rpm for 10 min, 4 °C, Centrifuge 5414 R; Eppendorf, Hamburg, Germany). Samples were rapidly centrifuged (1000 × g for 10 min, 4 °C), to gain plasma aliquots. An aliquot of around 100 µL was immediately acidified with HCl (0.1 N final) to preserve ghrelin acylation; the residual plasma was kept for other assays. Plasma aliquots were then frozen in dry ice before being stored at -80 °C until they were assayed. From brains, the hypothalamus was carefully dissected, immediately frozen in liquid nitrogen and stored at -80 °C until analysis. Subcutaneous and perigonadic adipose tissues were gently dissected and weighed.
### 2.2.4. Immunocytochemistry

In another set of experiments, mice were deeply anesthetized with an overdose of ketamine (100 mg/kg) and xylazine mix (20 mg/kg) after a short-term food restriction protocol (FR, *n* = 2; FRW, *n* = 2; AL, *n* = 2) and transcardially perfused with 100 mL of 0.9% saline followed by 100 mL of 4% paraformaldehyde in 0.1 mol/L phosphate-buffered solution (PBS, pH 7.4). Brains were removed, postfixed at 4 °C overnight in 4% paraformaldehyde in PBS and frozen in dry ice before processing. Immunocytochemical procedure was previously described in Bouret et al., 2004. Coronal hypothalamic sections (20 µm) were obtained with a Cryostat (Leica) and collected in PBS. Sections were incubated for 48 hours in a rabbit anti-AgRP (polyclonal rabbit antiserum against AgRP, 1:4000; Phoenix Pharmaceutical, Belmont, CA). The primary antibodies were localized with a biotinylated goat antirabbit IgG (1:600, 90 min, Vector Laboratories). Tyramide signal amplification (TSA) was accomplished by placing the sections in an avidin-biotin solution (Vectastain) for 1 h, followed by incubation in TSA solution for 20 min. according to the manufacturer's instruction (TSA-Indirect kit, New England Nuclear Life Science). Deposited biotin was detected with Alexa 488-conjugated streptavidin (Molecular Probes). Sections were mounted and coverslipped with buffered glycerol (pH 8.5). Quantification of the fiber density was performed in the paraventricular nucleus (PVN, main location of ARC AgRP projections) using a Leica SP-confocal microscope equipped with a 10× objective (numerical aperture, 0.40; working distance, 360 µm). The density of immunolabeling was evaluated with Image J (NIH, 1987, https://imagej.nih.gov/ij/download.html (accessed on 12 August 2021) for 3 sections per nucleus per animal. The procedure was carried out on each image plane in the stack, and the values for all image planes in a stack summed. The resulting value was an estimate of the AgRP fiber density in the PVN.

### 2.2.5. Blood Assays

All samples were analyzed in duplicates. Plasma leptin was measured using an ELISA kit (R&D Systems Quantikine Europe, Abingdon, UK). Intra- and interassay coefficients of variations were 4.4 and <4.7%, respectively. Plasma acyl-ghrelin and des-acyl ghrelin concentrations were evaluated by specific EIA (A05118 for the acylated form and A05117 for the des-acyl form; Bertin Bioreagent, Montigny le Bretonneux, France). Intra- and inter-assay coefficients of variations were 6.1% and 5.7% for AG and 5.5% and 4.8% for DAG, respectively.

### 2.2.6. RT-PCR Extraction and Analysis

Total hypothalamus RNAs were extracted in guanidinium thiocyanate and phenol-chloroform according to Chomczynski and Sacchi method [58]. Two micrograms of RNA were treated with DNase I (Roche Diagnostics, Penzberg, Germany) then reverse-transcribed using SuperScript III (Invitrogen, Thermo Scientific, Waltham, MA, USA) and real-time PCR was performed using SYBR Green mix (2X, Roche Diagnostics, Meylan, France). The glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as internal control and housekeeping gene to normalize gene expression. The purity of the PCR products was assessed by dissociation curves. The amount of target cDNA was calculated by the comparative threshold (Ct) method and expressed by means of the $2^{-\Delta\Delta C_t}$ method. PCR primers were described in Table 1.
<table><caption>**Table 1.** List of PCR primers. AgRP, agouti-related peptide; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; GHSR, ghrelin receptor, Lept-R, leptin receptor; MC3-R/MC4-R, melanocortin 3 and 4 receptors; NPY, neuropeptide Y; POMC, proopiomelanocortin; Y1-R/Y5-R, neuropeptide Y receptor 1/receptor 5.</caption><thead><tr><td>Gene</td><td>Reference</td><td>Forward Primer</td><td>Reverse Primer</td></tr></thead><tbody><tr><td>AgRP</td><td>U89486.1</td><td>5'-CCCAGAGTTCCCAGGTCTAAGTCT-3'</td><td>5'-CACCTCCGCCAAAGCTTCT-3'</td></tr><tr><td>NPY</td><td>NM_023456.2</td><td>5'-CCGCTCTGCGACACTACAT-3'</td><td>5'-TGTCTCAGGGCTGGATCT-3'</td></tr><tr><td>POMC</td><td>BC061215.1</td><td>5'-AGTGCCAGGACCTCACCA-3'</td><td>5'-CAGCGAGAGGTCGAGTTTG-3'</td></tr><tr><td>MC3R</td><td>NM_008561</td><td>5'-TCCGATGCTGCCTAACCTCT-3'</td><td>5'-GGATGTTTTCCATCAGACTGACG-3'</td></tr><tr><td>MC4R</td><td>NM_016977</td><td>5'-CCCGGACGGAGGATGCTAT-3'</td><td>5'-TCGCCACGATCACTAGAATGT-3'</td></tr><tr><td>Y1R</td><td>NM_010934</td><td>5'-TGATCTCCACCTGCGTCAAC-3'</td><td>5'-ATGGCTATGGTCTCGTAGTCAT-3'</td></tr><tr><td>Y2R</td><td>NM_008731</td><td>5'-GCCAGGGCACACTACTCCTA-3'</td><td>5'-CTACCCCTAGCAAGATGATGGA-3'</td></tr><tr><td>Y5R</td><td>NM_016708.3</td><td>5'-CATCTCAAGCAGAAGCGACC-3'</td><td>5'-CTCCATACTAGAGTCCTCGGG-3'</td></tr><tr><td>Lept-R</td><td>NM_010704.2</td><td>5'-GTCTTCGGGATGTGAATGTC-3'</td><td>5'-ACCTAAGGGTGGATCGGGTTT-3'</td></tr><tr><td>GHSR</td><td>NM_177330</td><td>5'-TGGAGATCGCGCAGATCAG-3'</td><td>5'-CCGGGAACTCTCATCCTTCAG-3'</td></tr><tr><td>GADPH</td><td>GU214026</td><td>5'-GAACATCATCCCTGCATCC-3'</td><td>5'-CCAGTGAGCTTCCCGTTCA-3'</td></tr></tbody></table>

## 2.3. Statistical Analysis

The clinical data were expressed as means ± SD and analyzed with jamovi (the jamovi project 2021, version 1.6.23). Analysis of normality of variances was tested by Shapiro-Wilk test. Mann-Whitney test was used to compare aAN and prAN populations. Spearman's rank correlation coefficient was used to analyze correlations between clinical dimensions (Shapiro-Wilk *p* < 0.05). For the cardiovascular data, raw data were analyzed through a homemade algorithm using Python. The preclinical data were expressed as means ± SEM. GraphPad Prism 5.01 (Abacus Concepts, Berkeley, CA, USA) was used to analyze the data and generate graphs. Analysis of normality was tested by Shapiro-Wilk test. Equality of variance was tested by Levene's test. Statistical analysis was performed using one-way ANOVA or two-way ANOVA when appropriate, followed by a Tukey post hoc test when the *p* value of the ANOVA was significant (*p* < 0.05). When equality of variance was violated, Games-Howell post hoc test was used. We eliminated three individual values related to mRNA expression of GHSR (group AL—short term; group FR—long term) and of leptin receptor (group FRW-LT recovery) using the z-score test which deciphers whether value is aberrant or not (Z-score = (value - mean) / SD)).

# 3. Results

## 3.1. Experiment 1. Preliminary Data: Clinical Investigation of Partial Recovery

### 3.1.1. Clinical Characteristics: Acute AN versus Partially Recovered AN Patients

We recruited 266 patients suffering from acute (aAN, *n* = 225) or partial recovery (prAN, *n* = 41) AN. Illness duration was similar between aAN and prAN (U = 3274, *p* = 0.199), as was the level of anxiety (U = 2969, *p* = 0.069), depression (U = 3485, *p* = 0.663), obsession (U = 3061, *p* = 0.117) and compulsion (U = 3410, *p* = 0.525). On the other hand, patients with aAN had a lower minimal BMI (U = 1956, *p* < 0.001), a lower total EDI-2 score (U = 3133, *p* = 0.007) especially with the EDI-2 sub-scores, drive for thinness (U = 3030, *p* = 0.004), bulimia (U = 3290, *p* = 0.017), body dissatisfaction (U = 2750, *p* < 0.001), perfectionism (U = 3407, *p* = 0.044), interoceptive awareness (U = 3143, *p* = 0.008) and emotional dysregulation (U = 3185, *p* = 0.011). Lastly, albumin (U = 1670, *p* < 0.001) and ferritin (U = 845, *p* < 0.001) were higher in patients with prAN (Table 2).
<table><caption>**Table 2.** Clinical characteristics of Anorexia Nervosa patients (AN) and partially recovered AN (prAN). Significant *p*-values are in bold.</caption><thead><tr><th rowspan="2"></th><th colspan="2">AN (n = 225)</th><th colspan="2">prAN (n = 41)</th><th colspan="2">AN vs. prAN</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>27.8</td><td>9.76</td><td>27.2</td><td>8.43</td><td>4558</td><td>0.904</td></tr><tr><td>BMI</td><td>15.1</td><td>1.19</td><td>20.5</td><td>1.81</td><td>0</td><td><strong>&lt;0.001</strong></td></tr><tr><td>BMI min</td><td>13.4</td><td>3.07</td><td>15</td><td>1.34</td><td>1956</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Illness duration</td><td>8.9</td><td>8.8</td><td>10.3</td><td>8.6</td><td>3274</td><td>0.199</td></tr><tr><td>EDI</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Drive for thinness</td><td>10.5</td><td>6.95</td><td>13.8</td><td>6.33</td><td>3030</td><td><strong>0.004</strong></td></tr><tr><td>Bulimia</td><td>3.65</td><td>5.35</td><td>5.92</td><td>6.37</td><td>3290</td><td><strong>0.017</strong></td></tr><tr><td>Body dissatisfaction</td><td>12.1</td><td>7.16</td><td>16.9</td><td>8.05</td><td>2750</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Ineffectiveness</td><td>11.6</td><td>7.57</td><td>12.7</td><td>6.52</td><td>3890</td><td>0.383</td></tr><tr><td>Perfectionism</td><td>7.52</td><td>4.76</td><td>9.15</td><td>4.51</td><td>3407</td><td><strong>0.044</strong></td></tr><tr><td>Interpersonal distrust</td><td>6.29</td><td>4.59</td><td>6</td><td>4.56</td><td>4016</td><td>0.564</td></tr><tr><td>Interoceptive awareness</td><td>10.5</td><td>6.9</td><td>13.7</td><td>7.05</td><td>3143</td><td><strong>0.008</strong></td></tr><tr><td>Maturity fears</td><td>6.87</td><td>5.95</td><td>7.35</td><td>6.22</td><td>4040</td><td>0.604</td></tr><tr><td>Ascetism</td><td>7.73</td><td>4.95</td><td>8.68</td><td>4.49</td><td>3666</td><td>0.161</td></tr><tr><td>Emotional dysregulation</td><td>5.44</td><td>5.49</td><td>8.25</td><td>6.72</td><td>3185</td><td><strong>0.011</strong></td></tr><tr><td>Social insecurity</td><td>7.84</td><td>5.23</td><td>8.4</td><td>4.26</td><td>3904</td><td>0.402</td></tr><tr><td>Total</td><td>89.62</td><td>46.57</td><td>110.97</td><td>41.14</td><td>3133</td><td><strong>0.007</strong></td></tr><tr><td>HADS Anxiety</td><td>12.5</td><td>4.33</td><td>13.9</td><td>3.45</td><td>2969</td><td>0.069</td></tr><tr><td>HADS Depression</td><td>8.59</td><td>4.53</td><td>8.29</td><td>4.34</td><td>3485</td><td>0.663</td></tr><tr><td>YBS Compulsion</td><td>7.84</td><td>5.37</td><td>8.74</td><td>4.27</td><td>3410</td><td>0.525</td></tr><tr><td>YBS Obsession</td><td>8.93</td><td>5.01</td><td>10.6</td><td>4.16</td><td>3061</td><td>0.117</td></tr><tr><td>Ferritin (ng/mL)</td><td>81</td><td>75.1</td><td>30.8</td><td>17</td><td>845</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Albumin (g/L)</td><td>47.7</td><td>6.21</td><td>43.7</td><td>8.73</td><td>1670</td><td><strong>&lt;0.001</strong></td></tr></tbody></table>

AN: Anorexia Nervosa; BMI: body mass index; EDI: Eating Disorder Inventory II; HADS: Hospital Anxiety and Depression scale; prAN: partially recovered Anorexia Nervosa; YBS: Yale Brown Scale.

### 3.1.2. Continuous Cardiac Monitoring, Past and Actual BMI and Illness Duration

Among the 266 patients, 68 had a 72 h continuous heart rate monitoring: 55 aAN and 13 prAN. This measure was used to assess the level of bradycardia and of physical activity. Comparison between the sample with continuous cardiac monitoring and the total sample did not show significant differences except for anxiety in the prAN group (10.83 ± 4.24 with monitoring versus 13.87 ± 3.45, *p* = 0.019) and illness duration in the AN group (10.45 ± 7.36 with monitoring versus 9.09 ± 8.76; *p* = 0.032) and aAN group (10.54 ± 7.62 with monitoring versus 8.88 ± 8.79; *p* = 0.034, Table 3).

<table><caption>**Table 3.** Representativity of the sample with cardiac monitoring. Comparison of clinical characteristics of all patients with anorexia nervosa (AN), acute AN (aAN) or partially recovered AN (prAN) and the sample with continuous cardiac monitoring. Significant *p*-values are in bold.</caption><thead><tr><th rowspan="2"></th><th colspan="2">AN (n = 266) vs. AN with Cardiac Monitoring (n = 68)</th><th colspan="2">aAN (n = 225) vs. aAN with Cardiac Monitoring (n = 55)</th><th colspan="2">prAN (n = 41) vs. prAN with Cardiac Monitoring (n = 13)</th></tr><tr><th>U</th><th>p</th><th>U</th><th>p</th><th>U</th><th>P</th></tr></thead><tbody><tr><td>Age</td><td>7656</td><td>0.108</td><td>5030</td><td>0.050</td><td>218</td><td>0.559</td></tr><tr><td>BMI</td><td>7786</td><td>0.110</td><td>5256</td><td>0.124</td><td>209</td><td>0.248</td></tr><tr><td>BMI min</td><td>6640</td><td>0.371</td><td>4258</td><td>0.338</td><td>265</td><td>0.984</td></tr><tr><td>Illness duration</td><td>6721</td><td><strong>0.032</strong></td><td>4517</td><td><strong>0.034</strong></td><td>220</td><td>0.759</td></tr><tr><td>EDI</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Drive for thinness</td><td>7989</td><td>0.908</td><td>5456</td><td>0.910</td><td>224</td><td>0.727</td></tr><tr><td>Bulimia</td><td>7936</td><td>0.889</td><td>5408</td><td>0.825</td><td>232</td><td>0.860</td></tr><tr><td>Body dissatisfaction</td><td>7483</td><td>0.373</td><td>4968</td><td>0.270</td><td>238</td><td>0.965</td></tr><tr><td>Ineffectiveness</td><td>7927</td><td>0.834</td><td>5408</td><td>0.834</td><td>237</td><td>0.957</td></tr><tr><td>Perfectionism</td><td>7698</td><td>0.574</td><td>5149</td><td>0.462</td><td>215</td><td>0.586</td></tr><tr><td>Interpersonal distrust</td><td>6987</td><td>0.098</td><td>4841</td><td>0.172</td><td>195</td><td>0.326</td></tr><tr><td>Interoceptive Awareness</td><td>7784</td><td>0.668</td><td>5447</td><td>0.896</td><td>157</td><td>0.071</td></tr><tr><td>Maturity fears</td><td>7664</td><td>0.539</td><td>5170</td><td>0.488</td><td>240</td><td>1.000</td></tr><tr><td>Ascetism</td><td>7816</td><td>0.704</td><td>5500</td><td>0.981</td><td>191</td><td>0.290</td></tr><tr><td>Emotional Dysregulation</td><td>7170</td><td>0.169</td><td>5120</td><td>0.425</td><td>176</td><td>0.167</td></tr><tr><td>Social insecurity</td><td>7430</td><td>0.331</td><td>4986</td><td>0.285</td><td>232</td><td>0.861</td></tr><tr><td>Total</td><td>7552</td><td>0.460</td><td>5200</td><td>0.560</td><td>211</td><td>0.529</td></tr><tr><td>HADS Anxiety</td><td>6356</td><td>0.104</td><td>4672</td><td>0.512</td><td>125</td><td><strong>0.019</strong></td></tr><tr><td>HADS Depression</td><td>7043</td><td>0.634</td><td>4890</td><td>0.865</td><td>192</td><td>0.418</td></tr><tr><td>YBS Compulsion</td><td>7126</td><td>0.697</td><td>4948</td><td>0.922</td><td>200</td><td>0.531</td></tr><tr><td>YBS Obsession</td><td>6751</td><td>0.310</td><td>4610</td><td>0.397</td><td>207</td><td>0.629</td></tr></tbody></table>

aAN: acute Anorexia Nervosa; AN: Anorexia Nervosa; BMI: body mass index; EDI: Eating Disorder Inventory II; HADS: Hospital Anxiety and Depression scale; prAN: partially recovered Anorexia Nervosa; YBS: Yale Brown Scale.
We found a negative correlation between the lowest BMI and high intensity physical activity (rho = -0.271, *p* = 0.029). Illness duration was negatively correlated to bradycardia (rho = -0.298, *p* = 0.015) and positively correlated to moderate physical activity (rho = 0.277, *p* = 0.024). All correlations for the AN group are reported in Table 4.

<table><caption>**Table 4.** Non parametric correlation for BMI minimum, current BMI, illness duration and heart rate ranges in all patients with (AN *n* = 68). Significant *p*-values are in bold. BMI: body mass index; EDI: Eating Disorder Inventory II; HADS: Hospital Anxiety and Depression scale; YBS: Yale Brown Scale.</caption><thead><tr><th></th><th></th><th>BMI min</th><th>Current BMI</th><th>Illness Duration</th><th>Bradycardia</th><th>Physiological</th><th>Moderate</th><th>High Intensity</th></tr></thead><tbody><tr><td rowspan="2">BMI min</td><td>Rho</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Current BMI</td><td>Rho</td><td>0.538</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td><strong>&lt;0.001</strong></td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Illness duration</td><td>Rho</td><td>-0.131</td><td>-0.040</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.293</td><td>0.746</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Bradycardia</td><td>Rho</td><td>0.079</td><td>-0.143</td><td>-0.298</td><td>—</td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.531</td><td>0.251</td><td><strong>0.015</strong></td><td>—</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Physiological</td><td>Rho</td><td>0.168</td><td>0.322</td><td>0.158</td><td>-0.642</td><td>—</td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.181</td><td><strong>0.008</strong></td><td>0.206</td><td><strong>&lt;0.001</strong></td><td>—</td><td></td><td></td></tr><tr><td rowspan="2">Moderate</td><td>Rho</td><td>-0.234</td><td>-0.032</td><td>0.277</td><td>-0.733</td><td>0.343</td><td>—</td><td></td></tr><tr><td><i>p</i>-value</td><td>0.061</td><td>0.800</td><td><strong>0.024</strong></td><td><strong>&lt;0.001</strong></td><td><strong>0.005</strong></td><td>—</td><td></td></tr><tr><td rowspan="2">High intensity</td><td>Rho</td><td>-0.271</td><td>-0.084</td><td>0.213</td><td>-0.557</td><td>-0.032</td><td>0.669</td><td>—</td></tr><tr><td><i>p</i>-value</td><td><strong>0.029</strong></td><td>0.502</td><td>0.087</td><td><strong>&lt;0.001</strong></td><td>0.797</td><td><strong>&lt;0.001</strong></td><td>—</td></tr></tbody></table>

We then analyzed the differences between patients who have reached the normal BMI range for a sustained period (prAN, *n* = 13) and acute patients (aAN, *n* = 55). Compared to the prAN group, the aAN group had more time with bradycardia (31.5% versus 13.8%; U = 200, *p* = 0.02), less time with a physiological heart rate (51.2% versus 65.4%; U = 183, *p* = 0.009) and an equivalent time at moderate (10.1% versus 11.1%; U = 261, *p* = 0.181) or high (U = 262, *p* = 0.186) physical activity (Figure 1).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Data obtained from 72 h cardiovascular recording allowing to determine the percentage of time spent in bradycardia, physiological, moderate and high physical activity in patients with acute anorexia nervosa (aAN, *n* = 55) and partially recovered AN (prAN, *n* = 13). bpm: beat per minute; HRmax: maximum heart rate (220-age). * *p* < 0.05; ** *p* < 0.01.</figcaption></figure>
In aAN, lifetime minimal BMI, but not current BMI, was negatively correlated to a high level (rho = -0.332, *p* = 0.016) and moderate level (rho = -0.279, *p* = 0.045) physical activity. Illness duration was negatively correlated to bradycardia (rho = -0.423, *p* = 0.002) and positively correlated to moderate (rho = 0.324, *p* = 0.018) and high (rho = 0.307, *p* = 0.025) intensity of physical activity (Table 5). In the prAN group, only a trend for a negative correlation was found between minimal BMI and high intensity of physical activity (rho = -0.525, *p* = 0.065, Table 6).

<table><caption>**Table 5.** Non parametric correlation for BMI minimum, current BMI, illness duration and heart rate ranges in acute AN patients (n = 55). Significant *p*-values are in bold. Trend values are in italics. BMI: body mass index; EDI: Eating Disorder Inventory II; HADS: Hospital Anxiety and Depression scale; YBS: Yale Brown Scale.</caption><thead><tr><th></th><th></th><th>BMI min</th><th>Current BMI</th><th>Illness Duration</th><th>Brady-Cardia</th><th>Physiological</th><th>Moderate</th><th>High Intensity</th></tr></thead><tbody><tr><td rowspan="2">BMI min</td><td>Rho</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Current BMI</td><td>Rho</td><td>0.620</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td><strong>0.001</strong></td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Illness duration</td><td>Rho</td><td>-0.186</td><td>-0.073</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.182</td><td>0.600</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Bradycardia</td><td>Rho</td><td>0.121</td><td>0.081</td><td>-0.423</td><td>—</td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.394</td><td>0.566</td><td><strong>0.002</strong></td><td>—</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Physiological</td><td>Rho</td><td>0.050</td><td>0.173</td><td>0.195</td><td>-0.670</td><td>—</td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.725</td><td>0.217</td><td>0.161</td><td><strong>&lt;0.001</strong></td><td>—</td><td></td><td></td></tr><tr><td rowspan="2">Moderate</td><td>Rho</td><td>-0.279</td><td>-0.188</td><td>0.324</td><td>-0.730</td><td>0.342</td><td>—</td><td></td></tr><tr><td><i>p</i>-value</td><td><strong>0.045</strong></td><td>0.177</td><td><strong>0.018</strong></td><td><strong>&lt;0.001</strong></td><td><strong>0.013</strong></td><td>—</td><td></td></tr><tr><td rowspan="2">High intensity</td><td>Rho</td><td>-0.332</td><td>-0.269</td><td>0.307</td><td>-0.535</td><td>-0.059</td><td>0.677</td><td>—</td></tr><tr><td><i>p</i>-value</td><td><strong>0.016</strong></td><td>0.052</td><td><strong>0.025</strong></td><td><strong>&lt;0.001</strong></td><td>0.675</td><td><strong>&lt;0.001</strong></td><td>—</td></tr></tbody></table>

<table><caption>**Table 6.** Non parametric correlation for BMI min, current BMI, illness duration and heart rate ranges in partially recovered AN patients (n = 13). Trend values are in italics. BMI: body mass index; EDI: Eating Disorder Inventory II; HADS: Hospital Anxiety and Depression scale; YBS: Yale Brown Scale.</caption><thead><tr><th></th><th></th><th>BMI min</th><th>Current BMI</th><th>Illness Duration</th><th>Brady-cardia</th><th>Physiological</th><th>Moderate</th><th>High Intensity</th></tr></thead><tbody><tr><td rowspan="2">BMI min</td><td>Rho</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Current BMI</td><td>Rho</td><td>-0.234</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.442</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Illness duration</td><td>Rho</td><td>0.306</td><td>-0.113</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.310</td><td>0.714</td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Bradycardia</td><td>Rho</td><td>0.479</td><td>-0.126</td><td>-0.415</td><td>—</td><td></td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.098</td><td>0.683</td><td>0.158</td><td>—</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Physiological</td><td>Rho</td><td>0.193</td><td>0.297</td><td>-0.116</td><td>-0.462</td><td>—</td><td></td><td></td></tr><tr><td><i>p</i>-value</td><td>0.528</td><td>0.325</td><td>0.707</td><td>0.115</td><td>—</td><td></td><td></td></tr><tr><td rowspan="2">Moderate</td><td>Rho</td><td>-0.303</td><td>-0.121</td><td>0.072</td><td>-0.363</td><td>0.044</td><td>—</td><td></td></tr><tr><td><i>p</i>-value</td><td>0.315</td><td>0.696</td><td>0.816</td><td>0.224</td><td>0.892</td><td>—</td><td></td></tr><tr><td rowspan="2">High intensity</td><td>Rho</td><td>-0.525</td><td>-0.099</td><td>-0.487</td><td>-0.549</td><td>-0.198</td><td>0.275</td><td>—</td></tr><tr><td><i>p</i>-value</td><td>0.065</td><td>0.751</td><td>0.091</td><td>0.055</td><td>0.517</td><td>0.363</td><td>—</td></tr></tbody></table>

## 3.2. Experiment 2. Preclinical Investigation

### 3.2.1. Body Weight, Food Intake, Physical Activity and Fat Mass Evolution

At D0, body weight was not significantly different between all groups (n = 6 per group): AL: 18.22 g ± 0.29; FR: 18.22 g ± 0.24; FRW: 18.03 g ± 0.12, F(2,15) = 0.21, *p* = 0.815. The day before nutritional recovery, at day 65, the body weight of FR and FRW mice was significantly decreased compared to AL mice, respectively (14.07 g ± 0.14 and 14.65 g ± 0.09 versus 21.42 g ± 0.37; F(2,15) = 290.2, *p* < 0.001; Figure 2A). At the end of the food restriction protocol (D65), the body weight gain calculated from D0 was -22.8% and -18.9% for FR and FRW mice, respectively, with a significant difference (F(2,15) = 280.9, *p* < 0.001) between AL and FR mice (*p* < 0.001, Figure 2B).
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Body weight evolution (A–C), food intake during nutritional recovery (D) and daily physical activity (E,F) in the three protocols: short-term, long-term, long-term and nutritional recovery. AL: *ad libitum*; FR: food restriction; FRW: food restriction and wheel. Data represent mean ± SEM (n = 6–8 per group); \* *p* < 0.05 AL vs. FR; \# *p* < 0.05 FR vs. FRW. Ø, *p* < 0.05 Day 66 vs. Day 69.</figcaption></figure>
The nutritional recovery led to a rapid gain in body weight reaching +55.9% and +48.8% for FR and FRW mice, respectively, at the end of the protocol (D86, Figure 2C) with a significant difference ($F(2,15) = 165, p < 0.001$) between AL and FR mice ($p < 0.001$, Figure 2C). However, despite a full recovery in body weight within 3 days of ad libitum feeding (D68), we noticed that FRW mice showed a transient body weight reduction between D68 and D71 which was significantly different compared to FR mice with an interaction time × group ($F(12,90) = 67,19, p < 0.001$, Figure 2A). Indeed, AL mice showed a significantly higher body weight than FR mice at D65, D66, D67 ($p < 0.005$), and FR mice gained significantly more weight than FRW mice at D69, D70, D71 ($p < 0.05$). Such difference can be attributed to the quantity of food ingested during all the nutritional recovery period (Figure 2D). Two-way ANOVA indicated a protocol effect ($F(5,10) = 7.9, p = 0.003$) and a group effect ($F(2,10) = 4.7, p = 0.04$) with a significant difference at D66 ($p = 0.0004$), D67 ($p = 0.03$), D68 ($p = 0.003$), D69 ($p = 0.002$) and D70 ($p = 0.01$) between AL and FR mice. Surprisingly, the binge-eating behavior in the FR and FRW groups followed different patterns, with a longer duration for the FR group at D69 ($p = 0.01$) and D70 ($p = 0.005$).

FRW mice exhibited 24 h running wheel activity which decreased along the long-term protocol (Figure 2E, $F(9,52) = 5.57, p < 0.0001$). However, the percentage of time spent in activity during the day increased progressively in the long-term protocol (Figure 2F). Along the two weeks of nutritional recovery, the FRW mice showed an increased activity at D69 compared to the first day of recovery (D66, $p < 0.05$) which could explain the transient decreases in food intake and body weight (Figure 2E). However, nutritional recovery was accompanied by a decrease in the percentage spent in diurnal running wheel activity to a complete cessation of this activity at D81 (Figure 2F).

Finally, in parallel with the changes in feeding behavior and physical activity, FR and FRW mice showed a differential recovery of their fat mass (Figure 3A,B). Indeed, for both subcutaneous (SCAT, Figure 3A) and perigonadic (PG, Figure 3B) fat mass weight, two-way ANOVA indicated a group effect (SCAT, $F(2,60) = 23.75, p < 0.0001$; PG, $F(2,61) = 36.88, p < 0.0001$), a protocol duration effect (SCAT, $F(2,60) = 8.94, p = 0.0004$; PG, $F(2,61) = 47.07, p < 0.0001$), a trend for interaction protocols × groups ($F(4,60) = 2.46, p = 0.054$) for the SCAT and an interaction protocols × groups ($F(4,61) = 3.53, p = 0.012$) for PG. In the short-term protocol, and FR mice had a significant decrease in SCAT ($p = 0.0004$) and PG ($p < 0.001$) compared to AL mice. This decrease was higher in FRW mice than FR mice only for PG (trend, $p = 0.09$). The long-term protocol exacerbated these changes in fat mass (AL versus FR: SCAT, $p = 0.0003$; PG, $p = 0.0002$) and a trend only for SCAT (FR versus FRW, $p = 0.054$). Finally, only FR mice fully recovered their PG fat mass after two weeks of refeeding ($p = 0.004$).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Adipose tissue weight. (A) Subcutaneous fat mass and (B) perigonadic fat mass weighed in the three protocols: short-term, long-term, long-term (LT) and nutritional recovery. AL: *ad libitum*; FR: food restriction; FRW: food restriction and wheel. Data represent mean ± SEM ($n = 6–8$ per group); ** $p < 0.01$ AL vs. FR; ## $p < 0.01$ FR vs. FRW.</figcaption></figure>
### 3.2.2. Ghrelin and Leptin Assays

Plasma ghrelin and leptin concentrations were evaluated in the three protocols as nutritional sensors (Table 7). Two-way ANOVA analyses exhibited a protocol effect for AG (F(2,43) = 23.78, *p* < 0.001) and DAG (F(2,43) = 6.41, *p* = 0.004) but not for leptin (F(2, 51) = 0.99; *p* = 0.3785), a group effect for AG (F(2,43) = 4.85, *p* = 0.013), DAG (F(2,43) = 6.95; *p* = 0.002) and leptin (F(2,51) = 22.28, *p* < 0.0001) and an interaction between protocols and groups for AG (F(4,43) = 5.19, *p* = 0.002) DAG (F(4,43) = 2.41, *p* = 0.064) but not for leptin (F(4, 51) = 1.55; *p* = 0.203).

In the short-term protocol, AG and DAG plasma concentrations were significantly higher in FR mice compared to AL mice (AG, *p* = 0.002; DAG, *p* = 0.013), but no significant difference seen between FR and FRW mice (AG, *p* = 0.133; DAG, *p* = 0.772). Finally, leptin plasma concentrations were significantly lower for FR mice (*p* = 0.029) compared to AL mice, but no difference was seen between FR and FRW mice (*p* = 0.811). In the long-term protocol, AG and DAG plasma concentrations were not significantly different (AG, AL versus FR: *p* = 0.936, FR versus FRW: *p* = 0.225; DAG, AL versus FR: *p* = 0.130, FR versus FRW, *p* = 0.910). However, leptin plasma concentrations were lower for FR mice compared to AL mice (*p* = 0.006) and similar for FR mice compared to FRW mice (*p* = 0.912). After nutritional recovery, mice exhibited no differences of AG and DAG plasma concentration (AG, AL versus FR: *p* = 0.902, FR versus FRW: *p* = 0.774; DAG, AL versus FR: *p* = 0.577, FR versus FRW, *p* = 0.449) whereas leptin plasma concentration remained lower for FR mice compared to AL mice (*p* = 0.029) and similar for FR mice compared to FRW mice (*p* = 0.106).

<table><caption>**Table 7.** Acyl ghrelin, desacyl ghrelin and leptin plasma levels in short-term, long-term and recovery mice protocols. Significant values are in bold. Data are expressed in mean +/− SEM or *p*-value (*p*). AL: *ad libitum*; FR: food restricted; FRW: food restricted and wheel. AL versus FR: * *p* < 0.05; ** *p* < 0.01.</caption><thead><tr><th></th><th>Acyl-Ghrelin</th><th>Desacyl-Ghrelin</th><th>Leptin</th></tr></thead><tbody><tr><td><strong>Short term</strong></td><td></td><td></td><td></td></tr><tr><td>AL</td><td>229+/-35.4</td><td>679+/-79.8</td><td>5724+/-866</td></tr><tr><td>FR</td><td>576+/-63.7 **</td><td>1583+/-160.5 *</td><td>3548+/-248.9 *</td></tr><tr><td>FRW</td><td>958+/-165.7</td><td>1763+/-252.7</td><td>3080+/-202.7</td></tr><tr><td>Group effect (p)</td><td><strong>&lt;0.001</strong></td><td><strong>&lt;0.001</strong></td><td><strong>0.007</strong></td></tr><tr><td><strong>Long term</strong></td><td></td><td></td><td></td></tr><tr><td>AL</td><td>166+/-42.5</td><td>1085+/-132.4</td><td>4668+/-419.8</td></tr><tr><td>FR</td><td>147+/-20.3</td><td>1459+/-108.4</td><td>3365+/-155.5 **</td></tr><tr><td>FRW</td><td>240+/-41.5</td><td>1534+/-132.2</td><td>3518+/-88.3</td></tr><tr><td>Group effect (p)</td><td>0.208</td><td><strong>0.045</strong></td><td><strong>0.004</strong></td></tr><tr><td><strong>Long term-Recovery</strong></td><td></td><td></td><td></td></tr><tr><td>AL</td><td>150+/-14.6</td><td>790+/-13.7</td><td>4594+/-380.7</td></tr><tr><td>FR</td><td>137+/-27.3</td><td>913+/-125.7</td><td>3639+/-95.4 *</td></tr><tr><td>FRW</td><td>115+/-14</td><td>728+/-40.7</td><td>2915+/-100.9</td></tr><tr><td>Group effect (p)</td><td>0.633</td><td>0.433</td><td><strong>&lt;0.001</strong></td></tr></tbody></table>

In a longitudinal perspective, AL mice displayed the same profile for AG, DAG and leptin plasma concentration in the three protocols (AG, short-term versus long-term: *p* = 0.451; long-term versus LT-rec: *p* = 0.959; short-term versus LT-rec: *p* = 0.399. DAG, short-term versus long-term: *p* = 0.062, long-term versus LT-rec; *p* = 0.146, short-term versus LT-rec: *p* = 0.415. leptin, short-term versus long-term: *p* = 0.398; long-term versus LT-rec: *p* = 0.995; short-term versus LT-rec: *p* = 0.398). For FR mice, AG plasma concentrations were significantly higher in short-term compared to long-term (*p* < 0.001) and LT-Rec (*p* < 0.001) protocols and similar between long-term and LT-rec (*p* = 0.948). DAG plasma concentrations were significantly higher in both short-term (*p* = 0.021) and long-term (*p* = 0.032) protocols compared to the LT-Rec protocol and similar between short-term and long-term (*p* = 0.802). The absence of a protocol effect on leptin plasma concentrations did not allow a post hoc analysis. For FRW mice, AG plasma concentrations were higher in the short-term than in the long-term protocol
(p = 0.008) and LT-rec protocol (p = 0.003) with a trend of to be lower between long-term and LT-Rec protocols (p = 0.057). DAG plasma concentrations exhibited no differences between short-term and long-term protocol (p = 0.711) but were lower in LT-rec than in the short-term protocol (p = 0.011) and in the long-term protocol (p = 0.002). Leptin plasma concentrations were significantly lower in the LT-rec protocol compared to long-term protocol (p = 0.012) but exhibited no differences for short-term vs. long-term (p = 0.072) and LT-rec (p = 0.685).

### 3.2.3. Immunocytochemical Detection of AgRP

Considering the metabolic changes induced by the short-term protocol [32,49], we further investigated whether a chronic food restriction might induce noticeable changes in the immunocytochemical AgRP labeling distribution in ARC and PVN of the hypothalamus. Surprisingly, in the ARC, the FRW mice displayed AgRP-immunoreactive cell bodies not seen in AL and FR mice (Figure 4), whereas the AgRP fiber density was not significantly different in the PVN (data not shown, *p* = 0.6).

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Immunohistochemical detection of AgRP (agouti-related peptide) in hypothalamic arcuate nucleus (ARC), fibers (arrows) and cell bodies (double arrows) and paraventricular nucleus (PVN) in AL, FR and FRW groups submitted to a short-term protocol. The three groups did not show any significant changes in the AgRP fibers density in the PVN. AL: *ad libitum*; FR: food restricted; FRW: food restricted with wheel; 3rd V: 3rd ventricle.</figcaption></figure>
3.2.4. mRNA Expression of AgRP/NPY and Melanocortin Systems in the Hypothalamus

The changes in the immunolabeling of AgRP in FRW mice, associated with the varia-
tions in the plasma concentrations of ghrelin and leptin as well as a distinguished fat mass
recovery compared to FR mice led us to investigate whether these changes might lead to a
differential regulation in the expression of neuropeptides regulated by these hormones and
involved energy metabolism (Figure 5).

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Impact of short-term, long-term (LT) and nutritional recovery (LT-Rec) on the mRNA expression of hypothalamic of agouti-related peptides (A, AgRP), neuropeptide Y (B, NPY) and pro-opiomelanocortin (C, POMC). AL: *ad libitum*; FR: food restricted; FRW: food restricted with wheel. Data represent mean ± SEM (n = 4–8 per group). \* *p* < 0.05 AL vs. FR; \# *p* < 0.05 FR vs. FRW; \*\* *p* < 0.01 AL vs. FR; \*\*\* *p* < 0.01 FR vs. FRW.3.2.5. mRNA expression receptors of leptin, ghrelin and AgRP/NPY and melanocortin systems in the hypothalamus.</figcaption></figure>
Two-way ANOVA analysis indicated a time effect for AgRP and POMC mRNA expression (Figure 5A,C; AgRP, F(2,46) = 8.86, *p* = 0.0006; POMC, F(2,45) = 8.45, *p* = 0.0008), a group effect for AgRP, NPY and POMC mRNA expression (Figure 5A–C; AgRP, F(2,46) = 5.83, *p* = 0.005; NPY, F(2,47) = 6.51, *p* = 0.003; POMC, F(2,45) = 13.54, *p* < 0.0001) and an interaction protocol × group for AgRP and POMC and a tendency for NPY (AgRP, F(4,46) = 4.46, *p* = 0.003; NPY, F(4,47) = 2.29, *p* = 0.074; POMC, F(4,45) = 4.41, *p* = 0.004).

In the short-term protocol, AgRP and NPY mRNA expressions were significantly higher in FR mice compared to AL mice (*p* < 0.05) and in FRW compared to FR mice (trend for AgRP, *p* = 0.09). POMC mRNA expression was significantly decreased only for FRW mice (*p* = 0.007) compared to FR mice. The long-term protocol exacerbated significantly the changes in AgRP and NPY mRNA levels especially for the FR mice compared to AL mice (AgRP, *p* = 0.003; NPY, *p* = 0.002) and FRW mice (AgRP, *p* = 0.05; NPY, *p* = 0.04). Indeed, compared to AL mice, the long-term protocol led to a 10-fold and a 3-fold higher AgRP and NPY mRNA levels in FR mice, whereas the short-term protocol led to a 3-fold and a 2-fold higher expression, respectively. Finally, the POMC mRNA levels were significantly lower, especially for FR mice compared to AL mice (*p* < 0.001). Nutritionally recovered mice presented a normalization of NPY and POMC mRNA levels. For AgRP mRNA levels, one-way ANOVA revealed a tendency with FRW mice expressing more AgRP than FR mice (*p* = 0.07).

When we considered the evolution of mRNA expression of each group in the three protocols, we noticed first that AL mice displayed the same profile for each period and for the three investigated neuropeptides (Figure 5). Secondly, along the three protocols, FR and FRW mice did not display similar changes in the neuropeptide mRNA expressions. For AgRP, significant differences were both noted for FR and FRW groups: mRNA AgRP levels of FR mice were significantly increased in the long-term protocol compared to the short-term (tendency, *p* = 0.08) and LT-Rec (*p* = 0.017) protocols, whereas for FRW mice, it was significantly lower in the LT-Rec protocol than in the short-term (*p* = 0.004) and long-term protocols (*p* = 0.003). For NPY, only FRW mice showed significant changes with a higher mRNA expression in the short-term than in the long-term (*p* = 0.017) and LT-Rec (*p* = 0.009) protocols. Finally, for POMC, once again, the FR and FRW mice did not show a similar evolution of their mRNA expression: FR mice displayed a significant decrease in their POMC mRNA expression between the short-term and long-term protocols (*p* = 0.04) and a significant increase between the long-term and the LT-Rec (*p* = 0.02) protocols, whereas FRW mice showed a POMC mRNA level significantly higher in the LT-Rec protocol than for the short-term (*p* = 0.001) and long-term protocols (*p* = 0.03).

AgRP/NPY and POMC neurons were targeted by ghrelin and leptin. Therefore, we aimed at evaluating whether the mRNA expression of these receptors was modified in the three protocols (Figure 6). A two-way ANOVA analysis indicated a time effect for GHSR, Y1-R and MC4-R mRNA expression (GHSR, F(2,37) = 8.67, *p* = 0.0008; Y1-R, F(2,40) = 10.21, *p* = 0.0003; MC4-R, F(2,42) = 3.1, *p* = 0.055), a group effect for GHSR, Lept-R, Y1-R, Y5-R and MC3-R mRNA expression (GHSR, F(2,37) = 6.87, *p* = 0.003; Lept-R, F(2,40) = 10.89, *p* = 0.0002, Y1-R, F(2,40) = 7.73, *p* = 0.001; Y5-R, F(2,41) = 3.46, *p* = 0.04; MC3-R, F(2,42) = 12.41, *p* < 0.0001) and an interaction protocol × group for GHSR, Lept-R, Y1-R, Y5-R, MC3-R and MC4-R mRNA expression (GHSR, F(4,37) = 3.78, *p* = 0.01; Lept-R, F(4,40) = 3.89, *p* = 0.009, Y1-R, F(4,40) = 8.62, *p* < 0.0001; Y5-R, F(4,41) = 4.06, *p* = 0.007; MC3-R, F(4,42) = 5.21, *p* = 0.002; MC4-R, F(4,42) = 6.07, *p* = 0.0006).
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Impact of short-term, long-term (LT) and nutritional recovery (LT-Rec) on the mRNA expression of hypothalamic of (A) ghrelin receptor (GHSR), (B) leptin receptor, (C) Y1-receptor, (D) Y5-receptor, (E) melanocortin 3-receptor, (F) melanocortin 4-receptor. AL: *ad libitum*; FR: food restricted; FRW: food restricted with wheel. Data represent mean ± SEM (*n* = 3–8 per group). \* *p* < 0.05 AL vs. FR; \# *p* < 0.05 FR vs. FRW; \*\* *p* < 0.01 AL vs. FR.</figcaption></figure>

In the short-term protocol, GHSR, Lept-R, Y1-R and Y5-R mRNA expressions were significantly lower in FR mice compared to AL mice (p < 0.05), and no significant difference was noted between FR and FRW mice. In contrast, MC3-R mRNA expression was significantly higher only for FR mice compared to AL (p = 0.024) and FRW mice (p = 0.002) with no significant changes for MC4-R mRNA expression. As previously mentioned, the long-term protocol induced no major changes compared to short-term: GHRS, AL versus FR (p = 0.03); Lept-R, AL versus FR (p = 0.005); Y1-R, AL versus FR (p < 0.0001); MC3-R, AL versus FR (p = 0.03) and FR versus FRW (p = 0.008); whereas the mRNA levels of Y5-R and MC4-R of FR mice were significantly higher compared to FRW mice (Y5-R, p = 0.016) and to AL (MC4-R, p = 0.015) and FRW mice (MC4-R, p = 0.003). The nutritional recovery led to a return to normal values for mRNA expression of GHSR. However, FR mice showed a differential recovery compared to FRW mice for the Lept-R, Y1-R and Y5-R mRNA levels which were higher than, respectively, in FRW mice (Lept-R, p = 0.006; Y1-R, p = 0.016; Y5-R, p = 0.006) and in AL mice (Y1-R, p = 0.021). Finally, the MC3-R and MC4-R mRNA expressions were significantly higher in FR mice compared to AL mice (MC3-R, p = 0.02; MC4-R, p = 0.009).

In AL mice, except for GHSR and Y1-R mRNA expression whose levels were significantly higher in the short-term (p < 0.05) than in long-term protocol, we noticed a similar
profile of expression for each time and for the four other investigated receptors. In the three protocols, FR and FRW mice did not display similar changes in the receptor mRNA expressions. For GHS-R, significant differences were only observed for FRW mice with a significant decrease in the mRNA expression in the long-term protocol compared to short-term (p = 0.004). For Lept-R, the mRNA levels were significantly higher both for FR and FRW compared to the LT-Rec protocol between short-term (FR, p = 0.0006; FRW, p = 0.002) and long-term (FR, p < 0.001; FRW, p = 0.003). Y1-R, FR and FRW did not show the same evolutions: the mRNA levels in FR mice were higher in the LT-Rec protocol between the short-term (FR, p = 0.003) and long-term (FR, p = 0.0004) protocols, whereas for FRW, the Y1-R mRNA levels in LT-Rec returned to values obtained in the short-term protocol (ST versus LT, p = 0.002; LT versus LT-Rec, p < 0.001). Finally, no significant changes were noted for FR mice in the three protocols concerning the mRNA expression of Y5-R, MC3-R and MC4-R. It was not the same for FRW. Indeed, the Y5-R mRNA levels in the LT-Rec were similar to values obtained in the short-term protocol (ST versus LT, p = 0.0001; LT versus LT-Rec, p = 0.002), whereas the mRNA levels of MC3-R and MC4-R were higher between the short-term (MC3-R, p < 0.0001; MC4-R, p = 0.0005), long-term (MC3-R, p < 0.0001; MC4-R, p = 0.0001) and LT-Rec protocols.

## 4. Discussion

Anorexia nervosa is characterized by a high percentage of relapse after weight recovery. To prevent worsening and disastrous outcomes, we need to understand the mechanisms sustaining AN symptomatology and to decipher the factors which could contribute to restore an appropriate homeostatic balance between physical activity food intake, and body weight.

In the present study, we combined a clinical evaluation of acutely ill and partially recovered patients with restrictive-type AN using a new ecological cardiac monitoring device to assess physical activity with preclinical mouse model mimicking several aspects of AN in which we could record running wheel activity as an index of physical activity. For the first time, we showed an effect of a history of undernutrition, uncovered by the minimal BMI, on the level of physical activity in AN. Furthermore, the preclinical model supports an important role of physical activity in the recovery process, and indicates the leptinergic system as one factor that can drive a reliable restoration of metabolic variables through the hypothalamic regulation of neuropeptides involved in feeding behavior.

### 4.1. Physical Activity in Acute AN and prAN

To examine indirectly the effect of weight gain on AN symptomatology, we compared, in a treatment-seeking sample, patients with aAN and prAN. We found comparable anxiety and depressive symptoms and a more severe global symptomatology severity assessed with EDI-2 in the prAN group. Some dimensions appeared to be more severe, especially those connected to body image and bulimia. It is, thus, possible that the prevalence of binge/purging type of AN is higher in the prAN group (not controlled for), as indicated by a higher minimal BMI in this group. It is also possible that this particular group may constitute the most severe patients with prAN, since others may be less susceptible to seek treatment in a specialized center. It is known that partial recovery is a common evolution in AN [13].

To specifically address physical activity, we presented here original data of continuous cardiac monitoring permitting an ecological assessment of physical activity through the remission process. Most of the knowledge about physical activity in AN was found through declarative data or accelerometry, which may underestimate physical activity. Original and objective data on this core part of symptom of AN are still lacking [33]. The objective assessment of physical activity in daily life remained difficult until the recent development of new devices. Tools covering just a part of the physical activity (accelerometry) or the patient him- or herself (auto- and hetero- questionnaires) may underestimate physical activity. In contrast, cardiac monitoring presents the advantage of capturing all physical
activities, including movement-free activity. These new approaches will hopeful allow a translational benefit with preclinical models to better understand, assess and treat excessive physical activity in AN (Figure 7). Our initial results described a link between both a history of low body weight and a high duration of illness, and a high level of physical activity. Variations in BMI are usually considered as a major event of in the recovery process and are illustrated in daily practice by weight history. We showed here that a history of severe weight loss, and not current BMI, is associated with higher levels of current physical activity. It could be considered such as a sequelae of the undernutrition state. Weight restoration is a cornerstone of treatment [59] and a strong prognostic factor for treatment efficacy [60]. However, persistence of symptoms after weight gain, illustrated here through treatment-seeking patients with prAN, has to be further investigated. The interaction between physical activity, caloric restriction and weight loss has been widely described and explored with various rodent models mimicking several physiological, metabolic and endocrine aspects of AN [32,50,51]. We further investigated such models here.

<figure><img src="image_12.png" /><figcaption>**Figure 7.** A translational approach of recovery in AN. Weight restoration is an important step of clinical management of patients suffering from AN. However, there is a high degree of inter-individual variability in the response to weight loss correction. Thus, some patients recover from their symptoms as soon as they start weight restoration (1) and others require significant weight gain to see their symptoms improving (2). Others still have symptoms despite full weight recovery: it is partial recovery (prAN). Our initial data on patients with prAN describe a long-term effect of undernutrition state on physical activity regulation. Beyond the recovery of a normal weight, it is necessary to describe objective biomarkers of remission. To identify them at the cellular or molecular levels, mouse models mimicking the effect of prolonged weight deficit on brain tissues after weight recovery is proving very useful. We demonstrated here the persistence of alteration of the central leptinergic pathway ten days after weight recovery. AgRP: agouti related-peptide; NPY: neuropeptide Y; POMC: pro-opiomelanocortin; Y1-R/Y5-R: neuropeptide Y receptor 1/receptor 5; MC3R/MC4-R: melanocortin 3 and 4 receptors; GHSR: ghrelin receptor; Lept-R: leptin receptor.</figcaption></figure>
## 4.2. Impact of Moderate Physical Activity during Chronic Food Restriction on Recovery

After a long period of chronic food restriction, FR and FRW mice recovered in three days, a body weight similar to AL mice by adopting a binge-eating-like behavior. However, the running wheel activity performed by FRW mice induced several changes compared to FR mice such as a failure to recover perigonadic fat mass and a transient decrease in food intake associated with a slight increase in physical activity and a transient reduction in body weight. The moderate physical activity adopted during the long-term protocol might be adaptive considering the low body weight of these food restricted mice and may have a protective function to avoid rapid weight gain during nutritional recovery. Indeed, using a rat model of AN (ABA model), Giles et al. (2016) demonstrated that a rapid regain in body weight was metabolically unhealthy as it resulted in an increased lipid accumulation in the liver such as we can observe in the refeeding syndrome described in patients with AN [61].

These differences between FR and FRW might be attributed to plasma leptin concentrations, which were lower only in FRW mice at the end of the nutritional recovery compared to the long-term protocol. These decreased leptin levels were possibly the consequence of perigonadic fat mass that remained low after recovery in FRW mice, despite a comparable body weight with FR mice. Similar data were obtained with another mouse model mimicking anorexia, the “separation-induced weight loss” model (SWL). However, SWL mice, which did not have access to a running wheel such as FRW mice, recovered a fat mass similar to AL mice [50]. Furthermore, a recent analysis of the evolution of plasma leptin levels in patients with AN at inclusion, at discharge and two months later, showed that low leptin levels at discharge were associated with poor outcome [29]. For more than a decade, rodent and human studies have pointed to a causal link between excessive activity and low circulating leptin levels [17,18,62,63] leading to a recent preliminary case series using a recombinant analog of human leptin to treat with success three seriously ill patients with AN [30]. Thus, leptin might be considered as a potential reliable indicator of recovery linked to a physical activity assessment. The evolution of plasma concentrations of ghrelin was not different between FR and FRW mice in our conditions. As a metabolic sensor, AG and DAG were more sensible to the nutritional state than to the presence of physical activity [4,28].

## 4.3. Nutritional Recovery Differentially Impacts Hypothalamic Activity

When we considered the primary hypothalamic targets of ghrelin and leptin, namely, the AgRP/NPY and POMC neurons, we also uncovered noticeable changes. FRW mice displayed after two weeks of food restriction a distinguishable phenotype from FR mice. They displayed an accumulation of AgRP in neuronal cell bodies as also described in the genetic model of AN, the anx/anx mice [64,65], but without any changes in the fibers density at the level of PVN, one of the main target of AgRP arcuate neurons. This could be related to the higher level of AgRP mRNA in the hypothalamus after the short-term protocol described here, but also in another study using rats submitted to a short-term food restriction associated with running wheels (ABA model; [66]). A recent elegant study pointed out that selective arcuate NPY/AgRP activation, using DREADD technology, significantly increased not only the food intake, but also acute moderate-intensity exercise [37]. Furthermore, considering MC3/MC4 receptors which are naturally antagonized by AgRP [67], their mRNA expression was strongly diminished during the food restriction period, only in FRW mice. These data may reflect the impact of physical activity rather than the nutritional state. The NPY mRNA expression followed the evolution of AgRP, with an adapted and stable increase for FRW mice, whereas in FR mice this augmentation was enhanced during the duration of the protocol. The orexigenic effect of NPY is relayed through Y1-R and Y5-R in the PVN and lateral hypothalamic area [68,69]. These receptors are also abundantly present on POMC neurons [70], where Y1-R stimulation leads to an inhibition of electrical activity [71] and Y5R activation inhibits α-MSH release [72]. Their mRNA expressions were decreased except for Y5-R, for which expression increased only in FR mice during the long-term food restriction, suggesting once again a differential impact
of physical activity on the regulation mechanisms of feeding. The down-regulation of Y1-R expression might be considered as a consequence of the upregulation of NPY and the associated feedback control. The mRNA expression of POMC mirrored the relative gene expression of orexigenic peptides, with a decrease which appeared earlier for FRW mice than for FR mice. This adapted decrease was also described in the ABA model [66].

Finally, the chronic food restriction was also accompanied by a reduced mRNA expression of GHSR and leptin receptor. To our knowledge, our study is the first to measure hypothalamic GHSR expression in an AN-like rodent model. Despite an initial increase in its expression in the FRW mice after two weeks of protocol, probably due to a feedback control induced by high circulating AG levels, GHSR mRNA expression decreased when the protocol was prolonged. We can hypothesize that a continuous release of AG and DAG might establish a desensitization of the hypothalamic GHSR. Concerning the leptin receptor, a similar decrease occurred in the relative gene expression of the leptin receptor in the SWL model during a long-term food restriction protocol [73]. This decrease might be paralleled with the total breakdown of fat masses.

The two-week of ad libitum refeeding led to a recovery of the mRNA expression of AgRP, NPY and POMC. However, this was not the case for the receptors. The upregulation of MC3/MC4-R mRNA expression in both FR and FRW mice may have reflected the changes in energy status since these two groups of mice gained more than 50% of their body weight during the duration of the nutritional recovery period. The mRNA expression Y1- and Y5-R of FRW mice returned to AL basal values, which was not the case for FR mice. Y1 and Y5 receptors are known to mediate the hyperphagic effect of NPY [74]. Indeed, the simultaneous ablation of both receptors in the germline Y1Y5 receptor double knockout (Y1Y5-/-) mice led to reductions in spontaneous and/or fasting-induced food intake. This indicates that food intake in mice requires the coordinated action of both the Y1 and the Y5 receptors [75]. Thus, here, we hypothesized that this increase observed in FR compared to FRW mice might have been a consequence of the hyperphagic behavior observed after refeeding, especially for FR mice. In parallel, the leptin receptor mRNA expression remained low only in FRW mice, whereas GHSR mRNA expression was stabilized. Once again, food restriction associated with moderate physical activity, as it was the case for FRW mice, induced a differential recovery compared to sole food restriction. The impact of physical activity was particularly salient on the leptin system and on related pathways. To link our preclinical data with the clinical reality, further studies are now needed to investigate, for example, the behavior of these mice faced to a food challenge, such as high caloric food.

### 4.4. Limitations

Our study of partial remission patients was based on treatment-seeking patients only. This sampling bias has to be considered as a limitation to address the effect of renutrition on AN symptom. These individuals may be a subgroup of partial recovery patients. The included patients may be among the most severe patients with partial recovery, as they seek help even though having normal weight. A longitudinal design could bring more accurate data to reinforce the effects observed in the present study.

Nevertheless, our sample remains informative. We presented for the first-time initial data obtained from continuous cardiac monitoring. To better distinguish the activity profiles of our AN patient cohorts, we need to develop precise algorithms.

Concerning the mouse model, the main limitation was the difficulty to gather all the data obtained in each protocol, because the brain analyses of peptides and receptors need to use three different groups of mice. There is currently no consensus about the nutritional recovery in mouse model mimicking the metabolic aspects of AN. To our knowledge, the analysis on peripheral and central aspects of nutritional recovery after a long-term protocol of food restriction associated or not with physical activity is lacking. Finally, we deliberately did not include another control group consisting of mice with ad libitum food access and a running wheel (ALW). Indeed, our previous works have underlined that AL
and ALW mice displayed similar behavior in anxiety tests and comparable changes in endocrine and metabolic parameters evolution during the protocols [32,49,52].

**5. Conclusions**

Taken together, our translational approach allowed us to study the role of physical activity in the modulation of feeding behavior in the context of nutritional recovery in AN. The present clinical data bring a new ecological measure of cardiac activity that reflects the physical profile of patients with AN. In combination with the undernutrition history and likely biomarkers, it will allow us to better characterize and treat aAN and prAN patients to reach a full recovery (Figure 7).

The preclinical data point out several potential central mechanisms, which can be key elements in the recovery processes. Further studies using molecular tools are now necessary to better apprehend the role of the leptin system in chronic food restriction and recovery. Our data reflecting the central role of leptin and related pathways in nutritional recovery are supported by recent papers on AN, since systemic leptin levels after weight restoration predict the short-term [29] and long-term outcomes of patients with AN [76]. The results obtained here allow us, considering a potential combination of biomarkers, to take into account in order to validate a definition of remission: the physical activity profile according to an undernutrition history, biological assays and functional metabolic neuroimaging. Indeed, on hypothalamic biosensors of nutritional status in mice might be paralleled with the structural and functional changes observed in the hypothalamus of patients with AN using proton magnetic resonance spectroscopy and diffusion tensor imaging [48]. The neurobiological and neuroendocrine data shed new light on the mechanisms of nutritional recovery and open new perspectives to prevent relapse in AN.

**Author Contributions:** Conceptualization: O.V., P.D., C.C., P.G. and V.T.; Methodology: O.V., P.D., I.A.K.N., D.A., O.L.T., N.C., C.R. and V.T.; Formal analysis: O.V., P.D. and V.T.; Writing—original draft preparation, O.V. and P.D.; Writing—review and editing: all authors; Project administration: O.V., C.C., P.G., V.T. and P.D.; Funding acquisition: P.D., P.G., O.V., C.C. and V.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This clinical study was funded by Fondation Sandrine Castellotti, and was also under the aegis of Fondation de l'Avenir pour la Recherche Médicale Appliquée, whose aims are to help the prevention of eating disorders, promote medical research and improve the care sector. This preclinical study was supported by grants from the Nord-Pas de Calais Région "ARCIR-GIRAM" (France) and from INSERM (France).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki. Ethical review and approval were waived for this study due to the fact that the data used in this study were being systematically collected for clinical assessments. All patients signed an informed consent stating that they agreed with the usage of their clinical data for the purpose of research (i.e., above clinical utility for themselves).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data that support the findings of this study are not available.

**Acknowledgments:** We acknowledge the animal facility staff of the UMR-S-1172 (JPARC, Lille) for their help and technical assistance. We warmly thank Aurore Pailler, Marion Deloulay, Maria-Manel Okka and Etienne Kimmel for managing devices of continuous cardiac monitoring. We also kindly thank Mathieu Anderhalt from Données Brutes® for the conception of the Python algorithm and Chloé Tezenas du Montcel for her statistical advices.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Duncan, L.; Yilmaz, Z.; Gaspar, H.; Walters, R.; Goldstein, J.; Anttila, V.; Bulik-Sullivan, B.; Ripke, S.; Eating Disorders Working Group of the Psychiatric Genomics Consortium; Thornton, L.; et al. Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa. *Am. J. Psychiatry* **2017**, *174*, 850–858. [CrossRef] [PubMed]
2. Schorr, M.; Miller, K.K. The endocrine manifestations of anorexia nervosa: Mechanisms and management. *Nat. Rev. Endocrinol.* **2017**, *13*, 174–186. [CrossRef] [PubMed]

3. Watson, H.J.; Yilmaz, Z.; Thornton, L.M.; Hübel, C.; Coleman, J.R.I.; Gaspar, H.A.; Bryois, J.; Hinney, A.; Leppä, V.M.; Mattheisen, M.; et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. *Nat. Genet.* **2019**, *51*, 1207–1214. [CrossRef] [PubMed]

4. Duriez, P.; Ramoz, N.; Gorwood, P.; Viltart, O.; Tolle, V. A Metabolic Perspective on Reward Abnormalities in Anorexia Nervosa. *Trends Endocrinol. Metab.* **2019**, *30*, 915–928. [CrossRef] [PubMed]

5. Viltart, O.; Duriez, P.; Tolle, V. Metabolic and neuroendocrine adaptations to undernutrition in anorexia nervosa: From a clinical to a basic research point of view. *Horm. Mol. Biol. Clin. Investig.* **2018**, *36*. [CrossRef]

6. Holtkamp, K.; Hebebrand, J.; Herpertz-Dahlmann, B. The contribution of anxiety and food restriction on physical activity levels in acute anorexia nervosa. *Int. J. Eat. Disord.* **2004**, *36*, 163–171. [CrossRef]

7. Duclos, M.; Ouerdani, A.; Mormede, P.; Konsman, J.P. Food restriction-induced hyperactivity: Addiction or adaptation to famine? *Psychoneuroendocrinology* **2013**, *38*, 884–897. [CrossRef]

8. Pleplé, A.; Lalanne, C.; Huas, C.; Mattar, L.; Hanachi, M.; Flament, M.F.; Carchon, I.; Jouen, F.; Berthoz, S.; Godart, N. Nutritional status and anxious and depressive symptoms in anorexia nervosa: A prospective study. *Sci. Rep.* **2021**, *11*, 771. [CrossRef]

9. Van Furth, E.F.; van der Meer, A.; Cowan, K. Top 10 research priorities for eating disorders. *Lancet Psychiatry* **2016**, *3*, 706–707. [CrossRef]

10. Strober, M.; Freeman, R.; Morrell, W. The Long-Term Course of Severe Anorexia Nervosa in Adolescents: Survival Analysis of Recovery, Relapse, and Outcome Predictors over 10–15 Years in a Prospective Study. *Int. J. Eat. Disord.* **1997**, *22*, 339–360. [CrossRef]

11. Treasure, J.; Zipfel, S.; Micali, N.; Wade, T.; Stice, E.; Claudino, A.; Schmidt, U.; Frank, G.K.; Bulik, C.M.; Wentz, E. Anorexia nervosa. *Nat. Rev. Dis. Primers* **2015**, *1*, 1–21. [CrossRef] [PubMed]

12. Muratore, A.F.; Attia, E. Current Therapeutic Approaches to Anorexia Nervosa: State of the Art. *Clin. Ther.* **2021**, *43*, 85–94. [CrossRef] [PubMed]

13. Eddy, K.T.; Dorer, D.J.; Franko, D.L.; Tahilani, K.; Thompson-Brenner, H.; Herzog, D.B. Diagnostic Crossover in Anorexia Nervosa and Bulimia Nervosa: Implications for DSM-V. *Am. J. Psychiatry* **2008**, *165*, 245–250. [CrossRef] [PubMed]

14. Kawakami, Y.; Arai, T.; Okuda, Y.; Yamashita, K. The Regulation of obese(ob) Gene Expression by Intracellular Fatty Acid Concentration in Adipocytes. *Intern. Med.* **1999**, *38*, 213–215. [CrossRef] [PubMed]

15. Fried, S.K.; Ricci, M.R.; Russell, C.D.; Laferrère, B. Regulation of Leptin Production in Humans. *J. Nutr.* **2000**, *130*, 3127S–3131S. [CrossRef]

16. Peelman, F.; Zabeau, L.; Moharana, K.; Savvides, S.N.; Tavernier, J. 20 YEARS OF LEPTIN: Insights into signaling assemblies of the leptin receptor. *J. Endocrinol.* **2014**, *223*, T9–T23. [CrossRef]

17. Fernandes, M.F.A.; Matthys, D.; Hryhorczuk, C.; Sharma, S.; Mogra, S.; Alquier, T.; Fulton, S. Leptin Suppresses the Rewarding Effects of Running via STAT3 Signaling in Dopamine Neurons. *Cell Metab.* **2015**, *22*, 741–749. [CrossRef]

18. Exner, C.; Hebebrand, J.; Remschmidt, H.; Wewetzer, C.; Ziegler, A.; Herpertz, S.; Schweiger, U.; Blum, W.F.; Preibisch, G.; Heldmaier, G.; et al. Leptin suppresses semi-starvation induced hyperactivity in rats: Implications for anorexia nervosa. *Mol. Psychiatry* **2000**, *5*, 476–481. [CrossRef] [PubMed]

19. Kojima, M.; Hosoda, H.; Date, Y.; Nakazato, M.; Matsuo, H.; Kangawa, K. Ghrelin is a growth-hormone-releasing acylated peptide from stomach. *Nature* **1999**, *402*, 656–660. [CrossRef]

20. Date, Y.; Kojima, M.; Hosoda, H.; Sawaguchi, A.; Mondal, M.S.; Suganuma, T.; Matsukura, S.; Kangawa, K.; Nakazato, M. Ghrelin, a Novel Growth Hormone-Releasing Acylated Peptide, Is Synthesized in a Distinct Endocrine Cell Type in the Gastrointestinal Tracts of Rats and Humans. *Endocrinology* **2000**, *141*, 4255–4261. [CrossRef]

21. Müller, T.; Nogueiras, R.; Andermann, M.; Andrews, Z.B.; Anker, S.; Argente, J.; Batterham, R.; Benoit, S.; Bowers, C.; Broglio, F.; et al. Ghrelin. *Mol. Metab.* **2015**, *4*, 437–460. [CrossRef]

22. Méquinion, M.; Langlet, F.; Zgheib, S.; Dickson, S.; Dehouck, B.; Chauveau, C.; Viltart, O. Ghrelin: Central and Peripheral Implications in Anorexia Nervosa. *Front. Endocrinol.* **2013**, *4*, 15. [CrossRef]

23. Labarthe, A.; Fiquet, O.; Hassouna, R.; Zizzari, P.; Lanfumey, L.; Ramoz, N.; Grouselle, D.; Epelbaum, J.; Tolle, V. Ghrelin-Derived Peptides: A Link between Appetite/Reward, GH Axis, and Psychiatric Disorders? *Front. Endocrinol.* **2014**, *5*, 163. [CrossRef]

24. Ralevski, E.; Shanabrough, M.; Newcomb, J.; Gandelman, E.; Hayden, R.; Horvath, T.L.; Petrakis, I. Ghrelin is Related to Personality Differences in Reward Sensitivity and Impulsivity. *Alcohol Alcohol.* **2018**, *53*, 52–56. [CrossRef]

25. Bake, T.; Edvardsson, C.E.; Cummings, C.J.; Dickson, S.L. Ghrelin's effects on food motivation in rats are not limited to palatable foods. *J. Neuroendocrinol.* **2019**, *31*, e12665. [CrossRef]

26. Hosoda, H.; Kojima, M.; Matsuo, H.; Kangawa, K. Ghrelin and Des-acyl Ghrelin: Two Major Forms of Rat Ghrelin Peptide in Gastrointestinal Tissue. *Biochem. Biophys. Res. Commun.* **2000**, *279*, 909–913. [CrossRef] [PubMed]

27. Asakawa, A.; Inui, A.; Fujimiya, M.; Sakamaki, R.; Shinfuku, N.; Ueta, Y.; Meguid, M.M.; Kasuga, M. Stomach regulates energy balance via acylated ghrelin and desacyl ghrelin. *Gut* **2005**, *54*, 18–24. [CrossRef] [PubMed]

28. Gortan Cappellari, G.; Barazzoni, R. Ghrelin forms in the modulation of energy balance and metabolism. *Eat. Weight. Disord.* **2019**, *24*, 997–1013. [CrossRef] [PubMed]
29. Dardennes, R.; Tolle, V.; Lavoisy, G.; Grouselle, D.; Alanbar, N.; Duriez, P.; Gorwood, P.; Ramoz, N.; Epelbaum, J. Lower leptin level at discharge in acute anorexia nervosa is associated with early weight-loss. *Eur. Eat. Disord. Rev.* **2021**, *29*, 634–644. [CrossRef]

30. Milos, G.; Antel, J.; Kaufmann, L.-K.; Barth, N.; Koller, A.; Tan, S.; Wiesing, U.; Hinney, A.; Libuda, L.; Wabitsch, M.; et al. Short-term metreleptin treatment of patients with anorexia nervosa: Rapid on-set of beneficial cognitive, emotional, and behavioral effects. *Transl. Psychiatry* **2020**, *10*, 303. [CrossRef]

31. Seidel, M.; Markmann Jensen, S.; Healy, D.; Dureja, A.; Watson, H.J.; Holst, B.; Bulik, C.M.; Sjögren, J.M. A Systematic Review and Meta-Analysis Finds Increased Blood Levels of All Forms of Ghrelin in Both Restricting and Binge-Eating/Purging Subtypes of Anorexia Nervosa. *Nutrients* **2021**, *13*, 709. [CrossRef]

32. Duriez, P.; Robichon, L.; Dardennes, R.; Lavoisy, G.; Grouselle, D.; Epelbaum, J.; Ramoz, N.; Gorwood, P.; Tolle, V.; Viltart, O. Unexpected Association of Desacyl-Ghrelin with Physical Activity and Chronic Food Restriction: A Translational Study on Anorexia Nervosa. *J. Clin. Med.* **2020**, *9*, 2782. [CrossRef]

33. Meyer, C.; Taranis, L.; Touyz, S. Excessive exercise in the eating disorders: A need for less activity from patients and more from researchers. *Eur. Eat. Disord. Rev.* **2008**, *16*, 81–83. [CrossRef]

34. Krashes, M.J.; Koda, S.; Ye, C.; Rogan, S.C.; Adams, A.C.; Cusher, D.S.; Maratos-Flier, E.; Roth, B.L.; Lowell, B.B. Rapid, reversible activation of AgRP neurons drives feeding behavior in mice. *J. Clin. Investig.* **2011**, *121*, 1424–1428. [CrossRef]

35. Betley, J.N.; Cao, Z.F.H.; Ritola, K.D.; Sternson, S.M. Parallel, Redundant Circuit Organization for Homeostatic Control of Feeding Behavior. *Cell* **2013**, *155*, 1337–1350. [CrossRef]

36. Zhan, C.; Zhou, J.; Feng, Q.; Zhang, J.-E.; Lin, S.; Bao, J.; Wu, P.; Luo, M. Acute and Long-Term Suppression of Feeding Behavior by POMC Neurons in the Brainstem and Hypothalamus, Respectively. *J. Neurosci.* **2013**, *33*, 3624–3632. [CrossRef]

37. Bunner, W.; Landry, T.; Laing, B.T.; Li, P.; Rao, Z.; Yuan, Y.; Huang, H. ARCAgRP/NPY Neuron Activity Is Required for Acute Exercise-Induced Food Intake in Un-Trained Mice. *Front. Physiol.* **2020**, *11*, 411. [CrossRef]

38. Williams, K.W.; Elmquist, J.K. From neuroanatomy to behavior: Central integration of peripheral signals regulating feeding behavior. *Nat. Neurosci.* **2012**, *15*, 1350–1355. [CrossRef]

39. Dietrich, M.O.; Bober, J.; Ferreira, J.G.; Tellez, L.A.; Mineur, Y.S.; Souza, D.O.; Gao, X.-B.; Picciotto, M.R.; Araújo, I.; Liu, Z.-W.; et al. AgRP neurons regulate development of dopamine neuronal plasticity and nonfood-associated behaviors. *Nat. Neurosci.* **2012**, *15*, 1108–1110. [CrossRef]

40. Dietrich, M.O.; Zimmer, M.R.; Bober, J.; Horvath, T.L. Hypothalamic Agrp Neurons Drive Stereotypic Behaviors beyond Feeding. *Cell* **2015**, *160*, 1222–1232. [CrossRef]

41. Denis, R.G.P.; Joly-Amado, A.; Webber, E.; Langlet, F.; Schaeffer, M.; Padilla, S.L.; Cansell, C.; Dehouck, B.; Castel, J.; Delbès, A.-S.; et al. Palatability Can Drive Feeding Independent of AgRP Neurons. *Cell Metab.* **2017**, *25*, 975. [CrossRef]

42. Miletta, M.C.; Iyilikci, O.; Shanabrough, M.; Šestan-Peša, M.; Cammisa, A.; Zeiss, C.J.; Dietrich, M.O.; Horvath, T.L. AgRP neurons control compulsive exercise and survival in an activity-based anorexia model. *Nat. Metab.* **2020**, *2*, 1204–1211. [CrossRef] [PubMed]

43. He, Z.; Gao, Y.; Alhadeff, A.L.; Castorena, C.M.; Huang, Y.; Lieu, L.; Afrin, S.; Sun, J.; Betley, J.N.; Guo, H.; et al. Cellular and synaptic reorganization of arcuate NPY/AgRP and POMC neurons after exercise. *Mol. Metab.* **2018**, *18*, 107–119. [CrossRef] [PubMed]

44. Li, H.E.; Rossi, M.A.; Watson, G.D.R.; Moore, H.G.; Cai, M.T.; Kim, N.; Vokt, K.A.; Lu, D.; Bartholomew, R.A.; Hughes, R.N.; et al. Hypothalamic-Extended Amygdala Circuit Regulates Temporal Discounting. *J. Neurosci.* **2021**, *41*, 1928–1940. [CrossRef] [PubMed]

45. Li, C.; Hou, Y.; Zhang, J.; Sui, G.; Du, X.; Licinio, J.; Wong, M.-L.; Yang, Y. AGRP neurons modulate fasting-induced anxiolytic effects. *Transl. Psychiatry* **2019**, *9*, 111. [CrossRef]

46. Schaeffer, M.; Langlet, F.; Lafont, C.; Molino, F.; Hodson, D.J.; Roux, T.; Lamarque, L.; Verdié, P.; Bourrier, E.; Dehouck, B.; et al. Rapid sensing of circulating ghrelin by hypothalamic appetite-modifying neurons. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 1512–1517. [CrossRef]

47. Wang, Q.; Liu, C.; Uchida, A.; Chuang, J.-C.; Walker, A.; Liu, T.; Osborne-Lawrence, S.; Mason, B.L.; Mosher, C.; Berglund, E.D.; et al. Arcuate AgRP neurons mediate orexigenic and glucoregulatory actions of ghrelin. *Mol. Metab.* **2014**, *3*, 64–72. [CrossRef]

48. Florent, V.; Baroncini, M.; Jissendi-Tchofo, P.; Lopes, R.; Vanhoutte, M.; Rasika, S.; Pruvo, J.-P.; Vignau, J.; Verdun, S.; Johansen, J.E.; et al. Hypothalamic Structural and Functional Imbalances in Anorexia Nervosa. *Neuroendocrinology* **2020**, *110*, 552–562. [CrossRef]

49. Méquinion, M.; Caron, E.; Zgheib, S.; Stievenard, A.; Zizzari, P.; Tolle, V.; Cortet, B.; Lucas, S.; Prevot, V.; Chauveau, C.; et al. Physical activity: Benefit or weakness in metabolic adaptations in a mouse model of chronic food restriction? *Am. J. Physiol. Endocrinol. Metab.* **2015**, *308*, E241–E255. [CrossRef] [PubMed]

50. Méquinion, M.; Chauveau, C.; Viltart, O. The Use of Animal Models to Decipher Physiological and Neurobiological Alterations of Anorexia Nervosa Patients. *Front. Endocrinol.* **2015**, *6*, 68. [CrossRef]

51. Scharner, S.; Stengel, A. Animal Models for Anorexia Nervosa—A Systematic Review. *Front. Hum. Neurosci.* **2020**, *14*, 596381. [CrossRef] [PubMed]

52. Duriez, P.; Eddarkaoui, S.; Blum, D.; Dickson, S.L.; Gorwood, P.; Tolle, V.; Viltart, O. Does physical activity associated with chronic food restriction alleviate anxiety like behaviour, in female mice? *Horm. Behav.* **2020**, *124*, 104807. [CrossRef]
53. Verhagen, L.A.W.; Egecioglu, E.; Luijendijk, M.C.M.; Hillebrand, J.J.G.; Adan, R.A.H.; Dickson, S.L. Acute and chronic suppression of the central ghrelin signaling system reveals a role in food anticipatory activity. *Eur. Neuropsychopharmacol.* **2011**, *21*, 384–392. [CrossRef] [PubMed]

54. American Psychiatric Association. *Diagnostic and Statistical Manual of Mental Disorders*, 5th ed.; American Psychiatric Association: Arlington, TX, USA, 2013; ISBN 978-0-89042-555-8.

55. Garner, D.M. *Eating Disorder Inventory—2nd Manual*; Psychological Assessment Resources: Odessa, FL, USA, 1991.

56. Zigmond, A.S.; Snaith, R.P. The Hospital Anxiety and Depression Scale. *Acta Psychiatr. Scand.* **1983**, *67*, 361–370. [CrossRef] [PubMed]

57. Goodman, W.K.; Price, L.H.; Rasmussen, S.A.; Mazure, C.; Fleischmann, R.L.; Hill, C.L.; Heninger, G.R.; Charney, D.S. The Yale-Brown Obsessive Compulsive Scale. I. Development, Use, and Reliability. *Arch. Gen. Psychiatry* **1989**, *46*, 1006–1011. [CrossRef] [PubMed]

58. Chomczynski, P.; Sacchi, N. Single-step method of RNA isolation by acid guanidinium thiocyanate–phenol–chloroform extraction. *Anal. Biochem.* **1987**, *162*, 156–159. [CrossRef]

59. Ehrlich, S.; Tam, F.I. Editorial: Refeeding in Anorexia Nervosa: Quo Vadis? *J. Am. Acad. Child. Adolesc. Psychiatry* **2021**, *60*, 566–567. [CrossRef]

60. Redgrave, G.W.; Schreyer, C.C.; Coughlin, J.W.; Fischer, L.K.; Pletch, A.; Guarda, A.S. Discharge Body Mass Index, Not Illness Chronicity, Predicts 6-Month Weight Outcome in Patients Hospitalized With Anorexia Nervosa. *Front. Psychiatry* **2021**, *12*, 641861. [CrossRef]

61. Giles, E.D.; Hagman, J.; Pan, Z.; MacLean, P.S.; Higgins, J.A. Weight restoration on a high carbohydrate refeeding diet promotes rapid weight regain and hepatic lipid accumulation in female anorexic rats. *Nutr. Metab.* **2016**, *13*, 18. [CrossRef]

62. Holtkamp, K.; Herpertz-Dahlmann, B.; Hebebrand, K.; Mika, C.; Kratzsch, J.; Hebebrand, J. Physical Activity and Restlessness Correlate with Leptin Levels in Patients with Adolescent Anorexia Nervosa. *Biol. Psychiatry* **2006**, *60*, 311–313. [CrossRef]

63. Van Elburg, A.A.; Kas, M.J.H.; Hillebrand, J.J.G.; Eijkemans, R.J.C.; Van Engeland, H. The impact of hyperactivity and leptin on recovery from anorexia nervosa. *J. Neural Transm.* **2007**, *114*, 1233–1237. [CrossRef]

64. Broberger, C.; Johansen, J.; Johansson, C.; Schalling, M.; Hokfelt, T. The neuropeptide Y/agouti gene-related protein (AGRP) brain circuitry in normal, anorectic, and monosodium glutamate-treated mice. *Proc. Natl. Acad. Sci. USA* **1998**, *95*, 15043–15048. [CrossRef]

65. Nilsson, I.; Lindfors, C.; Fetissov, S.O.; Hökfelt, T.; Johansen, J.E. Aberrant agouti-related protein system in the hypothalamus of theanx/anx mouse is associated with activation of microglia. *J. Comp. Neurol.* **2008**, *507*, 1128–1140. [CrossRef] [PubMed]

66. De Rijke, C.E.; Hillebrand, J.J.G.; Verhagen, L.A.W.; Roeling, T.A.P.; Adan, R.A.H. Hypothalamic neuropeptide expression following chronic food restriction in sedentary and wheel-running rats. *J. Mol. Endocrinol.* **2005**, *35*, 381–390. [CrossRef]

67. Ollmann, M.M.; Wilson, B.D.; Yang, Y.K.; Kerns, J.A.; Chen, Y.; Gantz, I.; Barsh, G.S. Antagonism of Central Melanocortin Receptors in Vitro and in Vivo by Agouti-Related Protein. *Science* **1997**, *278*, 135–138. [CrossRef] [PubMed]

68. Kopp, J.; Xu, Z.-Q.; Zhang, X.; Pedrazzini, T.; Herzog, H.; Kresse, A.; Wong, H.; Walsh, J.H.; Hökfelt, T. Expression of the neuropeptide Y Y1 receptor in the CNS of rat and of wild-type and Y1 receptor knock-out mice. Focus on immunohistochemical localization. *Neuroscience* **2002**, *111*, 443–532. [CrossRef]

69. Morin, S.M.; Gehlert, D.R. Distribution of NPY Y5-Like Immunoreactivity in the Rat Brain. *J. Mol. Neurosci.* **2006**, *29*, 109–114. [CrossRef]

70. Mercer, R.E.; Chee, M.J.S.; Colmers, W.F. The role of NPY in hypothalamic mediated food intake. *Front. Neuroendocr.* **2011**, *32*, 398–415. [CrossRef]

71. Acuna-Goycolea, C.; Tamamaki, N.; Yanagawa, Y.; Obata, K.; van den Pol, A.N. Mechanisms of Neuropeptide Y, Peptide YY, and Pancreatic Polypeptide Inhibition of Identified Green Fluorescent Protein-Expressing GABA Neurons in the Hypothalamic Neuroendocrine Arcuate Nucleus. *J. Neurosci.* **2005**, *25*, 7406–7419. [CrossRef]

72. Galas, L.; Tonon, M.-C.; Beaujean, D.; Fredriksson, R.; Larhammar, D.; Lihrmann, I.; Jégou, S.; Fournier, A.; Chartrel, N.; Vaudry, H. Neuropeptide Y Inhibits Spontaneous Alpha-Melanocyte-Stimulating Hormone (Alpha-MSH) Release via a Y(5) Receptor and Suppresses Thyrotropin-Releasing Hormone-Induced Alpha-MSH Secretion via a Y(1) Receptor in Frog Melanotrope Cells. *Endocrinology* **2002**, *143*, 1686–1694. [CrossRef]

73. Méquinion, M.; Le Thuc, O.; Zgheib, S.; Alexandre, D.; Chartrel, N.; Rovère, C.; Hardouin, P.; Viltart, O.; Chauveau, C. Long-Term Energy Deficit in Mice Causes Long-Lasting Hypothalamic Alterations after Recovery. *Neuroendocrinology* **2017**, *105*, 372–383. [CrossRef]

74. Iyengar, S.; Li, D.L.; Simmons, R.M. Characterization of neuropeptide Y-induced feeding in mice: Do Y1-Y6 receptor subtypes mediate feeding? *J. Pharmacol. Exp. Ther.* **1999**, *289*, 1031–1040. [PubMed]

75. Nguyen, A.D.; Mitchell, N.F.; Lin, S.; Macia, L.; Yulyaningsih, E.; Baldock, P.A.; Enriquez, R.F.; Zhang, L.; Shi, Y.-C.; Zolotukhin, S.; et al. Y1 and Y5 Receptors Are Both Required for the Regulation of Food Intake and Energy Homeostasis in Mice. *PLoS ONE* **2012**, *7*, e40191. [CrossRef] [PubMed]

76. Kim, Y.; Hersch, J.; Bodell, L.P.; Schebendach, J.; Hildebrandt, T.; Walsh, B.T.; Mayer, L.E.S. The association between leptin and weight maintenance outcome in anorexia nervosa. *Int. J. Eat. Disord.* **2021**, *54*, 527–534. [CrossRef] [PubMed]