Nagahide Takahashi1,2,3
Tomoko Nishimura2,3
Taeko Harada2,3
Akemi Okumura2,3
Toshiki Iwabuchi2,3
Md. Shafiur Rahman2,3Hitoshi Kuwabara4Shu Takagai5
Yoko Nomura6Nori Takei2,3
Kenji J. Tsuchiya2,3*- 1Department of Child and Adolescent Psychiatry, Nagoya University Graduate School of Medicine, Nagoya, Japan
- 2Research Center for Child Mental Development, Hamamatsu University School of Medicine, Hamamatsu, Japan
- 3United Graduate School of Child Development, Hamamatsu University School of Medicine, Hamamatsu, Japan
- 4Department of Psychiatry, Hamamatsu University School of Medicine, Hamamatsu, Japan
- 5Department of Child and Adolescent Psychiatry, Hamamatsu University School of Medicine, Hamamatsu, Japan
- 6Queens College and Graduate Center, City University of New York, New York, NY, United States
Introduction: Obesity is highly heritable, and recent evidence demonstrates that obesity is associated with cognitive deficits, specifically working memory. However, the relationship between genetic risks for obesity and working memory is not clear. In addition, whether the effect of these genetic risks on working memory in children is mediated by increased body mass index (BMI) has not been elucidated.
Methods: In order to test whether the polygenic risk score (PRS) for obesity in adulthood (adulthood-BMI-PRS) is associated with working memory at 8 years of age, and whether the effect is mediated by childhood BMI, in children from the general population, participants in the Hamamatsu Birth Cohort for Mothers and Children (HBC) study in Hamamatsu, Japan, underwent testing for association of adulthood-BMI-PRS with working memory. HBC data collection began in December 2007 and is ongoing. Adulthood-BMI-PRS values were generated using summary data from the recent genome-wide association study (GWAS) undertaken in Japan, and the significance of thresholds was calculated for each outcome. Outcomes measured included the working memory index (WMI) of Weschler Intelligence Scale-4 (WISC-IV) scores and the BMI at 8 years of age. Gene-set enrichment analysis was conducted to clarify the molecular basis common to adulthood-BMI and childhood-WMI. Mediation analysis was performed to assess whether childhood-BMI of children mediated the association between adulthood-BMI-PRS and working memory.
Results: A total of 734 participants (377 males, 357 females) were analyzed. Adulthood-BMI-PRS was associated with lower childhood-WMI (β[SE], −1.807 [0.668]; p = 0.010, corrected) of WISC-IV. Gene-set enrichment analyses found that regulation of neurotrophin Trk receptor signaling (β[SE], −2.020 [6.39]; p = 0.002, corrected), negative regulation of GTPase activity (β[SE], 2.001 [0.630]; p = 0.002, corrected), and regulation of gene expression epigenetic (β[SE], −2.119 [0.664]; p = 0.002, corrected) were enriched in BMI in adulthood and WMI in childhood. Mediation analysis showed that there is no mediation effect of childhood-BMI between the adulthood-BMI-PRS and working memory deficits in children.
Conclusion: Adulthood-BMI-PRS was associated with working memory among children in the general population. These genetic risks were not mediated by the childhood-BMI itself and were directly associated with working memory deficits.
Introduction
Accumulating evidence shows that obesity is associated with deficits in neurocognitive functioning, such as deficits in working memory in children (Smith et al., 2011; Cornier et al., 2013; Crone and Steinbeis, 2017). Obesity is associated with adiposity, cellular stress, and excessive inflammation, all of which can lead to insulin resistance and cerebral structural alteration (Unamuno et al., 2018). Numerous neuroimaging studies have found structural alterations in obese children (Maayan et al., 2011; Veit et al., 2014; Sharkey et al., 2015; Saute et al., 2018); however, due to the limited sample size and methodological difference, results regarding the relationship between obesity, cortical thickness, and cognitive deficits have been inconsistent (Maayan et al., 2011; Sharkey et al., 2015; Saute et al., 2018). A recent large neuroimaging study reported that higher body mass index (BMI) was associated with lower working memory, and this association was mediated by reduced prefrontal cortex thickness in children (Laurent et al., 2019).
Although obesity is a complex phenotype, the role of genetic factors in the development of obesity remains undisputed. Additionally, a recent genome-wide association study (GWAS) identified 85 loci associated with BMI in adults from the Japanese general population (Akiyama et al., 2017). Furthermore, a previous study showed that BMI was genetically correlated with general cognitive function in adults using LD score regression analysis (rg = 0.51) (Marioni et al., 2016). However, Riggs et al. found that alteration in working memory is antecedent to weight gain in children, suggesting the possibility that cortical thickness or subsequent cognitive dysfunctions cause obesity (Riggs et al., 2010; Groppe and Elsner, 2017). Together, it is reasonable to hypothesize that children with genetic risks for obesity might have a higher risk for neurocognitive deficits compared to those with low genetic risks for obesity and the genetic risks for obesity potentially directly affect working memory without mediation of BMI.
As such, we examined (1) whether polygenic risk scores (PRSs) for BMI in adulthood (adulthood-BMI-PRS) are associated with deficits in working memory, leveraged on our birth cohort, composed of representative samples of Japanese, and (2) whether childhood-BMI mediated the association between adulthood-BMI-PRS and working memory.
Materials and Methods
Participants
Participants included infants (n = 832; 426 boys, 406 girls) born in Japan between December 2007 and June 2011. The recruitment procedures are described in detail in our previous study (Takagai et al., 2016). The study procedures were approved by the ethical committee. Written informed consent was obtained from each mother for the participation of her infant. Participants with parents of non-Japanese descent were excluded from the study (n = 8). No other screening, such as neurodevelopmental disorders or psychiatric disorders was conducted for the analysis. Hamamatsu University School of Medicine and the University Hospital Ethics Committee accepted the study methods.
Measurement
Working memory was assessed, using the working memory index (WMI) of the Wechsler Intelligence Scale for Children-4 (WISC-IV) when the children became 8 years of age. Information pertaining to BMI was obtained on the same day when the WISC-IV assessments were carried out.
Genotyping, Quality Control, and Imputation
Genotyping was conducted using the Japonica array designed specifically for single-nucleotide polymorphism (SNP) genotyping for a Japanese population (Kawai et al., 2015). The quality controls retaining SNPs and subjects were as follows: missing data for SNP < 0.02, Pi-hat calculated by identity-by-descent analysis < 0.2, SNP Hardy-Weinberg equilibrium of p > 10–6, and minor allele frequency > 0.01. Genotyping imputation was performed using BEAGLE 5.0 (Browning et al., 2018) to the Japanese reference panel phase 3 of 1000 Genome Project. SNPs with an imputation INFO score < 0.8 were excluded. We also excluded SNPs located within the MHC region, because of high linkage disequilibrium (LD) in this region. The number of SNPs analyzed for PRS was 5,606,655.
Statistical Analysis
Polygenic Risk Score Analysis
PRS was generated by PRSice-2 (Choi and O’Reilly, 2019) using a recent BMI-GWAS in the Japanese adult population as a discovery cohort1 (Akiyama et al., 2017). Four main components calculated with PLINK 1.9 (Chang et al., 2015) were used to account for population stratification. The criterion for SNP clumping was pairwise linkage disequilibrium of r2< 0.1 within a 1-Mb window. PRSs were calculated with different p-value thresholds: 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5. Standardized PRS scores (mean = 0; standard deviation = 1) were used for the analyses. p-values for WMI were corrected using 10,000 permutation tests. Sex and small-for-gestational-age (SGA) were included as covariates. The statistical power of the PRS at each p-value threshold was calculated using the AVENGEME R-package (Dudbridge et al., 2018). Since SNP-based heritability (h2 SNP) estimated from the original BMI study was 0.29812, all generated PRSs showed adequate power between 90 and 100%.
Gene-set enrichment analyses were conducted using PRSet (Choi et al., 2021) to identify gene sets that contain SNPs associated with both BMI in adulthood and WMI in childhood. Gene ontology (GO) sets (c5: biological process) were obtained from the MSigDB database2 and used for the analyses by PRSice-2 (Choi and O’Reilly, 2019). The p-value threshold for PRSet was set at 1, since gene-set PRSs containing a small portion of SNPs may be unrepresentative of the entire gene set (Fanelli et al., 2020). The p-values for PRSet were corrected by 10,000 permutation tests.
Mediation Analysis
Mediation analysis was performed to assess whether childhood-BMI mediated the association between adulthood-BMI-PRS and working memory. Best-fit PRSs computed with highest R2, obtained from linear regression analyses, were used for the medication analysis. The R package “lavaan” was used, and the significance of indirect effect of childhood-BMI was assessed by 1,000 bootstrap at a 95% confidential interval (Rosseel, 2012). The data were tested for normal distribution by a Shapiro–Wilk test.
Results
Association Between Adulthood-Body Mass Index-Polygenic Risk Score and Childhood-Working Memory Index
Participant characteristics are summarized in Table 1. Genotyping quality control and identity-by-decent analysis were used to remove 98 individuals from the analysis, resulting in a total of 726 participants (373 males, 353 females) for further analysis. The adulthood-BMI-PRS was significantly associated with lower WMI of WISC-IV at various p-value thresholds (Table 2).
Table 1. Sample characteristics: participating children and their parents.
Table 2. Association between BMI PRS and WMI of WISC-IV.
Gene-Set Enrichment Analysis of Adulthood-Body Mass Index-Polygenic Risk Score and Childhood-Working Memory Index
Gene-set enrichment analysis identified that several gene ontologies, such as regulation of neurotrophin Trk receptor signaling (β[SE], −2.020 [6.39]; p = 0.002, corrected), negative regulation of GTPase activity (β[SE], 2.001 [0.630]; p = 0.002, corrected), and regulation of gene expression epigenetic (β[SE], −2.119 [0.664]; p = 0.002, corrected), were enriched in BMI in adulthood and WMI in children (Table 3).
Table 3. Top 20 gene sets significantly enriched for WMI and BMI.
Mediation Analysis of Adulthood- Body Mass Index-Polygenic Risk Score, Childhood-Body Mass Index, and Childhood-Working Memory Index
All data including adulthood-BMI-PRS, childhood-BMI, and childhood-WMI were normally distributed by a Shapiro–Wilk test (p > 0.05). Mediation analysis was conducted for BMI-PRS (best fit p-value threshold at 0.319), childhood-BMI, and childhood-WMI. Analysis revealed that although there was a significant total effect of BMI-PRS on WMI (β[SE], −1.874 [0.633]; p = 0.004), no indirect effect of childhood-BMI (β[SE], 0.110 [0.102]; p = .280) was observed, indicating that most of the effect originated from the direct effect of BMI-PRS (β[SE], −1.736 [0.638]; p = 0.006) (Figure 1). Furthermore, regression analysis showed that there is no association between adulthood-BMI-PRS and childhood-BMI in the present study (β[SE], −0.185 [0.130]; p = 0.156).
Figure 1. Mediation analysis of Adulthood-BMI-PRS, Childhood-BMI, and Childhood-Working memory. The solid line indicates the path that was statistically significant and the dashed lines indicate those that were estimated but not statistically significant. BMI, body mass index; PRS, polygenic risk score. ∗∗p < 0.01.
Discussion
We report that genetic risks for obesity are linked to working memory deficits in children. This finding partially supports the previous report (Laurent et al., 2019); however, we found that the effect of BMI-PRS on WMI was not mediated by childhood-BMI, indicating that childhood-BMI itself does not affect this cognitive domain directly, but the genes involved in adulthood-BMI are responsible for these functions in the brain. This finding is also consistent with a previous study showing strong genetic correlation between BMI and general cognitive functions in adults using LD score regression analysis (Marioni et al., 2016).
Multiple reasons could be considered for the lack of association between adulthood-BMI-PRS and childhood-BMI itself in our study population. First, as we targeted 8-year-old children, it might be too early to detect the effect of genes related to adulthood-BMI on BMI in children. Second, there is a possibility that genes related to adulthood-BMI are different between children and adults.
In accordance with the previous GWAS results (Akiyama et al., 2017), GO sets commonly enriched in adulthood-BMI and childhood-WMI have been reported to be involved in the brain cortical maturation (Anton-Fernandez et al., 2015; Budday et al., 2015; Lian and Sheen, 2015). For example, genes in oligodendrocyte progenitor proliferation or layer formation in cerebral cortex have been reported to be expressed in the brain from the early postnatal period.3 Among these GOs, it is noteworthy that regulation of gene expression epigenetic was identified, since a recent GWAS demonstrated that a SNP on the HDAC4 gene was associated with selective attention (Pinar et al., 2018). Taken together, it can be considered that these genes are involved in brain maturation and working memory directly, without mediation of body weight.
Limitations
There are a few limitations in our study. First, compared to the previous study (Laurent et al., 2019), the average childhood-BMI [mean (SD), 18.64 (3.9) in the previous study vs. 16.2 (2.4) in this study] was lower in our cohort; thus, there is a possibility that we did not target a population with a broad range of BMI and, thus, could not detect an effect of adulthood-BMI-PRS on childhood-BMI in this study. Incidence of children classified as overweight (i.e., 85–95%) in this study was 9.68%, compared to 13.4% in the previous study. Similarly, the incidence of children classified as obese (i.e., > 95%) was 6.98% in this study, compared to 15.4% in the previous study (Laurent et al., 2019). Second, as BMI does not always reflect obesity, more sophisticated measurement to evaluate obesity, such as body fat percentage, is needed (Nickerson et al., 2018).
Conclusion
In this study, the adulthood-BMI-PRS was associated with working memory among children in the general population. These genetic risks were directly associated with working memory deficits, and not mediated by children’s BMI. Future studies are warranted in order to replicate these findings.
Data Availability Statement
The data generated for this study is subject to the following licenses/restrictions: Privacy and Confidentiality of Participants. Requests to access these datasets should be directed to KT, tsuchiya@hama-med.ac.jp.
Ethics Statement
The studies involving human participants were reviewed and approved by Hamamatsu University School of Medicine and the University Hospital Ethics Committee. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
NaT had full access to all the data used in the study and takes responsibility for the integrity of the data and accuracy of the data analysis, and statistical analysis. NaT, NoT, and KT study concept and design. NaT and KT drafting the manuscript. NoT and KT study supervision. TH, TN, and AO administrative, technical, and material support. All authors contributed significantly to the study and the creation of this manuscript, acquisition, analysis, interpretation of data, and critical revision of the manuscript for important intellectual content.
Funding
This work was supported by grants from the Ministry of Education, Culture, Sports, Science, and Technology in Japan (Grant No. 19H03582 to KT).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to the individuals who participated in the study. We would like to thank Ms. Chikako Nakayasu, Ms. Yuko Anma, and Ms. Haruka Suzuki for data collection, and Ms. Noriko Kodera and Ms. Emi Higashimoto for administration. This is a short text to acknowledge the contributions of specific colleagues, institutions, or agencies that aided the efforts of the authors.
Footnotes
- ^ https://humandbs.biosciencedbc.jp/hum0014-v23
- ^ https://www.gsea-msigdb.org/gsea/msigdb/index.jsp
- ^ informatics.jax.org
References
Akiyama, M., Okada, Y., Kanai, M., Takahashi, A., Momozawa, Y., Ikeda, M., et al. (2017). Genome-wide association study identifies 112 new loci for body mass index in the Japanese population. Nat. Genet. 49, 1458–1467. doi: 10.1038/ng.3951
Anton-Fernandez, A., Leon-Espinosa, G., DeFelipe, J., and Munoz, A. (2015). Changes in the golgi apparatus of neocortical and hippocampal neurons in the hibernating hamster. Front. Neuroanat. 9:157. doi: 10.3389/fnana.2015.00157
Browning, B. L., Zhou, Y., and Browning, S. R. (2018). A one-penny imputed genome from next-generation reference panels. Am. J. Hum. Genet. 103, 338–348. doi: 10.1016/j.ajhg.2018.07.015
Budday, S., Steinmann, P., and Kuhl, E. (2015). Physical biology of human brain development. Front. Cell. Neurosci. 9:257. doi: 10.3389/fncel.2015.00257
Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7. doi: 10.1186/s13742-015-0047-8
Choi, S. W., Garcia-Gonzalez, J., Ruan, Y., Wu, H. M., Johnson, J., O’Reilly, P., et al. (2021). The Power of Pathway-Based Polygenic Risk Scores. doi: 10.21203/rs.3.rs-643696/v1
Choi, S. W., and O’Reilly, P. F. (2019). PRSice-2: polygenic risk score software for biobank-scale data. Gigascience 8:giz082. doi: 10.1093/gigascience/giz082
Cornier, M. A., McFadden, K. L., Thomas, E. A., Bechtell, J. L., Eichman, L. S., Bessesen, D. H., et al. (2013). Differences in the neuronal response to food in obesity-resistant as compared to obesity-prone individuals. Physiol. Behav. 110–111, 122–128. doi: 10.1016/j.physbeh.2013.01.002
Crone, E. A., and Steinbeis, N. (2017). Neural perspectives on cognitive control development during childhood and adolescence. Trends Cogn. Sci. 21, 205–215. doi: 10.1016/j.tics.2017.01.003
Dudbridge, F., Pashayan, N., and Yang, J. (2018). Predictive accuracy of combined genetic and environmental risk scores. Genet. Epidemiol. 42, 4–19. doi: 10.1002/gepi.22092
Fanelli, G., Benedetti, F., Kasper, S., Kautzky, A., Zohar, J., Souery, D., et al. (2020). Higher polygenic risk scores for schizophrenia may be suggestive of treatment non-response in major depressive disorder. medRxiv [Preprint] doi: 10.1101/2020.01.15.20017699
Groppe, K., and Elsner, B. (2017). Executive function and weight status in children: a one-year longitudinal perspective. Child Neuropsychol. 23, 129–147. doi: 10.1080/09297049.2015.1089981
Kawai, Y., Mimori, T., Kojima, K., Nariai, N., Danjoh, I., Saito, R., et al. (2015). Japonica array: improved genotype imputation by designing a population-specific SNP array with 1070 Japanese individuals. J. Hum. Genet. 60, 581–587. doi: 10.1038/jhg.2015.68
Laurent, J. S., Watts, R., Adise, S., Allgaier, N., Chaarani, B., Garavan, H., et al. (2019). Associations among body mass index, cortical thickness, and executive function in children. JAMA Pediatr. 174, 170–177. doi: 10.1001/jamapediatrics.2019.4708
Lian, G., and Sheen, V. L. (2015). Cytoskeletal proteins in cortical development and disease: actin associated proteins in periventricular heterotopia. Front. Cell. Neurosci. 9:99. doi: 10.3389/fncel.2015.00099
Maayan, L., Hoogendoorn, C., Sweat, V., and Convit, A. (2011). Disinhibited eating in obese adolescents is associated with orbitofrontal volume reductions and executive dysfunction. Obesity (Silver Spring) 19, 1382–1387. doi: 10.1038/oby.2011.15
Marioni, R. E., Yang, J., Dykiert, D., Mottus, R., Campbell, A., Group, C. C. W., et al. (2016). Assessing the genetic overlap between BMI and cognitive function. Mol. Psychiatry 21, 1477–1482. doi: 10.1038/mp.2015.205
Nickerson, B. S., Esco, M. R., Bishop, P. A., Fedewa, M. V., Snarr, R. L., Kliszczewicz, B. M., et al. (2018). Validity of BMI-based body fat equations in men and women: a 4-compartment model comparison. J. Strength Cond. Res. 32, 121–129. doi: 10.1519/JSC.0000000000001774
Pinar, A., Hawi, Z., Cummins, T., Johnson, B., Pauper, M., Tong, J., et al. (2018). Genome-wide association study reveals novel genetic locus associated with intra-individual variability in response time. Transl. Psychiatry 8:207. doi: 10.1038/s41398-018-0262-z
Riggs, N., Chou, C. P., Spruijt-Metz, D., and Pentz, M. A. (2010). Executive cognitive function as a correlate and predictor of child food intake and physical activity. Child Neuropsychol. 16, 279–292. doi: 10.1080/09297041003601488
Rosseel, Y. (2012). lavaan: an R package for structural equation modeling. J. Stat. Softw. 48:36. doi: 10.18637/jss.v048.i02
Saute, R. L., Soder, R. B., Alves Filho, J. O., Baldisserotto, M., and Franco, A. R. (2018). Increased brain cortical thickness associated with visceral fat in adolescents. Pediatr. Obes. 13, 74–77. doi: 10.1111/ijpo.12190
Sharkey, R. J., Karama, S., and Dagher, A. (2015). Overweight is not associated with cortical thickness alterations in children. Front. Neurosci. 9:24. doi: 10.3389/fnins.2015.00024
Smith, E., Hay, P., Campbell, L., and Trollor, J. N. (2011). A review of the association between obesity and cognitive function across the lifespan: implications for novel approaches to prevention and treatment. Obes. Rev. 12, 740–755. doi: 10.1111/j.1467-789X.2011.00920.x
Takagai, S., Tsuchiya, K. J., Itoh, H., Kanayama, N., Mori, N., Takei, N., et al. (2016). Cohort profile: Hamamatsu birth cohort for mothers and children (HBC Study). Int. J. Epidemiol. 45, 333–342. doi: 10.1093/ije/dyv290
Unamuno, X., Gomez-Ambrosi, J., Rodriguez, A., Becerril, S., Fruhbeck, G., and Catalan, V. (2018)., Adipokine dysregulation and adipose tissue inflammation in human obesity. Eur. J. Clin. Invest. 48:e12997. doi: 10.1111/eci.12997
Keywords: polygenic risk score, obesity, cognition, GWAS, child development
Citation: Takahashi N, Nishimura T, Harada T, Okumura A, Iwabuchi T, Rahman MS, Kuwabara H, Takagai S, Nomura Y, Takei N and Tsuchiya KJ (2021) Association Between Genetic Risks for Obesity and Working Memory in Children. Front. Neurosci. 15:749230. doi: 10.3389/fnins.2021.749230
Received: 29 July 2021; Accepted: 26 August 2021;
Published: 22 September 2021.
Edited by:
Stefano Berto, University of Texas Southwestern Medical Center, United StatesCopyright © 2021 Takahashi, Nishimura, Harada, Okumura, Iwabuchi, Rahman, Kuwabara, Takagai, Nomura, Takei and Tsuchiya. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kenji J. Tsuchiya, tsuchiya@hama-med.ac.jp