Obesity-associated deficits in inhibitory control are phenocopied to mice through gut microbiota changes in one-carbon and aromatic amino acids metabolic pathways

María Arnoriaga-Rodríguez, Jordi Mayneris-Perxachs, Oren Contreras-Rodríguez, Aurelijus Burokas, Juan-Antonio Ortega-Sanchez, Gerard Blasco, Claudia Coll, Carles Biarnés, Anna Castells-Nobau, Josep Puig, Josep Garre-Olmo, Rafel Ramos, Salvador Pedraza, Ramon Brugada, Joan C Vilanova, Joaquín Serena, Jordi Barretina, Jordi Gich, Vicente Pérez-Brocal, Andrés Moya, Xavier Fernández-Real, Lluis Ramio-Torrentà, Reinald Pamplona, Joaquim Sol, Mariona Jové, Wifredo Ricart, Manuel Portero-Otin, Rafael Maldonado, Jose Manuel Fernández-Real, María Arnoriaga-Rodríguez, Jordi Mayneris-Perxachs, Oren Contreras-Rodríguez, Aurelijus Burokas, Juan-Antonio Ortega-Sanchez, Gerard Blasco, Claudia Coll, Carles Biarnés, Anna Castells-Nobau, Josep Puig, Josep Garre-Olmo, Rafel Ramos, Salvador Pedraza, Ramon Brugada, Joan C Vilanova, Joaquín Serena, Jordi Barretina, Jordi Gich, Vicente Pérez-Brocal, Andrés Moya, Xavier Fernández-Real, Lluis Ramio-Torrentà, Reinald Pamplona, Joaquim Sol, Mariona Jové, Wifredo Ricart, Manuel Portero-Otin, Rafael Maldonado, Jose Manuel Fernández-Real

Abstract

Background: Inhibitory control (IC) is critical to keep long-term goals in everyday life. Bidirectional relationships between IC deficits and obesity are behind unhealthy eating and physical exercise habits.

Methods: We studied gut microbiome composition and functionality, and plasma and faecal metabolomics in association with cognitive tests evaluating inhibitory control (Stroop test) and brain structure in a discovery (n=156), both cross-sectionally and longitudinally, and in an independent replication cohort (n=970). Faecal microbiota transplantation (FMT) in mice evaluated the impact on reversal learning and medial prefrontal cortex (mPFC) transcriptomics.

Results: An interplay among IC, brain structure (in humans) and mPFC transcriptomics (in mice), plasma/faecal metabolomics and the gut metagenome was found. Obesity-dependent alterations in one-carbon metabolism, tryptophan and histidine pathways were associated with IC in the two independent cohorts. Bacterial functions linked to one-carbon metabolism (thyX,dut, exodeoxyribonuclease V), and the anterior cingulate cortex volume were associated with IC, cross-sectionally and longitudinally. FMT from individuals with obesity led to alterations in mice reversal learning. In an independent FMT experiment, human donor's bacterial functions related to IC deficits were associated with mPFC expression of one-carbon metabolism-related genes of recipient's mice.

Conclusion: These results highlight the importance of targeting obesity-related impulsive behaviour through the induction of gut microbiota shifts.

Keywords: intestinal microbiology; obesity.

Conflict of interest statement

Competing interests: None declared.

© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Figures

