Body mass index and childhood symptoms of depression, anxiety, and attention-deficit hyperactivity disorder: A within-family Mendelian randomization study
Abstract
Background:
Higher BMI in childhood is associated with emotional and behavioural problems, but these associations may not be causal. Results of previous genetic studies imply causal effects but may reflect influence of demography and the family environment.
Methods:
This study used data on 40,949 8-year-old children and their parents from the Norwegian Mother, Father and Child Cohort Study (MoBa) and Medical Birth Registry of Norway (MBRN). We investigated the impact of BMI on symptoms of depression, anxiety, and attention-deficit hyperactivity disorder (ADHD) at age 8. We applied within-family Mendelian randomization, which accounts for familial effects by controlling for parental genotype.
Results:
Within-family Mendelian randomization estimates using genetic variants associated with BMI in adults suggested that a child’s own BMI increased their depressive symptoms (per 5 kg/m2 increase in BMI, beta = 0.26 S.D., CI = −0.01,0.52, p=0.06) and ADHD symptoms (beta = 0.38 S.D., CI = 0.09,0.63, p=0.009). These estimates also suggested maternal BMI, or related factors, may independently affect a child’s depressive symptoms (per 5 kg/m2 increase in maternal BMI, beta = 0.11 S.D., CI:0.02,0.09, p=0.01). However, within-family Mendelian randomization using genetic variants associated with retrospectively-reported childhood body size did not support an impact of BMI on these outcomes. There was little evidence from any estimate that the parents’ BMI affected the child’s ADHD symptoms, or that the child’s or parents’ BMI affected the child’s anxiety symptoms.
Conclusions:
We found inconsistent evidence that a child’s BMI affected their depressive and ADHD symptoms, and little evidence that a child’s BMI affected their anxiety symptoms. There was limited evidence of an influence of parents’ BMI. Genetic studies in samples of unrelated individuals, or using genetic variants associated with adult BMI, may have overestimated the causal effects of a child’s own BMI.
Funding:
This research was funded by the Health Foundation. It is part of the HARVEST collaboration, supported by the Research Council of Norway. Individual co-author funding: the European Research Council, the South-Eastern Norway Regional Health Authority, the Research Council of Norway, Helse Vest, the Novo Nordisk Foundation, the University of Bergen, the South-Eastern Norway Regional Health Authority, the Trond Mohn Foundation, the Western Norway Regional Health Authority, the Norwegian Diabetes Association, the UK Medical Research Council. The Medical Research Council (MRC) and the University of Bristol support the MRC Integrative Epidemiology Unit.
Editor's evaluation
The manuscript uses genetic effects on BMI to test whether BMI affecxts childhood emotional and behavioural problems: symptoms of depression, anxiety, and attention-deficit and hyperactivity disorder (ADHD) at age 8. By using a within-family design in a large sample of children with genotyped parents in Norway, the study finds that previous estimates of the effect of BMI on childhood emotional and behavioural symptoms may have been overestimated due to confounding with the environment. Larger samples will be needed to determine whether there is a causal effect of BMI on childhood emotional or behavioural problems, and what size it is.
https://doi.org/10.7554/eLife.74320.sa0eLife digest
Some studies show that children with obesity are more likely to receive a diagnosis of depression, anxiety, or attention-deficit hyperactivity disorder (ADHD). But this does not necessarily mean obesity causes these conditions. Depression, anxiety, or ADHD could cause obesity. A child's environment, including family income or their parents' mental health, could also affect a child's weight and mental health. Understanding the nature of these relationships could help scientists develop better interventions for both obesity and mental health conditions.
Genetic studies may help scientists better understand the role of the environment in these conditions, but it's important to consider both the child's and their parents’ genetics in these analyses. This is because parents and children share not only genes, but also environmental conditions. For example, families that carry genetic variants associated with higher body weight might also have lower incomes, if parents have been affected by biases against heavier people in society and the workplace. Children in these families could have worse mental health because of effects of their parent’s weight, rather than their own weight. Looking at both child and adult genetics can help disentangle these processes.
Hughes et al. show that a child's own body mass index, a ratio of weight and height, is not strongly associated with the child’s mental health symptoms. They analysed genetic, weight, and health survey data from about 41,000 8-year-old children and their parents. The results suggest that a child's own BMI does not have a large effect on their anxiety symptoms. There was also no clear evidence that a child's BMI affected their symptoms of depression or ADHD.
These results contradict previous studies, which did not account for parental genetics. Hughes et al. suggest that, at least for eight-year-olds, factors linked with adult weight and which differ between families may be more critical to a child's mental health than a child’s own weight. For older children and adolescents, this may not be the case, and the individual’s own weight may be more important. As a result, policies designed to reduce obesity in mid-childhood are unlikely to greatly improve the mental health of children. On the other hand, policies targeting the environmental or societal factors contributing to higher body weights, bias against people with higher weights, and poor child mental health directly may be more beneficial.
Introduction
Children with high body mass index (BMI) have been found to have greater risk of emotional and behavioural problems, including symptoms and diagnoses of depression (Lindberg et al., 2020; Patalay and Hardman, 2019; Geoffroy et al., 2014; Quek et al., 2017) anxiety (Lindberg et al., 2020) and attention-deficit hyperactivity disorder (ADHD) (Cortese and Tessari, 2017; Griffiths et al., 2011). Prior to the COVID-19 pandemic, prevalence of childhood overweight and childhood obesity, respectively, was 21.3% and 5.7% in Europe (Garrido-Miguel et al., 2019) and 20.1% and 4.3% in Norway (Glavin et al., 2014). The estimated prevalence in Europe of mid-childhood emotional disorders was around 4% (Kovess-Masfety et al., 2016; Sadler et al., 2018) while the global prevalence of child and adolescent ADHD was estimated at 5% (Sayal et al., 2018). These rates may have increased considerably in the wake of the pandemic (Vizard et al., 2020). In this context, there is a clear need to understand the relationship between these factors, but it is not known if child body weight causes emotional or behavioural problems.
High BMI in childhood could affect emotional symptoms through social mechanisms, for example bullying victimization (Puhl et al., 2017). An impact on ADHD has been proposed via sleep disturbance and neurocognitive functioning (Vogel et al., 2015). However, even if children with high BMI are more likely than normal weight children to experience these symptoms, associations may not be causal. Aspects of the family environment may independently affect children’s BMI and their likelihood of developing emotional and behavioural symptoms, for example socioeconomic disadvantage (Russell et al., 2016) and parental mental health (Hope et al., 2019a; Hope et al., 2019b). Some studies have suggested that prenatal maternal obesity may confound associations of childhood BMI with emotional and behavioural symptoms (Sanchez et al., 2018), although the evidence is mixed (Li et al., 2020; Arafat and Minică, 2018). Reverse causality is also plausible: depressive, anxiety or ADHD symptoms could cause higher BMI, for instance via disordered eating patterns or decreased physical activity (Blaine, 2008; Martins-Silva et al., 2019). To avoid confounding and reverse causation, recent studies have applied Mendelian randomization (MR), a causal inference approach which uses genetic variants as instrumental variables for putative risk factors (Davies et al., 2018). Results, principally based on adult populations, are consistent with a causal influence of BMI on ADHD (Martins-Silva et al., 2019) and depression (Tyrrell et al., 2019). They are inconclusive for anxiety, reporting both positive (Walter et al., 2015) and negative (Millard et al., 2019) predicted causal effects of body weight.
However, although MR studies avoid classical confounding and reverse causation, they can be vulnerable to other sources of bias. Specifically, estimates from ‘classic’ MR studies – those conducted on samples of unrelated individuals – may be affected by demographic and familial factors (Davies et al., 2019; Morris et al., 2020). Bias can firstly arise from uncontrolled population stratification, where systematic differences in genotype between individuals from different ancestral clusters correlates with differences in environmental or cultural factors. This is an example of gene-environment correlation, which can lead to biased associations of genotypes and phenotypes. Secondly, indirect genetic effects may exist whereby parental genotype influences a child’s phenotype via environmental pathways, termed ‘dynastic effects’ or ‘genetic nurture’ (Kong et al., 2018). Thirdly, assortative mating in the parents’ generation, where parents are more (or less) similar to each other than would be expected by chance, can distort genotype-phenotype associations in the child’s generation. Recent work has suggested that these biases may be especially pronounced for complex social and behavioural phenotypes (Brumpton et al., 2020; Howe et al., 2022). Previously reported MR estimates of the effect of BMI on emotional and behavioural problems may therefore partly reflect demographic or familial biases rather than a causal influence of BMI. To investigate this, we used a ‘within-family’ Mendelian randomization (within-family MR) design. This approach uses the child’s, mother’s, and father’s genotype data as instruments for the BMI of the child, mother, and father. Within family Mendelian randomization estimates of the effect of the child’s BMI on the outcomes are robust to demographic and family-level biases. We compared within-family MR estimates with estimates from multivariable regression of the child’s outcomes on the child’s, mother’s and father’s reported BMI, and with estimates from ‘classic’ Mendelian randomization (classic MR), in which the child’s genotype data was used to instrument the child’s BMI without controlling for the parents’ genotype.
Methods
Study population
The Norwegian Mother, Father and Child Cohort Study (MoBa) is a population-based pregnancy cohort study over 114,500 children, 95,200 mothers, and 75,200 fathers conducted by the Norwegian Institute of Public Health (Magnus et al., 2016). Participants were recruited from all over Norway from 1999 to 2008, with 41% of all pregnant women invited consenting to participate. The first child was born in October 1999 and the last in July 2009. The cohort now includes over 114,500 children, 95,200 mothers, and 75,200 fathers (for more details see Appendix 1: MoBa study details). As of May 2022, genotype data which had passed quality control filters was available for 76,577 children, 53,358 fathers, and 77,634 mothers. This analysis was restricted to 40,949 mother-father-child ‘trios’ for whom genetic data were available for all three individuals, and at least one questionnaire had been completed.
The numbers of participants excluded are shown in a STROBE flow chart in Appendix 1—figure 1. From all records in MoBa (N=114,030 after removing consent withdrawals), participants were excluded if the parents had not completed any of the MoBa questionnaires used in imputation models. Of the 104,915 records remaining, there were 40,949 births for which genetic data were available and had passed QC filters for mother, father, and child (for details see Appendix 1: Genotyping and imputation, and Appendix 1: Genetic quality control). Missing values in phenotypic information for these participants were estimated using multiple imputation (details in Appendix 1: Multiple imputation). Related participants were retained, but all models were clustered by genetic family ID derived using KING software (Manichaikul et al., 2010). This genetic family ID groups first, second, and third-degree relatives (i.e. siblings in the parental generation and their children as well as nuclear families), in this way accounting for non-independence of observations.
Measures
Children’s BMI was calculated from height and weight values reported by mothers when the children were 8 years old. Maternal pre-pregnancy BMI was calculated from height and weight reported at ~17 weeks gestation. Father’s BMI was calculated from self-reported height and weight at ~17 weeks gestation. This information was missing from around 60% of fathers, and in these cases the mother’s report of the father’s height and weight was used instead (observed values of BMI from the two sources were correlated at 0.98). Values of height and weight more than 4 standard-deviations from the mean were treated as outliers and coded to missing.
Depressive, anxiety, and ADHD symptoms were reported by the mother when the child was 8 years old using validated measures. For depressive symptoms, the 13-item Short Mood and Feelings Questionnaire (SMFQ) was used, for anxiety symptoms the 5-item Short Screen for Child Anxiety Related Disorders (SCARED) (Birmaher et al., 1999) and for ADHD symptoms the Parent/Teacher Rating Scale for Disruptive Behaviour Disorders (RS-DBD) (total score and subdomain scores for inattention and hyperactivity) (Silva et al., 2005). Prorated summary scores were calculated for individuals with at least 80% of item-level information. Full details of all questions asked in MoBa are available at https://mobawiki.fhi.no/mobawiki/index.php/Questionnaires.
Blood samples were obtained from both parents during pregnancy and from mothers and children (umbilical cord) at birth. Details of genotyping and genetic quality control are described in Appendix 1: Genotyping and imputation and Appendix 1: Genetic quality control. Polygenic scores (PGS) for BMI were calculated using SNPs previously associated in GWAS with BMI at p<5.0 × 10–8 and weighted using the individual SNP-coefficients from the GWAS. We first constructed a PGS based on the largest existing GWAS of BMI in adults (Yengo et al., 2018). Since genetic influences on BMI in childhood and adulthood differ (Silventoinen et al., 2016) we also constructed a PGS based on a GWAS of body size in childhood as recalled by adult participants of UK Biobank (Richardson et al., 2020). These SNPs have been shown in external validation samples to predict BMI in childhood better than SNPs associated with adult BMI (Richardson et al., 2020; Brandkvist et al., 2021). From the full GWAS results, we excluded SNPs not available in MoBa, then used the TwoSampleMR package (Hemani et al., 2018b) to identify SNPs independently associated with BMI (with a clumping threshold of r=0.01, LD = 10,000 kb) at p<5.0 × 10–8. This left 954 SNPs associated with adult BMI, and 321 associated with childhood body size. Full details of SNPs included in both PGSs are provided in Supplementary file 1a and b. Equivalent PGSs were derived for depression and ADHD based on SNPs previously associated with these conditions at p<5.0 × 10–8 in GWAS (Wray et al., 2018; Demontis et al., 2019). This was not possible for anxiety, due to few known SNPs associated with these traits at p<0.05 × 10–8. Details of the SNPs in the depression and ADHD PGSs are provided in Supplementary file 1c and d.
Statistical analysis
Among trios with genetic data, multiple imputation by chained equations was performed in STATAv16 to estimate missing phenotypic information (details in Appendix 1: Multiple imputation of phenotypes). We used non-genetic linear regression, classic MR, and within-family MR to estimate the effects of the child’s BMI on the following outcomes: depressive, anxiety, and ADHD symptoms, and subdimensions of ADHD (inattention and hyperactivity). Non-genetic regression models were adjusted for child’s sex, year of birth, mother’s and father’s BMI, and likely confounders of observational associations: mother’s and father’s educational qualifications, mother’s and father’s depressive/anxiety symptoms (using selected items from the 25-item Hopkins Checklist Hesbacher et al., 1980) and ADHD symptoms (from the 6-item adult ADHD self-report scale Kessler et al., 2005), mother’s and father’s smoking status during pregnancy, and maternal parity at the child’s birth. For comparability, these models also included all covariates included in genetic models: genotyping centre, genotyping chip, and 20 principal components of ancestry for the child, mother, and father (for detailed information on principal components see Appendix 1: Genetic quality control). All MR models were conducted with two-stage least squares instrumental-variable regression using Stata’s ivregress, with F-statistics and R2 values obtained using ivreg2. Classic MR models, which do not account for parental genotype, used the child’s own PGS but not those of the parents to instrument the child’s BMI. Within-family MR models were multivariable MR models, in which we used PGSs for all members of a child-mother-child trio to instrument the BMI of all three individuals (model equations are provided in Appendix 1: Model equations). Classic and within-family MR models were adjusted for the child’s sex and year of birth, and the genotyping centre, genotyping chip, and the first 20 principal components of ancestry for the child, mother, and father. Given skew in outcomes variables, all models used robust standard errors (Stata’s vce option) and thus made no assumptions about the distribution of outcomes. We report two sets of results, in which either the adult BMI GWAS, or the childhood body size GWAS, was used to create the BMI PGS for the child, mother, and father. Z tests of difference were used to formally compare the classic MR and within-family MR estimates. To assess the extent of assortative mating in the parental generation based on phenotype data, we ran linear regression models of standardized paternal BMI, depressive symptoms, and ADHD symptoms on standardized maternal BMI, depressive symptoms, and ADHD symptoms. We then regressed paternal polygenic scores for BMI, depression, and ADHD on maternal polygenic scores for BMI, depression, and ADHD. All models investigating assortative mating adjusted for both parents’ principal ancestry components and genotyping covariates. We did not examine correlations with polygenic scores for anxiety, due to few known SNPs associated with these traits at p<0.05 × 10–8. All statistical tests were two-tailed.
Sensitivity analyses
To check sensitivity of results to outliers, all analyses were repeated using log-transformed versions of outcome measures (as all symptoms scales began at 0, we added 1 to scores before log-transforming). Genetic studies designed to assess causation can be biased by horizontal pleiotropy (Davies et al., 2018). This is when genetic variants in a polygenic score influence the outcome via pathways which do not involve the exposure. Pleiotropic effects can inflate estimated associations, or bias estimates towards the null. Methods have been developed to test for the presence of horizontal pleiotropy by comparing SNP-specific associations of exposures and outcomes, although these tests themselves rest on assumptions (Hemani et al., 2018a). We therefore performed additional robustness checks based on associations of individual SNPs included in the polygenic scores with BMI in the GWAS, and associations of the same SNPs with each outcome in MoBa. It was not computationally feasible to include individual SNPs in the imputation models, so SNP-outcome associations in MoBa were calculated using unimputed SNP data with imputed outcome data. For robustness checks of classic MR models, SNP-outcome associations were adjusted for the child’s sex and birth year, and the genotyping centre, genotyping chip, and ancestry principal components of the child, mother, and father. For robustness checks of within-family MR models, SNP-outcome associations were adjusted for the child’s sex and birth year, mother’s and father’s genotype, and the genotyping centre, genotyping chip, and principal components of the child, mother, and father. We conducted inverse-variance weighted, MR-Median, MR-Mode, and MR-Egger regression in STATAv16 with the MRRobust package (Spiller et al., 2019). A non-zero intercept from an MR-Egger model indicates presence of horizontal pleiotropy. We repeated main analyses without using imputed data in the sample of participants who had full genetic, exposure, outcome, and covariate data. To explore nonlinearities in associations of BMI with depression, anxiety, and ADHD symptoms, we ran non-genetic models with the child’s BMI divided into quintiles. Finally, MR models were run with additional adjustment for parental education. Attenuation of classic MR estimates in these models would be consistent with confounding by aspects of the family environment linked to parental education.
Results
This analysis was restricted to 40,949 mother-father-child ‘trios’ for whom genetic data were available for all three individuals, and at least one questionnaire had been completed. To assess whether participants included in the analytic sample (N=40,949) differed from the rest of the MoBa sample (N=72,742), we conducted t-tests and chi-squared tests for key characteristics at birth, BMI, and outcomes using unimputed data. There were modest differences, described in Appendix 1: Comparison of analytic sample and excluded participants. BMI did not differ for mothers, fathers or children, but children in the analytic sample had slightly lower depressive symptoms (mean SMFQ = 1.81 vs 1.91), anxiety symptoms (mean SCARED = 1.04 vs 1.00) and ADHD symptoms (mean RS-DBD ADHD = 8.4 vs 8.7). Descriptive characteristics of the full MoBa sample are in Appendix 1—table 1.
Descriptive statistics of the analytic sample after multiple imputation is presented in Table 1. The mean BMI for children was 16.3 (SD = 2.0), for mothers 24.0 (SD = 4.1), and for fathers 25.9 (SD = 3.2). Corresponding descriptive characteristics from unimputed data are included in Appendix 1—table 2. Both polygenic scores used to instrument BMI were strong instruments, even when used in within-family models. For the adult BMI PGS, conditional first-stage F-statistics for children, mothers, and fathers were 718.7, 1338.2, and 1272.5. The conditional R2 showed that the score explained 1.7%, 3.2%, and 3.0% of the variation in BMI for children, mothers, and fathers respectively. For the childhood body size PGS, conditional first-stage F-statistics were 919.8, 1071.8, and 960.2 for children, mothers, and fathers, with the scores explaining 2.2%, 2.6% and 2.3% of the variation in BMI. The correlation of the polygenic scores for adult BMI and for childhood body size was 0.38 for children, 0.36 for mothers and 0.37 for fathers.
Associations of BMI with depressive, anxiety, and ADHD symptoms at age 8
Depressive symptoms (SMFQ)
In adjusted non-genetic regression models (Figure 2, Appendix 1—table 3), children’s higher BMI at age 8 was associated with slightly higher depressive symptoms. Per 5 kg/m2 increase in BMI, SMFQ score was 0.05 standard deviations (SD) higher (95% CI: 0.01,0.09, p=0.02). Classic MR using the adult BMI PGS suggested that for each 5 kg/m2 increase in the child’s BMI, the child’s SMFQ score increased by 0.45 SD (95% CI: 0.26,0.64, p<0.001). Within-family MR using the adult BMI PGS also provided some evidence for an effect (beta: 0.26 SD, 95% CI: –0.01,0.52, p=0.06), but the within-family MR estimate was less precise (70% as precise as the classic MR estimate, and 15% as precise as the OLS estimates, from the ratio of standard errors). A z test for the difference (p=0.26) indicated that the within-family MR estimate was consistent with the classic MR estimate. Using the childhood body size PGS (Figure 3, Appendix 1—table 4) there was little evidence that a child’s own BMI affected their depressive symptoms from either classic MR (beta: 0.08 (95% CI: –0.07,0.22, p=0.29) or within-family MR (beta: 0.02 (95%CI: –0.20,0.23, p=0.88). In summary, evidence for an effect of childhood BMI on depressive symptoms was strongest using the genetic variants for adult BMI.
Anxiety symptoms (SCARED)
In non-genetic models (Figure 2, Appendix 1—table 3), each 5 kg/m2 increase in BMI was associated with a 0.07 SD lower (95% CI: −0.11,–0.03, p=0.001) SCARED score. Using the adult BMI PGS, there was little evidence for an effect from classic MR (beta: –0.06, 95% CI: –0.25,0.12, p=0.51), or within-family MR models (beta: 0.01, 95% CI: –0.25,0.29, p=0.96). Again, the within-family MR estimate was less precise than the classic MR estimate (68% as precise), or the OLS estimate (15% as precise), and the classic and within-family MR estimates were consistent (p=0.54). Using the childhood body size PGS (Figure 3, Appendix 1—table 4), MR estimates were similar (classic MR beta: –0.04, 95%: –0.18,0.11, P=0.62, within-family MR beta: 0.02, 95% CI: –0.18,0.22, P=0.83). In summary, there was little evidence from any genetic model that childhood BMI affects anxiety symptoms.
ADHD symptoms (RS-DBD)
In non-genetic models (Figure 2, Appendix 1—table 3) children’s BMI was negatively associated with ADHD symptoms after adjusting for confounders. Per 5 kg/m2 increase in BMI, ADHD symptoms from the RS-DBD were 0.07 SD lower (95% CI: −0.11,–0.03, p=0.001), with similar associations observed for the inattention or hyperactivity subscales (Figure 2, Appendix 1—table 3). Using the adult BMI PGS there was evidence from both classic and within-family MR models of a positive association of BMI and ADHD. In the classic MR model ADHD symptoms were 0.35 SD higher (95% CI: 0.17,0.53, p<0.001) per 5 kg/m2 increase in BMI; the within-family MR estimate, at 0.36 SD (CI: 0.09,0.63, p=0.009) was almost identical (p for difference = 0.95). A similar pattern was seen with the inattention and hyperactivity subscales (Figure 2, Appendix 1—table 3). The within-family MR estimate was again the least precise (65% as precise as the classic MR estimate, 14% as precise as the non-genetic estimate). Using the childhood body size PGS (Figure 3, Appendix 1—table 4) there was little evidence of an association from either classic MR (beta: –0.07, 95% CI: –0.21,0.07, p=0.35) or within-family MR models (beta: –0.03, 95% CI: –0.22,0.17, p=0.80). Thus, as for depressive symptoms, evidence for an effect of childhood BMI on ADHD symptoms was inconsistent and only detected using the adult BMI polygenic score.
Association of mother’s and father’s BMI with child’s symptoms
In non-genetic models which adjusted for the child’s BMI as well as covariates, the mother’s BMI was associated with slightly more depressive symptoms in the child (the child’s SMFQ score was 0.05 SD higher (95% CI: 0.03,0.07, p<0.001), per 5 kg/m2 increase in maternal BMI). Maternal BMI was also associated with more ADHD symptoms in the child: the child’s RS-DBD score was 0.04 S.D. higher (95% CI: 0.02,0.06, p<0.001) per 5 kg/m2 increase in maternal BMI, with similar associations for inattention and hyperactivity subscales. No such associations were seen with paternal BMI.
Within-family MR models also provide estimates for the effect of factors linked to maternal and paternal BMI on child outcomes, conditional on the child’s own BMI. However, compared to within-family MR estimates for the child’s own genotype, the interpretation of parental estimates differs. Like classic MR estimates for the child’s BMI, within-family MR estimates for each parent’s BMI will capture the causal effect of the parent’s BMI on the child’s outcome, but can also reflect residual population stratification and assortative mating in the parents’ generation or earlier. For an unbiased estimate of parental effects, we would need to account for grandparental genotype. Within-family MR models provided inconsistent consistent evidence that maternal BMI affected the child’s depressive symptoms: using the adult BMI PGS (Figure 2, Appendix 1—table 3), estimates suggested that higher maternal BMI increased depressive symptoms in the child (0.11 SD higher SMFQ score (95% CI: 0.02,0.19, p=0.01) per 5 kg/m2 increase in maternal BMI), but within-family MR models using the childhood body size PGS did not (Figure 3, Appendix 1—table 4). There was little evidence from within-family MR of other maternal or paternal effects on the child’s emotional or behavioural outcomes.
In the parents’ generation, phenotypes were associated within parental pairs, consistent with assortative mating on these traits (Appendix 1—table 5). Adjusted for ancestry and other genetic covariates, maternal and paternal BMI were positively associated (beta: 0.23, 95% CI: 0.22,0.25, p<0.001), as were maternal and paternal depressive symptoms (beta: 0.18, 95% CI: 0.16,0.20, p<0.001), and maternal and paternal ADHD symptoms (beta: 0.11, 95% CI: 0.09,0.13, p<0.001). Consistent with cross-trait assortative mating, there was an association of mother’s BMI with father’s ADHD symptoms (beta: 0.03, 95% CI: 0.02,0.05, p<0.001) and mother’s ADHD symptoms with father’s depressive symptoms (beta: 0.05,95% CI: 0.05,0.06, p<0.001). Phenotypic associations can reflect the influence of one partner on another as well as selection into partnerships, but regression models of paternal polygenic scores on maternal polygenic scores also pointed to a degree of assortative mating. Adjusted for ancestry and genotyping covariates, there were small associations between parents’ BMI polygenic scores (beta: 0.01, 95% CI: 0.00,0.02, p=0.02 for the adult BMI PGS, and beta: 0.01, 95% CI: 0.00,0.02, p=0.008 for the childhood body size PGS), and of the mother’s childhood body size PGS with the father’s ADHD PGS (beta: 0.01, 95% CI: 0.00,0.02, p=0.03). We did not detect associations with pairs of other polygenic scores, which may be due to insufficient statistical power.
Sensitivity analyses
Analyses using log-transformed versions of the outcomes (Appendix 1—Tables 6 and 7) were consistent with main results. Robustness checks based on comparing associations of individual SNPs with BMI in the GWAS and with children’s outcomes in MoBa (Appendix 1—Tables 8 and 9) were consistent with the main results. MR-Egger models found little evidence of horizontal pleiotropy, although MR-Egger estimates were imprecise (Appendix 1—Tables 8 and 9). Results of analyses using the complete-case sample were qualitatively similar to results using imputed data (Appendix 1—Tables 10 and 11). In non-genetic models where the child’s BMI was divided into quintiles (Appendix 1—table 12), there was little evidence of nonlinear associations. With additional adjustment for parental education, point estimates for depressive and ADHD symptoms in classic MR models were closer to the null, but confidence intervals substantially overlapped (Appendix 1—Tables 13 and 14).
Discussion
In a large cohort of Norwegian 8-year-olds, higher childhood BMI was phenotypically associated with slightly more depressive symptoms, but fewer anxiety symptoms and ADHD symptoms. Genetic analyses using the adult BMI PGS suggested that higher BMI in childhood increased symptoms of both depression and ADHD. This was clearest in classic MR models, but also suggested by within-family MR models, whose precision is lower but which account for parental genotype. Compared to associations from non-genetic models, effect sizes for depression and ADHD from genetic models based on the adult BMI PGS were larger. However, these estimates were less precise, and confidence intervals for the classic MR and within-family MR estimates substantially overlapped for all outcomes. The childhood body size PGS explained more variation in children’s BMI than the adult BMI PGS did, consistent with other studies (Richardson et al., 2020; Brandkvist et al., 2021), while the adult BMI PGS explained more variation in maternal and paternal BMI. Genetic analyses which used the childhood body size SNPs provided little evidence that the child’s BMI affected their depressive or ADHD symptoms outcomes. This suggests that genetic variation associated with adult BMI has a greater impact on these outcomes than genetic variation associated with recalled childhood body size. This is consistent with the moderate correlation observed between the two polygenic scores, indicating that they capture both overlapping and unique variation. Our results may therefore reflect differences in how each set of SNPs relate to traits other than childhood BMI which are relevant to a child’s depressive and ADHD symptoms. Nevertheless, within-family MR estimates using the childhood body size PGS were still consistent with small effects of the child’s BMI on all outcomes, with upper confidence limits around a 0.2 standard-deviation increase in each outcome per 5 kg/m2 increase in BMI. There was little evidence that maternal or paternal BMI affected a child’s ADHD or anxiety symptoms. In within-family MR models using the adult BMI PGS, but not the childhood body size PGS, maternal BMI was positively associated with children’s depressive symptoms. This is consistent with a causal impact of the mother’s recent BMI but not their BMI in childhood, but it may also reflect family-level biases from previous generations.
The positive association between BMI and depressive symptoms in non-genetic models accords with previous observational studies (Lindberg et al., 2020; Patalay and Hardman, 2019; Quek et al., 2017; Geoffroy et al., 2014). The inverse association between BMI and anxiety symptoms in non-genetic models contrasts with the results of a recent study, in which Swedish 6–17 year olds receiving treatment for obesity had a greater likelihood of a diagnosis or prescription for anxiety disorder compared to controls (Lindberg et al., 2020). The discrepancy may reflect confounding (we adjusted for more factors, including parental BMI), age of the participants (children in our study were younger) or differences in the outcome or exposure, since we considered anxiety symptoms rather than diagnosis, and a continuous BMI measure rather than obesity. However, anxiety symptoms in our sample were not raised in the top BMI quintile. Another difference concerns the population: children receiving obesity treatment may be more likely than other children with obesity to experience anxiety symptoms or to receive a diagnosis. The inverse association between BMI and ADHD symptoms in non-genetic models contrasts with previous reports of positive or null associations with obesity, which typically adjusted for fewer confounders (Cortese and Tessari, 2017; Nigg et al., 2016). Since previous studies have found more evidence of an association in adults than children, and often considered ADHD diagnoses rather than symptoms, the discrepancy may also point to age-varying associations, or to different influences on likelihood of diagnosis compared to parent-reported symptoms (Nigg et al., 2016; Cortese and Tessari, 2017).
For depressive symptoms and ADHD, classic and within-family MR estimates using the adult BMI PGS were larger than estimates from non-genetic models. Horizontal pleiotropy, which we could not rule out, could have inflated MR estimates. It could also help explain the discrepancy in results using the adult BMI and childhood body size polygenic scores, if SNPs in the adult BMI polygenic score have a greater impact on depressive or ADHD symptoms via pathways independent of childhood BMI. We found little evidence of pleiotropy using MR-Egger estimators, but the power to detect pleiotropy with this method is low. Additionally, classic MR estimates may be inflated by demographic and familial factors, but within-family MR estimates for effects of a child’s own BMI are robust to these factors. For depressive symptoms, the within-family MR estimate was closer to the non-genetic estimate than the classic MR estimate, which may reflect bias in the classic MR estimate due to demographic and familial factors. At the same time, the within-family MR estimate was imprecise, and confidence limits consistent with a substantial effect of children’s BMI on depressive symptoms. For ADHD, point estimates from the classic MR and within-family MR models using the adult BMI PGS were very similar, and both statistically distinguishable from the null. These results therefore accord with a recent study which accounted for family-level biases by using dizygotic twin pairs, obtaining between-family and within-family estimates for the effect of BMI on ADHD symptoms (Liu et al., 2021). Using a PGS of SNPs associated with adult BMI, within-family analysis found a 0.07 S.D. increase in ADHD symptoms at age 8 per S.D. increase in BMI PGS, which was consistent with the between-family estimate. The between-family estimate was attenuated by adjustment for parental education, suggesting an influence of family-level processes. In our study, classic MR estimates for depressive and ADHD symptoms which adjusted for parental education were consistent with the main results, with largely overlapping confidence intervals, although point estimates were closer to the null. Thus, our results are also consistent with an influence of demographic or family-level effects, and with earlier evidence that such processes impact the relationship between BMI and ADHD (Chen et al., 2014; Geuijen et al., 2019).
Several sources of genetic familial bias may have influenced classic MR estimates of the impact of the child’s own BMI. Firstly, frequencies of BMI-associated variants may differ between sub-populations in a similar manner to environmental influences on emotional or behavioural functioning (population stratification). Such gene-environment correlation can inflate estimates from classic MR models, but are unlikely to affect within-family MR models, where ancestry is fully controlled for via parental genotypes. Although we included principal components of ancestry in all models, residual population stratification may nevertheless have influenced the classic MR results. Secondly, there may be indirect effects of parental BMI via the family environment (dynastic effects, or genetic nurture). This could explain the association of maternal BMI with children’s depressive symptoms in the within-family MR model using the adult BMI PGS. In observational studies, maternal pre-pregnancy obesity is linked with children’s risk of emotional disorders and ADHD (Sanchez et al., 2018). Although mechanisms are not well understood, an in utero effect on children’s neurodevelopment of metabolic correlates of obesity has been proposed (Edlow, 2017). Our within-family MR results suggest that previously reported associations of maternal BMI with a child’s ADHD are not causal, but are consistent with an effect on the child’s depressive symptoms. This could reflect an impact of maternal BMI later in the child’s life. A well-documented ‘wage penalty’ exists for high BMI (Howe et al., 2020), especially for women (Bozoyan and Wolbring, 2018) reflecting social consequences of obesity being a stigmatized condition (Giel et al., 2010). High BMI in adulthood is also linked to worse mental health, with stronger associations for women again pointing to gendered social processes (Rubino et al., 2020). Maternal BMI may therefore influence children’s emotional and behavioural problems via economic consequences, or via maternal mental health, throughout childhood. However, while our results are consistent with an influence of maternal BMI on child’s depressive symptoms, these results should be interpreted with caution. In contrast to estimated effects for the child’s BMI, where controlling for parental genotype is likely to eliminate familial biases, estimated maternal and paternal effects from within-family MR models may have been impacted by familial biases in previous generations. Adjustment for grandparental genotype would be required to obtain similarly unbiased estimates for the parents. Thirdly, people with high BMI may be more likely to partner with people with emotional or behavioural conditions (cross-trait assortative mating). Over generations, this would induce an association of not only the phenotypes but of associated genetic variants. We found some genomic evidence of assortative mating for BMI, and cross-trait assortative mating between BMI and ADHD, but not between other traits. However, associations between polygenic scores, which only capture some of the genetic variation associated with these phenotypes, may not capture the full extent of genetic assortment on these traits.
Despite a high participation rate, MoBa is not perfectly representative, and selection biases linked to participation could have affected our results. The current analyses were restricted to families with complete genetic data and at least some relevant questionnaire data. These families were found to have slightly more years of education than the wider MoBa sample, and the children to score slightly lower for depressive, anxiety, and ADHD symptoms. Reflecting the requirement of genetic data for fathers, single mothers were under-represented. Analyses were restricted to individuals of European ancestry, with polygenic scores based on results of GWAS which were also restricted to individuals of European ancestry. Consequently, our results may not be generalisable to other populations. Outcomes were based on mother-reported symptoms of depression, anxiety disorders and ADHD, and estimates based on diagnoses may have differed. However, a child’s sociodemographic characteristics can influence their likelihood of diagnosis independently of symptoms (Thompson et al., 2021), indicating that such an approach is not always preferable. BMI measurements were based on reported height and weight, so reporting bias may have influenced relationships. In many families, fathers’ BMI was based on height and weight reported by the mothers. However, these measures were very highly correlated with father’s self-reports, so additional measurement error is unlikely to have greatly affected our results for father’s BMI. Due to attrition, a substantial proportion of values for the child’s BMI and outcomes were imputed, and we cannot be sure that observations were missing at random conditional on variables included in imputation models. Effects of parental BMI may be time-varying, for example a parent’s own BMI during childhood could influence their child independent of the parent’s later BMI. We could not explore these effects because information on parent’s childhood BMI was not available. Within-family MR may still be affected by horizontal pleiotropy, and recent genetic work points to genetic overlap between BMI and psychiatric disorders including major depression (Bahrami et al., 2020). While robustness checks found little evidence of pleiotropy, these methods rely on assumptions. Moreover, MR-Egger is known to give imprecise estimates (Burgess and Thompson, 2017), and confidence intervals from MR-Egger models were wide. Thus, pleiotropy cannot be ruled out. The Mendelian randomization methods employed here assume any causal impact of BMI is linear – that a kg/m2 increase in BMI will have the same impact regardless of the child’s initial BMI. There is substantial evidence for a ‘J-shaped’ phenotypic association of BMI with common mental disorders, consistent with an impact of both high and low BMI on risk of depression or anxiety (McCrea et al., 2012; Gaysina et al., 2011; Geoffroy et al., 2014). Genetic methods exist for exposures with nonlinear effects but require much larger samples (Sun et al., 2019). If there exist nonlinear effects of BMI on mental health, rather than vice versa, our results may underestimate the effects of high BMI. Finally, the effects of BMI on emotional and behavioural functioning likely differ by age, and relationships may be substantially different for older children or adolescents. In particular, depressive symptoms do not tend to occur until the teenage years (Kwong et al., 2019) and observational associations of BMI and ADHD become clearer with age (Nigg et al., 2016). Work in larger samples of related individuals will be needed to precisely estimate the influence of a child’s BMI on their emotional and behavioural outcomes. In response to a reviewer’s request, we conducted post-hoc power calculations to estimate the minimum effect on the outcomes of the child’s BMI which could be detected with 80% power in a dataset of this size (40,949 trios). Simulations indicated that, using the adult BMI PGS, an effect on each outcome of 0.15 S.D. per 5 kg/m2 could be detected using classic MR, and an effect on each outcome of 0.22 S.D. per 5 kg/m2 using within-family MR. Using the childhood body size PGS, equivalent detectable effects were 0.12 S.D. and 0.16 S.D. per 5 kg/m2. MoBa is currently the largest individual study in which this approach can be applied, but new data are becoming available which will allow analyses of this kind within and across studies, such as through the Within Family Consortium https://www.withinfamilyconsortium.com/home/. Meanwhile, studies with extensive intergenerational information will be needed to fully explore mechanisms linking child outcomes to maternal BMI.
Conclusion
Our results suggest that genetic variation associated with BMI in adulthood affects a child’s depressive and ADHD symptoms, but genetic variation associated with recalled childhood body size does not substantially affect these outcomes. There was little evidence that BMI affects anxiety. However, our estimates were imprecise, and these differences may be due to estimation error. There was little evidence that parental BMI affects a child’s ADHD or anxiety symptoms, but factors associated with maternal BMI may independently influence a child’s depressive symptoms. Genetic studies using unrelated individuals, or polygenic scores for adult BMI, may have overestimated the causal effects of a child’s own BMI.
Appendix 1
Additional methods
1. MoBa study details
This study is based on the Norwegian Mother, Father and Child Cohort Study (MoBa) and uses data from the Medical Birth Registry of Norway (MBRN). The Medical Birth Registry (MBRN) is a national health registry containing information about all births in Norway. The current analysis is based on version 12 of the quality-assured data files released for research in January 2019. The establishment of MoBa and initial data collection was based on a license from the Norwegian Data Protection Agency and approval from The Regional Committees for Medical and Health Research Ethics. The MoBa cohort is now based on regulations related to the Norwegian Health Registry Act. The current study was approved by The Regional Committees for Medical and Health Research Ethics (2016/1702).
The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. We are grateful to all the participating families in Norway who take part in this on-going cohort study.
2. Genotyping
Genotyping of MoBa has been conducted through multiple research projects, spanning several years. The research projects (HARVEST, SELECTIONpreDISPOSED, and NORMENT) provided genotype data to MoBa Genetics. In total, 238,001 MoBa samples were sent to be genotyped in 24 genotyping batches. This was carried out at 3 centres (1. Genomics Core Facility, Trondheim, Norway, 2. ERASMUS MC, Rotterdam, Netherlands, and 3. deCODE Genetics, Reykjavik, Iceland) using six genotyping arrays. The 24 batches had varying selection criteria; this included a batch of ADHD child cases and their parents, and another of matched control children and their parents. Detailed information on batch selection and the genotyping process are described elsewhere (Corfield et al., 2022).
3. Genetic quality control
Pre-imputation QC, phasing, imputation, and post-imputation QC were carried out according to the MoBaPsychGen pipeline, which includes QC on both single nucleotide polymorphism (SNP) and individual level, and whose full details are provided elsewhere (Corfield et al., 2022). Output of the pipeline included 207,569 unique individuals and 6,981,748 SNPs after post-imputation QC; unique individuals comprised 76,577 children, 53,358 fathers, and 77,634 mothers. Phasing and imputation were performed using the publicly available Haplotype Reference Consortium release 1.1 panel as a reference. Information from the Medical Birth Registry of Norway and MoBa questionnaires were used to identify biological sex, year of birth, reported parent-offspring (PO) relationships, and, in the offspring generation, multiple births. Ancestry outliers were identified based on principal component analysis with the 1000 Genomes phase 1 unrelated data (1,083 individuals) (Auton et al., 2015). Approximately 95% of the participants were identified as having European ancestry. Ancestry outliers were removed based on visual inspection, using pairwise plots for the first seven principal components. Similarly, principal component analysis was used to identify substructure within each subpopulation of all MoBa batches. PCs were first estimated in founders only, and non-founders projected into the PC space of the founders. Outliers were removed based on visual inspection, using pairwise plots for the first ten principal components.
Relatedness was inferred using KING version 2.2.5 (Manichaikul et al., 2010). In the presence of admixture, KING accurately infers MZ twin or duplicate pairs (kinship coefficient >0.3540), first-degree (PO, FS, DZ twin pairs; kinship coefficient range 0.1770–0.3540), second-degree (HS, GO, AUNN; kinship coefficient range 0.0884–0.1770), and third-degree (first cousins; kinship coefficient range 0.0442–0.0884) relationships. After within-family and between-family relationships were confirmed by genetic data, a check for Mendelian errors (ME) was performed in PLINK. The ME check included families with one or two parents present in the data. Families with >5% errors and SNPs with >1% errors were removed. The remaining ME were set to missing.
4. Multiple imputation
Multiple imputation by chained equations was performed in STATAv16 to estimate missing phenotypic information for the 40,949 trios with complete genetic data. 100 imputed datasets were produced and analysis across these datasets conducted with STATA’s mi estimate commands. The imputation model included all BMI variables used in the main analyses (child’s BMI at age 8, mother’s pre-pregnancy BMI and father’s BMI as reported at 17 weeks gestation), the child’s sex and year of birth, and other phenotypic covariates used in non-genetic models, including mother’s and father’s smoking status reported at 17 weeks gestation, mother’s and father’s depressive/anxiety symptoms [using selected items from the 25-item Hopkins Checklist (Hesbacher et al., 1980)], and ADHD symptoms [from the 6-item adult ADHD self-report scale (Kessler et al., 2005)], maternal parity at the child’s birth, and family socioeconomic characteristics, including parental educational qualifications and categorical variables of income and subjective financial strain. Variables from the birth registry file were also included as auxiliary variables: the mother’s marital status, the age of the mother and father, and the child’s birthweight and length. Approximately normally-distributed continuous variables including BMI were imputed using truncated regression, specifying as upper and lower limits the smallest and largest values observed in the full MoBa sample. Ordered categorical variables were imputed with ordered logistic regression. There was no missingness in genetic information within the analytic sample. Polygenic scores for adult BMI, childhood body size, depression, ADHD, and educational attainment were included on the right-hand side of the imputation equations, along with indicators for genotyping centre and chip and the 20 principal components of ancestry for all individuals. Continuous variables which were not normally-distributed were imputed with predictive mean matching, specifying knn(10). This included child’s depressive and anxiety symptoms at age 8 (SMFQ and SCARED summary scores), mother’s and father’s depressive/anxiety symptoms at 17 weeks gestation (summary scores based on items from the 25-item Hopkins Checklist (Hesbacher et al., 1980), and mother’s and father’s ADHD symptoms from the 6-item adult ADHD self-report scale (Kessler et al., 2005). To facilitate analysis of ADHD inattention and hyperactivity subscales, the two subscales were imputed, again with predictive mean matching, and the full scale calculated post-imputation with mi passive. An earlier measure of the child’s ADHD symptoms at 5 yrs, based on questions from the Short-Form Conners Parent Rating Scale (Kumar and Steer, 2003), was included as an auxiliary variable. The percentage of imputed data in the analytic sample for each variable was: mother’s education 3.4%, father’s education 2.2%, mother’s smoking 2.0%, father’s smoking 0.7%, mother’s depressive symptoms 3.2%, father’s depressive symptoms 7.1%, mother’s ADHD symptoms 40.9%, father’s ADHD symptoms 60.1%, mother’s BMI: 4.0%, father’s BMI: 3.8%, child’s BMI at age 8: 60.6%, child’s depressive symptoms at age 8: 54.2%, child’s anxiety symptoms at age 8: 54.1%, child’s ADHD symptoms (inattention and hyperactivity) both 54.1%.
5. Polygenic score construction
Beginning with the full GWAS results for each phenotype (Yengo et al., 2018; Richardson et al., 2020; Wray et al., 2018; Demontis et al., 2019), R was used to subset SNPs included in full GWAS results to SNPs also available in the quality controlled MoBa data. From SNPs available in both, independent genome-wide significant associations were identified by clumping in MRBase, specifying r=0.01 and p<5.0 × 10–8 (Hemani et al., 2018b). First subsetting to SNPs available in MoBa and then clumping within these avoids the need for an additional step identifying proxy SNPs. This left 954 SNPs associated with adult BMI, and 321 associated with childhood body size.
Dosage data for MoBa participants for the genetic variants relevant to each polygenic score were extracted from the final quality controlled -bfiles using PLINK, and.raw files imported to STATA. SNPs were harmonized by comparing effect alleles in the GWAS and reference alleles in MoBa.
Polygenic scores were calculated as the sum of the number of effect alleles (0, 1, or 2) multiplied in each case by the harmonized SNP-coefficient from the GWAS. There was a very small amount of missingness in individual SNPs due to quality control measures. Among SNPs included in the adult BMI PGS, the median proportion of individuals missing the SNP was 0.001 and the maximum 0.04. For SNPs included in the childhood body size PGS, equivalent values were 0.002 and 0.04, for SNPs included in the depression and the ADHD PGS these were 0.003 and 0.01. Where individual SNP information was missing for an individual, the sample mean for the number of effect alleles (between 0 and 2) was added for that SNP.
6. Model equations
In within-family MR models, we used polygenic scores for all members of a child-mother-child trio to instrument the BMI of all three individuals. Within-families MR models were adjusted for child’s sex, the first 20 principal components of ancestry for the child and both parents, and the genotyping chip and centre of the child and both parents.
In Stata, this was specified in the form: ivregress 2sls outcome (c_bmi m_bmi f_bmi = c_pgs m_pgs f_pgs) sex c_PC* m_PC* f_PC* c_genotyping_centre* m_ genotyping_centre* f_ genotyping_centre* c_genotyping_chip* m_ genotyping_chip* f_ genotyping_chip*
This can be represented with the following equations:
The outcomes (depressive, anxiety, and ADHD symptoms) can be expressed as:
yi = β0 + β1 Xic + β2Xim + β3Xif + βcC +ei
The three exposures are offspring’s, mother’s, and father’s BMI:
Xic = γoc + γ1cZic + γ2cZim + γ3cZif + γ4cC+ vic
Xim=γom + γ1mZic + γ2mZim + γ3mZif + γ4mC + vim
Xif=γoc + γ1fZic + γ2fZim + γ3fZif + γ4fC + vif where: yi = outcome
Xic = child’s BMI
Xim = mother’s BMI
Xif = father’s BMI
Zic = child’s polygenic score
Zim = mother’s polygenic score
Zif = father’s polygenic score
C=covariates including principal components for offspring, mother and father and genotyping centre and chip ei = error term for the outcome equation vic = error term for the child exposure (BMI) equation vim = error term for the mother exposure (BMI) equation vif = error term for the father exposure (BMI) equation
Additional Results
1. Comparison of analytic sample and excluded participants
To assess if participants included in the analytic sample (N=40,949) differed from others in the birth registry file (N=72,742) (Appendix 1—figure 1), we conducted t-tests and chi-squared tests for key characteristics at birth, BMI, and outcomes using unimputed data. Reflecting the large size of the sample population, several differences reached statistical significance. Children in the analytic sample did not differ from those excluded on sex but were slightly older (mean year of birth: 2005.4 vs 2004.5). They had slightly higher birthweight (mean = 3.6 kg vs 3.4 kg) and slightly younger fathers (mean paternal 32.6 vs 32.8 years). Mothers and fathers included in the analytic sample had slightly higher educational qualifications (e.g., 24.1% and 24.3% of mothers and fathers respectively had 4 years of college education, against 21.0% and 21.5% for those not included). At the time of the child’s birth, mothers in the analytic sample had fewer existing children (e.g., 46.8% vs 42.4% had none), were more likely to be married or cohabiting (97.4% vs 94.4%), and less likely to have smoked in pregnancy (e.g., daily smoking: 4.5% vs 6.3%). Fathers in the analytic sample were also more likely to have stopped smoking during the pregnancy (20.9% vs 15.8%). There was little of a difference in BMI for mothers, fathers, or children. Children in the analytic sample had slightly lower depressive symptoms (mean SMFQ = 1.81 vs 1.91),anxiety symptoms (mean SCARED: 1.04 vs 1.00),and ADHD symptoms (mean RS-DBD ADHD: 8.4 vs 8.7). Descriptive characteristics of the full MoBa sample are in Appendix 1—table 1.
Data availability
The consent given by the participants does not open for storage of data on an individual level in repositories or journals. Researchers who want access to data sets for replication should submit an application to datatilgang@fhi.no. Access to data sets requires approval from The Regional Committee for Medical and Health Research Ethics in Norway and an agreement with MoBa.
References
-
Fetal origins of mental disorders? an answer based on Mendelian randomizationTwin Research and Human Genetics 21:485–494.https://doi.org/10.1017/thg.2018.65
-
Psychometric properties of the screen for child anxiety related emotional disorders (scared): a replication studyJournal of the American Academy of Child and Adolescent Psychiatry 38:1230–1236.https://doi.org/10.1097/00004583-199910000-00011
-
Does depression cause obesity?Journal of Health Psychology 13:1190–1197.https://doi.org/10.1177/1359105308095977
-
The weight wage penalty: a mechanism approach to discriminationEuropean Sociological Review 34:254–267.https://doi.org/10.1093/esr/jcy009
-
Interpreting findings from Mendelian randomization using the MR-egger methodEuropean Journal of Epidemiology 32:377–389.https://doi.org/10.1007/s10654-017-0255-x
-
Attention-Deficit/Hyperactivity disorder (ADHD) and obesity: update 2016Current Psychiatry Reports 19:4.https://doi.org/10.1007/s11920-017-0754-1
-
Within family Mendelian randomization studiesHuman Molecular Genetics 28:R170–R179.https://doi.org/10.1093/hmg/ddz204
-
Maternal obesity and neurodevelopmental and psychiatric disorders in offspringPrenatal Diagnosis 37:95–110.https://doi.org/10.1002/pd.4932
-
Overweight in family members of probands with ADHDEuropean Child & Adolescent Psychiatry 28:1659–1669.https://doi.org/10.1007/s00787-019-01331-7
-
Weight bias in work settings-a qualitative reviewObesity Facts 3:33–40.https://doi.org/10.1159/000276992
-
Is obesity associated with emotional and behavioural problems in children? findings from the millennium cohort studyInternational Journal of Pediatric Obesity 6:e423–e432.https://doi.org/10.3109/17477166.2010.526221
-
Evaluating the potential role of pleiotropy in mendelian randomization studiesHuman Molecular Genetics 27:R195–R208.https://doi.org/10.1093/hmg/ddy163
-
Psychiatric illness in family practiceThe Journal of Clinical Psychiatry 41:6–10.
-
Effects of body mass index on relationship status, social contact and socio-economic position: Mendelian randomization and within-sibling study in UK BiobankInternational Journal of Epidemiology 49:1173–1184.https://doi.org/10.1093/ije/dyz240
-
Comparing the prevalence of mental health problems in children 6-11 across europeSocial Psychiatry and Psychiatric Epidemiology 51:1093–1103.https://doi.org/10.1007/s00127-016-1253-0
-
Factorial validity of the conners’ parent rating scale-revised: short form with psychiatric outpatientsJournal of Personality Assessment 80:252–259.https://doi.org/10.1207/S15327752JPA8003_04
-
Identifying critical points of trajectories of depressive symptoms from childhood to young adulthoodJournal of Youth and Adolescence 48:815–827.https://doi.org/10.1007/s10964-018-0976-5
-
Are there causal relationships between attention-deficit/hyperactivity disorder and body mass index? evidence from multiple genetically informed designsInternational Journal of Epidemiology 50:496–509.https://doi.org/10.1093/ije/dyaa214
-
Cohort profile update: the Norwegian mother and child cohort study (MobA)International Journal of Epidemiology 45:382–388.https://doi.org/10.1093/ije/dyw029
-
Assessing causality in the association between attention-deficit/hyperactivity disorder and obesity: A mendelian randomization studyInternational Journal of Obesity 43:2500–2508.https://doi.org/10.1038/s41366-019-0346-8
-
Body mass index and common mental disorders: exploring the shape of the association and its moderation by age, gender and educationInternational Journal of Obesity 36:414–421.https://doi.org/10.1038/ijo.2011.65
-
The association between socioeconomic disadvantage and attention deficit/hyperactivity disorder (ADHD): a systematic reviewChild Psychiatry and Human Development 47:440–458.https://doi.org/10.1007/s10578-015-0578-3
-
ADHD in children and young people: prevalence, care pathways, and service provisionThe Lancet. Psychiatry 5:175–186.https://doi.org/10.1016/S2215-0366(17)30167-0
-
A rating scale for disruptive behavior disorders, based on the DSM-IV item poolThe Psychiatric Quarterly 76:327–339.https://doi.org/10.1007/s11126-005-4966-x
-
Genetic and environmental effects on body mass index from infancy to the onset of adulthood: an individual-based pooled analysis of 45 twin cohorts participating in the Collaborative project of development of anthropometrical measures in twins (codatwins) studyThe American Journal of Clinical Nutrition 104:371–379.https://doi.org/10.3945/ajcn.116.130252
-
Software application profile: mrrobust—A tool for performing two-sample summary Mendelian randomization analysesInternational Journal of Epidemiology 48:684–690.https://doi.org/10.1093/ije/dyy195
-
Social characteristics as predictors of ADHD labeling across the life courseSociety and Mental Health 11:91–112.https://doi.org/10.1177/2156869320916535
-
Using genetics to understand the causal influence of higher BMI on depressionInternational Journal of Epidemiology 48:834–848.https://doi.org/10.1093/ije/dyy223
-
Circadian rhythm disruption as a link between attention-deficit/hyperactivity disorder and obesity?Journal of Psychosomatic Research 79:443–450.https://doi.org/10.1016/j.jpsychores.2015.10.002
Article and author information
Author details
Funding
Medical Research Council (MC_UU_00011/1)
- Amanda M Hughes
- Eleanor Sanderson
- Tim Morris
- George Davey Smith
- Laura D Howe
- Neil M Davies
University of Bristol
- Amanda M Hughes
- Eleanor Sanderson
- Tim Morris
- George Davey Smith
- Laura D Howe
- Neil M Davies
European Research Council (894675)
- Ziada Ayorech
South-Eastern Norway Regional Health Authority (2019097)
- Ziada Ayorech
Research Council of Norway (295989)
- Neil M Davies
Research Council of Norway (288083)
- Alexandra Havdahl
South-Eastern Norway Regional Health Authority (2018059)
- Alexandra Havdahl
Research Council of Norway (274611)
- Ted Reichborn-Kjennerud
Research Council of Norway (223273)
- Ole A Andreassen
European Research Council (847776 CoMorMent)
- Ole A Andreassen
University of Bergen
- Øyvind Helgeland
South-Eastern Norway Regional Health Authority (2020022)
- Alexandra Havdahl
Helse Vest (912250)
- Stefan Johansson
Helse Vest (F-12144)
- Stefan Johansson
Novo Nordisk (NNF19OC0057445)
- Stefan Johansson
Research Council of Norway (315599)
- Stefan Johansson
Medical Research Council (MR/M020894/1)
- Laura D Howe
European Research Council (293574)
- Pål Njølstad
Trond Mohn stiftelse
- Pål Njølstad
Research Council of Norway (240413)
- Pål Njølstad
Western Norway Regional Health Authority
- Pål Njølstad
Novo Nordisk (54741)
- Pål Njølstad
Norwegian Diabetes Association
- Pål Njølstad
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. We are grateful to all the participating families in Norway who take part in this on-going cohort study. We thank the Norwegian Institute of Public Health (NIPH) for generating high-quality genomic data. This work was performed on the TSD (Tjeneste for Sensitive Data) facilities, owned by the University of Oslo, operated and developed by the TSD service group at the University of Oslo, IT-Department (USIT). (tsd-drift@usit.uio.no). The analyses were performed on resources provided by Sigma2 - the National Infrastructure for High Performance Computing and Data Storage in Norway. This research is part of the HARVEST collaboration, supported by the Research Council of Norway (#229624). We also thank the NORMENT Centre for providing genotype data, funded by the Research Council of Norway (#223273) and deCODE Genetics, South East Norway Health Authority and KG Jebsen Stiftelsen. We further thank the Center for Diabetes Research, the University of Bergen for providing genotype data and performing quality control and imputation of the data funded by the ERC AdG project SELECTionPREDISPOSED, Stiftelsen Kristian Gerhard Jebsen, Trond Mohn Foundation, the Research Council of Norway, the Novo Nordisk Foundation, the University of Bergen, and the Western Norway Health Authorities (Helse Vest). Funding: This research was funded by a project entitled ‘social and economic consequences of health: causal inference methods and longitudinal, intergenerational data’, which is part of the Health Foundation’s Social and Economic Value of Health Programme (Grant ID: 807293). The Health Foundation is an independent charity committed to bringing about better health and health care for people in the UK. This research is part of the HARVEST collaboration, supported by the Research Council of Norway (#229624). Individual co-authors area also supported by specific sources of funding. ZA is supported by a Marie Skłodowska-Curie Fellowship from the European Union (894675) and the South-Eastern Norway Regional Health Authority (2019097). TR is supported by the Research Council of Norway (274611 PI: Reichborn-Kjennerud). OAA is funded by the Research Council of Norway (223273) and EU H2020 RIA (847776 CoMorMent). ØH is supported by the University of Bergen, Norway. SJ was supported by Helse Vest’s Open Research Grant (grants #912250 and F-12144), the Novo Nordisk Foundation (grant NNF19OC0057445) and the Research Council of Norway (grant #315599). PN is supported by the European Research Council (AdG SELECTionPREDISPOSED #293574), the Trond Mohn Foundation (Mohn Center for Diabetes Precision Medicine), the Research Council of Norway (FRIPRO grant #240413), the Western Norway Regional Health Authority (Strategic Fund “Personalized Medicine for Children and Adults”), the Novo Nordisk Foundation (grant #54741), and the Norwegian Diabetes Association. AH was supported by the South-Eastern Norway Regional Health Authority (2018059, 2020022) and the Research Council of Norway (288083). LDH is supported by a Career Development Award from the UK Medical Research Council (MR/M020894/1). NMD was supported via a Research Council of Norway grant (295989). The Medical Research Council (MRC) and the University of Bristol support the MRC Integrative Epidemiology Unit [MC_UU_00011/1] (AMH, ES, LDH, NMD, GDS, TM). The funders had no role in the design or execution of this analysis, interpretation of results, or the decision to publish.
Ethics
The establishment of MoBa and initial data collection was based on a license from the Norwegian Data Protection Agency and The Regional Committees (REC) for Medical and Health Research Ethics. The REC South East Norway, one of four in Norway, was the ethical committee that evaluated the ethics of this study. Approval from the REC was granted (2016/1702). Informed consent was obtained from each MoBa participant upon recruitment, which included consent to link to the Medical Birth Registry of Norway (MBRN). The MoBa cohort is now based on regulations related to the Norwegian Health Registry Act.
Copyright
© 2022, Hughes et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,240
- views
-
- 167
- downloads
-
- 10
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.