# The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice

Yong Fan, René Klinkby Støving, Samar Berreira Ibraim, Tuulia Hyötyläinen, Florence Thirion, Tulika Arora, Liwei Lyu, Evelina Stankevic, Tue Haldor Hansen, Pierre Déchelotte, et al.

► To cite this version:

Yong Fan, René Klinkby Støving, Samar Berreira Ibraim, Tuulia Hyötyläinen, Florence Thirion, et al. The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice. Nature Microbiology, 2023, 8 (5), pp.787-802. 10.1038/s41564-023-01355-5. hal-04103775

HAL Id: hal-04103775

https://normandie-univ.hal.science/hal-04103775

Submitted on 29 May 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice

Yong Fan, René Klinkby Støving, Samar Berreira Ibraim, Tuulia Hyötyläinen, Florence Thirion, Tulika Arora, Liwei Lyu, Evelina Stankevic, Tue Haldor Hansen, Pierre Déchelotte, et al.

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

## To cite this version:

Yong Fan, René Klinkby Støving, Samar Berreira Ibraim, Tuulia Hyötyläinen, Florence Thirion, et al. The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice. Nature Microbiology, 2023, 8 (5), pp.787-802. 10.1038/s41564-023-01355-5. hal-04103775

HAL Id: hal-04103775

https://normandie-univ.hal.science/hal-04103775v1

Submitted on 29 May 2024

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

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

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

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

# The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice

Received: 23 June 2022

Accepted: 3 March 2023

Published online: 17 April 2023

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

Check for updates

A list of authors and their affiliations appears at the end of the paper

Anorexia nervosa (AN) is an eating disorder with a high mortality. About 95% of cases are women and it has a population prevalence of about 1%, but evidence-based treatment is lacking. The pathogenesis of AN probably involves genetics and various environmental factors, and an altered gut microbiota has been observed in individuals with AN using amplicon sequencing and relatively small cohorts. Here we investigated whether a disrupted gut microbiota contributes to AN pathogenesis. Shotgun metagenomics and metabolomics were performed on faecal and serum samples, respectively, from a cohort of 77 females with AN and 70 healthy females. Multiple bacterial taxa (for example, *Clostridium* species) were altered in AN and correlated with estimates of eating behaviour and mental health. The gut virome was also altered in AN including a reduction in viral-bacterial interactions. Bacterial functional modules associated with the degradation of neurotransmitters were enriched in AN and various structural variants in bacteria were linked to metabolic features of AN. Serum metabolomics revealed an increase in metabolites associated with reduced food intake (for example, indole-3-propionic acid). Causal inference analyses implied that serum bacterial metabolites are potentially mediating the impact of an altered gut microbiota on AN behaviour. Further, we performed faecal microbiota transplantation from AN cases to germ-free mice under energy-restricted feeding to mirror AN eating behaviour. We found that the reduced weight gain and induced hypothalamic and adipose tissue gene expression were related to aberrant energy metabolism and eating behaviour. Our 'omics' and mechanistic studies imply that a disruptive gut microbiome may contribute to AN pathogenesis.

Anorexia nervosa (AN) is a serious mental health condition and eating disorder characterized by distorted body image, obsessive thoughts about food, ritualistic patterns of behaviour including reduced food intake, loss of body weight, raised physical activity and emotional rigidity¹. AN primarily affects women in about 95% of cases and has a population prevalence of about 1%². It can be classified into two subtypes, the

common restricting (AN-RS) type and the less prevalent binge-eating or purging (AN-BP) type¹. The evidence base for treatment is lacking³ and although specialized multidisciplinary treatment can reduce mortality⁴, less than half of AN cases achieve complete remission⁵. The aggregate mortality rate is estimated to be 5.6% per decade, much higher than in the general population⁶.

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

e-mail: oluf@sund.ku.dk
<figure><img src="image_3.png" /></figure>

# The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice

Received: 23 June 2022

Accepted: 3 March 2023

Published online: 17 April 2023

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

Check for updates

A list of authors and their affiliations appears at the end of the paper

Anorexia nervosa (AN) is an eating disorder with a high mortality. About 95% of cases are women and it has a population prevalence of about 1%, but evidence-based treatment is lacking. The pathogenesis of AN probably involves genetics and various environmental factors, and an altered gut microbiota has been observed in individuals with AN using amplicon sequencing and relatively small cohorts. Here we investigated whether a disrupted gut microbiota contributes to AN pathogenesis. Shotgun metagenomics and metabolomics were performed on faecal and serum samples, respectively, from a cohort of 77 females with AN and 70 healthy females. Multiple bacterial taxa (for example, *Clostridium* species) were altered in AN and correlated with estimates of eating behaviour and mental health. The gut virome was also altered in AN including a reduction in viral-bacterial interactions. Bacterial functional modules associated with the degradation of neurotransmitters were enriched in AN and various structural variants in bacteria were linked to metabolic features of AN. Serum metabolomics revealed an increase in metabolites associated with reduced food intake (for example, indole-3-propionic acid). Causal inference analyses implied that serum bacterial metabolites are potentially mediating the impact of an altered gut microbiota on AN behaviour. Further, we performed faecal microbiota transplantation from AN cases to germ-free mice under energy-restricted feeding to mirror AN eating behaviour. We found that the reduced weight gain and induced hypothalamic and adipose tissue gene expression were related to aberrant energy metabolism and eating behaviour. Our 'omics' and mechanistic studies imply that a disruptive gut microbiome may contribute to AN pathogenesis.

Anorexia nervosa (AN) is a serious mental health condition and eating disorder characterized by distorted body image, obsessive thoughts about food, ritualistic patterns of behaviour including reduced food intake, loss of body weight, raised physical activity and emotional rigidity¹. AN primarily affects women in about 95% of cases and has a population prevalence of about 1%². It can be classified into two subtypes, the

common restricting (AN-RS) type and the less prevalent binge-eating or purging (AN-BP) type¹. The evidence base for treatment is lacking³ and although specialized multidisciplinary treatment can reduce mortality⁴, less than half of AN cases achieve complete remission⁵. The aggregate mortality rate is estimated to be 5.6% per decade, much higher than in the general population⁶.

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

e-mail: oluf@sund.ku.dk
Despite research to determine the aetiology of AN, it remains a syndrome, that is, a collection of symptoms without a well-defined unifying cause. Twin studies have reported heritability estimates of 50–60%⁷ and genome-wide association studies have identified eight genomic loci showing correlations with psychiatric disorders, physical activity, and metabolic and anthropometric traits. This is independent of common variants associated with body-mass index⁸,⁹. At the pathophysiological level, AN is characterized by multiple endocrine changes¹⁰ and perturbed signalling of neurotransmitters in various parts of the brain¹¹.

The human digestive tract contains complex assemblies of microorganisms that can impact host metabolism, immunity and neurobiology via metabolites and other pathways¹². This may include the gut-microbiota-brain axis, which can affect brain functions including regulation of appetite, behaviour and emotions¹³. For example, the bacterial metabolite caseinolytic peptidase B (ClpB), predominantly produced by enterobacteria, is an antigenic mimic of α-melanocyte-stimulating hormone, which can exert anorexigenic effects¹⁴,¹⁵.

It has been hypothesized that an aberrant gut microbiota may be involved in the pathogenesis of AN. Several small studies that used amplicon sequencing to characterize the gut microbiota at the genus level in AN have been published¹⁶⁻¹⁹, showing dysbiosis of gut bacterial microbiota (see Supplementary Note 1). Moreover, in a mouse model of anorexia, changes in the gut microbiota have been shown to be associated with changes in eating behaviour and expression of hypothalamic neuropeptides²⁰.

Here we explored the hypothesis that a perturbed intestinal gut microbiota and serum metabolome contribute to the complex pathogenesis of AN. To do so, we performed shotgun metagenomics on faecal samples from 77 female AN cases and 70 age-matched female controls allowing for in-depth analyses of the gut bacterial and archaeal microbiota at taxonomic, functional and genetics levels, as well as analyses of the viral gut microbiota. We also characterized the serum metabolome, which was analysed with the gut metagenome data in relation to individual markers of eating- and psychological behaviour. Causal mechanisms were explored in silico using bidirectional mediation analyses and in vivo through faecal microbiota transplantation (FMT) of gut microbiota from AN cases to female germ-free littermates. Our findings lend support to the hypothesis that a disrupted AN gut microbiota and associated bacterial metabolites contribute to AN pathogenesis (Extended Data Fig. 1).

## Results

### Phenotypes of women with AN and controls

Summarized statistics of clinical characteristics of the 77 enrolled women with AN and 70 age-matched control women considered to be of healthy weight (HC) are given in Supplementary Table 1. The validated Eating Disorder Inventory-3 (EDI-3) questionnaire was used to estimate levels of specific eating behaviour²¹ and a detailed description of the EDI-3 subscale is shown in Supplementary Note 2. As expected, women with AN were much leaner, had lower fasting serum concentrations of glucose and insulin, higher insulin sensitivity as estimated by homoeo-static model assessment for insulin resistance (HOMA-IR) and lower serum C-reactive protein. Detailed baseline characteristics of study participants are shown in Supplementary Table 2. Additionally, within AN cases, AN-RS type cases were characterized by higher values of serum insulin and lower insulin sensitivity than AN-BP individuals (Supplementary Table 1). Comparing stool samples from AN and HC, there was no significant difference in bacterial cell counts between AN and HC or within subtypes of AN ($P_{Wilcoxon}$ > 0.05, Supplementary Table 1).

### Gut microbiota composition is altered in AN

At the phylum level, AN microbiota samples were characterized by a reduction in Bacteroidota and Actinobacteriota (Extended Data Fig. 2a). At the family level, Bacteroidaceae was dominant in both groups (Extended Data Fig. 2b). Among the top 20 most abundant families,

the abundance of Christensenellales CAG-138 was higher in AN as a new observation for this cohort, whereas the abundance of Ruminococcaceae and Lachnospiraceae was higher in HC. Among the 89 identified bacterial families, Christensenellaceae was the most significantly enriched in AN (Extended Data Fig. 2c). We observed higher β-diversity of AN microbiome at the genus level (Extended Data Fig. 2d), with *Bacteroides* being the dominant phylotype in both groups (Extended Data Fig. 2e). Among the top 30 genera, *Faecalibacterium*, *Agathobacter*, *Gemmiger*, unclassified Lachnospiraceae G, *Ruminococcus 2*, *Roseburia*, *Dysosmobacter* – *Oscillibacter*, *Coprococcus*, *Oscillospirales* 4 CAG-103 and *Eisenbergiella* were more abundant in HC, while unclassified Christensenellales CAG-138 was more abundant in AN (Extended Data Fig. 2e). Additionally, among the 225 identified bacterial genera, *Lactobacillus* was the most significantly enriched in AN (Extended Data Fig. 2f). Despite the difference in major genera between AN and HC, we found that richness of metagenomic species pangenomes (MSP, hereafter called species) was similar between the two groups (Extended Data Fig. 2g). In enterotype analyses²², we found a higher prevalence of the Ruminococcacea-enterotype (R-enterotype) in AN compared with HC, and a higher prevalence of the same enterotype in AN-BP than in AN-RS subtype (Extended Data Fig. 2h).

At the species level (Supplementary Table 3), we observed that the AN gut microbiota is characterized by higher β-diversity (Fig. 1a). Within AN subgroups, the AN-BP subtype had a more heterogeneous bacterial community at the species level than the AN-RS subtype (Fig. 1b). In the comparison between in- and out-patient AN cases, we found that β-diversity of in-patients was higher than that of out-patients at the species level (Supplementary Fig. 1 and Data). Recent body weight change within 4 weeks was not associated with gut bacterial composition (Supplementary Table 4). The species that were significantly different in abundance distribution between AN and HC after deconfounding interferences of multiple medications (selective serotonin re-uptake inhibitors, antipsychotics and benzodiazepines, specified in Supplementary Table 2) are shown in Fig. 1c. Among the depleted species in AN were *Roseburia intestinalis* and *Roseburia inulinivorans*, species that have a high capacity for digesting plant polysaccharides and are considered to be part of the health-related gut microbiota²³. In a co-abundance undirected network analysis (Extended Data Fig. 3), we identified a bacterial community consisting of *Eisenbergiella*, butyrate-producing bacterium SS3/4 - (*Clostridium*) sp. CAG:81, *Faecalibacterium prausnitzii* 3, (*Oscillibacter*) sp. ER4/Firmicutes bacterium CAG: 129_59_24, *Oscillibacter* sp. 57_20, (*Clostridium*) sp. 2789STDY5834924, unclassified Lachnospiraceae and unclassified *Dysosmobacter* – *Oscillibacter*, which was more abundant in HC. A community that was highly enriched in AN comprised *Erysipelatoclostridium ramosum*, *Enterocloster boltae*, (*Clostridium*) *innocuum* and *Blautia* sp. CAG:257.

### Associations between absolute bacterial abundance and bioclinical variables in AN

We investigated numerical covariations between the absolute abundance of bacteria at genus and species levels and bioclinical variables in the combined HC and AN cohort. We used a linear regression model adjusting for confounders including age, smoking and multiple drug intake (Methods and Supplementary Note 3, Fig. 2 and Data).

Interestingly, some bacterial taxa were linked to eating disorder scores and psychological conditions after adjusting for multiple confounding factors including age, body mass index (BMI), smoking history and medications. At the species level, we found that *Clostridium* species were positively correlated with eating disorder scores (Fig. 1d), indicating a potential role of these species in the regulation of eating behaviour and neuropsychiatric symptoms²⁴. Moreover, among the bacterial species that were inversely correlated with eating disorder scores, we found that the absolute abundances of *Lactococcus acidophilus*²⁵ and *Faecalibacterium prausnitzii*²⁶, both of which are associated with depressive symptoms, were related to a score for interpersonal
Despite research to determine the aetiology of AN, it remains a syndrome, that is, a collection of symptoms without a well-defined unifying cause. Twin studies have reported heritability estimates of 50–60%⁷ and genome-wide association studies have identified eight genomic loci showing correlations with psychiatric disorders, physical activity, and metabolic and anthropometric traits. This is independent of common variants associated with body-mass index⁸,⁹. At the pathophysiological level, AN is characterized by multiple endocrine changes¹⁰ and perturbed signalling of neurotransmitters in various parts of the brain¹¹.

The human digestive tract contains complex assemblies of microorganisms that can impact host metabolism, immunity and neurobiology via metabolites and other pathways¹². This may include the gut-microbiota-brain axis, which can affect brain functions including regulation of appetite, behaviour and emotions¹³. For example, the bacterial metabolite caseinolytic peptidase B (ClpB), predominantly produced by enterobacteria, is an antigenic mimic of α-melanocyte-stimulating hormone, which can exert anorexigenic effects¹⁴,¹⁵.

It has been hypothesized that an aberrant gut microbiota may be involved in the pathogenesis of AN. Several small studies that used amplicon sequencing to characterize the gut microbiota at the genus level in AN have been published¹⁶⁻¹⁹, showing dysbiosis of gut bacterial microbiota (see Supplementary Note 1). Moreover, in a mouse model of anorexia, changes in the gut microbiota have been shown to be associated with changes in eating behaviour and expression of hypothalamic neuropeptides²⁰.

Here we explored the hypothesis that a perturbed intestinal gut microbiota and serum metabolome contribute to the complex pathogenesis of AN. To do so, we performed shotgun metagenomics on faecal samples from 77 female AN cases and 70 age-matched female controls allowing for in-depth analyses of the gut bacterial and archaeal microbiota at taxonomic, functional and genetics levels, as well as analyses of the viral gut microbiota. We also characterized the serum metabolome, which was analysed with the gut metagenome data in relation to individual markers of eating- and psychological behaviour. Causal mechanisms were explored in silico using bidirectional mediation analyses and in vivo through faecal microbiota transplantation (FMT) of gut microbiota from AN cases to female germ-free littermates. Our findings lend support to the hypothesis that a disrupted AN gut microbiota and associated bacterial metabolites contribute to AN pathogenesis (Extended Data Fig. 1).

## Results

### Phenotypes of women with AN and controls

Summarized statistics of clinical characteristics of the 77 enrolled women with AN and 70 age-matched control women considered to be of healthy weight (HC) are given in Supplementary Table 1. The validated Eating Disorder Inventory-3 (EDI-3) questionnaire was used to estimate levels of specific eating behaviour²¹ and a detailed description of the EDI-3 subscale is shown in Supplementary Note 2. As expected, women with AN were much leaner, had lower fasting serum concentrations of glucose and insulin, higher insulin sensitivity as estimated by homoeo-static model assessment for insulin resistance (HOMA-IR) and lower serum C-reactive protein. Detailed baseline characteristics of study participants are shown in Supplementary Table 2. Additionally, within AN cases, AN-RS type cases were characterized by higher values of serum insulin and lower insulin sensitivity than AN-BP individuals (Supplementary Table 1). Comparing stool samples from AN and HC, there was no significant difference in bacterial cell counts between AN and HC or within subtypes of AN ($P_{Wilcoxon}$ > 0.05, Supplementary Table 1).

### Gut microbiota composition is altered in AN

At the phylum level, AN microbiota samples were characterized by a reduction in Bacteroidota and Actinobacteriota (Extended Data Fig. 2a). At the family level, Bacteroidaceae was dominant in both groups (Extended Data Fig. 2b). Among the top 20 most abundant families,

the abundance of Christensenellales CAG-138 was higher in AN as a new observation for this cohort, whereas the abundance of Ruminococcaceae and Lachnospiraceae was higher in HC. Among the 89 identified bacterial families, Christensenellaceae was the most significantly enriched in AN (Extended Data Fig. 2c). We observed higher β-diversity of AN microbiome at the genus level (Extended Data Fig. 2d), with *Bacteroides* being the dominant phylotype in both groups (Extended Data Fig. 2e). Among the top 30 genera, *Faecalibacterium*, *Agathobacter*, *Gemmiger*, unclassified Lachnospiraceae G, *Ruminococcus 2*, *Roseburia*, *Dysosmobacter* – *Oscillibacter*, *Coprococcus*, *Oscillospirales* 4 CAG-103 and *Eisenbergiella* were more abundant in HC, while unclassified Christensenellales CAG-138 was more abundant in AN (Extended Data Fig. 2e). Additionally, among the 225 identified bacterial genera, *Lactobacillus* was the most significantly enriched in AN (Extended Data Fig. 2f). Despite the difference in major genera between AN and HC, we found that richness of metagenomic species pangenomes (MSP, hereafter called species) was similar between the two groups (Extended Data Fig. 2g). In enterotype analyses²², we found a higher prevalence of the Ruminococcacea-enterotype (R-enterotype) in AN compared with HC, and a higher prevalence of the same enterotype in AN-BP than in AN-RS subtype (Extended Data Fig. 2h).

At the species level (Supplementary Table 3), we observed that the AN gut microbiota is characterized by higher β-diversity (Fig. 1a). Within AN subgroups, the AN-BP subtype had a more heterogeneous bacterial community at the species level than the AN-RS subtype (Fig. 1b). In the comparison between in- and out-patient AN cases, we found that β-diversity of in-patients was higher than that of out-patients at the species level (Supplementary Fig. 1 and Data). Recent body weight change within 4 weeks was not associated with gut bacterial composition (Supplementary Table 4). The species that were significantly different in abundance distribution between AN and HC after deconfounding interferences of multiple medications (selective serotonin re-uptake inhibitors, antipsychotics and benzodiazepines, specified in Supplementary Table 2) are shown in Fig. 1c. Among the depleted species in AN were *Roseburia intestinalis* and *Roseburia inulinivorans*, species that have a high capacity for digesting plant polysaccharides and are considered to be part of the health-related gut microbiota²³. In a co-abundance undirected network analysis (Extended Data Fig. 3), we identified a bacterial community consisting of *Eisenbergiella*, butyrate-producing bacterium SS3/4 - (*Clostridium*) sp. CAG:81, *Faecalibacterium prausnitzii* 3, (*Oscillibacter*) sp. ER4/Firmicutes bacterium CAG: 129_59_24, *Oscillibacter* sp. 57_20, (*Clostridium*) sp. 2789STDY5834924, unclassified Lachnospiraceae and unclassified *Dysosmobacter* – *Oscillibacter*, which was more abundant in HC. A community that was highly enriched in AN comprised *Erysipelatoclostridium ramosum*, *Enterocloster boltae*, (*Clostridium*) *innocuum* and *Blautia* sp. CAG:257.

### Associations between absolute bacterial abundance and bioclinical variables in AN

We investigated numerical covariations between the absolute abundance of bacteria at genus and species levels and bioclinical variables in the combined HC and AN cohort. We used a linear regression model adjusting for confounders including age, smoking and multiple drug intake (Methods and Supplementary Note 3, Fig. 2 and Data).

Interestingly, some bacterial taxa were linked to eating disorder scores and psychological conditions after adjusting for multiple confounding factors including age, body mass index (BMI), smoking history and medications. At the species level, we found that *Clostridium* species were positively correlated with eating disorder scores (Fig. 1d), indicating a potential role of these species in the regulation of eating behaviour and neuropsychiatric symptoms²⁴. Moreover, among the bacterial species that were inversely correlated with eating disorder scores, we found that the absolute abundances of *Lactococcus acidophilus*²⁵ and *Faecalibacterium prausnitzii*²⁶, both of which are associated with depressive symptoms, were related to a score for interpersonal
<figure><img src="image_5.png" /><figcaption>**Fig. 1** | Alterations in gut bacterial species in AN cases compared to healthy controls, and associations with eating disorder scores. **a,b** Box plot (line, median; box, interquartile range (IQR); whiskers, 1.5× IQR) of β-diversity of AN (n = 77) and HC (n = 70) gut microbiota (a) and of two AN subtypes (AN-RS n = 56, AN-BP n = 21) and HC gut microbiota (b) at bacterial species level (Canberra distance). Statistical significance of differences between two groups was determined by Wilcoxon rank-sum test (two-sided). **c** Significantly contrasted bacterial species between AN and HC. Differences in abundance were detected using the metadecounfDR pipeline where covariates including age, BMI, smoking and</figcaption></figure>

multiple drug intake were corrected. Cliff's delta values give estimates of effect size. For each contrasted MSP, prevalence in the whole cohort, HC, AN, and P<sub>adj</sub> are given next to the MSP annotation. **d** Heat map showing that gut bacterial species are linked to eating disorder scores in AN cases, using a linear regression model where age, BMI, smoking and multiple drug intake were defined as covariates and adjusted for. Variables in specific eating disorder scale are marked in blue, and general psychological scale is marked in red. Right panel to the heat map indicates the direction of each variable. For each MSP, prevalence in AN is given next to the MSP annotation. +, P<sub>adj</sub> < 0.05 by Benjamini-Hochberg method (see for exact P values).