Figure 1
Figure 1
A microbiota taxonomic and functional signature is associated with inhibitory control. (A) Volcano plot of differential bacterial abundance (pFDR <0.05) associated with the Stroop Color Word Test (SCWT) as calculated by DESeq2 from shotgun metagenomic sequencing in the IRONMET cohort (n=114), adjusting for age, sex and education years. Fold change (FC) associated with a unit change in the SCWT and Benjamini-Hochberg adjusted p values (pFDR) are plotted for each taxon. Significantly different taxa are coloured according to phylum. In the same graph, the violin plots for the SCWT scores in patients with and without obesity are also shown. Differences between groups were analysed by a Wilcoxon tests. (B) Volcano plot of differential bacterial abundance (pFDR <0.1) associated with Iowa Gambling Test (IGT) as calculated by DESeq2 from shotgun metagenomic sequencing in an independent cohort (n=24), adjusting for age, sex and education years. (C) Manhattan-like plot of significantly expressed KEGG metagenome functions associated with the SCWT (pFDR <0.020) identified from DESeq2 analysis in the IRONMET cohort (n=114) adjusted for age, sex and educations. The −log10(pFDR) values are multiplied by the FC sign to take into account the direction of the association. Bars are coloured according to the pFDR. (D) Manhattan-like plot of significantly expressed KEGG metagenome functions associated with the IGT (pFDR <0.05) identified from DESeq2 analysis in an independent cohort (n=24) adjusted for age, sex and educations.
Figure 2
Figure 2
The microbiota taxonomic and functional signature linked to inhibitory control is modulated by obesity. (A, B) Volcano plot of differential expressed (pFDR <0.1) bacterial abundance and (C, D) bacterial functions associated with the Stroop Color and Word Test (SCWT) as calculated by DESeq2 from shotgun metagenomic sequencing in the patients without and with obesity from the IRONMET cohort, respectively, controlling for age, sex and education years. Fold change associated with a unit change in the SCWT and Benjamini-Hochberg adjusted p values (pFDR) are plotted for each taxon. Significantly different taxa are coloured according to phylum.
Figure 3
Figure 3
Plasma and faecal metabolomics linked to inhibitory control in the Ironmet and Imageomics cohorts. Bar plots of normalised variable importance measure (VIM) for the metabolites associated with the Stroop Color and Word Test (SCWT) in (A, B) plasma and (C–E) faecal samples identified by HPLC-ESI-MS/MS in positive mode (n=130), negative mode (n=130) and NMR (n=156), respectively, in the Ironmet cohort. Bar plots of VIM for the metabolites associated with the SCWT in plasma samples of the Imageomics cohort (n=970) identified by HPCL-ESI-MS/MS in (F) positive and (G) negative mode. In all cases, metabolites were identified using a multiple random forest-based machine learning variable selection strategy using the Boruta algorithm with 5000 trees and 500 iterations. All metabolites were identified based on exact mass, retention time and MS/MS spectrum, except those with (*) that were only identified based on exact mass and retention time. Unidentified metabolites are shown as exact mass at retention time. 2-PPA, 2-phenylpropanoic acid; 4-HPLA, 4-hydroxyphenyllactic acid; 5-HIAA, 5-hydroxyindoleacetic acid; DMSO2, dimethyl sulfone; FA, fatty acid; IPAM, indolepropionamide; MA, methylamine; TMA, trimethylamine.
Figure 4
Figure 4
Plasma and faecal metabolomics linked to inhibitory control in the Ironmet and Imageomics cohorts according to obesity status. Bar plots of normalised variable importance measure (VIM) for the metabolites associated with the Stroop Color and Word Test (SCWT) in (A–D) plasma and (E–J) faecal samples identified by HPLC-ESI-MS/MS in positive mode (n=130), negative mode (n=130) and NMR (n=156), respectively, in the Ironmet cohort in patients with and without obesity. Bar plots of VIM for the metabolites associated with the SCWT in plasma samples of the Imageomics cohort (n=970) identified by HPCL-ESI-MS/MS in (K) positive and (L, M) negative mode in patients with and without obesity. The above colour bar indicates the sign of the association among the metabolites and the SCWT, with red indicating negative correlation and green positive correlation. In all cases, metabolites were identified using a multiple random forest-based machine learning variable selection strategy using the Boruta algorithm with 5000 trees and 500 iterations. All metabolites were identified based on exact mass, retention time and MS/MS spectrum, except those with (*) that were only identified based on exact mass and retention time. Unidentified metabolites are shown as exact mass at retention time. 4-HPLA, 4-hydroxyphenyllactic acid; BA1, Bile acid1: 4,4-dimethyl-5-α-cholesta-8,14-dien-3β-ol; EPA, eicosapentaenoic acid; FA, fatty acid; IPAM, indolepropionamide; MA, methylamine; TMA, trimethylamine.
Figure 5
Figure 5
Main metabolic pathways involved in the associations among metagenomics, metabolomics and the Stroop Color and Word Test (SCWT). (A) Overview of the main catabolic pathways of tryptophan and tyrosine. Tryptophan and tyrosine are the precursors for the synthesis of the neurotransmitters serotonin and dopamine, respectively. The gut microbiota can also metabolise tryptophan and tyrosine to indoles and hydroxyphenolic acids, respectively. Dietary tryptophan is mostly metabolised via the Kynurenine pathway, which is activated by inflammation. (B) Overview of the folate-mediated one-carbon metabolism. The folate cycle (green) is required for the synthesis of DNA (pink and blue) as well as methylation reaction (DNA, proteins and lipids) through the methionine cycle (orange). Histidine (purple), choline and betaine are two sources of 1C units feeding into the one-carbon metabolism. Bacterial pathways have been shaded in red. Metabolites involved in the one-carbon metabolism and significantly associated with the SCWT are highlighted in bold in a yellow box. Bacterial functions participating in the one-carbon metabolism and significantly associated with the SCWT are highlighted in bold in a red box. AICAR, 5-aminoimidazole-4-carboxamide 1-β-D-ribofuranoside; bgtB, arginine/lysine/histidine/glutamine transport system substrate-binding and permease protein; bluB, 5,6-dimethylbenzimidazole synthase; btuB, vitamin B12 transporter; DHF, dihydrofolate; DMA, dimethylamine; DMSO2, dimethylsulfone; dut, dUTP pyrophosphatase; FAICAR, 5-formamidoimidazole-4-carboxamide ribotide; FGAR, 5’-phosphoribosyl-N-formylglycineamide; FIGlu, N-Formimino-glutamate; folK, 2-amino-4-hydroxy-6-hydroxymethyldihydropteridine diphosphokinase; GAR, 5’-phosphoribosylglycineamide; hisB, imidazoleglycerol-phosphate dehydratase; IGP, imidazole glycerol-phosphate; IMP, inosine 5’-monophosphate; MA, methylamine; nudB, dihydroneopterin triphosphate diphosphatase; pdxA, 4-hydroxythreonine-4-phosphate dehydrogenase; proV, glycine betaine/proline transport system ATP-binding protein; proW, glycine betaine/proline transport system permease protein; proX, glycine betaine/proline transport system substrate-binding protein; purT, phosphoribosylglycinamide formyltransferase 2; ribBA, 3,4-dihydroxy 2-butanone 4-phosphate synthase/GTP cyclohydrolase II; ribE, riboflavin synthase; SAH, S-adenosylhomocysteine; SAM, S-adenosylmethionine; THF, tetrahydrofolate; thyX, thymidylate synthase; TMA, trimethylamine.
Figure 6
Figure 6
The gut microbiota is associated with the brain area involved in inhibitory control. (A, B) The anterior cingulate cortex (ACC) volume was positively associated with the Stroop Color and Word Test (SCWT) in the Ironmet cohort (n=95) after controlling for age, sex, education years and total intracranial volume. (C) Volcano plots of differential bacterial abundance and (D) KEGG metagenome functions associated with the ACC volume as calculated by DESeq2 controlling for previous covariates. Fold change (FC) associated with a unit change in the corresponding volumes and Benjamini-Hochberg adjusted p values (pFDR) are plotted for each taxon. (E) Manhattan-like plot of significantly expressed KEGG metagenome functions associated with the ACC volume highlighting those bacterial functions also associated with the SCWT in blue. The −log10(pFDR) values are multiplied by the FC sign to take into account the direction of the association. Bars are coloured according to the pFDR.
Figure 7
Figure 7
Faecal microbiota transplantation (FMT) mice studies. (A) Experimental design for the first FMT study. The microbiota from human donors without obesity (body mass index (BMI) 2, n=11) and with obesity (BMI ≥30 kg/m2, n=11) was delivered to recipient mice pretreated with antibiotics for 14 days. Reversal learning tests (RLT) were performed after 18 days. Violin plot for the recipient’s mice RLT scores at 18 days based on human donor obesity status. (B) Volcano plot of differential human donor bacterial abundance associated with the recipient’s mice RLT at day 18, identified from DESeq2 analysis controlling for donor’s age, sex and education years. Fold change (FC) associated with a unit change in the corresponding memory test and Benjamini-Hochberg adjusted p values (pFDR) are plotted for each taxon. Significantly different taxa are coloured according to phylum. (C) Manhattan-like plot showing only the significantly expressed KEGG metagenome functions associated with the recipient’s mice RLT (pFDR <0.1) that were also associated with the Stroop Color and Word Test (SCWT) in humans. The −log10(pFDR) values are multiplied by the FC sign to take into account the direction of the association. Bars are coloured according to the pFDR. (D) Experimental design for the second FMT study. The microbiota from human donors with low (n=11) and high (n=11) SCWT scores was delivered to recipient mice pretreated with antibiotics for 14 days. RNA sequencing of the medial prefrontal cortex (mPFC) was performed after 4 weeks. Violin plots for the SCWT according to the human donor scores. (E) Volcano plots of recipient’s mice differential mPFC genes associated (pFDR <0.1) with the human’s donor metagenome functions dUTP pyrophosphatase and (F) exodeoxyribonuclease V controlling for donor’s age, sex and education years. FC associated with a unit change in the expression of the corresponding bacterial function and the Benjamini-Hochberg adjusted p values (pFDR) are plotted for each gene. (G) Bar plot of the normalised variable importance measure (VIM) for the recipient’s mice mPFC genes associated with the human donor’s SCWT identified by machine learning using multiple random forest-based variable selection strategy with the Boruta algorithm with 5000 trees and 500 iterations. (H) Volcano plot of differential human donor’s KEGG metagenome functions associated with recipient’s mice Ms4a4a and (I) Slc16a12 genes. FC associated with a unit change in the expression of both genes and Benjamini-Hochberg adjusted p values (pFDR) are plotted for each metagenome function.

References

    1. Harvey PD. Domains of cognition and their assessment. Dialogues Clin Neurosci 2019;21:227–37. 10.31887/DCNS.2019.21.3/pharvey
    1. Yang Y, Shields GS, Guo C, et al. . Executive function performance in obesity and overweight individuals: a meta-analysis and review. Neurosci Biobehav Rev 2018;84:225–44. 10.1016/j.neubiorev.2017.11.020
    1. Jansen A, Houben K, Roefs A. A cognitive profile of obesity and its translation into new interventions. Front Psychol 2015;6:1807. 10.3389/fpsyg.2015.01807
    1. Arnoriaga-Rodríguez M, Fernández-Real JM. Microbiota impacts on chronic inflammation and metabolic syndrome - related cognitive dysfunction. Rev Endocr Metab Disord 2019;20:473–80. 10.1007/s11154-019-09537-5
    1. Puig J, Biarnes C, Pedraza S, et al. . The aging Imageomics study: rationale, design and baseline characteristics of the study population. Mech Ageing Dev 2020;189:111257. 10.1016/j.mad.2020.111257
    1. Vioque J, Navarrete-Muñoz E-M, Gimenez-Monzó D, et al. . Reproducibility and validity of a food frequency questionnaire among pregnant women in a Mediterranean area. Nutr J 2013;12:26. 10.1186/1475-2891-12-26
    1. Yu G, Wang L-G, Han Y, et al. . clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 2012;16:284–7. 10.1089/omi.2011.0118
    1. Kamburov A, Stelzl U, Lehrach H, et al. . The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res 2013;41:D793–800. 10.1093/nar/gks1055
    1. Kursa MB, Rudnicki WR. Feature Selection with the Boruta Package. J Stat Softw 2010;36:1–13. 10.18637/jss.v036.i11
    1. Degenhardt F, Seifert S, Szymczak S. Evaluation of variable selection methods for random forests and omics data sets. Brief Bioinform 2019;20:492–503. 10.1093/bib/bbx124
    1. Morris MS, Jacques PF, Rosenberg IH, et al. . Circulating unmetabolized folic acid and 5-methyltetrahydrofolate in relation to anemia, macrocytosis, and cognitive test performance in American seniors. Am J Clin Nutr 2010;91:1733–44. 10.3945/ajcn.2009.28671
    1. Johnson S, Wozniak DF, Imai S. CA1 Nampt knockdown recapitulates hippocampal cognitive phenotypes in old mice which nicotinamide mononucleotide improves. NPJ Aging Mech Dis 2018;4:10. 10.1038/s41514-018-0029-z
    1. Murphy FC, Smith KA, Cowen PJ, et al. . The effects of tryptophan depletion on cognitive and affective processing in healthy volunteers. Psychopharmacology 2002;163:42–53. 10.1007/s00213-002-1128-9
    1. Walderhaug E, Lunde H, Nordvik JE, et al. . Lowering of serotonin by rapid tryptophan depletion increases impulsiveness in normal individuals. Psychopharmacology 2002;164:385–91. 10.1007/s00213-002-1238-4
    1. Evers EAT, van der Veen FM, Jolles J, et al. . Acute tryptophan depletion improves performance and modulates the BOLD response during a Stroop task in healthy females. Neuroimage 2006;32:248–55. 10.1016/j.neuroimage.2006.03.026
    1. Scholes KE, Harrison BJ, O'Neill BV, et al. . Acute serotonin and dopamine depletion improves attentional control: findings from the stroop task. Neuropsychopharmacology 2007;32:1600–10. 10.1038/sj.npp.1301262
    1. Colzato LS, van den Wildenberg WPM, van Wouwe NC, et al. . Dopamine and inhibitory action control: evidence from spontaneous eye blink rates. Exp Brain Res 2009;196:467–74. 10.1007/s00221-009-1862-x
    1. Roesch-Ely D, Scheffel H, Weiland S, et al. . Differential dopaminergic modulation of executive control in healthy subjects. Psychopharmacology 2005;178:420–30. 10.1007/s00213-004-2027-z
    1. Meyer-Lindenberg A, Miletich RS, Kohn PD, et al. . Reduced prefrontal activity predicts exaggerated striatal dopaminergic function in schizophrenia. Nat Neurosci 2002;5:267–71. 10.1038/nn804
    1. Beloborodova N, Bairamov I, Olenin A, et al. . Effect of phenolic acids of microbial origin on production of reactive oxygen species in mitochondria and neutrophils. J Biomed Sci 2012;19:89. 10.1186/1423-0127-19-89
    1. Cussotto S, Delgado I, Anesi A, et al. . Tryptophan metabolic pathways are altered in obesity and are associated with systemic inflammation. Front Immunol 2020;11:557. 10.3389/fimmu.2020.00557
    1. Shane B. Folate and vitamin B12 metabolism: overview and interaction with riboflavin, vitamin B6, and polymorphisms. Food Nutr Bull 2008;29:S5–16. 10.1177/15648265080292S103
    1. He X, Slupsky CM. Metabolic fingerprint of dimethyl sulfone (DMSO2) in microbial-mammalian co-metabolism. J Proteome Res 2014;13:5281–92. 10.1021/pr500629t
    1. Ducker GS, Rabinowitz JD. One-Carbon metabolism in health and disease. Cell Metab 2017;25:27–42. 10.1016/j.cmet.2016.08.009
    1. Panula P, Nuutinen S. The histaminergic network in the brain: basic organization and role in disease. Nat Rev Neurosci 2013;14:472–87. 10.1038/nrn3526
    1. Zhu H, Wang N, Yao L, et al. . Moderate UV exposure enhances learning and memory by promoting a novel glutamate biosynthetic pathway in the brain. Cell 2018;173:1716–27. 10.1016/j.cell.2018.04.014
    1. Nouchi R, Takeuchi H, Taki Y, et al. . Neuroanatomical bases of effortful control: evidence from a large sample of young healthy adults using voxel-based morphometry. Sci Rep 2016;6:31231. 10.1038/srep31231
    1. Xue G, Xue F, Droutman V, et al. . Common neural mechanisms underlying reversal learning by reward and punishment. PLoS One 2013;8:e82169. 10.1371/journal.pone.0082169
    1. Izquierdo A, Jentsch JD. Reversal learning as a measure of impulsive and compulsive behavior in addictions. Psychopharmacology 2012;219:607–20. 10.1007/s00213-011-2579-7
    1. Batterink L, Yokum S, Stice E. Body mass correlates inversely with inhibitory control in response to food among adolescent girls: an fMRI study. Neuroimage 2010;52:1696–703. 10.1016/j.neuroimage.2010.05.059
    1. Weygandt M, Mai K, Dommes E, et al. . Impulse control in the dorsolateral prefrontal cortex counteracts post-diet weight regain in obesity. Neuroimage 2015;109:318–27. 10.1016/j.neuroimage.2014.12.073
    1. Cho J, Yu N-K, Choi J-H, et al. . Multiple repressive mechanisms in the hippocampus during memory formation. Science 2015;350:82–7. 10.1126/science.aac7368
    1. Gulmez Karaca K, Brito DVC, Oliveira AMM. MeCP2: a critical regulator of chromatin in neurodevelopment and adult brain function. Int J Mol Sci 2019;20. 10.3390/ijms20184577. [Epub ahead of print: 16 Sep 2019].
    1. Kim B, Choi Y, Kim H-S, et al. . Methyl-CpG binding protein 2 in Alzheimer dementia. Int Neurourol J 2019;23:S72–81. 10.5213/inj.1938196.098
    1. Silva CS, Eira J, Ribeiro CA, et al. . Transthyretin neuroprotection in Alzheimer's disease is dependent on proteolysis. Neurobiol Aging 2017;59:10–14. 10.1016/j.neurobiolaging.2017.07.002
    1. Pjetri E, Zeisel SH. Deletion of one allele of MTHFD1 (methylenetetrahydrofolate dehydrogenase 1) impairs learning in mice. Behav Brain Res 2017;332:71–4. 10.1016/j.bbr.2017.05.051
    1. Glasgow SD, Ruthazer ES, Kennedy TE. Guiding synaptic plasticity: novel roles for netrin-1 in synaptic plasticity and memory formation in the adult brain. J Physiol 2020. 10.1113/JP278704. [Epub ahead of print: 03 Feb 2020].
    1. Joubert BR, den Dekker HT, Felix JF, et al. . Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns. Nat Commun 2016;7:10577. 10.1038/ncomms10577
    1. Lev-Vachnish Y, Cadury S, Rotter-Maskowitz A, et al. . L-Lactate promotes adult hippocampal neurogenesis. Front Neurosci 2019;13:403. 10.3389/fnins.2019.00403
    1. Baroncelli L, Molinaro A, Cacciante F, et al. . A mouse model for creatine transporter deficiency reveals early onset cognitive impairment and neuropathology associated with brain aging. Hum Mol Genet 2016;25:4186–200. 10.1093/hmg/ddw252
    1. Skelton MR, Schaefer TL, Graham DL, et al. . Creatine transporter (CrT; Slc6a8) knockout mice as a model of human CrT deficiency. PLoS One 2011;6:e16187. 10.1371/journal.pone.0016187
    1. Deming Y, Filipello F, Cignarella F, et al. . The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk. Sci Transl Med 2019;11:eaau2291. 10.1126/scitranslmed.aau2291

Source: PubMed

3
Subscribe