alienation (Fig. 1d). Additionally, the absolute abundance of *Parasutterella* correlated positively with body dissatisfaction and absolute abundance of *Bifidobacterium* correlated with a marker of perfectionism. Despite a similar absolute abundance of *Brachyspira* in AN and HC, this genus was positively correlated with markers for ‘drive for thinness’ in AN (Extended Data Fig. 4). As shown in Extended Data Fig. 5, we found no difference in circulating levels of anorexigenic ClpB<sup>14</sup> between AN and HC groups (Supplementary Note 4).

## Predicted bacterial growth rates are altered in AN

We estimated growth dynamics of the bacterial gut microbiota from the metagenomic data by calculating the peak-to-trough ratio (PTR)<sup>27</sup> for 50 bacterial species. Thirty-five of these were present in more than 20 samples. The median PTR values differed markedly between AN and HC (P<sub>Wilcoxon</sub> = 2.0 × 10<sup>-4</sup>, Extended Data Fig. 6), which might be related to the severe reduction in food intake in AN patients. Six bacteria were predicted to have significantly lower growth rate in AN (P<sub>Wilcoxon</sub> < 0.05, Extended Data Fig. 6). These were *Akkermansia muciniphila*, *Alistipes finegoldii*, *Coprococcus catus*, *Eubacterium siraeum*, *Odoribacter splanchicus* and *butyrate-producing bacterium* SS3/4.

## Lactococcus phages and weakened trans-kingdom interactions in AN

We observed higher viral richness (Chao1, Fig. 2a) and diversity (Shannon, Fig. 2b) in AN faecal samples compared with HC. Recent body weight change within 4 weeks was not associated with gut viral composition (Supplementary Table 4). After deconfounding for covariates (age, smoking and drug intake), we identified 31 viral species that

were enriched or decreased in AN (Fig. 2c). Of major interest, 25 of the 30 increased viral species in AN were *Lactococcus* phages with known *Lactococcus lactis* hosts—bacteria that have been extensively used in the production of fermented food products.

Analysis of viral-bacterial correlations within the AN and HC gut microbiota revealed a remarkable decrease in the number of viral-bacterial interactions in AN (219 for HC versus 84 for AN, P<sub>Fisher's exact test</sub> = 8.8 × 10<sup>-15</sup>, Fig. 2d). This was primarily driven by weakened interactions between viral species and short-chain fatty acid bacterial producers, such as *Roseburia inulinivorans*, *Faecalibacterium prausnitzii* and *Roseburia hominis* (Supplementary Fig. 3 and Data). We did not observe any interplay between *Lactococcus* phages and *Lactococcus* bacteria in the trans-kingdom analysis (Supplementary Figs. 3 and 4 and Data).

In analyses of AN subgroups, principal coordinate analysis (PCoA) on the Canberra distance showed no remarkable alterations in gut viral composition (P<sub>PERMANOVA</sub> = 0.571, Supplementary Fig. 5 and Data). However, when comparing the number of viral-bacterial correlations in AN-RS and AN-BP gut microbiota, we found a reduction in the number and a much lower ratio of positive correlations in AN-RS (164 for AN-BP versus 44 for AN-RS, P<sub>Fisher's exact test</sub> < 2.2 × 10<sup>-16</sup>) (Fig. 2d and Supplementary Fig. 4). This suggests that the gut microbiota in AN-RS has weakened gut viral-bacterial interactions.

## Predicted gut microbiota functions correlate with eating behaviours and metabolism

Using Gut Metabolic Modules (GMMs)<sup>28</sup> and Gut Brain Modules (GBMs)<sup>29</sup> to predict gut bacterial functional potentials, we identified 159 functional modules. Notably, the abundance of GBMs for serotonin
a

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

c

<figure><img src="image_7.png" /><figcaption>Cliff's delta</figcaption></figure>

b

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

d

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

**Fig. 1 | Alterations in gut bacterial species in AN cases compared to healthy controls, and associations with eating disorder scores.** a,b Box plot (line, median; box, interquartile range (IQR); whiskers, 1.5× IQR) of β-diversity of AN (n = 77) and HC (n = 70) gut microbiota (a) and of two AN subtypes (AN-RS n = 56, AN-BP n = 21) and HC gut microbiota (b) at bacterial species level (Canberra distance). Statistical significance of differences between two groups was determined by Wilcoxon rank-sum test (two-sided). c, Significantly contrasted bacterial species between AN and HC. Differences in abundance were detected using the metadecounfDR pipeline where covariates including age, BMI, smoking and

multiple drug intake were corrected. Cliff's delta values give estimates of effect size. For each contrasted MSP, prevalence in the whole cohort, HC, AN, and P<sub>adj</sub> are given next to the MSP annotation. **d.** Heat map showing that gut bacterial species are linked to eating disorder scores in AN cases, using a linear regression model where age, BMI, smoking and multiple drug intake were defined as covariates and adjusted for. Variables in specific eating disorder scale are marked in blue, and general psychological scale is marked in red. Right panel to the heat map indicates the direction of each variable. For each MSP, prevalence in AN is given next to the MSP annotation. +, P<sub>adj</sub> < 0.05 by Benjamini-Hochberg method (see for exact P values).

alienation (Fig. 1d). Additionally, the absolute abundance of *Parasutterella* correlated positively with body dissatisfaction and absolute abundance of *Bifidobacterium* correlated with a marker of perfectionism. Despite a similar absolute abundance of *Brachyspira* in AN and HC, this genus was positively correlated with markers for ‘drive for thinness’ in AN (Extended Data Fig. 4). As shown in Extended Data Fig. 5, we found no difference in circulating levels of anorexigenic ClpB<sup>14</sup> between AN and HC groups (Supplementary Note 4).

Predicted bacterial growth rates are altered in AN

We estimated growth dynamics of the bacterial gut microbiota from the metagenomic data by calculating the peak-to-trough ratio (PTR)²⁷ for 50 bacterial species. Thirty-five of these were present in more than 20 samples. The median PTR values differed markedly between AN and HC (P<sub>Wilcoxon</sub> = 2.0 × 10⁻⁴, Extended Data Fig. 6), which might be related to the severe reduction in food intake in AN patients. Six bacteria were predicted to have significantly lower growth rate in AN (P<sub>Wilcoxon</sub> < 0.05, Extended Data Fig. 6). These were *Akkermansia muciniphila*, *Alistipes finegoldii*, *Coprococcus catus*, *Eubacterium siraeum*, *Odoribacter splanchicus* and *butyrate-producing bacterium* SS3/4.

Lactococcus phages and weakened trans-kingdom
interactions in AN

We observed higher viral richness (Chao1, Fig. 2a) and diversity (Shannon, Fig. 2b) in AN faecal samples compared with HC. Recent body weight change within 4 weeks was not associated with gut viral composition (Supplementary Table 4). After deconfounding for covariates (age, smoking and drug intake), we identified 31 viral species that

were enriched or decreased in AN (Fig. 2c). Of major interest, 25 of the 30 increased viral species in AN were *Lactococcus* phages with known *Lactococcus lactis* hosts—bacteria that have been extensively used in the production of fermented food products.

Analysis of viral-bacterial correlations within the AN and GC gut microbiota revealed a remarkable decrease in the number of viral-bacterial interactions in AN (219 for HC versus 84 for AN, P<sub>Fisher's exact test</sub> = 8.8 × 10⁻¹⁵, Fig. 2d). This was primarily driven by weakened interactions between viral species and short-chain fatty acid bacterial producers, such as *Roseburia inulinivorans*, *Faecalibacterium prausnitzii* and *Roseburia hominis* (Supplementary Fig. 3 and Data). We did not observe any interplay between *Lactococcus* phages and *Lactococcus* bacteria in the trans-kingdom analysis (Supplementary Figs. 3 and 4 and Data).

In analyses of AN subgroups, principal coordinate analysis (PCoA) on the Canberra distance showed no remarkable alterations in gut viral composition (P<sub>PERMANOVA</sub> = 0.571, Supplementary Fig. 5 and Data). However, when comparing the number of viral-bacterial correlations in AN-RS and AN-BP gut microbiota, we found a reduction in the number and a much lower ratio of positive correlations in AN-RS (164 for AN-BP versus 44 for AN-RS, P<sub>Fisher's exact test</sub> < 2.2 × 10⁻¹⁶) (Fig. 2d and Supplementary Fig. 4). This suggests that the gut microbiota in AN-RS has weakened gut viral-bacterial interactions.

Predicted gut microbiota functions correlate with eating
behaviours and metabolism

Using Gut Metabolic Modules (GMMs)²⁸ and Gut Brain Modules (GBMs)²⁹ to predict gut bacterial functional potentials, we identified 159 functional modules. Notably, the abundance of GBMs for serotonin
<figure><img src="image_6.png" /><figcaption>**Fig. 2 | The viral gut microbiota differs between AN cases and controls.**
**a,b,** Box plot (line, median; box, IQR; whiskers, 1.5× IQR) of changes in Chao1 richness (a) and Shannon diversity (b) of the viral gut microbiota between AN (n = 77) and HC (n = 70) at viral species level. Significance was examined by two-sided Wilcoxon rank-sum test (a,b). **c,** Cliff's delta values of contrasted gut viral species between AN and HC with P<0.05 by Benjamini-Hochberg correction</figcaption></figure>

(given next to the viral annotation). The differential species were identified by metadefoundR pipeline where impacts of cofactors including age, smoking and multiple drug intake were deconfounded. **d,** Difference in number of transkingdom ecologic correlations between the viral and bacterial gut microbiota in AN (n = 77) compared to HC (n = 70), and between two AN subtypes (AN-RS n = 56, AN-BP n = 21) using the SparCC algorithm.

biosynthesis and the degradation of dopamine, glutamate and tryptophan, which are metabolites with effects on mood and appetite, were enriched in AN (Fig. 3a). Conversely, the abundance of the glutamate synthesis II and vitamin K2 pathways were higher in HC (Fig. 3a)³⁰. In addition, we found that serotonin synthesis and glutamate degradation pathways were inversely correlated with circulating concentrations of glucose and insulin, or insulin sensitivity (Fig. 3b). While we observed differences in GBMs, we did not identify differences in GMMs between AN and HC (Supplementary Table 5).

Bacterial genomes may have structural variants (SVs) that potentially interfere with functional genes impacting the interplay between microbes and their host³¹. Therefore, differences in the presence or abundance of SVs between otherwise identical bacterial strains may underlie critical phenotypic and functional differences³¹,³². Here we profiled SVs in all samples and identified 5,056 deletion SVs and 2,423 variable SVs across 56 bacterial species (Fig. 4a,b). For some species, we observed marked differences in copy number variation. We identified 87 deletion SVs and 18 variable SVs in *Bacteroides uniformis* across 134 individuals, 78 deletions and 15 variable SVs in *Faecalibacterium prausnitzii* in 124 individuals, and 110 deletion SVs and 55 variable SVs in *Parabacteroides distasonis* in 121 individuals. For the archaeal microbiota, we only identified SVs in *Methanobrevibacter smithii* in 13 individuals (Fig. 4a, and Supplementary Fig. 6 and Data). To explore potential differences in bacterial genetics, we further computed the Canberra distance of bacterial SV profiles between all 147 samples (Fig. 4c). AN and HC samples were significantly different in the β-diversity of SV composition (P< wilcoxon < 2.2 × 10⁻¹⁶). Collectively, these results suggest the composition of bacterial SVs differs between the two groups.

Exploring relationships between gut bacterial SVs and markers of eating behaviour, we found that in AN cases, bacterial SVs were significantly associated with eating disorder scores after deconfounding for multiple covariates (Fig. 4d). As a noteworthy example, in the association analysis between SVs in the *B. uniformis* genome and host eating scores, we found that a 10-kbp deletion was directly associated with markers of bulimia and self-denial, indicating that this deletion SV may be involved in the regulation of eating disorder and psychological traits in AN (Fig. 4e). Indeed, gene analysis showed that the beginning of this specific genomic region in *B. uniformis* encodes a thiamine-monophosphate kinase (Fig. 4f and Supplementary Table 6), which is the distal enzyme involved in the thiamine (vitamin B1) biosynthesis pathway. Thiamine deficiency has been associated with mental health including memory loss, anxiety, depression, irritability, insomnia as well as appetite loss and gastrointestinal complaints³³. About one third of AN cases may suffer from thiamine deficiency³⁴. We found that AN cases lacking this bacterial genomic region had higher scores for bulimia (a key characteristic for the AN-BP subtype) and self-denial (Fig. 4g), a pattern also suggested by our correlation analyses (Fig. 4e). Another example linking bacterial genetics to metabolism-related bioclinical variables is given in (Extended Data Fig. 7) and Supplementary Note 5.

## Serum metabolites are associated with markers of appetite and food intake regulation

We performed untargeted metabolomics profiling of serum from AN cases and controls. This revealed a serum metabolome profile consisting of 28 polar metabolites and 35 microbiota-related metabolites, which was significantly different between AN and HC (Fig. 5a), while
<figure><img src="image_10.png" /><figcaption>**Fig. 2 | The viral gut microbiota differs between AN cases and controls.**
**a,b,** Box plot (line, median; box, IQR; whiskers, 1.5× IQR) of changes in Chao1 richness (a) and Shannon diversity (b) of the viral gut microbiota between AN (n = 77) and HC (n = 70) at viral species level. Significance was examined by two-sided Wilcoxon rank-sum test (a,b). **c,** Cliff's delta values of contrasted gut viral species between AN and HC with P<0.05 by Benjamini-Hochberg correction</figcaption></figure>

(given next to the viral annotation). The differential species were identified by metadefoundR pipeline where impacts of cofactors including age, smoking and multiple drug intake were deconfounded. **d,** Difference in number of trans-kingdom ecologic correlations between the viral and bacterial gut microbiota in AN (n = 77) compared to HC (n = 70), and between two AN subtypes (AN-RS n = 56, AN-BP n = 21) using the SparCC algorithm.

biosynthesis and the degradation of dopamine, glutamate and tryptophan, which are metabolites with effects on mood and appetite, were enriched in AN (Fig. 3a). Conversely, the abundance of the glutamate synthesis II and vitamin K2 pathways were higher in HC (Fig. 3a)³⁰. In addition, we found that serotonin synthesis and glutamate degradation pathways were inversely correlated with circulating concentrations of glucose and insulin, or insulin sensitivity (Fig. 3b). While we observed differences in GBMs, we did not identify differences in GMMs between AN and HC (Supplementary Table 5).

Bacterial genomes may have structural variants (SVs) that potentially interfere with functional genes impacting the interplay between microbes and their host³¹. Therefore, differences in the presence or abundance of SVs between otherwise identical bacterial strains may underlie critical phenotypic and functional differences³¹,³². Here we profiled SVs in all samples and identified 5,056 deletion SVs and 2,423 variable SVs across 56 bacterial species (Fig. 4a,b). For some species, we observed marked differences in copy number variation. We identified 87 deletion SVs and 18 variable SVs in *Bacteroides uniformis* across 134 individuals, 78 deletions and 15 variable SVs in *Faecalibacterium prausnitzii* in 124 individuals, and 110 deletion SVs and 55 variable SVs in *Parabacteroides distasonis* in 121 individuals. For the archaeal microbiota, we only identified SVs in *Methanobrevibacter smithii* in 13 individuals (Fig. 4a, and Supplementary Fig. 6 and Data). To explore potential differences in bacterial genetics, we further computed the Canberra distance of bacterial SV profiles between all 147 samples (Fig. 4c). AN and HC samples were significantly different in the β-diversity of SV composition (P< wilcoxon < 2.2 × 10⁻¹⁶). Collectively, these results suggest the composition of bacterial SVs differs between the two groups.

Exploring relationships between gut bacterial SVs and markers of eating behaviour, we found that in AN cases, bacterial SVs were significantly associated with eating disorder scores after deconfounding for multiple covariates (Fig. 4d). As a noteworthy example, in the association analysis between SVs in the *B. uniformis* genome and host eating scores, we found that a 10-kbp deletion was directly associated with markers of bulimia and self-denial, indicating that this deletion SV may be involved in the regulation of eating disorder and psychological traits in AN (Fig. 4e). Indeed, gene analysis showed that the beginning of this specific genomic region in *B. uniformis* encodes a thiamine-monophosphate kinase (Fig. 4f and Supplementary Table 6), which is the distal enzyme involved in the thiamine (vitamin B1) biosynthesis pathway. Thiamine deficiency has been associated with mental health including memory loss, anxiety, depression, irritability, insomnia as well as appetite loss and gastrointestinal complaints³³. About one third of AN cases may suffer from thiamine deficiency³⁴. We found that AN cases lacking this bacterial genomic region had higher scores for bulimia (a key characteristic for the AN-BP subtype) and self-denial (Fig. 4g), a pattern also suggested by our correlation analyses (Fig. 4e). Another example linking bacterial genetics to metabolism-related bioclinical variables is given in (Extended Data Fig. 7) and Supplementary Note 5.

## Serum metabolites are associated with markers of appetite and food intake regulation

We performed untargeted metabolomics profiling of serum from AN cases and controls. This revealed a serum metabolome profile consisting of 28 polar metabolites and 35 microbiota-related metabolites, which was significantly different between AN and HC (Fig. 5a), while
<figure><img src="image_7.png" /><figcaption>**Fig. 3 | Predicted functional potentials of the bacterial gut microbiome in AN cases and healthy control participants.** a, Cliff's delta effect size of contrasted functional modules between AN (n = 77) and HC (n = 70) using the metadeconfoundR pipeline where interferences from covariates including age, BMI, smoking and multiple drug intake were corrected. Gold bars, functional modules more abundant in AN; blue bars, functional modules more abundant in HC.</figcaption></figure>

For each contrasted module, *P* value after Benjamini-Hochberg correction is given next to the module annotation. b, Heat map of the associations between clinical variables and functional potentials of gut bacteriome by linear regression model where impacts of covariates including age, smoking and multiple drug intake were deconfounded. + indicates *P* < 0.05 by Benjamini-Hochberg correction (see for exact *P* values).

only slightly altered between the two AN subtypes (Extended Data Fig. 8a). We identified 25 serum metabolites with concentration differences between cases and controls after adjusting for confounders (adjusted *P* value (*P*<sub>adj</sub>) < 0.05) (Fig. 5b).

Serum concentrations of both primary (cholic acid (CA), glyco-
cholic acid (GCA)) and secondary bile acids (glycohyocholic acid
(GHCA), 7-oxo-hyocholic acid (7-oxo-HCA), glycohyodeoxycholic
acid (GHDCA), 7-oxo-deoxycholic acid (7-oxo-DCA), ω/α-muricholic
acid (ω/α-MCA), ursodeoxycholic acid (UDCA)) were higher in AN,
indicating a potential role of the gut microbiota in AN-related changes
in secondary bile acid synthesis and metabolism, and satiety regula-
tion³⁵. Moreover, serum concentrations of indole-3-acetic acid and
indole-3-propionic acid (IPA), two tryptophan metabolites, were
higher in AN compared with HC. Interestingly, IPA is associated with

the secretion of glucagon-like peptide 1, which can stimulate satiety
and slow gastric emptying<sup>36,37</sup>. Dysregulation of valine, saturated and
long-chain unsaturated fatty acids were also observed in the AN group
(Fig. 5b and Supplementary Note 6).

**Associations between serum metabolites, gut microbiota, and markers of appetite and mental health**

To explore the role of serum metabolites in the interaction between gut
microbiota and host phenotypes, we constructed bidirectional media-
tion models. For direction 1, we hypothesized that serum metabolites
(including microbiota-related metabolites) as variables mediate the
causal effect of gut microbial features (bacterial species, gut brain
modules and bacterial genetics) on host phenotypes (Eating Disorder
Inventory-3 (EDI-3) scores and metabolic traits). For direction 2, we
<figure><img src="image_11.png" /><figcaption>**Fig. 3 | Predicted functional potentials of the bacterial gut microbiome in AN cases and healthy control participants.** a, Cliff's delta effect size of contrasted functional modules between AN (n = 77) and HC (n = 70) using the metadeconfoundR pipeline where interferences from covariates including age, BMI, smoking and multiple drug intake were corrected. Gold bars, functional modules more abundant in AN; blue bars, functional modules more abundant in</figcaption></figure>

HC. For each contrasted module, *P* value after Benjamini-Hochberg correction is given next to the module annotation. **b**, Heat map of the associations between clinical variables and functional potentials of gut bacteriome by linear regression model where impacts of covariates including age, smoking and multiple drug intake were deconfounded. + indicates *P* < 0.05 by Benjamini-Hochberg correction (see for exact *P* values).

only slightly altered between the two AN subtypes (Extended Data Fig. 8a). We identified 25 serum metabolites with concentration differences between cases and controls after adjusting for confounders (adjusted *P* value (*P*<sub>adj</sub>) < 0.05) (Fig. 5b).

Serum concentrations of both primary (cholic acid (CA), glyco-
cholic acid (GCA)) and secondary bile acids (glycohyocholic acid
(GHCA), 7-oxo-hyocholic acid (7-oxo-HCA), glycohyodeoxycholic
acid (GHDCA), 7-oxo-deoxycholic acid (7-oxo-DCA), ω/α-muricholic
acid (ω/α-MCA), ursodeoxycholic acid (UDCA)) were higher in AN,
indicating a potential role of the gut microbiota in AN-related changes
in secondary bile acid synthesis and metabolism, and satiety regula-
tion³⁵. Moreover, serum concentrations of indole-3-acetic acid and
indole-3-propionic acid (IPA), two tryptophan metabolites, were
higher in AN compared with HC. Interestingly, IPA is associated with

the secretion of glucagon-like peptide 1, which can stimulate satiety
and slow gastric emptying<sup>36,37</sup>. Dysregulation of valine, saturated and
long-chain unsaturated fatty acids were also observed in the AN group
(Fig. 5b and Supplementary Note 6).

**Associations between serum metabolites, gut microbiota, and markers of appetite and mental health**

To explore the role of serum metabolites in the interaction between gut
microbiota and host phenotypes, we constructed bidirectional media-
tion models. For direction 1, we hypothesized that serum metabolites
(including microbiota-related metabolites) as variables mediate the
causal effect of gut microbial features (bacterial species, gut brain
modules and bacterial genetics) on host phenotypes (Eating Disorder
Inventory-3 (EDI-3) scores and metabolic traits). For direction 2, we
treated phenotypes as mediators and alterations in metabolites as outcomes of changes in microbial features (Fig. 5c). This in silico bidirectional analysis enabled us to quantify the extent to which a hypothesized mediator (a metabolite) participates in the interaction between a cause (microbial features) and its effect (host phenotypic traits).

We first performed the causal inference of bacterial features on EDI-3 questionnaire scores in AN cases (Extended Data Fig. 8b). Inferred causal relationships in direction 1 consisted of 13 microbial features as initiators, 11 metabolites as mediators and 7 EDI-3 scores as outcomes (Fig. 5d). As a noteworthy example, we identified 'drive for thinness' as an outcome of the AN-enriched bacterial species *C. paraputrificum* (Fig. 1c), which was linked to changes in serum IPA levels, also enriched in the AN group. *C. paraputrificum* is a producer of multiple tryptophan catabolites, including IPA, indoleacrylic acid, indoleacetic acid and tryptamine, all of which are involved in regulating appetite and mental health<sup>24</sup> (Fig. 5e).

In another example, indoxyl sulfate, which was enriched in AN patients, was identified as a mediator of AN-enriched (*Clostridium*) sp. CAG:269 in the score for body dissatisfaction (Fig. 5e). Indoxyl sulfate is a cardio- and uraemic toxin, and has been shown to induce anxiety- or depression-like behaviours in humans and animal models<sup>38,39</sup>. *Clostridium* is a genus of indole-producing bacteria that can encode tryptophanase<sup>40</sup>, the critical enzyme converting tryptophan to indole, pyruvate and ammonia<sup>41</sup>.

Next, we analysed causal inference between microbial features and host metabolic traits across the whole cohort (Extended Data Fig. 8c). The causal network in direction 1 consisted of 14 microbial features, 10 metabolites and 5 metabolic traits as causal treatments, mediators and outcomes, respectively (Extended Data Fig. 8d). Notably, we found that the serotonin synthesis module causally affected host BMI via the secondary bile acid glycoursodeoxycholic acid, which is upregulated by serotonin<sup>42</sup> (Extended Data Fig. 8e). Finally, consistent with our previous findings, serum leucine mediated the impact of *B. vulgatus* on glucose homoeostasis<sup>43</sup> (Extended Data Fig. 8e). We observed no unidirectional causal relationship between changes in eating behaviour and psychological status, gut microbial features and metabolites (Supplementary Table 7).

## Gut microbiota induces reduced weight gain and altered energy metabolism in mice

To investigate potential causal relationships between an altered gut microbiota in AN and relevant phenotypes, we transplanted faecal microbiota from three randomly chosen AN-RS cases (to achieve uniformity as AN-BP has a more heterogeneous phenotype than the AN-RS subtype) and three age-matched HC participants to three independent litters of female germ-free (GF) mice (Extended Data Fig. 9a). To minimize variations in genetic background, we included littermates as control mice. In each litter study, 8, 6 and 6 littermates, respectively, were randomly assigned as recipients of AN or HC faecal microbiota. Following stool transplantation, recipient mice were singly housed and received a 30% calorie restricted diet for 3 weeks to mimic reduced food intake in human AN (Extended Data Fig. 9b). Ad libitum chow diet feeding did not generate any phenotypic alterations in GF mice<sup>44</sup>, an observation consistent with a previous report on kwashiorkor<sup>45</sup> (Extended Data Fig. 10a).

After 21 d, GF mice transplanted with stools from AN cases showed a larger initial decrease in body weight and a slower weight gain over time compared with mice that received HC FMT (Fig. 6a and Extended Data Fig. 10b; see Supplementary Note 7 for discussion).

We performed hypothalamic gene expression analysis following FMT. AN-transplanted and HC-transplanted recipients differed in expression of several hypothalamic genes involved in the control of eating behaviour and energy expenditure (Fig. 6b). This included increased expression of the appetite suppressors *Bdnf*<sup>46</sup> and *Cartpt*<sup>47</sup>, and the receptor for serotonin, *Htr1b* (involved in the downstream regulation of serotonin), in AN FMT recipients. Expression of *Snca*, which encodes the neuronal protein alpha-synuclein associated with several neurodegenerative diseases, was higher in AN-transplanted mice<sup>48</sup>. In addition, we analysed messenger RNA (mRNA) levels of genes encoding proteins regulating adipose tissue thermogenesis. We found that abundances of *Ucp1*, *Elovl3* and *Pgc1α* mRNA were increased in inguinal fat of AN-transplanted mice, indicating enhanced adipose tissue thermogenesis in this group of mice (Fig. 6c).

16S ribosomal ribonucleic acid (rRNA) gene amplicon sequencing of human donor stools and GF mouse recipient caecal contents identified 85 overlapping amplicon sequence variants (ASVs; Fig. 6d), of which 45 (53%) altered ASVs in donors were transferred to recipients (Fig. 6e). Serum metabolome profiling detected 31 conserved metabolites between donors and recipients, and alterations of 19 of these (61%) were transferred from human to mouse (Extended Data Fig. 10c). We identified three ASVs: ASV_021, ASV_229 and ASV_002 that were annotated as *Bacteroides* at the genus level. The relative abundance of these ASVs and the expression of browning genes, including *Ucp1* (Extended Data Fig. 10d), were strongly positively correlated, indicating a potential role of these ASVs in body weight loss or lower weight gain through enhanced adipose browning. The inverse correlation between the relative abundance of ASV_122, annotated as genus *Akkermansia*, and the hypothalamic appetite-suppressing gene *Htr1b* is also notable and suggests a potential role of ASV_122 in appetite regulation (Extended Data Fig. 10d). Taken together, alterations in hypothalamic and adipose tissue gene expression and changes in body weight over time in mice suggest that a disrupted gut microbiota in human AN may contribute to some of the elements in the complex pathogenesis of AN.

## Discussion

Using a combination of sequencing and metabolomics to characterize the gut microbiome and metabolome in humans and mice, we show that the bacterial and viral components of the microbiome and the serum metabolome are altered in those with AN compared with healthy individuals. We also find that bacterial species SVs are different between AN and healthy controls, and in silico causal inference analyses imply that bacterial metabolites mediate some effects of an altered gut microbiota on AN behaviour. Finally, GF mice transplanted with AN stools on an energy-restricted diet initially lose more weight and have slower weight gain over time compared with mice transplanted with stools from healthy individuals. This was associated with higher expression of appetite suppressor genes in the hypothalamus and higher expression of thermogenesis-related genes in adipose tissue of AN-transplanted mice.

**Fig. 4 | Structural variations in the bacterial and archaeal gut microbiota in AN cases and healthy controls.** a, Number of SVs of each bacterial or archaeal species in 147 (77 cases and 70 controls) study participants. For each species, the number of deletion and variable SVs are given. b, Pie chart showing the total identified SVs numbers. c, Box plot (line, median; box, IQR; whiskers, 1.5× IQR) of β-diversity (Canberra distance) of SV-based genetic composition in AN (*n* = 77) and HC (*n* = 70) bacteriome. *P* value was determined by two-sided Wilcoxon rank-sum test. d, Chord diagram showing significant associations between eating disorder scores and bacterial SVs after adjusting for age, BMI, smoking and multiple drug intake. e, Heat map showing the associations between SVs

of *Bacteroides uniformis* and EDI-3 scores using linear regression model where impacts of age, BMI, smoking and multiple drug intake were deconfounded. + indicates *P* < 0.05 corrected by Benjamini-Hochberg (see for exact *P* values). In d and e, variables of eating disorder scale are coloured in blue, general psychological scale are in red, and bacterial SVs are in black. f, The deletion rate of the 10-kbp deletion SV harbouring thiamine-monophosphate kinase in *B. uniformis* genome in the AN group. g, Box plot (line, median; box, IQR; whiskers, 1.5× IQR) showing the EDI-3 scores in anorexia individuals with (*n* = 49) and without (*n* = 28) the 10-kbp deletion. Significance was determined by Wilcoxon rank-sum test (two-sided).
treated phenotypes as mediators and alterations in metabolites as outcomes of changes in microbial features (Fig. 5c). This in silico bidirectional analysis enabled us to quantify the extent to which a hypothesized mediator (a metabolite) participates in the interaction between a cause (microbial features) and its effect (host phenotypic traits).

We first performed the causal inference of bacterial features on EDI-3 questionnaire scores in AN cases (Extended Data Fig. 8b). Inferred causal relationships in direction 1 consisted of 13 microbial features as initiators, 11 metabolites as mediators and 7 EDI-3 scores as outcomes (Fig. 5d). As a noteworthy example, we identified 'drive for thinness' as an outcome of the AN-enriched bacterial species *C. paraputrificum* (Fig. 1c), which was linked to changes in serum IPA levels, also enriched in the AN group. *C. paraputrificum* is a producer of multiple tryptophan catabolites, including IPA, indoleacrylic acid, indoleacetic acid and tryptamine, all of which are involved in regulating appetite and mental health<sup>24</sup> (Fig. 5e).

In another example, indoxyl sulfate, which was enriched in AN patients, was identified as a mediator of AN-enriched (*Clostridium*) sp. CAG:269 in the score for body dissatisfaction (Fig. 5e). Indoxyl sulfate is a cardio- and uraemic toxin, and has been shown to induce anxiety- or depression-like behaviours in humans and animal models<sup>38,39</sup>. *Clostridium* is a genus of indole-producing bacteria that can encode tryptophanase<sup>40</sup>, the critical enzyme converting tryptophan to indole, pyruvate and ammonia<sup>41</sup>.

Next, we analysed causal inference between microbial features and host metabolic traits across the whole cohort (Extended Data Fig. 8c). The causal network in direction 1 consisted of 14 microbial features, 10 metabolites and 5 metabolic traits as causal treatments, mediators and outcomes, respectively (Extended Data Fig. 8d). Notably, we found that the serotonin synthesis module causally affected host BMI via the secondary bile acid glycoursodeoxycholic acid, which is upregulated by serotonin<sup>42</sup> (Extended Data Fig. 8e). Finally, consistent with our previous findings, serum leucine mediated the impact of *B. vulgatus* on glucose homoeostasis<sup>43</sup> (Extended Data Fig. 8e). We observed no unidirectional causal relationship between changes in eating behaviour and psychological status, gut microbial features and metabolites (Supplementary Table 7).

## Gut microbiota induces reduced weight gain and altered energy metabolism in mice

To investigate potential causal relationships between an altered gut microbiota in AN and relevant phenotypes, we transplanted faecal microbiota from three randomly chosen AN-RS cases (to achieve uniformity as AN-BP has a more heterogeneous phenotype than the AN-RS subtype) and three age-matched HC participants to three independent litters of female germ-free (GF) mice (Extended Data Fig. 9a). To minimize variations in genetic background, we included littermates as control mice. In each litter study, 8, 6 and 6 littermates, respectively, were randomly assigned as recipients of AN or HC faecal microbiota. Following stool transplantation, recipient mice were singly housed and received a 30% calorie restricted diet for 3 weeks to mimic reduced food intake in human AN (Extended Data Fig. 9b). Ad libitum chow diet feeding did not generate any phenotypic alterations in GF mice<sup>44</sup>, an observation consistent with a previous report on kwashiorkor<sup>45</sup> (Extended Data Fig. 10a).

After 21 d, GF mice transplanted with stools from AN cases showed a larger initial decrease in body weight and a slower weight gain over time compared with mice that received HC FMT (Fig. 6a and Extended Data Fig. 10b; see Supplementary Note 7 for discussion).

We performed hypothalamic gene expression analysis following FMT. AN-transplanted and HC-transplanted recipients differed in expression of several hypothalamic genes involved in the control of eating behaviour and energy expenditure (Fig. 6b). This included increased expression of the appetite suppressors *Bdnf*<sup>46</sup> and *Cartpt*<sup>47</sup>, and the receptor for serotonin, *Htr1b* (involved in the downstream regulation of serotonin), in AN FMT recipients. Expression of *Snca*, which encodes the neuronal protein alpha-synuclein associated with several neurodegenerative diseases, was higher in AN-transplanted mice<sup>48</sup>. In addition, we analysed messenger RNA (mRNA) levels of genes encoding proteins regulating adipose tissue thermogenesis. We found that abundances of *Ucp1*, *Elovl3* and *Pgc1α* mRNA were increased in inguinal fat of AN-transplanted mice, indicating enhanced adipose tissue thermogenesis in this group of mice (Fig. 6c).

16S ribosomal ribonucleic acid (rRNA) gene amplicon sequencing of human donor stools and GF mouse recipient caecal contents identified 85 overlapping amplicon sequence variants (ASVs; Fig. 6d), of which 45 (53%) altered ASVs in donors were transferred to recipients (Fig. 6e). Serum metabolome profiling detected 31 conserved metabolites between donors and recipients, and alterations of 19 of these (61%) were transferred from human to mouse (Extended Data Fig. 10c). We identified three ASVs: ASV_021, ASV_229 and ASV_002 that were annotated as *Bacteroides* at the genus level. The relative abundance of these ASVs and the expression of browning genes, including *Ucp1* (Extended Data Fig. 10d), were strongly positively correlated, indicating a potential role of these ASVs in body weight loss or lower weight gain through enhanced adipose browning. The inverse correlation between the relative abundance of ASV_122, annotated as genus *Akkermansia*, and the hypothalamic appetite-suppressing gene *Htr1b* is also notable and suggests a potential role of ASV_122 in appetite regulation (Extended Data Fig. 10d). Taken together, alterations in hypothalamic and adipose tissue gene expression and changes in body weight over time in mice suggest that a disrupted gut microbiota in human AN may contribute to some of the elements in the complex pathogenesis of AN.

## Discussion

Using a combination of sequencing and metabolomics to characterize the gut microbiome and metabolome in humans and mice, we show that the bacterial and viral components of the microbiome and the serum metabolome are altered in those with AN compared with healthy individuals. We also find that bacterial species SVs are different between AN and healthy controls, and in silico causal inference analyses imply that bacterial metabolites mediate some effects of an altered gut microbiota on AN behaviour. Finally, GF mice transplanted with AN stools on an energy-restricted diet initially lose more weight and have slower weight gain over time compared with mice transplanted with stools from healthy individuals. This was associated with higher expression of appetite suppressor genes in the hypothalamus and higher expression of thermogenesis-related genes in adipose tissue of AN-transplanted mice.

**Fig. 4 | Structural variations in the bacterial and archaeal gut microbiota in AN cases and healthy controls.** a, Number of SVs of each bacterial or archaeal species in 147 (77 cases and 70 controls) study participants. For each species, the number of deletion and variable SVs are given. b, Pie chart showing the total identified SVs numbers. c, Box plot (line, median; box, IQR; whiskers, 1.5× IQR) of β-diversity (Canberra distance) of SV-based genetic composition in AN (*n* = 77) and HC (*n* = 70) bacteriome. *P* value was determined by two-sided Wilcoxon rank-sum test. d, Chord diagram showing significant associations between eating disorder scores and bacterial SVs after adjusting for age, BMI, smoking and multiple drug intake. e, Heat map showing the associations between SVs

of *Bacteroides uniformis* and EDI-3 scores using linear regression model where impacts of age, BMI, smoking and multiple drug intake were deconfounded. + indicates *P* < 0.05 corrected by Benjamini-Hochberg (see for exact *P* values). In d and e, variables of eating disorder scale are coloured in blue, general psychological scale are in red, and bacterial SVs are in black. f, The deletion rate of the 10-kbp deletion SV harbouring thiamine-monophosphate kinase in *B. uniformis* genome in the AN group. g, Box plot (line, median; box, IQR; whiskers, 1.5× IQR) showing the EDI-3 scores in anorexia individuals with (*n* = 49) and without (*n* = 28) the 10-kbp deletion. Significance was determined by Wilcoxon rank-sum test (two-sided).
In both FMT experiments and in silico inference analyses, we observed changes in circulating levels of glycine-chenodeoxycholic acid, indole-3-propionic acid, taurine-α-muricholic acid and

taurine-hyodeoxycholic acid. We propose that these metabolites may act as potential mediators of some of the AN phenotypes. For example, indole-3-propionic acid is a tryptophan metabolite, which is implicated

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

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

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

Bacteroides uniformis 3591_3597, 3599_3610

Bacteroides uniformis 2749_2750 and 12 segments

Bacteroides uniformis 2765_2769, 3625_3626

Bacteroides uniformis 3886_3887 and 2 segments

Bacteroides uniformis 32956_2957 and 12 segments

Bacteroides uniformis 1707_1708 and 2 segments

Bacteroides uniformis 3203_3204, 3204_3205

Bacteroides uniformis 3880_3881

Bacteroides uniformis 217_218

Bacteroides uniformis 605_608 and 5 segments

Bacteroides uniformis 702_703 and 2 segments

Bacteroides uniformis 829_832

Bacteroides uniformis 2550_2553 and 5 segments

Bacteroides uniformis 2584_2587, 3613_3615

Bacteroides uniformis 2588_2589

Bacteroides uniformis 2942_2944

Bacteroides uniformis 3065_3066

Bacteroides uniformis 33114_3115 and 10 segments

Bacteroides uniformis 3212_3222

Bacteroides uniformis 3300_3303, 3303_3304

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

<figure><img src="image_12.png" /></figure>
In both FMT experiments and in silico inference analyses, we observed changes in circulating levels of glycine-chenodeoxycholic acid, indole-3-propionic acid, taurine-α-muricholic acid and

taurine-hyodeoxycholic acid. We propose that these metabolites may act as potential mediators of some of the AN phenotypes. For example, indole-3-propionic acid is a tryptophan metabolite, which is implicated

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

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

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

Bacteroides uniformis 3591_3597, 3599_3610

Bacteroides uniformis 2749_2750 and 12 segments

Bacteroides uniformis 2765_2769, 3625_3626

Bacteroides uniformis 3886_3887 and 2 segments

Bacteroides uniformis 32956_2957 and 12 segments

Bacteroides uniformis 1707_1708 and 2 segments

Bacteroides uniformis 3203_3204, 3204_3205

Bacteroides uniformis 3880_3881

Bacteroides uniformis 217_218

Bacteroides uniformis 605_608 and 5 segments

Bacteroides uniformis 702_703 and 2 segments

Bacteroides uniformis 829_832

Bacteroides uniformis 2550_2553 and 5 segments

Bacteroides uniformis 2584_2587, 3613_3615

Bacteroides uniformis 2588_2589

Bacteroides uniformis 2942_2944

Bacteroides uniformis 3065_3066

Bacteroides uniformis 33114_3115 and 10 segments

Bacteroides uniformis 3212_3222

Bacteroides uniformis 3300_3303, 3303_3304

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

<figure><img src="image_16.png" /></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 5 | Serum metabolites differ between AN cases and controls and may be mediating the impact of gut microbial features on eating disorder traits.**</figcaption></figure>

a, Principal component analysis (PCA) of the serum metabolome profile of AN cases and HC participants. b, Cliff's delta values of contrasted metabolites between AN (n = 77) and HC (n = 70) after adjusting for age, BMI, smoking and multiple drug intake. Gold lollipops are metabolites enriched in AN, and blue lollipops show serum metabolites enriched in HC. c, Workflow for the bidirectional mediation analysis for gut microbial features, serum metabolites and host phenotypes.

d, Sankey diagram showing the inferred causal relationship network of direction

1 where gut microbial features including bacterial species, gut brain/metabolic modules and bacterial genetics were treated as causal factors, metabolites are mediators, and EDI-3 scores are outcomes. e, Examples of inferred causal relationships between microbial features, metabolites and EDI-3 scores. Direction 1 means microbial features → eating disorder scores mediated by metabolites, illustrated with a black line; direction 2 means microbial features → metabolites mediated by EDI-3 scores, illustrated with a dashed red line. The proportions of mediation effects are shown at the centre of ring charts. FFA, free fatty acid.

in serotonin activity⁴⁹, and hyodeoxycholic acid is a 6α-hydroxylated bile acid, also called muricholic acid, which reduces body weight gain⁵⁰. This bile acid is involved in satiety regulation⁵¹. Serotonin activity as well as appetite regulation could be implicated in development and/or maintenance of the AN syndrome. Future studies will need to explore the individual and combined effects of these metabolites on energy metabolism. Still, many more gut bacterial species and derived metabolites may mediate the observed AN traits in humans and mice, as reported in an activity-based anorexia model²⁰,⁵².

Analyses of deletion and variable SVs in gut bacterial species indicated that bacterial genetics may influence AN-relevant behavioural traits and pathophysiology. Of special interest is a 10-kbp deletion in the *B. uniformis* genome that was associated with estimates of bulimia and self-denial. We predicted that this deletion results in the loss of a thiamine-monophosphate kinase encoding gene, which may result in a relative deficiency of microbiota-produced thiamine. This is of interest in the context of AN pathology since various mental and intestinal health issues are known to be linked to thiamine deficiency³³.
<figure><img src="image_17.png" /><figcaption>**Fig. 5 | Serum metabolites differ between AN cases and controls and may be mediating the impact of gut microbial features on eating disorder traits.**</figcaption></figure>

a, Principal component analysis (PCA) of the serum metabolome profile of AN cases and HC participants. b, Cliff's delta values of contrasted metabolites between AN (n = 77) and HC (n = 70) after adjusting for age, BMI, smoking and multiple drug intake. Gold lollipops are metabolites enriched in AN, and blue lollipops show serum metabolites enriched in HC. c, Workflow for the bidirectional mediation analysis for gut microbial features, serum metabolites and host phenotypes.

d, Sankey diagram showing the inferred causal relationship network of direction

1 where gut microbial features including bacterial species, gut brain/metabolic modules and bacterial genetics were treated as causal factors, metabolites are mediators, and EDI-3 scores are outcomes. e, Examples of inferred causal relationships between microbial features, metabolites and EDI-3 scores. Direction 1 means microbial features → eating disorder scores mediated by metabolites, illustrated with a black line; direction 2 means microbial features → metabolites mediated by EDI-3 scores, illustrated with a dashed red line. The proportions of mediation effects are shown at the centre of ring charts. FFA, free fatty acid.

in serotonin activity⁴⁹, and hyodeoxycholic acid is a 6α-hydroxylated bile acid, also called muricholic acid, which reduces body weight gain⁵⁰. This bile acid is involved in satiety regulation⁵¹. Serotonin activity as well as appetite regulation could be implicated in development and/or maintenance of the AN syndrome. Future studies will need to explore the individual and combined effects of these metabolites on energy metabolism. Still, many more gut bacterial species and derived metabolites may mediate the observed AN traits in humans and mice, as reported in an activity-based anorexia model²⁰,⁵².

Analyses of deletion and variable SVs in gut bacterial species indicated that bacterial genetics may influence AN-relevant behavioural traits and pathophysiology. Of special interest is a 10-kbp deletion in the *B. uniformis* genome that was associated with estimates of bulimia and self-denial. We predicted that this deletion results in the loss of a thiamine-monophosphate kinase encoding gene, which may result in a relative deficiency of microbiota-produced thiamine. This is of interest in the context of AN pathology since various mental and intestinal health issues are known to be linked to thiamine deficiency³³.
<figure><img src="image_14.png" /><figcaption>**Fig. 6** | FMT using samples from AN donors induces AN-relevant phenotypes in GF mice. **a**, Body weight (BW) change compared to the body weight at day 0 after energy-restricted diet (AN-T *n* = 10, HC-T, *n* = 10 examined over 3 independent experiments). Significance was calculated by two-way analysis of variance (ANOVA), followed by Benjamini-Hochberg post hoc test. **b,c**, mRNA levels of the indicated mice genes in hypothalamus (b) and inguinal white adipose tissue (c) in the faecal microbiota mouse recipients (AN-T *n* = 10, HC-T, *n* = 10, examined over 3 independent experiments). Significance between the two</figcaption></figure>

groups was tested using unpaired two-tailed Student's t-test. Data are presented as mean ± s.e.m. (a-c). **d**, Venn diagram of the identified and transferred ASVs between human donors and GF mouse recipients. **e**, Left: heat map of the 84 conserved ASVs in human donors. Middle: differences in the 84 conserved ASVs derived from the caecal content between AN-T and HC-T GF mouse recipients (AN-T *n* = 10, HC-T, *n* = 10, examined over 3 independent experiments). Transferred microbial alterations are marked in blue. Right: taxonomic information of ASVs.
<figure><img src="image_18.png" /><figcaption>**Fig. 6** | FMT using samples from AN donors induces AN-relevant phenotypes in GF mice. **a**, Body weight (BW) change compared to the body weight at day 0 after energy-restricted diet (AN-T *n* = 10, HC-T, *n* = 10 examined over 3 independent experiments). Significance was calculated by two-way analysis of variance (ANOVA), followed by Benjamini-Hochberg post hoc test. **b,c**, mRNA levels of the indicated mice genes in hypothalamus (b) and inguinal white adipose tissue (c) in the faecal microbiota mouse recipients (AN-T *n* = 10, HC-T, *n* = 10, examined over 3 independent experiments). Significance between the two</figcaption></figure>

groups was tested using unpaired two-tailed Student's t-test. Data are presented as mean ± s.e.m. (a–c). **d**, Venn diagram of the identified and transferred ASVs between human donors and GF mouse recipients. **e**, Left: heat map of the 84 conserved ASVs in human donors. Middle: differences in the 84 conserved ASVs derived from the caecal content between AN-T and HC-T GF mouse recipients (AN-T *n* = 10, HC-T, *n* = 10, examined over 3 independent experiments). Transferred microbial alterations are marked in blue. Right: taxonomic information of ASVs.
Our studies of the viral gut microbiota in AN showed a partial uncoupling of the ecological interactions between viral species and short-chain producing bacterial species with impact on brain biology⁵³. In addition, AN samples were enriched in *Lactococcus* phages with known bacterial *L. lactis* hosts. The AN-associated enrichment in *Lactococcus* phages may result in disruption in food fermentation and suggest the possible use of fermented food items in future treatment of AN. The reason for this enrichment of *Lactococcus* phages is unsettled but our finding may justify testing a multistrain probiotic containing *L. lactis* in AN adolescents⁵⁴.

Our study has limitations: (1) the use of a cross-sectional Danish AN cohort limits generalizability of findings to other ethnicities; (2) the same limitation holds for the use of AN patients in treatment in a specialized centre, as these patients may not be representative of milder forms of AN; and (3) although effects of multiple covariates were deconfounded, we had no information on diet and physical activity—behaviours that impact the gut microbiota.

In conclusion, the present multi-omics study uncovers profound and complex disruptions of the gut microbiota in individuals with AN, with functional implications and altered serum metabolites. These compounds may act via the blood circulation or via gut-microbiota-brain neuronal signalling pathways affecting brain regulation of appetite, emotions and behaviour. FMT from human AN donors to GF mice under energy-restricted feeding resulted in lower body weight gain and a number of changes in expression of hypothalamic and adipose tissue genes involved in controlling behaviour and energy homeostasis. The combination of multi-omics and in vivo experiments complement our causal inference analyses to allow the identification of specific bacterial metabolites that potentially mediate human host AN traits. Our findings lend support to the hypothesis that a severely disrupted intestinal microbiota contributes to some of the stages in the pathogenesis of AN.

# Methods

## Study participants

The AN patients were all diagnosed by an experienced consultant psychiatrist and they met the DSM-5 criteria for the restricting or binge/purging subtypes of AN¹. Since 90–95% of individuals with diagnosed AN are females, we decided to include only women in the present study and since ethnicity may influence gut microbiota, we only included Danish Caucasian women with AN cases, recruited from three specialized centres in Denmark from 1 September 2014 to 31 July 2016. The centres were: Center for Eating Disorders (Odense University Hospital), Child and Adolescent Psychiatric Unit (Aarhus University Hospital) and Unit for Psychiatric Research (Aalborg University Hospital). Exclusion criteria comprised antibiotic or antifungal treatment within the previous 3 months, any acute or chronic somatic diseases or infections. All the included patients were treated in specialized centres, and they were interviewed by an experienced and specialized psychologist or psychiatrist at the start of their treatment. The validated Eating Disorder Inventory (EDI, details given in Supplementary Note 2) was used for the interview and as a questionnaire filled out by trained health professional specialists²¹. The exclusion criteria for the age-matched healthy control women were BMI below 18.5 or above 25 kg m⁻², regular medication of any kind apart from birth control pills, and antibiotics within the last 3 months. The control participants were recruited via public advertisement and via direct contact to health staff, medical students and their relatives. BMI and other clinical characteristics of the healthy controls are listed in Supplementary Table 2.

The study protocol was registered at ClinicalTrials.gov (NCT02217384) and the study was conducted in accordance with the Helsinki declaration and approved by the Regional Scientific Ethical Committee for Southern Denmark (file no 42053 S-20140040). All participants involved in this study provided written informed consent.

In-patients. All the below-described measurements were conducted during routine treatment in the specialized unit for somatic

and psychological stabilization of patients with severe AN. Safe and effective weight restoration of 2.0–3.0% per week is the goal of the treatment in the inpatient unit. Care was given by a multidisciplinary team and is in accordance with international guidelines. Individually customized meals were given under the supervision of trained nurses or dietitians at scheduled times. If the meals could not be consumed within the predefined timeframes (15 min for a snack and 30 min for a main meal), supplemental nutrition drinks were added either orally or via a duodenal tube. To account for individual preferences, a choice of three different meals was offered. The macronutrient content was consistent and within recommended energy percent ranges: 40–50% carbohydrate (maximum of 10% sugar), 30–40% fat and 20–25% protein. The daily energy intake was individualized according to the weight course. If a patient failed to reach 2% of weekly weight gain, the energy content of the daily menu was increased. All meals were followed by a supervised rest varying from 30 to 60 min in a seated position. Between the rests, light physical activity such as a walk was allowed; however, no forms of exercise training were allowed. For patients with an urge for excessive exercise or otherwise a lack of compliance with behaviour rules, behaviour supervision was extended and if needed, to 24 h a day.

Out-patients. AN patients paid regular visits to outpatient units where care was given by a multidisciplinary team involving medical doctors, nurses, psychologists and health behaviour educators.

They were given an individual diet plan, which had the same macronutrient composition as mentioned above for in-patients. As with the in-patients, all out-patients were also enrolled in cognitive behavioural therapeutic treatment courses.

Height was measured on a wall-mounted stadiometer and weight was measured in the morning before breakfast on a calibrated platform scale. BMI was calculated as weight divided by the square of height (=kg/m²).

## Biochemical analyses

Blood samples were taken in the morning after an overnight fast. The blood samples were collected on ice and processed to obtain serum and plasma, and subsequently stored at -80 °C. Serum concentrations of sodium, potassium, albumin and creatinine were measured by enzymatic assays on a Roche/Hitachi cobas c system. Serum concentrations of total cholesterol, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol were determined using the phosphotungstic acid magnesium chloride precipitation method. Serum immunoreactive insulin levels were measured using an enzyme-linked immunosorbent assay, while serum concentration of alanine aminotransferase was analysed with an enzymatic coulometric method including pyridoxal phosphate activation.

Analysis of plasma ClpB was performed as previously described⁵⁵.

## Faecal sample collection and DNA extraction

Stools were collected according to International Human Microbiome Standards (IHMS) guidelines (SOP 03 V1) in kits by AN cases and HC at home and immediately frozen at -20 °C until they were transported on dry ice and frozen 4–24 h later at -80 °C in plastic tubes at the biobanks. DNA extraction from aliquots of faecal samples was performed following IHMS SOP P7 V2⁵⁶.

## Bacterial cell counting

For bacterial cell counting (Supplementary Fig. 7), 0.08–0.12 g of frozen (-80 °C) faecal samples were diluted 15 times in pH 7.2 Dulbecco's phosphate-buffered saline (DPBS) (Sigma-Aldrich), mechanically homogenized using tissue lyser (40 min, 12.5 agitations per second; QIAGEN) and fixed with 2% paraformaldehyde (10 min, room temperature; Biotum). Then the samples were diluted 120 times in filtered staining buffer (1 mM EDTA, 0.01% Tween20, pH 7.2 DPBS, 1% BSA (Sigma-Aldrich)). To minimize clumps, the samples were filtered through a cell strainer (pore size 5 µm; pluriSelect), pre-wet in the
Our studies of the viral gut microbiota in AN showed a partial uncoupling of the ecological interactions between viral species and short-chain producing bacterial species with impact on brain biology⁵³. In addition, AN samples were enriched in *Lactococcus* phages with known bacterial *L. lactis* hosts. The AN-associated enrichment in *Lactococcus* phages may result in disruption in food fermentation and suggest the possible use of fermented food items in future treatment of AN. The reason for this enrichment of *Lactococcus* phages is unsettled but our finding may justify testing a multistrain probiotic containing *L. lactis* in AN adolescents⁵⁴.

Our study has limitations: (1) the use of a cross-sectional Danish AN cohort limits generalizability of findings to other ethnicities; (2) the same limitation holds for the use of AN patients in treatment in a specialized centre, as these patients may not be representative of milder forms of AN; and (3) although effects of multiple covariates were deconfounded, we had no information on diet and physical activity—behaviours that impact the gut microbiota.

In conclusion, the present multi-omics study uncovers profound and complex disruptions of the gut microbiota in individuals with AN, with functional implications and altered serum metabolites. These compounds may act via the blood circulation or via gut-microbiota-brain neuronal signalling pathways affecting brain regulation of appetite, emotions and behaviour. FMT from human AN donors to GF mice under energy-restricted feeding resulted in lower body weight gain and a number of changes in expression of hypothalamic and adipose tissue genes involved in controlling behaviour and energy homeostasis. The combination of multi-omics and in vivo experiments complement our causal inference analyses to allow the identification of specific bacterial metabolites that potentially mediate human host AN traits. Our findings lend support to the hypothesis that a severely disrupted intestinal microbiota contributes to some of the stages in the pathogenesis of AN.

# Methods

## Study participants

The AN patients were all diagnosed by an experienced consultant psychiatrist and they met the DSM-5 criteria for the restricting or binge/purging subtypes of AN¹. Since 90–95% of individuals with diagnosed AN are females, we decided to include only women in the present study and since ethnicity may influence gut microbiota, we only included Danish Caucasian women with AN cases, recruited from three specialized centres in Denmark from 1 September 2014 to 31 July 2016. The centres were: Center for Eating Disorders (Odense University Hospital), Child and Adolescent Psychiatric Unit (Aarhus University Hospital) and Unit for Psychiatric Research (Aalborg University Hospital). Exclusion criteria comprised antibiotic or antifungal treatment within the previous 3 months, any acute or chronic somatic diseases or infections. All the included patients were treated in specialized centres, and they were interviewed by an experienced and specialized psychologist or psychiatrist at the start of their treatment. The validated Eating Disorder Inventory (EDI, details given in Supplementary Note 2) was used for the interview and as a questionnaire filled out by trained health professional specialists²¹. The exclusion criteria for the age-matched healthy control women were BMI below 18.5 or above 25 kg m⁻², regular medication of any kind apart from birth control pills, and antibiotics within the last 3 months. The control participants were recruited via public advertisement and via direct contact to health staff, medical students and their relatives. BMI and other clinical characteristics of the healthy controls are listed in Supplementary Table 2.

The study protocol was registered at ClinicalTrials.gov (NCT02217384) and the study was conducted in accordance with the Helsinki declaration and approved by the Regional Scientific Ethical Committee for Southern Denmark (file no 42053 S-20140040). All participants involved in this study provided written informed consent.

In-patients. All the below-described measurements were conducted during routine treatment in the specialized unit for somatic

and psychological stabilization of patients with severe AN. Safe and effective weight restoration of 2.0–3.0% per week is the goal of the treatment in the inpatient unit. Care was given by a multidisciplinary team and is in accordance with international guidelines. Individually customized meals were given under the supervision of trained nurses or dietitians at scheduled times. If the meals could not be consumed within the predefined timeframes (15 min for a snack and 30 min for a main meal), supplemental nutrition drinks were added either orally or via a duodenal tube. To account for individual preferences, a choice of three different meals was offered. The macronutrient content was consistent and within recommended energy percent ranges: 40–50% carbohydrate (maximum of 10% sugar), 30–40% fat and 20–25% protein. The daily energy intake was individualized according to the weight course. If a patient failed to reach 2% of weekly weight gain, the energy content of the daily menu was increased. All meals were followed by a supervised rest varying from 30 to 60 min in a seated position. Between the rests, light physical activity such as a walk was allowed; however, no forms of exercise training were allowed. For patients with an urge for excessive exercise or otherwise a lack of compliance with behaviour rules, behaviour supervision was extended and if needed, to 24 h a day.

Out-patients. AN patients paid regular visits to outpatient units where care was given by a multidisciplinary team involving medical doctors, nurses, psychologists and health behaviour educators.

They were given an individual diet plan, which had the same macronutrient composition as mentioned above for in-patients. As with the in-patients, all out-patients were also enrolled in cognitive behavioural therapeutic treatment courses.

Height was measured on a wall-mounted stadiometer and weight was measured in the morning before breakfast on a calibrated platform scale. BMI was calculated as weight divided by the square of height (=kg/m²).

## Biochemical analyses

Blood samples were taken in the morning after an overnight fast. The blood samples were collected on ice and processed to obtain serum and plasma, and subsequently stored at -80 °C. Serum concentrations of sodium, potassium, albumin and creatinine were measured by enzymatic assays on a Roche/Hitachi cobas c system. Serum concentrations of total cholesterol, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol were determined using the phosphotungstic acid magnesium chloride precipitation method. Serum immunoreactive insulin levels were measured using an enzyme-linked immunosorbent assay, while serum concentration of alanine aminotransferase was analysed with an enzymatic coulometric method including pyridoxal phosphate activation.

Analysis of plasma ClpB was performed as previously described⁵⁵.

## Faecal sample collection and DNA extraction

Stools were collected according to International Human Microbiome Standards (IHMS) guidelines (SOP 03 V1) in kits by AN cases and HC at home and immediately frozen at -20 °C until they were transported on dry ice and frozen 4–24 h later at -80 °C in plastic tubes at the biobanks. DNA extraction from aliquots of faecal samples was performed following IHMS SOP P7 V2⁵⁶.

## Bacterial cell counting

For bacterial cell counting (Supplementary Fig. 7), 0.08–0.12 g of frozen (-80 °C) faecal samples were diluted 15 times in pH 7.2 Dulbecco's phosphate-buffered saline (DPBS) (Sigma-Aldrich), mechanically homogenized using tissue lyser (40 min, 12.5 agitations per second; QIAGEN) and fixed with 2% paraformaldehyde (10 min, room temperature; Biotum). Then the samples were diluted 120 times in filtered staining buffer (1 mM EDTA, 0.01% Tween20, pH 7.2 DPBS, 1% BSA (Sigma-Aldrich)). To minimize clumps, the samples were filtered through a cell strainer (pore size 5 µm; pluriSelect), pre-wet in the
staining buffer. Next, the bacterial cell suspension was stained with SYBR Green I (1:200,000; Thermo Fisher) in DMSO (Sigma-Aldrich) and incubated in the dark for 30 min. For accurate determination of bacterial cell counts, a known concentration of 123count eBeads (Invitrogen) was added to the samples before the analysis. Measurements were performed using a BD Fortessa LSRII flow cytometer (BD Biosciences) and data were acquired using BD FACSDiVaTM software. A threshold value of 200 was applied on the FITC (530/30 nm) channel. Fluorescence intensity at green (530/30 nm, FITC), blue (450/50 nm, Pacific Blue), yellow (575/26 nm, PE) and red (695/40 nm, PerCP-Cy5-5) fluorescence channels as well as forward- and side-scattered (FSC and SSC) light intensities were collected. Measurements were performed at a pre-set flow rate of 0.5 µl s⁻¹. Data were processed in R using flow-core package (v1.11.20)⁵⁷ in R (v4.1.2). Fixed gating strategy separated the microbial fluorescent events from the faecal sample background.

## Shotgun sequencing

DNA was quantitated using Qubit fluorometric quantitation (Thermo Fisher) and qualified using DNA size profiling on a fragment analyser (Agilent). High molecular weight DNA (>10 kbp; 3 µg) was used to build the library. Shearing of DNA into fragments of approximately 150 bp was performed using an ultrasonicator (Covaris) and DNA fragment library construction was performed using the Ion Plus Fragment Library and Ion Xpress Barcode Adapters kits (Thermo Fisher). Purified and amplified DNA fragment libraries were sequenced using the Ion Proton Sequencer (Thermo Fisher), with a minimum of 20 million high-quality reads of 150 bp (on average) generated per library.

## Construction of a gene count table

To construct the gene count table, METEOR software was used⁵⁸: first, reads were filtered for low quality by AlienTrimmer⁵⁹. After removal of low-quality reads and human DNA reads, 75.7% ± 2.7% high-quality metagenomics sequencing reads of faecal DNA were mapped onto the Integrated Gut Catalog 2 (IGC2)⁶⁰, comprising 10.4 million genes, using Bowtie2 (ref. 61). Reads mapped to a unique gene in the catalogue were attributed to their corresponding genes. Then, reads that mapped with the same alignment score to multiple genes in the catalogue were attributed according to the ratio of their unique mapping counts to the captured genes. The resulting count table was further processed using the R package MetaOMineR v1.31⁶². It was downsized at 14 million mapped reads to take into account differences in sequencing depth and mapping rate across samples. Then the downsized matrix was normalized for gene length and transformed into a frequency matrix by fragments per kilobase of transcript per million fragments mapped normalization. Gene count was computed as the number of genes present (abundance strictly positive) in the frequency matrix.

## Profiling and annotation of MSPs and gut enterotypes

IGC2 was previously organized into 1,990 MSPs with MSPminer⁶³ using a publicly available updated MSP dataset⁶⁴. Relative abundance of MSP was computed as the mean abundance of its 100 'marker' genes (that is, the genes that correlate the most altogether). If less than 10% of 'marker' genes were seen in a sample, the abundance of the MSP was set to 0. This approach was used in the MetaHIT⁶² and Metacardis⁶⁵ consortia. For the 4 MSPs with less than 100 core genes, all available core genes were used.

Abundances at higher taxonomical ranks were computed as the sum of the MSP that belong to a given taxa. MSP count was assessed as the number of MSPs present in a sample (that is, whose abundance is strictly positive). Enterotypes profiling was performed as previously demonstrated⁶⁶.

## Estimating functional modules of gut bacteriome

Genes from the IGC2 catalogue were mapped with diamond⁶⁷ onto KEGG orthologues (KO) from the KEGG database⁶⁸ (v8.9). Each gene was

assigned to the best-ranked KO among hits with e-value < 10 × 10⁻⁵ and bit score >60. Then we assessed presence and abundance of GMMs²⁸ and GBMs²⁹ in a metagenomic sample by the pipeline implemented in the R package omixerRpm (v0.3.2) as previously described²⁸,²⁹.

## Estimation of dynamic growth rate of bacteria from metagenomic samples

We used the computational pipeline to infer gut bacterial growth dynamics from metagenomic samples as previously described²⁷. Sequencing reads were mapped to a database that contains complete genomic references of 2,991 strains belonging to 1,509 microbial species. For each bacterial species, a reference strain with prevalence of 100% across the samples was selected. A coverage map was then assembled on the basis of aligned reads to the reference genome. Genomic segments were binned into 10-kbp regions and coverage of the resulting bins was calculated and smoothed. The location of the origin and terminus of replication was predicted by fits of the same strain across multiple samples. Lastly, PTRs were calculated for each bacterial species in every sample as the smoothed sequencing coverage of the representative strain at the predicted peak location, divided by that at the predicted trough location.

## Studies of bacterial structural variations

Before the SVs classification, the pipeline with iterative coverage-based read assignment algorithm was performed to reassign the ambiguous reads to the most likely reference with high accuracy³¹,³². The reference genomes provided in the proGenomes database (http://progenomes1.embl.de/) were concatenated and then divided into genomic 1-kbp bins and applied for the detection of highly variable genomic segments. The SGV-Finder pipeline³¹ was used to detect the SVs that are either (1) with deletion percentage of the genomic segment across the population of <25% (variable SVs, vSVs), (2) with deletion percentage between 25% and 75% (deletion SVs, dSVs; the absence or presence of this particular genomic segment was kept) or (3) with deletion percentage of >75% (this genomic segment was excluded from the analysis). All bacterial species with SV calling were present in at least 10% of the total samples and were used for subsequent analysis.

## Mediation analysis

The R package 'mediation⁶⁹ (v4.5.0) was used to infer causal relationships between gut microbial features, polar and microbiota-related metabolites, and metabolic traits and eating disorder scores. To reduce the testing numbers, we only kept the candidate groups consisting of variables that were strongly associated with each other; that is, for a candidate microbial feature-metabolite-phenotypic variable causal group, the association between gut microbial feature and serum metabolite was significant ($P_{adj} < 0.1$); the association between metabolite and phenotypic variable was significant ($P_{adj} < 0.1$); and the association between microbial feature and phenotypic variable was significant ($P_{adj} < 0.1$). After performing mediation analysis, only candidate groups with significance in direction 1 were kept for Sankey diagram visualization.

## Profiling and analysis of viral gut microbiota

We profiled the viral gut microbiota using MiCoP⁷⁰, as this method is optimized to call viruses directly from bulk metagenomics sequencing reads and compute relative abundance within the virome dataset. As a reference dataset, MiCoP draws upon the NCBI's RefSeq Viral database⁷¹. We identified a total of 209 viral species with prevalence of > and =10% and relative abundance of > and =0.01% for 147 (77 AN versus 70 HC) individuals included in the dataset. Richness, alpha and beta diversity were calculated with the R package 'fossil'⁷² and 'vegan'⁷³. Two-tailed Wilcoxon's rank-sum test was used to determine statistically significant differences in richness and alpha diversity indices between groups. Permutational multivariate analysis of variance (PERMANOVA) at *n* = 999
staining buffer. Next, the bacterial cell suspension was stained with SYBR Green I (1:200,000; Thermo Fisher) in DMSO (Sigma-Aldrich) and incubated in the dark for 30 min. For accurate determination of bacterial cell counts, a known concentration of 123count eBeads (Invitrogen) was added to the samples before the analysis. Measurements were performed using a BD Fortessa LSRII flow cytometer (BD Biosciences) and data were acquired using BD FACSDiVaTM software. A threshold value of 200 was applied on the FITC (530/30 nm) channel. Fluorescence intensity at green (530/30 nm, FITC), blue (450/50 nm, Pacific Blue), yellow (575/26 nm, PE) and red (695/40 nm, PerCP-Cy5-5) fluorescence channels as well as forward- and side-scattered (FSC and SSC) light intensities were collected. Measurements were performed at a pre-set flow rate of 0.5 µl s⁻¹. Data were processed in R using flow-core package (v1.11.20)⁵⁷ in R (v4.1.2). Fixed gating strategy separated the microbial fluorescent events from the faecal sample background.

## Shotgun sequencing

DNA was quantitated using Qubit fluorometric quantitation (Thermo Fisher) and qualified using DNA size profiling on a fragment analyser (Agilent). High molecular weight DNA (>10 kbp; 3 µg) was used to build the library. Shearing of DNA into fragments of approximately 150 bp was performed using an ultrasonicator (Covaris) and DNA fragment library construction was performed using the Ion Plus Fragment Library and Ion Xpress Barcode Adapters kits (Thermo Fisher). Purified and amplified DNA fragment libraries were sequenced using the Ion Proton Sequencer (Thermo Fisher), with a minimum of 20 million high-quality reads of 150 bp (on average) generated per library.

## Construction of a gene count table

To construct the gene count table, METEOR software was used⁵⁸: first, reads were filtered for low quality by AlienTrimmer⁵⁹. After removal of low-quality reads and human DNA reads, 75.7% ± 2.7% high-quality metagenomics sequencing reads of faecal DNA were mapped onto the Integrated Gut Catalog 2 (IGC2)⁶⁰, comprising 10.4 million genes, using Bowtie2 (ref. 61). Reads mapped to a unique gene in the catalogue were attributed to their corresponding genes. Then, reads that mapped with the same alignment score to multiple genes in the catalogue were attributed according to the ratio of their unique mapping counts to the captured genes. The resulting count table was further processed using the R package MetaOMineR v1.31⁶². It was downsized at 14 million mapped reads to take into account differences in sequencing depth and mapping rate across samples. Then the downsized matrix was normalized for gene length and transformed into a frequency matrix by fragments per kilobase of transcript per million fragments mapped normalization. Gene count was computed as the number of genes present (abundance strictly positive) in the frequency matrix.

## Profiling and annotation of MSPs and gut enterotypes

IGC2 was previously organized into 1,990 MSPs with MSPminer⁶³ using a publicly available updated MSP dataset⁶⁴. Relative abundance of MSP was computed as the mean abundance of its 100 'marker' genes (that is, the genes that correlate the most altogether). If less than 10% of 'marker' genes were seen in a sample, the abundance of the MSP was set to 0. This approach was used in the MetaHIT⁶² and Metacardis⁶⁵ consortia. For the 4 MSPs with less than 100 core genes, all available core genes were used.

Abundances at higher taxonomical ranks were computed as the sum of the MSP that belong to a given taxa. MSP count was assessed as the number of MSPs present in a sample (that is, whose abundance is strictly positive). Enterotypes profiling was performed as previously demonstrated⁶⁶.

## Estimating functional modules of gut bacteriome

Genes from the IGC2 catalogue were mapped with diamond⁶⁷ onto KEGG orthologues (KO) from the KEGG database⁶⁸ (v8.9). Each gene was

assigned to the best-ranked KO among hits with e-value < 10 × 10⁻⁵ and bit score >60. Then we assessed presence and abundance of GMMs²⁸ and GBMs²⁹ in a metagenomic sample by the pipeline implemented in the R package omixerRpm (v0.3.2) as previously described²⁸,²⁹.

## Estimation of dynamic growth rate of bacteria from metagenomic samples

We used the computational pipeline to infer gut bacterial growth dynamics from metagenomic samples as previously described²⁷. Sequencing reads were mapped to a database that contains complete genomic references of 2,991 strains belonging to 1,509 microbial species. For each bacterial species, a reference strain with prevalence of 100% across the samples was selected. A coverage map was then assembled on the basis of aligned reads to the reference genome. Genomic segments were binned into 10-kbp regions and coverage of the resulting bins was calculated and smoothed. The location of the origin and terminus of replication was predicted by fits of the same strain across multiple samples. Lastly, PTRs were calculated for each bacterial species in every sample as the smoothed sequencing coverage of the representative strain at the predicted peak location, divided by that at the predicted trough location.

## Studies of bacterial structural variations

Before the SVs classification, the pipeline with iterative coverage-based read assignment algorithm was performed to reassign the ambiguous reads to the most likely reference with high accuracy³¹,³². The reference genomes provided in the proGenomes database (http://progenomes1.embl.de/) were concatenated and then divided into genomic 1-kbp bins and applied for the detection of highly variable genomic segments. The SGV-Finder pipeline³¹ was used to detect the SVs that are either (1) with deletion percentage of the genomic segment across the population of <25% (variable SVs, vSVs), (2) with deletion percentage between 25% and 75% (deletion SVs, dSVs; the absence or presence of this particular genomic segment was kept) or (3) with deletion percentage of >75% (this genomic segment was excluded from the analysis). All bacterial species with SV calling were present in at least 10% of the total samples and were used for subsequent analysis.

## Mediation analysis

The R package 'mediation⁶⁹ (v4.5.0) was used to infer causal relationships between gut microbial features, polar and microbiota-related metabolites, and metabolic traits and eating disorder scores. To reduce the testing numbers, we only kept the candidate groups consisting of variables that were strongly associated with each other; that is, for a candidate microbial feature-metabolite-phenotypic variable causal group, the association between gut microbial feature and serum metabolite was significant ($P_{adj} < 0.1$); the association between metabolite and phenotypic variable was significant ($P_{adj} < 0.1$); and the association between microbial feature and phenotypic variable was significant ($P_{adj} < 0.1$). After performing mediation analysis, only candidate groups with significance in direction 1 were kept for Sankey diagram visualization.

## Profiling and analysis of viral gut microbiota

We profiled the viral gut microbiota using MiCoP⁷⁰, as this method is optimized to call viruses directly from bulk metagenomics sequencing reads and compute relative abundance within the virome dataset. As a reference dataset, MiCoP draws upon the NCBI's RefSeq Viral database⁷¹. We identified a total of 209 viral species with prevalence of > and =10% and relative abundance of > and =0.01% for 147 (77 AN versus 70 HC) individuals included in the dataset. Richness, alpha and beta diversity were calculated with the R package 'fossil'⁷² and 'vegan'⁷³. Two-tailed Wilcoxon's rank-sum test was used to determine statistically significant differences in richness and alpha diversity indices between groups. Permutational multivariate analysis of variance (PERMANOVA) at *n* = 999
was performed for Canberra distance. The viral-bacterial interactions in both AN-RS and AN-BP microbiome data were computed using the Sparse Correlations for Compositional (SparCC)⁷⁴ algorithm. Before the SparCC analysis, the AN bacterial and viral microbiota datasets were subset to AN-RS and AN-BP datasets, which were then separately submitted for SparCC analysis.

## Analysis of serum polar metabolites by gas chromatography-time-of-flight mass spectrometry

The metabolites listed as gut microbiota-related metabolites were based on literature mining⁷⁵,⁷⁶. Serum samples were randomized and sample preparation was done as described previously⁴,⁷⁷. Briefly, 400 µl of methanol (MeOH) containing internal standards (heptadecanoic acid, deuterium-labelled DL-valine, deuterium-labelled succinic acid and deuterium-labelled glutamic acid, 1 µg ml⁻¹) was added to 30 µl of the serum samples, which were then vortex mixed and incubated on ice for 30 min. Samples were then centrifuged (9,400 × g, 3 min) and 350 µl of the supernatant was collected after centrifugation. The solvent was evaporated to dryness, 25 µl of MOX reagent was added and the sample was incubated for 60 min at 45 °C. N-Methyl-N-(trimethylsilyl)trifluoroacetamide (25 µl) was added and after 60 min incubation at 45 °C, 25 µl of the retention index standard mixture (n-alkanes, 10 µg ml⁻¹) was added.

The analyses were done using an Agilent 7890B gas chromatograph coupled to an Agilent 7200 quadrupole time-of-flight mass spectrometer. The following parameters were used: injection volume was 1 µl with 100:1 split on PTV at 70 °C, heating to 300 °C at 120 °C min⁻¹; column: Zebron ZB-SemiVolatiles with length of 20 m, inner diameter of 0.18 mm, film thickness of 0.18 µm, with initial helium flow of 1.2 ml min⁻¹, increasing to 2.4 ml min⁻¹ after 16 min. Oven temperature programme: 50 °C (5 min), then to 270 °C at 20 °C min⁻¹ and then to 300 at 40 °C min⁻¹ (5 min). EI source: 250 °C, 70 eV electron energy, 35 µA emission, solvent delay 3 min. Mass range 55 to 650 amu, acquisition rate 5 spectra per second, acquisition time 200 ms per spectrum. Quad at 150 °C, 1.5 ml min⁻¹ N₂ collision flow, aux-2 temperature 280 °C.

Calibration curves were constructed using alanine, citric acid, fumaric acid, glutamic acid, glycine, lactic acid, malic acid, 2-hydroxybutyric acid, 3-hydroxybutyric acid, linoleic acid, oleic acid, palmitic acid, stearic acid, cholesterol, fructose, glutamine, indole-3-propionic acid, isoleucine, leucine, proline, succinic acid, valine, asparagine, aspartic acid, arachidonic acid, glycerol-3-phosphate, lysine, methionine, ornithine, phenylalanine, serine and threonine purchased from Sigma-Aldrich at a concentration range of 0.1–80 µg ml⁻¹. An aliquot of each sample was collected, pooled and used as quality-control samples, together with National Institute of Standards and Technology (NIST) CRM1950 serum sample, an in-house pooled serum sample. The relative standard deviation of the concentrations was on average 16% for the pooled quality-control samples and 10% for the NIST samples.

## Analysis of serum bile acids and serum semipolar metabolites

The sample preparation procedure was performed as described previously⁷⁸. The plate was preconditioned with 450 µl acetonitrile before the addition of 100 µl of sample and 10 µl of polyfluoroalkyl substances (PFAS) and bile acids (BAs) internal standard mixture (200 ng ml⁻¹ and 1,000 ng ml⁻¹, respectively). Thereafter, 450 µl of acetonitrile containing 1% formic acid were added to each well and the samples extracted using a 10" vacuum manifold. The eluate was evaporated to dryness under nitrogen gas flow and reconstituted to 80 µl of MeOH/2 mM aqueous ammonium ethanoate.

Chromatographic separation was carried out using an Acquity UPLC BEH C18 column (100 mm × 2.1 mm inner diameter, 1.7 µm particle size), fitted with a C18 precolumn (Waters). Mobile phase A consisted of H₂O:MeOH (v/v 70:30) and mobile phase B of MeOH with both phases containing 2 mM ammonium acetate as an ionization agent. The

flow rate was set at 0.4 ml min⁻¹ with the elution gradient as follows: 0–1.5 min, mobile phase B was increased from 5% to 30%; 1.5–4.5 min, mobile phase B was increased to 70%; 4.5–7.5 min, mobile phase B was increased to 100% and held for 5.5 min. A post-time of 5 min was used to regain the initial conditions for the next analysis. The total run time per sample was 18 min. The dual electrospray ionization source settings were as follows: capillary voltage was 4.5 kV, nozzle voltage 1,500 V, N₂ pressure in the nebulizer was 21 psi and the N₂ flow rate and temperature as sheath gas were 11 l min⁻¹ and 379 °C, respectively. To obtain accurate mass spectra in the mass spectrometry (MS) scan, the *m/z* range was set to 100–1,700 in negative ion mode. MassHunter B.06.01 software (Agilent) was used for all data acquisition.

Identification of compounds was done with an in-house spectral library using MS (and retention time), tandem mass spectrometry information. Quantitation was based on a matrix-matched calibration curve spiked with native compounds. The calibration curve consisted of concentrations ranging from 0–1,600 ng ml⁻¹ for BAs. The relative standard deviation for the BAs was on average 17.8% for the quality-control samples and 19.4% for the NIST samples.

## Animal experiments

Animal protocols were approved by the Science Ethics Committees of the Capital Region of Copenhagen, Denmark. Female germ-free Swiss Webster mice were bred and maintained in flexible film gnotobiotic isolators until the start of experiments at the Department of Experimental Medicine, University of Copenhagen. The mice were fed autoclaved chow diet (7% simple sugars, 3% fat, 50% polysaccharide, 15% protein (w/w), energy 3.5 kcal g⁻¹) and water ad libitum under a 12 h light/12 h dark cycle (lights on at 7:30 a.m.) and constant temperature (21–22 °C) and humidity (55 ± 5%).

Faecal samples from randomly selected subsets of three patients with AN (females aged 20, 22 and 20 yr with BMI of 10.3, 11.5 and 11.7 kg m⁻², respectively) and three HC (females aged 20, 22 and 21 yr with BMI of 22.6, 21.1 and 21.2 kg m⁻², respectively) who were representatives of cases and controls were used to colonize 6-week-old female GF littermates. Briefly, 250 mg of faecal samples were suspended with 5 ml of LYBHI media (supplemented with 0.05% Cysteine and 0.2% Hemin as reducing agents) diluted in 20% glycerol (20 ml g⁻¹ of faeces) in an anaerobic cryovial; these inoculum samples were then vortexed for 5 min, followed by 5 min standing to precipitate particles. The faecal slurries were then transferred to 1 ml cryovials and immediately frozen at –80 °C. At day one, both groups of mice were housed in autoclaved individually ventilated cages where they received the first dose (200 µl) of faecal slurries. The mice were then given autoclaved chow diet and water ad libitum for 2 d and their food intake was recorded. At day 3, mice were gavaged with a second dose of faecal material from the same matched AN and HC donors as before. Thereafter, mice in both groups were singly housed and subjected to 30% energy-restricted autoclaved chow diet (amount of given food was set at 70% of ad libitum food intake for each mouse) for 3 weeks where water was given ad libitum. Both the anorexia-transplanted (AN-T) and the normal control-transplanted (HC-T) mice were weighed every 5 d after the start of energy-restricted dieting.

At the end of the study, mice were anaesthetized with isoflurane and blood from the vena cava was collected in tubes containing EDTA. Blood samples were centrifuged for 6 min at 4,032 g at 4 °C. Plasma was isolated and stored at –80 °C for subsequent biochemical testing. Inguinal subcutaneous white adipose tissue, the caecum and the hypothalamus of each mouse were precisely dissected and collected for quantitative PCR analysis. No animal or data points were excluded from the present study.

Total RNA was extracted from tissues using Trizol reagent (Invitrogen) according to the manufacturer's instructions, followed by concentration measurement. One µg of RNA was transcribed to complementary DNA using the Reverse Transcription System (Promega).
was performed for Canberra distance. The viral-bacterial interactions in both AN-RS and AN-BP microbiome data were computed using the Sparse Correlations for Compositional (SparCC)⁷⁴ algorithm. Before the SparCC analysis, the AN bacterial and viral microbiota datasets were subset to AN-RS and AN-BP datasets, which were then separately submitted for SparCC analysis.

## Analysis of serum polar metabolites by gas chromatography-time-of-flight mass spectrometry

The metabolites listed as gut microbiota-related metabolites were based on literature mining⁷⁵,⁷⁶. Serum samples were randomized and sample preparation was done as described previously⁴,⁷⁷. Briefly, 400 µl of methanol (MeOH) containing internal standards (heptadecanoic acid, deuterium-labelled DL-valine, deuterium-labelled succinic acid and deuterium-labelled glutamic acid, 1 µg ml⁻¹) was added to 30 µl of the serum samples, which were then vortex mixed and incubated on ice for 30 min. Samples were then centrifuged (9,400 × g, 3 min) and 350 µl of the supernatant was collected after centrifugation. The solvent was evaporated to dryness, 25 µl of MOX reagent was added and the sample was incubated for 60 min at 45 °C. N-Methyl-N-(trimethylsilyl)trifluoroacetamide (25 µl) was added and after 60 min incubation at 45 °C, 25 µl of the retention index standard mixture (n-alkanes, 10 µg ml⁻¹) was added.

The analyses were done using an Agilent 7890B gas chromatograph coupled to an Agilent 7200 quadrupole time-of-flight mass spectrometer. The following parameters were used: injection volume was 1 µl with 100:1 split on PTV at 70 °C, heating to 300 °C at 120 °C min⁻¹; column: Zebron ZB-SemiVolatiles with length of 20 m, inner diameter of 0.18 mm, film thickness of 0.18 µm, with initial helium flow of 1.2 ml min⁻¹, increasing to 2.4 ml min⁻¹ after 16 min. Oven temperature programme: 50 °C (5 min), then to 270 °C at 20 °C min⁻¹ and then to 300 at 40 °C min⁻¹ (5 min). EI source: 250 °C, 70 eV electron energy, 35 µA emission, solvent delay 3 min. Mass range 55 to 650 amu, acquisition rate 5 spectra per second, acquisition time 200 ms per spectrum. Quad at 150 °C, 1.5 ml min⁻¹ N₂ collision flow, aux-2 temperature 280 °C.

Calibration curves were constructed using alanine, citric acid, fumaric acid, glutamic acid, glycine, lactic acid, malic acid, 2-hydroxybutyric acid, 3-hydroxybutyric acid, linoleic acid, oleic acid, palmitic acid, stearic acid, cholesterol, fructose, glutamine, indole-3-propionic acid, isoleucine, leucine, proline, succinic acid, valine, asparagine, aspartic acid, arachidonic acid, glycerol-3-phosphate, lysine, methionine, ornithine, phenylalanine, serine and threonine purchased from Sigma-Aldrich at a concentration range of 0.1–80 µg ml⁻¹. An aliquot of each sample was collected, pooled and used as quality-control samples, together with National Institute of Standards and Technology (NIST) CRM1950 serum sample, an in-house pooled serum sample. The relative standard deviation of the concentrations was on average 16% for the pooled quality-control samples and 10% for the NIST samples.

## Analysis of serum bile acids and serum semipolar metabolites

The sample preparation procedure was performed as described previously⁷⁸. The plate was preconditioned with 450 µl acetonitrile before the addition of 100 µl of sample and 10 µl of polyfluoroalkyl substances (PFAS) and bile acids (BAs) internal standard mixture (200 ng ml⁻¹ and 1,000 ng ml⁻¹, respectively). Thereafter, 450 µl of acetonitrile containing 1% formic acid were added to each well and the samples extracted using a 10" vacuum manifold. The eluate was evaporated to dryness under nitrogen gas flow and reconstituted to 80 µl of MeOH/2 mM aqueous ammonium ethanoate.

Chromatographic separation was carried out using an Acquity UPLC BEH C18 column (100 mm × 2.1 mm inner diameter, 1.7 µm particle size), fitted with a C18 precolumn (Waters). Mobile phase A consisted of H₂O:MeOH (v/v 70:30) and mobile phase B of MeOH with both phases containing 2 mM ammonium acetate as an ionization agent. The

flow rate was set at 0.4 ml min⁻¹ with the elution gradient as follows: 0–1.5 min, mobile phase B was increased from 5% to 30%; 1.5–4.5 min, mobile phase B was increased to 70%; 4.5–7.5 min, mobile phase B was increased to 100% and held for 5.5 min. A post-time of 5 min was used to regain the initial conditions for the next analysis. The total run time per sample was 18 min. The dual electrospray ionization source settings were as follows: capillary voltage was 4.5 kV, nozzle voltage 1,500 V, N₂ pressure in the nebulizer was 21 psi and the N₂ flow rate and temperature as sheath gas were 11 l min⁻¹ and 379 °C, respectively. To obtain accurate mass spectra in the mass spectrometry (MS) scan, the *m/z* range was set to 100–1,700 in negative ion mode. MassHunter B.06.01 software (Agilent) was used for all data acquisition.

Identification of compounds was done with an in-house spectral library using MS (and retention time), tandem mass spectrometry information. Quantitation was based on a matrix-matched calibration curve spiked with native compounds. The calibration curve consisted of concentrations ranging from 0–1,600 ng ml⁻¹ for BAs. The relative standard deviation for the BAs was on average 17.8% for the quality-control samples and 19.4% for the NIST samples.

## Animal experiments

Animal protocols were approved by the Science Ethics Committees of the Capital Region of Copenhagen, Denmark. Female germ-free Swiss Webster mice were bred and maintained in flexible film gnotobiotic isolators until the start of experiments at the Department of Experimental Medicine, University of Copenhagen. The mice were fed autoclaved chow diet (7% simple sugars, 3% fat, 50% polysaccharide, 15% protein (w/w), energy 3.5 kcal g⁻¹) and water ad libitum under a 12 h light/12 h dark cycle (lights on at 7:30 a.m.) and constant temperature (21–22 °C) and humidity (55 ± 5%).

Faecal samples from randomly selected subsets of three patients with AN (females aged 20, 22 and 20 yr with BMI of 10.3, 11.5 and 11.7 kg m⁻², respectively) and three HC (females aged 20, 22 and 21 yr with BMI of 22.6, 21.1 and 21.2 kg m⁻², respectively) who were representatives of cases and controls were used to colonize 6-week-old female GF littermates. Briefly, 250 mg of faecal samples were suspended with 5 ml of LYBHI media (supplemented with 0.05% Cysteine and 0.2% Hemin as reducing agents) diluted in 20% glycerol (20 ml g⁻¹ of faeces) in an anaerobic cryovial; these inoculum samples were then vortexed for 5 min, followed by 5 min standing to precipitate particles. The faecal slurries were then transferred to 1 ml cryovials and immediately frozen at –80 °C. At day one, both groups of mice were housed in autoclaved individually ventilated cages where they received the first dose (200 µl) of faecal slurries. The mice were then given autoclaved chow diet and water ad libitum for 2 d and their food intake was recorded. At day 3, mice were gavaged with a second dose of faecal material from the same matched AN and HC donors as before. Thereafter, mice in both groups were singly housed and subjected to 30% energy-restricted autoclaved chow diet (amount of given food was set at 70% of ad libitum food intake for each mouse) for 3 weeks where water was given ad libitum. Both the anorexia-transplanted (AN-T) and the normal control-transplanted (HC-T) mice were weighed every 5 d after the start of energy-restricted dieting.

At the end of the study, mice were anaesthetized with isoflurane and blood from the vena cava was collected in tubes containing EDTA. Blood samples were centrifuged for 6 min at 4,032 g at 4 °C. Plasma was isolated and stored at –80 °C for subsequent biochemical testing. Inguinal subcutaneous white adipose tissue, the caecum and the hypothalamus of each mouse were precisely dissected and collected for quantitative PCR analysis. No animal or data points were excluded from the present study.

Total RNA was extracted from tissues using Trizol reagent (Invitrogen) according to the manufacturer's instructions, followed by concentration measurement. One µg of RNA was transcribed to complementary DNA using the Reverse Transcription System (Promega).
Real-time PCR was performed using the LC480 detection system (Roche Diagnostics) and SYBR Green I Supermix (Takara). All qPCRs were run on the thermal cycles at 95 °C for 10 min, followed by 45 cycles of 0.01 s at 95 °C and of 20 s at 60 °C. Data were normalized to the housekeeping Rpl36 gene for adipose tissue and Rplp0 gene⁷⁹ for hypothalamus and analysed according to the delta-delta CT method. Sequences of oligonucleotides used in this study are provided in Supplementary Table 8.

## DNA extraction, 16S rRNA sequencing and data analyses in mice experiments

Microbial DNAs were isolated and purified from stool samples (-250 mg) of human donors and mouse recipients by using NucleoSpin soil mini kit (MACHEREY-NAGEL). The DNA was then amplified using Phusion High-Fidelity PCR master mix (New England Biolabs) by PCR targeting the V3-V4 region of the 16S rRNA gene (primer sequences provided in Supplementary Table 8). The following PCR programme was used: 98 °C for 30 s, 25× (98 °C for 10 s, 55 °C for 20 s, 72 °C for 20 s), 72 °C for 5 min. Amplification was verified by running the products on an agarose gel. Indices were added in a subsequent PCR using an Illumina Nextera kit with the following PCR programme: 98 °C for 30 s, 8× (98 °C for 10 s, 55 °C for 20 s, 72 °C for 20 s), 72 °C for 5 min. Attachment of indices was verified by running the products on an agarose gel. Products from the nested PCR were pooled on the basis of band intensity and the resulting library cleaned with magnetic beads. The DNA concentration of pooled libraries was measured fluorometrically. Sequencing was done on an Illumina MiSeq desktop sequencer using the MiSeq reagent kit V3 (Illumina) for 2 × 300 bp paired-end sequencing. Paired-end reads were subsequently trimmed, merged and analysed using the DADA2 (v1.16.0) pipeline⁸⁰.

## Statistical analysis

No data were excluded before the statistical analysis in the present study. No allocation and randomization were included as the study is observational. This study includes all available samples (n = 147) of patients with anorexia and healthy individuals. Although no statistical methods were used to predetermine sample sizes, our sample sizes are similar to those reported in previous publications⁴³,⁸¹. Samples were randomly distributed across metagenomics and metabolomics batches. Investigators were blinded to group allocation during data collection in metagenomic, biochemical and metabolomics analyses. All analyses of human samples were performed using R (v4.1.2). Gene expression analyses and body weight comparisons for animal studies were performed using GraphPad Prism (v9.3.0).

### (1) Differential analysis

We carried out the differential analysis using the metadecounfoundR pipeline implemented in the R package metadecounfoundR (v0.1.8; see https://github.com/TillBirkner/metadecounfoundR or https://doi.org/10.5281/zenodo.4721078) where we assessed the extent to which the observed differences between AN and HC participants in microbiome or metabolome analyses are confounded by covariates including age, BMI, smoking and medication. This pipeline initially used univariate statistics to find associations between microbiome features and disease status, followed by nested linear model comparison post hoc testing to check for the confounding effects of potential covariates and finally, returning a status label.

### (2) Association analysis

For the association and mediation analysis between omics features and eating behaviour and psychological traits within the AN group, we first checked the normality of continuous variables with Shapiro-Wilk normality test, finding most variables not to be normally distributed. Therefore, before association analysis, we standardized the continuous variables using empirical normal quantile transformation to follow a standard normal

distribution (N - (0, 1)). Then we implemented a linear regression model to assess the associations between omics features and eating behaviour and psychological traits using the following formula where confounding factors were added as covariates.

$$ 
\begin{array}{l}
\text{Eating behavior – and psychological traits ~ omics features} \\
\text{(for example, metabolite/msp/SV)} \\
\text{+Age + BMI + Smoking + Medication}
\end{array}
 $$

Medication included selective serotonin re-uptake inhibitors, antipsychotics and benzodiazepines.

### (3) For the association analysis between omics features and host metabolic traits, a normality check and data standardization were also performed before the linear regression analysis. In the linear regression model, the above-mentioned confounding factors were included as covariates, except for BMI as the extremely low BMI is the most remarkable phenotypic change for AN patients compared with HC individuals.

$$ 
\begin{array}{l}
\text{Metabolic traits(for example, BMI/plasma glucose) ~ omics features} \\
\text{(for example, metabolite/msp/SV) + Age + Smoking + Medication}
\end{array}
 $$

Differences in gut microbial diversity (gene richness, species count, taxonomic composition) between AN and HC were calculated using Wilcoxon tests, and Kruskal-Wallis test was used for assessing the significance of differences between multiple groups. Unless otherwise stated, all P values were corrected using the Benjamini-Hochberg method and P$_{adj}$ < 0.05 was considered statistically significant.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

Anonymized clinical data that are stored in Sharepoint via Odense Patient Data Explorative (file no OP_153) can be accessed by contacting rene.stoeving@rsyd.dk or can be found in Supplementary Table 2. Raw shotgun sequencing data and 16s rRNA gene amplicon sequencing data that support the findings of this study have been deposited in the European Nucleotide Archive with accession numbers PRJEB51776 and PRJEB60103, respectively. Metabolomics data are available from Metabolomics Workbench repository under the link https://doi.org/10.21228/M8KT5B. The KEGG Database is available at https://www.genome.jp/kegg/. Source data are provided with this paper.

## Code availability

Codes associated with the data analysis and visualization are available at https://github.com/fjw536/AnorexiaGutMicrobiome.git.

## References

1. American Psychiatric Association *Diagnostic and Statistical Manual of Mental Disorders* 5th edn (American Psychiatric Publishing, Inc., 2013).

2. Smink, F. R., Van Hoeken, D. & Hoek, H. W. Epidemiology of eating disorders: incidence, prevalence and mortality rates. *Curr. Psychiatry Rep.* **14**, 406–414 (2012).

3. Bulik, C. M. The challenges of treating anorexia nervosa. *Lancet* **383**, 105–106 (2014).

4. Winkler, L. A.-D., Bilenberg, N., Hørder, K. & Støving, R. K. Does specialization of treatment influence mortality in eating disorders?—A comparison of two retrospective cohorts. *Psychiatry Res.* **230**, 165–171 (2015).
Real-time PCR was performed using the LC480 detection system (Roche Diagnostics) and SYBR Green I Supermix (Takara). All qPCRs were run on the thermal cycles at 95 °C for 10 min, followed by 45 cycles of 0.01 s at 95 °C and of 20 s at 60 °C. Data were normalized to the housekeeping Rpl36 gene for adipose tissue and Rplp0 gene⁷⁹ for hypothalamus and analysed according to the delta-delta CT method. Sequences of oligonucleotides used in this study are provided in Supplementary Table 8.

## DNA extraction, 16S rRNA sequencing and data analyses in mice experiments

Microbial DNAs were isolated and purified from stool samples (-250 mg) of human donors and mouse recipients by using NucleoSpin soil mini kit (MACHEREY-NAGEL). The DNA was then amplified using Phusion High-Fidelity PCR master mix (New England Biolabs) by PCR targeting the V3-V4 region of the 16S rRNA gene (primer sequences provided in Supplementary Table 8). The following PCR programme was used: 98 °C for 30 s, 25× (98 °C for 10 s, 55 °C for 20 s, 72 °C for 20 s), 72 °C for 5 min. Amplification was verified by running the products on an agarose gel. Indices were added in a subsequent PCR using an Illumina Nextera kit with the following PCR programme: 98 °C for 30 s, 8× (98 °C for 10 s, 55 °C for 20 s, 72 °C for 20 s), 72 °C for 5 min. Attachment of indices was verified by running the products on an agarose gel. Products from the nested PCR were pooled on the basis of band intensity and the resulting library cleaned with magnetic beads. The DNA concentration of pooled libraries was measured fluorometrically. Sequencing was done on an Illumina MiSeq desktop sequencer using the MiSeq reagent kit V3 (Illumina) for 2 × 300 bp paired-end sequencing. Paired-end reads were subsequently trimmed, merged and analysed using the DADA2 (v1.16.0) pipeline⁸⁰.

## Statistical analysis

No data were excluded before the statistical analysis in the present study. No allocation and randomization were included as the study is observational. This study includes all available samples (n = 147) of patients with anorexia and healthy individuals. Although no statistical methods were used to predetermine sample sizes, our sample sizes are similar to those reported in previous publications⁴³,⁸¹. Samples were randomly distributed across metagenomics and metabolomics batches. Investigators were blinded to group allocation during data collection in metagenomic, biochemical and metabolomics analyses. All analyses of human samples were performed using R (v4.1.2). Gene expression analyses and body weight comparisons for animal studies were performed using GraphPad Prism (v9.3.0).

### (1) Differential analysis

We carried out the differential analysis using the metadecounfoundR pipeline implemented in the R package metadecounfoundR (v0.1.8; see https://github.com/TillBirkner/metadecounfoundR or https://doi.org/10.5281/zenodo.4721078) where we assessed the extent to which the observed differences between AN and HC participants in microbiome or metabolome analyses are confounded by covariates including age, BMI, smoking and medication. This pipeline initially used univariate statistics to find associations between microbiome features and disease status, followed by nested linear model comparison post hoc testing to check for the confounding effects of potential covariates and finally, returning a status label.

### (2) Association analysis

For the association and mediation analysis between omics features and eating behaviour and psychological traits within the AN group, we first checked the normality of continuous variables with Shapiro-Wilk normality test, finding most variables not to be normally distributed. Therefore, before association analysis, we standardized the continuous variables using empirical normal quantile transformation to follow a standard normal

distribution (N - (0, 1)). Then we implemented a linear regression model to assess the associations between omics features and eating behaviour and psychological traits using the following formula where confounding factors were added as covariates.

$$ 
\begin{array}{l}
\text{Eating behavior – and psychological traits ~ omics features} \\
\text{(for example, metabolite/msp/SV)} \\
\text{+Age + BMI + Smoking + Medication}
\end{array}
 $$

Medication included selective serotonin re-uptake inhibitors, antipsychotics and benzodiazepines.

### (3) For the association analysis between omics features and host metabolic traits, a normality check and data standardization were also performed before the linear regression analysis. In the linear regression model, the above-mentioned confounding factors were included as covariates, except for BMI as the extremely low BMI is the most remarkable phenotypic change for AN patients compared with HC individuals.

$$ 
\begin{array}{l}
\text{Metabolic traits(for example, BMI/plasma glucose) ~ omics features} \\
\text{(for example, metabolite/msp/SV) + Age + Smoking + Medication}
\end{array}
 $$

Differences in gut microbial diversity (gene richness, species count, taxonomic composition) between AN and HC were calculated using Wilcoxon tests, and Kruskal-Wallis test was used for assessing the significance of differences between multiple groups. Unless otherwise stated, all *P* values were corrected using the Benjamini-Hochberg method and *P*<sub>adj</sub> < 0.05 was considered statistically significant.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

Anonymized clinical data that are stored in Sharepoint via Odense Patient Data Explorative (file no OP_153) can be accessed by contacting rene.stoeving@rsyd.dk or can be found in Supplementary Table 2. Raw shotgun sequencing data and 16s rRNA gene amplicon sequencing data that support the findings of this study have been deposited in the European Nucleotide Archive with accession numbers PRJEB51776 and PRJEB60103, respectively. Metabolomics data are available from Metabolomics Workbench repository under the link https://doi.org/10.21228/M8KT5B. The KEGG Database is available at https://www.genome.jp/kegg/. Source data are provided with this paper.

## Code availability

Codes associated with the data analysis and visualization are available at https://github.com/fjw536/AnorexiaGutMicrobiome.git.

## References

1. American Psychiatric Association *Diagnostic and Statistical Manual of Mental Disorders* 5th edn (American Psychiatric Publishing, Inc., 2013).

2. Smink, F. R., Van Hoeken, D. & Hoek, H. W. Epidemiology of eating disorders: incidence, prevalence and mortality rates. *Curr. Psychiatry Rep.* **14**, 406–414 (2012).

3. Bulik, C. M. The challenges of treating anorexia nervosa. *Lancet* **383**, 105–106 (2014).

4. Winkler, L. A.-D., Bilenberg, N., Hørder, K. & Støving, R. K. Does specialization of treatment influence mortality in eating disorders?—A comparison of two retrospective cohorts. *Psychiatry Res.* **230**, 165–171 (2015).
5. Støving, R. K. et al. Purging behavior in anorexia nervosa and eating disorder not otherwise specified: a retrospective cohort study. *Psychiatry Res.* **198**, 253–258 (2012).

6. Sullivan, P. F. Mortality in anorexia nervosa. *Am. J. Psychiatry* **152**, 1073–1074 (1995).

7. Trace, S. E., Baker, J. H., Peñas-Lledó, E. & Bulik, C. M. The genetics of eating disorders. *Annu. Rev. Clin. Psychol.* **9**, 589–620 (2013).

8. Paolacci, S. et al. Genetic contributions to the etiology of anorexia nervosa: new perspectives in molecular diagnosis and treatment. *Mol. Genet. Genom. Med.* **8**, e1244 (2020).

9. Watson, H. J. et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. *Nat. Genet.* **51**, 1207–1214 (2019).

10. Støving, R. K. Mechanisms in endocrinology: anorexia nervosa and endocrinology: a clinical update. *Eur. J. Endocrinol.* **180**, R9–R27 (2019).

11. Bailer, U. F. et al. Interaction between serotonin transporter and dopamine D2/D3 receptor radioligand measures is associated with harm avoidant symptoms in anorexia and bulimia nervosa. *Psychiatry Res.* **211**, 160–168 (2013).

12. Fan, Y. & Pedersen, O. Gut microbiota in human metabolic health and disease. *Nat. Rev. Microbiol.* **19**, 55–71 (2021).

13. Morais, L. H., Schreiber, H. L. & Mazmanian, S. K. The gut microbiota-brain axis in behaviour and brain disorders. *Nat. Rev. Microbiol.* **19**, 241–255 (2021).

14. Tennoune, N. et al. Bacterial ClpB heat-shock protein, an antigen-mimetic of the anorexigenic peptide α-MSH, at the origin of eating disorders. *Transl. Psychiatry* **4**, e458 (2014).

15. Sobrino Crespo, C., Perianes Cachero, A., Puebla Jimenez, L., Barrios, V. & Arilla Ferreiro, E. Peptides and food intake. *Front. Endocrinol.* **5**, 58 (2014).

16. Morita, C. et al. Gut dysbiosis in patients with anorexia nervosa. *PLoS ONE* **10**, e0145274 (2015).

17. Kleiman, S. C. et al. The intestinal microbiota in acute anorexia nervosa and during renourishment: relationship to depression, anxiety, and eating disorder psychopathology. *Psychosom. Med.* **77**, 969 (2015).

18. Armougom, F., Henry, M., Vialettes, B., Raccah, D. & Raoult, D. Monitoring bacterial community of human gut microbiota reveals an increase in *Lactobacillus* in obese patients and Methanogens in anorexic patients. *PLoS ONE* **4**, e7125 (2009).

19. Mack, I. et al. Weight gain in anorexia nervosa does not ameliorate the faecal microbiota, branched chain fatty acid profiles, and gastrointestinal complaints. *Sci. Rep.* **6**, 26752 (2016).

20. Breton, J. et al. Gut microbiota alteration in a mouse model of anorexia nervosa. *Clin. Nutr.* **40**, 181–189 (2021).

21. Clausen, L., Rosenvinge, J. H., Friborg, O. & Rokkedal, K. Validating the Eating Disorder Inventory-3 (EDI-3): a comparison between 561 female eating disorders patients and 878 females from the general population. *J. Psychopathol. Behav. Assess.* **33**, 101–110 (2011).

22. Costea, P. I. et al. Enterotypes in the landscape of gut microbial community composition. *Nat. Microbiol.* **3**, 8–16 (2018).

23. La Rosa, S. L. et al. The human gut Firmicute Roseburia intestinalis is a primary degrader of dietary β-mannans. *Nat. Commun.* **10**, 905 (2019).

24. Roager, H. M. & Licht, T. R. Microbial tryptophan catabolites in health and disease. *Nat. Commun.* **9**, 3294 (2018).

25. Smith, C. J. et al. Probiotics normalize the gut-brain-microbiota axis in immunodeficient mice. *Am. J. Physiol. Gastrointest. Liver Physiol.* **307**, G793–G802 (2014).

26. Hao, Z., Wang, W., Guo, R. & Liu, H. Faecalibacterium prausnitzii (ATCC 27766) has preventive and therapeutic effects on chronic unpredictable mild stress-induced depression-like and anxiety-like behavior in rats. *Psychoneuroendocrinology* **104**, 132–142 (2019).

27. Korem, T. et al. Growth dynamics of gut microbiota in health and disease inferred from single metagenomic samples. *Science* **349**, 1101–1106 (2015).

28. Vieira-Silva, S. et al. Species-function relationships shape ecological properties of the human gut microbiome. *Nat. Microbiol.* **1**, 16088 (2016).

29. Valles-Colomer, M. et al. The neuroactive potential of the human gut microbiota in quality of life and depression. *Nat. Microbiol.* **4**, 623–632 (2019).

30. Maresz, K. Growing evidence of a proven mechanism shows vitamin K2 Can impact health conditions beyond bone and cardiovascular. *Integr. Med.* **20**, 34–38 (2021).

31. Zeevi, D. et al. Structural variation in the gut microbiome associates with host health. *Nature* **568**, 43–48 (2019).

32. Wang, D. et al. Characterization of gut microbial structural variations as determinants of human bile acid metabolism. *Cell Host Microbe* **29**, 1802–1814.e5 (2021).

33. Dhir, S., Tarasenko, M., Napoli, E. & Giulivi, C. Neurological, psychiatric, and biochemical aspects of thiamine deficiency in children and adults. *Front. Psychiatry* **4**, 207 (2019).

34. Winston, A., Jamieson, C., Madira, W., Gatward, N. & Palmer, R. L. Prevalence of thiamin deficiency in anorexia nervosa. *Int. J. Eat. Disord.* **28**, 451–454 (2000).

35. Perino, A. et al. Central anorexigenic actions of bile acids are mediated by TGR5. *Nat. Metab.* **3**, 595–603 (2021).

36. Germain, N. et al. Constitutional thinness and lean anorexia nervosa display opposite concentrations of peptide YY, glucagon-like peptide 1, ghrelin, and leptin. *Am. J. Clin. Nutr.* **85**, 967–971 (2007).

37. Kamal, N. et al. Delayed gastrointestinal transit times in anorexia nervosa and bulimia nervosa. *Gastroenterology* **101**, 1320–1324 (1991).

38. Sun, C.-Y. et al. Indoxyl sulfate caused behavioral abnormality and neurodegeneration in mice with unilateral nephrectomy. *Aging* **13**, 6681 (2021).

39. Brydges, C. R. et al. Indoxyl sulfate, a gut microbiome-derived uremic toxin, is associated with psychic anxiety and its functional magnetic resonance imaging-based neurologic signature. *Sci. Rep.* **11**, 21011 (2021).

40. Lee, J.-H. & Lee, J. Indole as an intercellular signal in microbial communities. *FEMS Microbiol. Rev.* **34**, 426–444 (2010).

41. Newton, W. A. & Snell, E. E. Formation and interrelationships of tryptophanase and tryptophan synthetases in *Escherichia coli*. *J. Bacteriol.* **89**, 355–364 (1965).

42. Yanchuk, P. et al. Role of serotonin in the regulation of respiration and bile secretory function of the liver. *Fiziol. Zh.* **61**, 102–110 (2015).

43. Pedersen, H. K. et al. Human gut microbes impact host serum metabolome and insulin sensitivity. *Nature* **535**, 376–381 (2016).

44. Glenny, E. M. et al. Gut microbial communities from patients with anorexia nervosa do not influence body weight in recipient germ-free mice. *Gut Microbes* **13**, 1897216 (2021).

45. Smith, M. I. et al. Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. *Science* **339**, 548–554 (2013).

46. Unger, T. J., Calderon, G. A., Bradley, L. C., Sena-Esteves, M. & Rios, M. Selective deletion of Bdnf in the ventromedial and dorsomedial hypothalamus of adult mice results in hyperphagic behavior and obesity. *J. Neurosci.* **27**, 14265–14274 (2007).

47. Yang, S.-C., Shieh, K.-R. & Li, H.-Y. Cocaine-and amphetamine-regulated transcript in the nucleus accumbens participates in the regulation of feeding behavior in rats. *Neuroscience* **133**, 841–851 (2005).

48. Hashimoto, M. & Masliah, E. Alpha-synuclein in Lewy body disease and Alzheimer's disease. *Brain Pathol.* **9**, 707–720 (1999).

49. Konopelski, P. et al. Indole-3-propionic acid, a tryptophan-derived bacterial metabolite, increases blood pressure via cardiac and vascular mechanisms in rats. *Am. J. Physiol. Regul.* **321**, R969–R981 (2021).
5. Støving, R. K. et al. Purging behavior in anorexia nervosa and eating disorder not otherwise specified: a retrospective cohort study. *Psychiatry Res.* **198**, 253–258 (2012).

6. Sullivan, P. F. Mortality in anorexia nervosa. *Am. J. Psychiatry* **152**, 1073–1074 (1995).

7. Trace, S. E., Baker, J. H., Peñas-Lledó, E. & Bulik, C. M. The genetics of eating disorders. *Annu. Rev. Clin. Psychol.* **9**, 589–620 (2013).

8. Paolacci, S. et al. Genetic contributions to the etiology of anorexia nervosa: new perspectives in molecular diagnosis and treatment. *Mol. Genet. Genom. Med.* **8**, e1244 (2020).

9. Watson, H. J. et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. *Nat. Genet.* **51**, 1207–1214 (2019).

10. Støving, R. K. Mechanisms in endocrinology: anorexia nervosa and endocrinology: a clinical update. *Eur. J. Endocrinol.* **180**, R9–R27 (2019).

11. Bailer, U. F. et al. Interaction between serotonin transporter and dopamine D2/D3 receptor radioligand measures is associated with harm avoidant symptoms in anorexia and bulimia nervosa. *Psychiatry Res.* **211**, 160–168 (2013).

12. Fan, Y. & Pedersen, O. Gut microbiota in human metabolic health and disease. *Nat. Rev. Microbiol.* **19**, 55–71 (2021).

13. Morais, L. H., Schreiber, H. L. & Mazmanian, S. K. The gut microbiota-brain axis in behaviour and brain disorders. *Nat. Rev. Microbiol.* **19**, 241–255 (2021).

14. Tennoune, N. et al. Bacterial ClpB heat-shock protein, an antigen-mimetic of the anorexigenic peptide α-MSH, at the origin of eating disorders. *Transl. Psychiatry* **4**, e458 (2014).

15. Sobrino Crespo, C., Perianes Cachero, A., Puebla Jimenez, L., Barrios, V. & Arilla Ferreiro, E. Peptides and food intake. *Front. Endocrinol.* **5**, 58 (2014).

16. Morita, C. et al. Gut dysbiosis in patients with anorexia nervosa. *PLoS ONE* **10**, e0145274 (2015).

17. Kleiman, S. C. et al. The intestinal microbiota in acute anorexia nervosa and during renourishment: relationship to depression, anxiety, and eating disorder psychopathology. *Psychosom. Med.* **77**, 969 (2015).

18. Armougom, F., Henry, M., Vialettes, B., Raccah, D. & Raoult, D. Monitoring bacterial community of human gut microbiota reveals an increase in *Lactobacillus* in obese patients and Methanogens in anorexic patients. *PLoS ONE* **4**, e7125 (2009).

19. Mack, I. et al. Weight gain in anorexia nervosa does not ameliorate the faecal microbiota, branched chain fatty acid profiles, and gastrointestinal complaints. *Sci. Rep.* **6**, 26752 (2016).

20. Breton, J. et al. Gut microbiota alteration in a mouse model of anorexia nervosa. *Clin. Nutr.* **40**, 181–189 (2021).

21. Clausen, L., Rosenvinge, J. H., Friborg, O. & Rokkedal, K. Validating the Eating Disorder Inventory-3 (EDI-3): a comparison between 561 female eating disorders patients and 878 females from the general population. *J. Psychopathol. Behav. Assess.* **33**, 101–110 (2011).

22. Costea, P. I. et al. Enterotypes in the landscape of gut microbial community composition. *Nat. Microbiol.* **3**, 8–16 (2018).

23. La Rosa, S. L. et al. The human gut Firmicute Roseburia intestinalis is a primary degrader of dietary β-mannans. *Nat. Commun.* **10**, 905 (2019).

24. Roager, H. M. & Licht, T. R. Microbial tryptophan catabolites in health and disease. *Nat. Commun.* **9**, 3294 (2018).

25. Smith, C. J. et al. Probiotics normalize the gut-brain-microbiota axis in immunodeficient mice. *Am. J. Physiol. Gastrointest. Liver Physiol.* **307**, G793–G802 (2014).

26. Hao, Z., Wang, W., Guo, R. & Liu, H. Faecalibacterium prausnitzii (ATCC 27766) has preventive and therapeutic effects on chronic unpredictable mild stress-induced depression-like and anxiety-like behavior in rats. *Psychoneuroendocrinology* **104**, 132–142 (2019).

27. Korem, T. et al. Growth dynamics of gut microbiota in health and disease inferred from single metagenomic samples. *Science* **349**, 1101–1106 (2015).

28. Vieira-Silva, S. et al. Species-function relationships shape ecological properties of the human gut microbiome. *Nat. Microbiol.* **1**, 16088 (2016).

29. Valles-Colomer, M. et al. The neuroactive potential of the human gut microbiota in quality of life and depression. *Nat. Microbiol.* **4**, 623–632 (2019).

30. Maresz, K. Growing evidence of a proven mechanism shows vitamin K2 Can impact health conditions beyond bone and cardiovascular. *Integr. Med.* **20**, 34–38 (2021).

31. Zeevi, D. et al. Structural variation in the gut microbiome associates with host health. *Nature* **568**, 43–48 (2019).

32. Wang, D. et al. Characterization of gut microbial structural variations as determinants of human bile acid metabolism. *Cell Host Microbe* **29**, 1802–1814.e5 (2021).

33. Dhir, S., Tarasenko, M., Napoli, E. & Giulivi, C. Neurological, psychiatric, and biochemical aspects of thiamine deficiency in children and adults. *Front. Psychiatry* **4**, 207 (2019).

34. Winston, A., Jamieson, C., Madira, W., Gatward, N. & Palmer, R. L. Prevalence of thiamin deficiency in anorexia nervosa. *Int. J. Eat. Disord.* **28**, 451–454 (2000).

35. Perino, A. et al. Central anorexigenic actions of bile acids are mediated by TGR5. *Nat. Metab.* **3**, 595–603 (2021).

36. Germain, N. et al. Constitutional thinness and lean anorexia nervosa display opposite concentrations of peptide YY, glucagon-like peptide 1, ghrelin, and leptin. *Am. J. Clin. Nutr.* **85**, 967–971 (2007).

37. Kamal, N. et al. Delayed gastrointestinal transit times in anorexia nervosa and bulimia nervosa. *Gastroenterology* **101**, 1320–1324 (1991).

38. Sun, C.-Y. et al. Indoxyl sulfate caused behavioral abnormality and neurodegeneration in mice with unilateral nephrectomy. *Aging* **13**, 6681 (2021).

39. Brydges, C. R. et al. Indoxyl sulfate, a gut microbiome-derived uremic toxin, is associated with psychic anxiety and its functional magnetic resonance imaging-based neurologic signature. *Sci. Rep.* **11**, 21011 (2021).

40. Lee, J.-H. & Lee, J. Indole as an intercellular signal in microbial communities. *FEMS Microbiol. Rev.* **34**, 426–444 (2010).

41. Newton, W. A. & Snell, E. E. Formation and interrelationships of tryptophanase and tryptophan synthetases in *Escherichia coli*. *J. Bacteriol.* **89**, 355–364 (1965).

42. Yanchuk, P. et al. Role of serotonin in the regulation of respiration and bile secretory function of the liver. *Fiziol. Zh.* **61**, 102–110 (2015).

43. Pedersen, H. K. et al. Human gut microbes impact host serum metabolome and insulin sensitivity. *Nature* **535**, 376–381 (2016).

44. Glenny, E. M. et al. Gut microbial communities from patients with anorexia nervosa do not influence body weight in recipient germ-free mice. *Gut Microbes* **13**, 1897216 (2021).

45. Smith, M. I. et al. Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. *Science* **339**, 548–554 (2013).

46. Unger, T. J., Calderon, G. A., Bradley, L. C., Sena-Esteves, M. & Rios, M. Selective deletion of Bdnf in the ventromedial and dorsomedial hypothalamus of adult mice results in hyperphagic behavior and obesity. *J. Neurosci.* **27**, 14265–14274 (2007).

47. Yang, S.-C., Shieh, K.-R. & Li, H.-Y. Cocaine-and amphetamine-regulated transcript in the nucleus accumbens participates in the regulation of feeding behavior in rats. *Neuroscience* **133**, 841–851 (2005).

48. Hashimoto, M. & Masliah, E. Alpha-synuclein in Lewy body disease and Alzheimer's disease. *Brain Pathol.* **9**, 707–720 (1999).

49. Konopelski, P. et al. Indole-3-propionic acid, a tryptophan-derived bacterial metabolite, increases blood pressure via cardiac and vascular mechanisms in rats. *Am. J. Physiol. Regul.* **321**, R969–R981 (2021).
50. Makki, K. et al. 6α-hydroxylated bile acids mediate TGR5 signalling to improve glucose metabolism upon dietary fiber supplementation in mice. *Gut* **72**, 314–324 (2023).

51. Higuchi, S. et al. Bile acid composition regulates GPR119-dependent intestinal lipid sensing and food intake regulation in mice. *Gut* **69**, 1620–1628 (2020).

52. Breton, J. et al. Characterizing the metabolic perturbations induced by activity-based anorexia in the C57Bl/6 mouse using 1H NMR spectroscopy. *Clin. Nutr.* **39**, 2428–2434 (2020).

53. Braniste, V. et al. The gut microbiota influences blood-brain barrier permeability in mice. *Sci. Transl. Med.* **6**, 263ra158 (2014).

54. Gröbner, E. M. et al. The effects of probiotics administration on the gut microbiome in adolescents with anorexia nervosa—a study protocol for a longitudinal, double-blind, randomized, placebo-controlled trial. *Eur. Eat. Disord. Rev.* **30**, 61–74 (2022).

55. Breton, J. et al. Elevated plasma concentrations of bacterial ClpB protein in patients with eating disorders. *Int. J. Eat. Disord.* **49**, 805–808 (2016).

56. Costea, P. I. et al. Towards standards for human fecal sample processing in metagenomic studies. *Nat. Biotechnol.* **35**, 1069–1076 (2017).

57. Hahne, F. et al. flowCore: a Bioconductor package for high throughput flow cytometry. *BMC Bioinformatics* **10**, 106 (2009).

58. Pons, N. et al. METEOR, a platform for quantitative metagenomic profiling of complex ecosystems. *JOBIM http://www.jobim2010. fr/sites/default/files/presentations/27Pons.pdf* (2010).

59. Criscuolo, A. & Brisse, S. AlienTrimmer: a tool to quickly and accurately trim off multiple short contaminant sequences from high-throughput sequencing reads. *Genomics* **102**, 500–506 (2013).

60. Wen, C. et al. Quantitative metagenomics reveals unique gut microbiome biomarkers in ankylosing spondylitis. *Genome Biol.* **18**, 142 (2017).

61. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **9**, 357–359 (2012).

62. Le Chatelier, E. et al. Richness of human gut microbiome correlates with metabolic markers. *Nature* **500**, 541–546 (2013).

63. Plaza Oñate, F. et al. MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data. *Bioinformatics* **35**, 1544–1552 (2019).

64. Plaza Oñate, F. P. et al. *Updated Metagenomic Species Pan-genomes (MSPs) of the Human Gastrointestinal Microbiota* (Recherche Data Gouv, 2021); https://doi.org/10.15454/FLANUP

65. Fromentin, S. et al. Microbiome and metabolome features of the cardiometabolic disease spectrum. *Nat. Med.* **28**, 303–314 (2022).

66. Vieira-Silva, S. et al. Statin therapy is associated with lower prevalence of gut microbiota dysbiosis. *Nature* **581**, 310–315 (2020).

67. Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods* **12**, 59–60 (2015).

68. Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. *Nucleic Acids Res.* **44**, D457–D462 (2016).

69. Tingley, D., Yamamoto, T., Hirose, K., Keele, L. & Imai, K. mediation: R package for causal mediation analysis. *J. Stat. Softw.* **59**, 1–38 (2014).

70. LaPierre, N. et al. MiCoP: microbial community profiling method for detecting viral and fungal organisms in metagenomic samples. *BMC Genomics* **20**, 423 (2019).

71. Brister, J. R., Ako-Adjei, D., Bao, Y. & Blinkova, O. NCBI viral genomes resource. *Nucleic Acids Res.* **43**, D571–D577 (2015).

72. Vavrek, M. J. Fossil: palaeoecological and palaeogeographical analysis tools. *Palaeontol. Electronica* **14**, 1–16 (2011).

73. Dixon, P. VEGAN, a package of R functions for community ecology. *J. Veg. Sci.* **14**, 927–930 (2003).

74. Friedman, J. & Alm, E. J. Inferring correlation networks from genomic survey data. *PLoS Comput. Biol.* **8**, e1002687 (2012).

75. Nicholson, J. K. et al. Host-gut microbiota metabolic interactions. *Science* **336**, 1262–1267 (2012).

76. Krautkramer, K. A., Fan, J. & Bäckhed, F. Gut microbial metabolites as multi-kingdom intermediates. *Nat. Rev. Microbiol.* **19**, 77–94 (2021).

77. Castillo, S., Mattila, I., Miettinen, J., Oresic, M. & Hyötyläinen, T. Data analysis tool for comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry. *Anal. Chem.* **83**, 3058–3067 (2011).

78. Salihović, S. et al. Simultaneous determination of perfluoroalkyl substances and bile acids in human serum using ultra-high-performance liquid chromatography-tandem mass spectrometry. *Anal. Bioanal. Chem.* **412**, 2251–2259 (2020).

79. Li, B. et al. Identification of optimal reference genes for RT-qPCR in the rat hypothalamus and intestine for the study of obesity. *Int. J. Obes.* **38**, 192–197 (2014).

80. Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. *Nat. Methods* **13**, 581–583 (2016).

81. Thirion, F. et al. The gut microbiota in multiple sclerosis varies with disease activity. *Genome Med.* **15**, 1 (2023).

## Acknowledgements

Y.F. received the Marie Skłodowska-Curie Individual Fellowship (797267) dedicated to the present study. The Novo Nordisk Foundation Center for Basic Metabolic Research is an independent research institution at the University of Copenhagen, partially funded by an unrestricted donation from the Novo Nordisk Foundation. R.K.S. was supported by Odense University Hospital Research Fund (R15-A800).

## Author contributions

O.P. conceived and designed the study concept and detailed the protocol and supervised all phases of the study protocol including paper writing. Y.F., R.K.S., S.B.I. and F.T. performed the data analysis, interpreted the results and contributed to drafting the paper. T.H., T.S. and O.R. supported the metabolome profiling. Y.F. designed and undertook all animal experiments with support from T.A. and L.L. L.L. and E.S. performed bacterial cell counting. R.K.S. and T.H.H. performed phenotyping of study participants. P.D, N.P., N.G., B.Q., F.L., H.R., G.F., S.V.-S., J.R., L.C., G.K.T, F.B, M.O. and S.D.E. contributed to technical project development and supervision and revised the drafts of the paper. All authors contributed to the review and editing of the paper.

## Competing interests

F.B. is a shareholder in Implexion Pharma. The other authors declare no competing interests.

## Additional information

Extended data is available for this paper at
https://doi.org/10.1038/s41564-023-01355-5.

Supplementary information The online version contains supplementary material available at
https://doi.org/10.1038/s41564-023-01355-5.

Correspondence and requests for materials should be addressed to Oluf Pedersen.

Peer review information *Nature Microbiology* thanks Tom Hildebrandt and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Reprints and permissions information is available at
www.nature.com/reprints.
50. Makki, K. et al. 6α-hydroxylated bile acids mediate TGR5 signalling to improve glucose metabolism upon dietary fiber supplementation in mice. *Gut* **72**, 314–324 (2023).

51. Higuchi, S. et al. Bile acid composition regulates GPR119-dependent intestinal lipid sensing and food intake regulation in mice. *Gut* **69**, 1620–1628 (2020).

52. Breton, J. et al. Characterizing the metabolic perturbations induced by activity-based anorexia in the C57Bl/6 mouse using 1H NMR spectroscopy. *Clin. Nutr.* **39**, 2428–2434 (2020).

53. Braniste, V. et al. The gut microbiota influences blood-brain barrier permeability in mice. *Sci. Transl. Med.* **6**, 263ra158 (2014).

54. Gröbner, E. M. et al. The effects of probiotics administration on the gut microbiome in adolescents with anorexia nervosa—a study protocol for a longitudinal, double-blind, randomized, placebo-controlled trial. *Eur. Eat. Disord. Rev.* **30**, 61–74 (2022).

55. Breton, J. et al. Elevated plasma concentrations of bacterial ClpB protein in patients with eating disorders. *Int. J. Eat. Disord.* **49**, 805–808 (2016).

56. Costea, P. I. et al. Towards standards for human fecal sample processing in metagenomic studies. *Nat. Biotechnol.* **35**, 1069–1076 (2017).

57. Hahne, F. et al. flowCore: a Bioconductor package for high throughput flow cytometry. *BMC Bioinformatics* **10**, 106 (2009).

58. Pons, N. et al. METEOR, a platform for quantitative metagenomic profiling of complex ecosystems. *JOBIM http://www.jobim2010. fr/sites/default/files/presentations/27Pons.pdf* (2010).

59. Criscuolo, A. & Brisse, S. AlienTrimmer: a tool to quickly and accurately trim off multiple short contaminant sequences from high-throughput sequencing reads. *Genomics* **102**, 500–506 (2013).

60. Wen, C. et al. Quantitative metagenomics reveals unique gut microbiome biomarkers in ankylosing spondylitis. *Genome Biol.* **18**, 142 (2017).

61. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **9**, 357–359 (2012).

62. Le Chatelier, E. et al. Richness of human gut microbiome correlates with metabolic markers. *Nature* **500**, 541–546 (2013).

63. Plaza Oñate, F. et al. MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data. *Bioinformatics* **35**, 1544–1552 (2019).

64. Plaza Oñate, F. P. et al. *Updated Metagenomic Species Pan-genomes (MSPs) of the Human Gastrointestinal Microbiota* (Recherche Data Gouv, 2021); https://doi.org/10.15454/FLANUP

65. Fromentin, S. et al. Microbiome and metabolome features of the cardiometabolic disease spectrum. *Nat. Med.* **28**, 303–314 (2022).

66. Vieira-Silva, S. et al. Statin therapy is associated with lower prevalence of gut microbiota dysbiosis. *Nature* **581**, 310–315 (2020).

67. Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods* **12**, 59–60 (2015).

68. Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. *Nucleic Acids Res.* **44**, D457–D462 (2016).

69. Tingley, D., Yamamoto, T., Hirose, K., Keele, L. & Imai, K. mediation: R package for causal mediation analysis. *J. Stat. Softw.* **59**, 1–38 (2014).

70. LaPierre, N. et al. MiCoP: microbial community profiling method for detecting viral and fungal organisms in metagenomic samples. *BMC Genomics* **20**, 423 (2019).

71. Brister, J. R., Ako-Adjei, D., Bao, Y. & Blinkova, O. NCBI viral genomes resource. *Nucleic Acids Res.* **43**, D571–D577 (2015).

72. Vavrek, M. J. Fossil: palaeoecological and palaeogeographical analysis tools. *Palaeontol. Electronica* **14**, 1–16 (2011).

73. Dixon, P. VEGAN, a package of R functions for community ecology. *J. Veg. Sci.* **14**, 927–930 (2003).

74. Friedman, J. & Alm, E. J. Inferring correlation networks from genomic survey data. *PLoS Comput. Biol.* **8**, e1002687 (2012).

75. Nicholson, J. K. et al. Host-gut microbiota metabolic interactions. *Science* **336**, 1262–1267 (2012).

76. Krautkramer, K. A., Fan, J. & Bäckhed, F. Gut microbial metabolites as multi-kingdom intermediates. *Nat. Rev. Microbiol.* **19**, 77–94 (2021).

77. Castillo, S., Mattila, I., Miettinen, J., Oresic, M. & Hyötyläinen, T. Data analysis tool for comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry. *Anal. Chem.* **83**, 3058–3067 (2011).

78. Salihović, S. et al. Simultaneous determination of perfluoroalkyl substances and bile acids in human serum using ultra-high-performance liquid chromatography-tandem mass spectrometry. *Anal. Bioanal. Chem.* **412**, 2251–2259 (2020).

79. Li, B. et al. Identification of optimal reference genes for RT-qPCR in the rat hypothalamus and intestine for the study of obesity. *Int. J. Obes.* **38**, 192–197 (2014).

80. Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. *Nat. Methods* **13**, 581–583 (2016).

81. Thirion, F. et al. The gut microbiota in multiple sclerosis varies with disease activity. *Genome Med.* **15**, 1 (2023).

## Acknowledgements

Y.F. received the Marie Skłodowska-Curie Individual Fellowship (797267) dedicated to the present study. The Novo Nordisk Foundation Center for Basic Metabolic Research is an independent research institution at the University of Copenhagen, partially funded by an unrestricted donation from the Novo Nordisk Foundation. R.K.S. was supported by Odense University Hospital Research Fund (R15-A800).

## Author contributions

O.P. conceived and designed the study concept and detailed the protocol and supervised all phases of the study protocol including paper writing. Y.F., R.K.S., S.B.I. and F.T. performed the data analysis, interpreted the results and contributed to drafting the paper. T.H., T.S. and O.R. supported the metabolome profiling. Y.F. designed and undertook all animal experiments with support from T.A. and L.L. L.L. and E.S. performed bacterial cell counting. R.K.S. and T.H.H. performed phenotyping of study participants. P.D, N.P., N.G., B.Q., F.L., H.R., G.F., S.V.-S., J.R., L.C., G.K.T, F.B, M.O. and S.D.E. contributed to technical project development and supervision and revised the drafts of the paper. All authors contributed to the review and editing of the paper.

## Competing interests

F.B. is a shareholder in Implexion Pharma. The other authors declare no competing interests.

## Additional information

Extended data is available for this paper at
https://doi.org/10.1038/s41564-023-01355-5.

Supplementary information The online version contains supplementary material available at
https://doi.org/10.1038/s41564-023-01355-5.

Correspondence and requests for materials should be addressed to Oluf Pedersen.

Peer review information *Nature Microbiology* thanks Tom Hildebrandt and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Reprints and permissions information is available at
www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this

article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023

Yong Fan Ⓓ1,20, wRené Klinkby Støving$^{2,20}$, Samar Berreira Ibraim$^{3,20}$, Tuulia Hyötyläinen$^{4}$, Florence Thirion Ⓓ3, Tulika Arora Ⓓ1, Liwei Lyu Ⓓ1,5, Evelina Stankevic Ⓓ1, Tue Haldor Hansen Ⓓ1, Pierre Déchelotte Ⓓ6, Tim Sinioja$^{4}$, Oddny Ragnarsdottir$^{4}$, Nicolas Pons$^{3}$, Nathalie Galleron$^{3}$, Benoît Quinquis$^{3}$, Florence Levenez$^{3}$, Hugo Roume Ⓓ3, Gwen Falony$^{7,8,9,10}$, Sara Vieira-Silva Ⓓ7,8,9,10, Jeroen Raes Ⓓ7,8, Loa Clausen Ⓓ11,12, Gry Kjaersdam Telléus$^{13,14}$, Fredrik Bäckhed$^{1,15,16}$, Matej Oresic Ⓓ17,18, S. Dusko Ehrlich$^{3,19}$ & Oluf Pedersen Ⓓ1,5✉

<footnote>$^{1}$Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Science, University of Copenhagen, Copenhagen, Denmark. $^{2}$Center for Eating Disorders, Odense University Hospital, and Research Unit for Medical Endocrinology, Mental Health Services in the Region of Southern Denmark, Open Patient data Explorative Network (OPEN) and Clinical Institute, University of Southern Denmark, Odense, Denmark. $^{3}$Université Paris-Saclay, INRAE, MGP, Jouy-en-Josas, France. $^{4}$School of Science and Technology, Örebro University, Örebro, Sweden.</footnote>

<footnote>$^{5}$Department of Medicine, University of Copenhagen and Herlev-Gentofte University Hospital, Copenhagen, Denmark. $^{6}$INSERM U1073 Research Unit and TargEDys, Rouen University, Rouen, France. $^{7}$Laboratory of Molecular bacteriology, Department of Microbiology and Immunology, Rega Institute Ku Leuven, Leuven, Belgium. $^{8}$Center for Microbiology, VIB, Leuven, Belgium. $^{9}$Institute of Medical Microbiology and Hygiene and Research Center for Immunotherapy (FZI), University Medical Center of the Johannes Gutenberg-University Mainz, Mainz, Germany. $^{10}$Institute of Molecular Biology (IMB), Mainz, Germany. $^{11}$Department of Child and Adolescent Psychiatry, Aarhus University Hospital, Aarhus, Denmark.</footnote>

<footnote>$^{12}$Department of Clinical Medicine, Faculty of Health, Aarhus University, Aarhus, Denmark. $^{13}$Unit for Psychiatric Research, Aalborg University Hospital, Aalborg, Denmark. $^{14}$Department of Communication and Psychology, The Faculty of Social Sciences and Humanities, Aalborg University, Aalborg, Denmark. $^{15}$The Wallenberg Laboratory, Department of Molecular and Clinical Medicine, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden. $^{16}$Department of Clinical Physiology, Sahlgrenska University Hospital, Region Västra Götaland, Gothenburg, Sweden. $^{17}$School of Medical Sciences, Örebro University, Örebro, Sweden. $^{18}$Turku Bioscience Centre, University of Turku and Åbo Akademi University, Turku, Finland. $^{19}$Department of Clinical and Movement Neurosciences, University College London, London, UK.</footnote>

<footnote>$^{20}$These authors contributed equally: Yong Fan, René Støving, Samar Berreira Ibraim. ✉e-mail: oluf@sund.ku.dk</footnote>
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this

article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023

Yong Fan Ⓓ1,20, wRené Klinkby Støving$^{2,20}$, Samar Berreira Ibraim$^{3,20}$, Tuulia Hyötyläinen$^{4}$, Florence Thirion Ⓓ3, Tulika Arora Ⓓ1, Liwei Lyu Ⓓ1,5, Evelina Stankevic Ⓓ1, Tue Haldor Hansen Ⓓ1, Pierre Déchelotte Ⓓ6, Tim Sinioja$^{4}$, Oddny Ragnarsdottir$^{4}$, Nicolas Pons$^{3}$, Nathalie Galleron$^{3}$, Benoît Quinquis$^{3}$, Florence Levenez$^{3}$, Hugo Roume Ⓓ3, Gwen Falony$^{7,8,9,10}$, Sara Vieira-Silva Ⓓ7,8,9,10, Jeroen Raes Ⓓ7,8, Loa Clausen Ⓓ11,12, Gry Kjaersdam Telléus$^{13,14}$, Fredrik Bäckhed$^{1,15,16}$, Matej Oresic Ⓓ17,18, S. Dusko Ehrlich$^{3,19}$ & Oluf Pedersen Ⓓ1,5✉

<footnote>$^{1}$Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Science, University of Copenhagen, Copenhagen, Denmark. $^{2}$Center for Eating Disorders, Odense University Hospital, and Research Unit for Medical Endocrinology, Mental Health Services in the Region of Southern Denmark, Open Patient data Explorative Network (OPEN) and Clinical Institute, University of Southern Denmark, Odense, Denmark. $^{3}$Université Paris-Saclay, INRAE, MGP, Jouy-en-Josas, France. $^{4}$School of Science and Technology, Örebro University, Örebro, Sweden.</footnote>

<footnote>$^{5}$Department of Medicine, University of Copenhagen and Herlev-Gentofte University Hospital, Copenhagen, Denmark. $^{6}$INSERM U1073 Research Unit and TargEDys, Rouen University, Rouen, France. $^{7}$Laboratory of Molecular bacteriology, Department of Microbiology and Immunology, Rega Institute Ku Leuven, Leuven, Belgium. $^{8}$Center for Microbiology, VIB, Leuven, Belgium. $^{9}$Institute of Medical Microbiology and Hygiene and Research Center for Immunotherapy (FZI), University Medical Center of the Johannes Gutenberg-University Mainz, Mainz, Germany. $^{10}$Institute of Molecular Biology (IMB), Mainz, Germany. $^{11}$Department of Child and Adolescent Psychiatry, Aarhus University Hospital, Aarhus, Denmark.</footnote>

<footnote>$^{12}$Department of Clinical Medicine, Faculty of Health, Aarhus University, Aarhus, Denmark. $^{13}$Unit for Psychiatric Research, Aalborg University Hospital, Aalborg, Denmark. $^{14}$Department of Communication and Psychology, The Faculty of Social Sciences and Humanities, Aalborg University, Aalborg, Denmark. $^{15}$The Wallenberg Laboratory, Department of Molecular and Clinical Medicine, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden. $^{16}$Department of Clinical Physiology, Sahlgrenska University Hospital, Region Västra Götaland, Gothenburg, Sweden. $^{17}$School of Medical Sciences, Örebro University, Örebro, Sweden. $^{18}$Turku Bioscience Centre, University of Turku and Åbo Akademi University, Turku, Finland. $^{19}$Department of Clinical and Movement Neurosciences, University College London, London, UK.</footnote>

<footnote>$^{20}$These authors contributed equally: Yong Fan, René Støving, Samar Berreira Ibraim. ✉e-mail: oluf@sund.ku.dk</footnote>
# Gut microbial alterations and their contributions to anorexia

## ① Phenotype information

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

✓ Bioclinical traits

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

✓ Scores of eating behavior and psychological traits

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

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

Anorexia
n = 77

Healthy control
n = 70

Age-matched

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

## ④ Serum metabolome

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

63 polar and microbiota-related metabolites

## ② Bacterial gut microbiota

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

### Taxonomic profile

(1,990 Metagenomic Species Pan-genomes)

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

### Bacterial genetics

(2,423 variable structural variants, 5,056 deletion structural variants)

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

### Growth dynamics

(50 peak-to-trough ratios)

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

### Functional profiles

(103 gut metabolic modules, 56 gut brain modules)

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

The gut microbiome

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

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

## ③ Viral gut microbiota

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

209 phages or virus

## ⑤ Gut microbiota transplantation

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

Female germ-free littermate breeding

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

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

Anorexia
microbiota

Fecal microbiota
transplantation

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

Healthy
control
microbiota

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

Differential bacterial species

27 amplicon sequence variants

18 amplicon sequence variants

## ⑥ Differential gene expression in hypothalamus and adipose tissue

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

Hypothalamic appetite-suppressing gene markers

Anorexia-
transplanted

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

Healthy control-
transplanted

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

Adipose
thermogenic genes

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

Differential gene expression High Low

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

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

Extended Data Fig. 1 | Graphical abstract of the study workflow and findings. Workflow was created with BioRender.com.
# Gut microbial alterations and their contributions to anorexia

<figure><img src="image_19.png" /><figcaption>Extended Data Fig. 1 | Graphical abstract of the study workflow and findings. Workflow was created with BioRender.com.</figcaption></figure>
<figure><img src="image_40.png" /></figure>

**Extended Data Fig. 2** | Taxonomic differences at phylum, family and genus level and differences in enterotypes in cases with AN compared with healthy women. **a,b,e** Box plot (line, median; box, interquartile range (IQR); whiskers, 1.5× IQR) of comparison between AN (gold, *n* = 77) and HC (blue, *n* = 70) of relative abundance of the 12 bacterial phyla detected in at least 10% of individuals (*a*), the 20 most abundant bacterial families (*b*), and the top 30 genera (*e*). Features are sorted by decreasing mean abundance. Zero values are set to 1e-10. Features colored in blue are enriched in HC group, and gold are enriched in AN group. Significance was determined by two-sided Wilcoxon rank-sum test, followed by drug-deconfounding and multiple testing correction by Benjamini-Hochberg method on all features (**a,b**, and **e**, see Source Data for exact *p* values).

**c,f,** Absolute values of Cliff's Delta effect size of families (**c**) and genera (**f**) contrasted between AN and HC after drug-deconfounding (adjusted *p*-value ≤ 0.1). Gold barplots indicate features more abundant in AN; blue barplots indicate features more abundant in HC (**c,f**). **d,g,** Box plot (line, median; box, IQR; whiskers, 1.5× IQR) showing β-diversity (Canberra distance) of gut bacteriome at genus level (**d**) and richness of Metagenomic Species Pan-genomes (MSPs) (**g**) between AN (gold, *n* = 77) and HC (blue, *n* = 70). Significance was determined by two-sided Wilcoxon rank-sum test (**d,g**). **h,** Upper panel demonstrates enterotype prevalence in AN (*n* = 77) and HC (*n* = 70), lower panel shows enterotype prevalence in HC and AN patients split into AN-RS (*n* = 56) and AN-BP (*n* = 21) groups.

Contingency test Chisq *p* = 0.0021
<figure><img src="image_20.png" /></figure>

**Extended Data Fig. 2** | Taxonomic differences at phylum, family and genus level and differences in enterotypes in cases with AN compared with healthy women. **a,b,e** Box plot (line, median; box, interquartile range (IQR); whiskers, 1.5× IQR) of comparison between AN (gold, *n* = 77) and HC (blue, *n* = 70) of relative abundance of the 12 bacterial phyla detected in at least 10% of individuals (*a*), the 20 most abundant bacterial families (*b*), and the top 30 genera (*e*). Features are sorted by decreasing mean abundance. Zero values are set to 1e-10. Features colored in blue are enriched in HC group, and gold are enriched in AN group. Significance was determined by two-sided Wilcoxon rank-sum test, followed by drug-deconfounding and multiple testing correction by Benjamini-Hochberg method on all features (**a,b**, and **e**, see Source Data for exact *p* values).

**c,f,** Absolute values of Cliff's Delta effect size of families (**c**) and genera (**f**) contrasted between AN and HC after drug-deconfounding (adjusted *p*-value ≤ 0.1). Gold barplots indicate features more abundant in AN; blue barplots indicate features more abundant in HC (**c,f**). **d,g,** Box plot (line, median; box, IQR; whiskers, 1.5× IQR) showing β-diversity (Canberra distance) of gut bacteriome at genus level (**d**) and richness of Metagenomic Species Pan-genomes (MSPs) (**g**) between AN (gold, *n* = 77) and HC (blue, *n* = 70). Significance was determined by two-sided Wilcoxon rank-sum test (**d,g**). **h,** Upper panel demonstrates enterotype prevalence in AN (*n* = 77) and HC (*n* = 70), lower panel shows enterotype prevalence in HC and AN patients split into AN-RS (*n* = 56) and AN-BP (*n* = 21) groups.
**Anorexia-enriched**

**Control-enriched**

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

**Extended Data Fig. 3** | Co-occurrence network deducted from bacterial species enriched in AN cases (left) and healthy controls (right) after drug-deconfounding. Node size and node colour represent the mean abundance and the genus of a given MSP, respectively. Genera represented by only one MSP are

coloured in grey. Red and blue lines indicate positive and negative correlations, respectively. Line thickness represents the absolute correlation coefficient. Only correlations with an absolute coefficient above 0.4 are shown; MSP without any correlation above the threshold are hidden.
**Anorexia-enriched**

**Control-enriched**

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

**Extended Data Fig. 3** | Co-occurrence network deducted from bacterial species enriched in AN cases (left) and healthy controls (right) after drug-deconfounding. Node size and node colour represent the mean abundance and the genus of a given MSP, respectively. Genera represented by only one MSP are

coloured in grey. Red and blue lines indicate positive and negative correlations, respectively. Line thickness represents the absolute correlation coefficient. Only correlations with an absolute coefficient above 0.4 are shown; MSP without any correlation above the threshold are hidden.
<figure><img src="image_42.png" /><figcaption>Extended Data Fig. 4 | Heatmap showing that gut bacterial genera are linked to eating disorder scores in AN cases after drug de-confounding. Variables in specific eating disorder scale are marked in blue, and in general psychological scale are marked in red. +, adjusted p < 0.05 (see Source Data for exact p values).</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Extended Data Fig. 4 | Heatmap showing that gut bacterial genera are linked to eating disorder scores in AN cases after drug de-confounding. Variables in specific eating disorder scale are marked in blue, and in general psychological scale are marked in red. +, adjusted p < 0.05 (see Source Data for exact p values).</figcaption></figure>
a

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

**Extended Data Fig. 5** | Fasting plasma concentration of Caseinolytic protease B (ClpB) in control subjects and AN subtypes. **a**, Relative abundance of Enterobacteriaceae family in HC (*n* = 70) and AN (*n* = 77) groups. **b**, Log-scale transformed fasting plasma ClpB concentration between HC (*n* = 70) and AN (*n* = 77) groups. **c**, Log-scale transformed fasting plasma ClpB concentration

b

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

between restrictive AN (AN-RS *n* = 56) and binge-purge AN (AN-BP *n* = 21) subtypes. Significance was calculated by two-sided Wilcoxon rank-sum test between two groups (a-c). Box plots indicate median and interquartile range (IQR) and whiskers represent 1.5× IQR (a-c).

c

<figure><img src="image_45.png" /></figure>
a

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

**Extended Data Fig. 5** | Fasting plasma concentration of Caseinolytic protease B (ClpB) in control subjects and AN subtypes. a, Relative abundance of Enterobacteriaceae family in HC (n = 70) and AN (n = 77) groups. b, Log-scale transformed fasting plasma ClpB concentration between HC (n = 70) and AN (n = 77) groups. c, Log-scale transformed fasting plasma ClpB concentration

b

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

between restrictive AN (AN-RS *n* = 56) and binge-purge AN (AN-BP *n* = 21) subtypes. Significance was calculated by two-sided Wilcoxon rank-sum test between two groups (a-c). Box plots indicate median and interquartile range (IQR) and whiskers represent 1.5× IQR (a-c).

c

<figure><img src="image_25.png" /></figure>
<figure><img src="image_46.png" /><figcaption>**Extended Data Fig. 6** | Altered dynamic growth rates of bacterial species between AN (n = 77) and control (n = 70) groups. Box plots indicate median and interquartile range (IQR), whiskers indicate 1.5× IQR. Significance was determined by two-sided Wilcoxon rank-sum test.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>Extended Data Fig. 6 | Altered dynamic growth rates of bacterial species between AN (n = 77) and control (n = 70) groups. Box plots indicate median and interquartile range (IQR), whiskers indicate 1.5× IQR. Significance was determined by two-sided Wilcoxon rank-sum test.</figcaption></figure>
<caption>a</caption>

Alistipes putredinis: 936_937
Rho=-0.31, p = 2.1e-02

<figure><img src="image_47.png" /><figcaption>b</figcaption></figure>

Alistipes putredinis: 936_937
Rho=-0.26, p = 5.9e-02

<figure><img src="image_48.png" /><figcaption>c</figcaption></figure>

Alistipes putredinis: 936_937
Rho=-0.33, p = 9.3e-03

<figure><img src="image_49.png" /><figcaption>d</figcaption></figure>

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

Extended Data Fig. 7 | SVs in bacterial gut microbiota from AN cases and associations with anorexia-relevant traits. a-c, Scatterplot showing, HOMA-IR, fasting plasma insulin and glucose in individuals harboring a 1-kbp variation in the *A. putredinis* genome (n = 147). Significance determined by two-sided spearman correlation test corrected by false discovery rate using Benjamini and

Hochberg method. Error band is linear regression line with 95% confidence band (a-c). Dots representing AN and HC individuals were colored in red and blue, respectively. d, Upper panel, standardized variability (y axis) along a genomic region of *A. putredinis* (x axis). Lower panel, locations (blue bar) of the gene of interest.
<caption>a</caption>

Alistipes putredinis: 936_937
Rho=-0.31, p = 2.1e-02

<figure><img src="image_27.png" /><figcaption>b</figcaption></figure>

Alistipes putredinis: 936_937
Rho=-0.26, p = 5.9e-02

<figure><img src="image_28.png" /><figcaption>c</figcaption></figure>

Alistipes putredinis: 936_937
Rho=-0.33, p = 9.3e-03

<figure><img src="image_29.png" /><figcaption>d</figcaption></figure>

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

Extended Data Fig. 7 | SVs in bacterial gut microbiota from AN cases and associations with anorexia-relevant traits. a-c, Scatterplot showing, HOMA-IR, fasting plasma insulin and glucose in individuals harboring a 1-kbp variation in the *A. putredinis* genome (n = 147). Significance determined by two-sided spearman correlation test corrected by false discovery rate using Benjamini and

Hochberg method. Error band is linear regression line with 95% confidence band (a-c). Dots representing AN and HC individuals were colored in red and blue, respectively. d, Upper panel, standardized variability (y axis) along a genomic region of *A. putredinis* (x axis). Lower panel, locations (blue bar) of the gene of interest.
<figure><img src="image_51.png" /></figure>

**Extended Data Fig. 8 | In silico analysis using bidirectional mediation inference.** a, Principal component analysis (PCA) plot of the serum metabolome of AN subtypes, AN-BP (binge/purge anorexia), AN-RS (restrictive anorexia). b, Summary number of inferred mediation relationship for direction 1 (gut microbial features → eating disorder scores mediated by serum metabolites), direction 2 (gut microbial features → serum metabolites mediated by eating disorder scores). c, Summary number of inferred mediation relationship for direction 1 (gut microbial features → phenotypes mediated by serum metabolites), direction 2 (gut microbial features → serum metabolites mediated by phenotypes) for the whole cohort. d, Sankey diagram showing the inferred causal relationship network of direction 1 where gut microbial features including bacterial species, gut brain and metabolic modules, and bacterial genetics

were treated as causal factors, serum metabolites are mediators, and metabolic traits are outcomes. e, Examples of inferred causal relationships between gut microbial features, metabolites, and host metabolic traits. Direction 1 that means microbial features → metabolic traits mediated by serum metabolites is illustrated with a black line while direction 2 that means microbial features → serum metabolites mediated by metabolic traits is illustrated with a stipulated red line. The proportion of mediation effect are shown at the center of ring charts. GDCA, glycodeoxycholic acid; GHCA, glycohyocholic acid; GHDCA, glycohyodeoxycholic acid; GUDCA, glycoursodeoxycholic acid; HCA, hyocholic acid; HOMA-IR homeostatic model assessment of insulin resistance; P-, plasma; S-, serum.
<figure><img src="image_31.png" /></figure>

**Extended Data Fig. 8** | *In silico* analysis using bidirectional mediation inference. **a**, Principal component analysis (PCA) plot of the serum metabolome of AN subtypes, AN-BP (binge/purge anorexia), AN-RS (restrictive anorexia). **b**, Summary number of inferred mediation relationship for direction 1 (gut microbial features → eating disorder scores mediated by serum metabolites), direction 2 (gut microbial features → serum metabolites mediated by eating disorder scores). **c**, Summary number of inferred mediation relationship for direction 1 (gut microbial features → phenotypes mediated by serum metabolites), direction 2 (gut microbial features → serum metabolites mediated by phenotypes) for the whole cohort. **d**, Sankey diagram showing the inferred causal relationship network of direction 1 where gut microbial features including bacterial species, gut brain and metabolic modules, and bacterial genetics

were treated as causal factors, serum metabolites are mediators, and metabolic traits are outcomes. **e**, Examples of inferred causal relationships between gut microbial features, metabolites, and host metabolic traits. Direction 1 that means microbial features → metabolic traits mediated by serum metabolites is illustrated with a black line while direction 2 that means microbial features → serum metabolites mediated by metabolic traits is illustrated with a stipulated red line. The proportion of mediation effect are shown at the center of ring charts. GDCA, glycodeoxycholic acid; GHCA, glycohyocholic acid; GHDCA, glycohyodeoxycholic acid; GUDCA, glycoursodeoxycholic acid; HCA, hyocholic acid; HOMA-IR homeostatic model assessment of insulin resistance; P-, plasma; S-, serum.
a

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

b

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

**Extended Data Fig. 9** | Workflow diagram of fecal microbiota transplantation from human donors to germ-free mice littermates. **a**, Workflow for the preparation of fecal microbiota slurry. Stool (250 mg) from both anorexia (AN) and control (HC) was cut on dry ice, transferred to anaerobic chamber, and resuspended with 5 ml of LYBHI media diluted in 20% glycerol. The resuspended fecal slurries were aliquoted in cryotubes and refrozen back quickly at -80 °C until further use. All stool samples of matched AN and HC donors were prepared on the same day and frozen aliquots were stored frozen for gavage to mice.

**b**, Experimental scheme for the GF mice transplantation study. In each independent litter study, GF female littermates at age of six weeks old were

taken out of breeding isolator and were randomly assigned to receive 200 µl of fecal slurries from three AN cases or three HC subjects. Both groups of mice were housed in autoclaved individually ventilated cages and were given autoclaved chow diet and water *ad libitum* for two days. After two days, mice were gavaged with a second dose of fecal material from the same matched AN and HC donors as before. Thereafter, mice in both groups were single housed and subjected to 30% calorie restricted chow diet for three weeks. Water was given *ad libitum* during this period. Both the anorexia-transplanted (AN-T) and the normal control-transplanted (HC-T) mice were weighed every five days after the start of energy-restricted diet. Created with Biorender.com.
a

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

b

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

**Extended Data Fig. 9** | Workflow diagram of fecal microbiota transplantation from human donors to germ-free mice littermates. **a**, Workflow for the preparation of fecal microbiota slurry. Stool (250 mg) from both anorexia (AN) and control (HC) was cut on dry ice, transferred to anaerobic chamber, and resuspended with 5 ml of LYBHI media diluted in 20% glycerol. The resuspended fecal slurries were aliquoted in cryotubes and refrozen back quickly at -80 °C until further use. All stool samples of matched AN and HC donors were prepared on the same day and frozen aliquots were stored frozen for gavage to mice.

**b**, Experimental scheme for the GF mice transplantation study. In each independent litter study, GF female littermates at age of six weeks old were

taken out of breeding isolator and were randomly assigned to receive 200 µl of fecal slurries from three AN cases or three HC subjects. Both groups of mice were housed in autoclaved individually ventilated cages and were given autoclaved chow diet and water *ad libitum* for two days. After two days, mice were gavaged with a second dose of fecal material from the same matched AN and HC donors as before. Thereafter, mice in both groups were single housed and subjected to 30% calorie restricted chow diet for three weeks. Water was given *ad libitum* during this period. Both the anorexia-transplanted (AN-T) and the normal control-transplanted (HC-T) mice were weighed every five days after the start of energy-restricted diet. Created with Biorender.com.