Introduction
Barker’s hypothesis suggested that early in development, especially during life in the intrauterine life, adverse effects can lead to permanent changes in physiology and metabolism, resulting in an increased risk of disease in adulthood [Reference Bulmer, J, M, K and A1]. Birth weight is considered to be an important predictor of neonatal and infant survival and is associated with the risk of noncommunicable, chronic illnesses in the offspring [Reference Martin-Gronert and Ozanne2], especially metabolic diseases [Reference Knop, Geng, Gorny, Ding, Li and Ley3, Reference Warrington, Beaumont, Horikoshi, Day, Helgeland and Laurin4] and cardiovascular diseases [Reference Knop, Geng, Gorny, Ding, Li and Ley3,Reference Warrington, Beaumont, Horikoshi, Day, Helgeland and Laurin4]. The development of brain can also be influenced by birth weight, which might lead to the unhealthily physical and mental development, even psychological illness [5–8]. For instance, some literature reported that birth weight was associated with adult cortical surface area hippocampal volume, which was mainly responsible for long-period memory storage conversion and orientation [Reference Raznahan, Greenstein, Lee, Clasen and Giedd5,Reference Walhovd, Fjell, Brown, Kuperman, Chung and Hagler6]. Wald et al. indicated the effects of birth weight variation within the normal range on children’s neuropsychological function, including language, executive functioning, and theory of mind [Reference Wade, Browne, Madigan, Plamondon and Jenkins7].
Previous studies also observed relationship between birth weight and depression. Warrington et al. conducted a study about maternal and fetal genetic effects on birth weight, and demonstrated genetic correlation between birth weight and mental disorder through linkage disequilibrium (LD) score regression, such as major depressive disorder and bipolar disorder (BD) [Reference Warrington, Beaumont, Horikoshi, Day, Helgeland and Laurin4]. However, previous literature showed that the association between depression and birth weight remains controversial [Reference Lyall, Inskip, Mackay, Deary, McIntosh and Hotopf9–12]. For example, an epidemiological study did not find compelling association between low birth weight and suspected major depression from early adolescence to early adulthood [Reference Levine11]. Another systematic review and meta-analysis displayed a weak association between low birth weight and later depression or psychological distress [Reference Wojcik, Lee, Colman, Hardy and Hotopf12].
Genetic factors contribute greatly to the variations of birth weight, brain volume and depression. Resent genome-wide association studies (GWAS) meta-analyses of own birth weight identified 146 independent single-nucleotide polymorphisms (SNPs) at genome-wide significance [Reference Warrington, Beaumont, Horikoshi, Day, Helgeland and Laurin4]. Zhao et al. performed a GWAS on 101 brain volume phenotypes using UK Biobank including 19,629 participants, identifying 365 independent genetic variations [Reference Zhao, Luo, Li, Li, Zhang and Shan13]. Twin and family estimates suggest that approximately 30–40% of the unipolar depression can be explained by genetic effects [Reference Sullivan, Neale and Kendler14]. Clarifying the genetic relationships between birth weight, depression and brain volumes can provide important clues for understanding the potential impact of birth weight on brain volumes and depression.
Polygenic risk score (PRS) analysis is a powerful approach that aggregates the effects of genetic variants from GWAS summary data to identify the relationship between traits and diseases. From the Suffolk County Mental Project, an inception cohort study of first-admission patients with psychosis, Jonas et al. indicated that schizophrenia (SZ) PRS was the strongest predictors of diagnostic shifts from affective to nonaffective psychosis and could predict whether patients who appear to have a mood disorder with psychotic at first admission will eventually be diagnosed with a SZ spectrum disorder [Reference Jonas, Lencz, Li, Malhotra, Perlman and Fochtmann15]. Another PRS study provided evidence that BD with manic psychosis genetically was more similar to SZ than any other tested BD subgroup through comparing SZ PRSs of different psychosis subtypes of BD [Reference Markota, Coombes, Larrabee, McElroy, Bond and Veldic16].
Utilizing the UK Biobank cohort, we first tested the associations between birth weight, depression, and brain volumes. The PRSs of birth weight were then calculated from the individual-level genotype data of UK Biobank cohort. Linear and logistic regression models were then performed to evaluate the genetic impact of birth weight on depression and brain volume. Finally, the mediation model was used to estimate the role of brain volume in the association between birth weight and depression. All analyses were conducted sex stratified to assess sex-specific role in the associations.
Aims of the Study
This study attempted to explore: (a) associations between birth weight on brain volumes; (b) associations between brain volumes on depression; and (c) associations between birth weight on depression, and whether brain volumes plays a mediating role in it.
Methods
UK Biobank cohort
Study individuals were from the UK Biobank health resource, which has performed a large prospective population-based cohort study, comprising linked health, hospital-record and genetic data of 502,656 participants aged 40–69 years in 2006 and 2010. All participants agreed to use their anonymous data and samples for any health-related research, to reconnect for further substudies [Reference Cathie, John, Naomi, Valerie, Paul and John17]. The present study was conducted under UK Biobank approvals, to access participants’ health-related records, including birth weight from screenshot question or verbal interview within Assessment Center, self-reported depression status and brain volume from magnetic resonance imaging assessments.
UK Biobank phenotypes of birth weight, depression, and brain
Participants were asked to enter their own birth-weight through verbal interview with 277,000 participants completed (UK Biobank field ID: 20022). In order to obtain a comprehensive and accurate control group of depressed and nondepressed patients, Patient Health Questionnaire (PHQ-9) and another strict criterion based on composite international diagnostic interview short-form (CIDI-SF) [Reference Davis, Cullen, Adams, Brailean, Breen and Coleman18,Reference Manea, Gilbody and McMillan19] were used to define the inclusion and exclusion criteria (full details in Appendix in the Supplementary Material). The following indicators, volume of gray matter (UK Biobank code: 25006), white matter (UK Biobank code: 25008), gray and white matter (UK Biobank code: 25010), left and right thalamus (UK Biobank code: 25011:25012), left and right accumbens (UK Biobank code: 25023:25024), gray matter in left frontal pole and gray matter in right frontal pole (UK Biobank code: 25782:25783), and gray matter in left VIIIa cerebellum and gray matter in right VIIIa cerebellum (UK Biobank code: 25909 and 25911) were selected as the measurement of brain structure. The total intracranial volume was defined as the sum of total gray matter, white matter, and cerebrospinal fluid. Birth weight and brain volume were mean-centered and normalized to one standard deviation (SD) before further analysis.
UK Biobank genotyping, imputation, and quality control
A total of 488,377 participants were genotyped undertaking either the UK Bileve array or the UK Biobank axiom array under favorable quality control. Details of the array design, genotyping, and quality control procedures have been described elsewhere [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp20]. These genotypes were subsequently imputed to the Haplotype Reference Consortium (HRC) reference panel [Reference Mccarthy, Das, Kretzschmar, Delaneau, Wood and Teumer21] (version 1.1) and UK10K and 1000 Genomes project reference panels [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp20]. We applied additional quality control filters to select high-quality SNPs: SNPs with high linkage disequilibrium (r 2 > 0.5) were pruned. We removed the participants, who reported inconsistencies between self-reported gender and genetic gender, who were genotyped but not imputed, and who withdraw their consents. Additionally, individuals were restricted to only “white British” based on self-reported ethnicity (UK Biobank field ID: 21000). UK Biobank participants was used to obtain relatively independent SNPs through King software.
GWAS datasets of birth weight and PRS calculation
The GWAS summary statistics of birth weight [Reference Knop, Geng, Gorny, Ding, Li and Ley3] was from a recent meta-analysis of early growth genetics (EGG) consortium, combining data from EGG consortium and the fully released data of UK Biobank (released May 2017). Briefly, data from EGG consortium were imputed against the 1000 Genomes Project (Phase 1 v3) reference panel and data from UK Biobank were imputed against the HRC reference panel. The association between each genetic variant and birth weight measure was tested using linear regression in EGG consortium and a linear mixed model in UK Biobank, with adjustment for gestational age. We selected the SNP with genome-wide association testing p value <5.0 × 10−8 to calculate PRSs for all individuals. The GWAS summary statistics of European-only meta-analysis of own birth weight including 298,142 individuals was used in the present study. More details about the GWAS of birth weight were described in the previous study [Reference Knop, Geng, Gorny, Ding, Li and Ley3].
The PRS of birth weight was calculated as the sum of allele, weighted by their effect sizes, defined by $ {\mathrm{PRS}}_i=\sum \limits_{n-1}^t{\beta}_n{\mathrm{SNP}}_{ni} $, where n (n = 1, 2, 3, …, t) and i (i = 1, 2, 3, …, k) denote the number of genetic markers and the sample size, respectively. βn represents the effect value of risk allele of the nth significant SNP related to birth weight obtained from the previous study [Reference Knop, Geng, Gorny, Ding, Li and Ley3]. SNPni represents the risk allele dosage of the nth SNP of the ith individual. At last, 626 SNPs associated with birth weight driven from EGG consortium were used for PRS calculation by using PLINK’s“--score” command [Reference Purcell, Neale, Todd-Brown, Thomas, Ferreira and Bender22]. We conducted data centralization and standardization of the PRSs of birth weight.
Statistical analysis
We applied a widely used approach for mediation analysis: the causal steps approach, which was described by Kenny et al. [Reference Baron and Kenny23] (The detailed information of the causal steps approach was provided in the Supplementary Material). First, logistic regression analyses were conducted to estimate the associations between birth weight and depression, where birth weight phenotype or related PRS was used as an exposure variable, and depression was used as an outcome variable. Second, linear regression analyses were applied for evaluating the associations between birth weight and brain volumes, where birth weight phenotype or related PRS was used as exposure variable and brain volumes were used as outcome variables. Third, logistic regression analyses were conducted to estimate the associations between brain volumes and depression, where brain volumes were used as exposure variables, and depression was used as an outcome variable. Fourth, logistic regression analyses were also performed for depression and birth weight, where each of brain volumes was used as a covariate to explore the role of brain volume in the association between birth weight and depression. For each regression analysis, sex, age, and 10 principle components of population structure, Townsend deprivation and alcohol use were used as covariates, and total intracranial volume was only used as a covariate for the third and fourth steps. All analyses were conducted sex stratified to assess sex-specific role in the associations through R 3.5.3.
Sensitive analyses
We used the R package “mediation,” which relied on the causal steps approach [Reference Baron and Kenny23] for the mediation analysis, and conducted the similar analyses with the same data. We used the mediate function to estimate the indirect effect (or average causal mediation effect (ACME)) and direct effect (or average direct effect (ADE)) of birth weight to depression. The mediation model needs two statistical models, the one is the linear regression model of the exposure variable (birth weight) and mediate variable (brain volume), the other one is the logistic regression model of the exposure variable (birth weight), mediate variable (brain volume), and outcome variable (depression). The default number of simulation (1000) was used, which type is the quasi-Bayesian Monte Carlo method based on normal approximation [Reference Imai, Keele and Tingley24]. From statistics level, only when the p value of ACME result was significant, suggesting that the brain volumes play the mediating role in the association between birth weight and depression. All analyses were conducted sex stratified to assess sex-specific role.
Result
Descriptive characteristics of study samples
Table 1 shows the descript characteristics of study samples.
Age was described as mean (standard deviation [SD]).
Abbreviation: PRS, polygenic risk scores.
Associations between birth weight and depression
We observed associations between birth weight and depression (odds ratio [OR] = 0.968, 95% confidence interval [CI] = 0.957–0.979, p = 2.29 × 10−6), as well as birth weight PRS and depression (OR = 0.984, 95% CI = 0.975–0.993, p = 0.002). In women, depression was associated with birth weight phenotype (OR = 0.962, 95% CI = 0.948–0.975, p = 4.87 × 10−6) and PRS (OR = 0.986, 95% CI = 0.974–0.997, p = 0.035). In men, depression was only associated with birth weight PRS (OR = 0.980, 95% CI = 0.961–0.998, p = 0.028) (Table 2).
Logistic regression was used to test the association of depression and birth weight (PRS and phenotype). Birth weight PRS indicates the polygenic scores for birth weight. Birth weight phenotype means the phenotype of birth weight.
Abbreviation: PRS, polygenic risk scores.
Associations between birth weight and brain volumes
As shown in Figure 1, we observed the associations between birth weight phenotype and the brain volumes, such as gray matter (B = 0.131, p = 3.51 × 10−74), gray and white matter (B = 0.137, p = 5.70 × 10−85), and thalamus (B = 0.136, p = 3.09 × 10−72 for left and B = 0.141, p = 1.86 × 10−79 for right). More additional details were provided in Table S1 in the Supplementary Material. In men and women, we obtained the similar results that birth weight was significant associated with each of the brain volumes. Detailed information was provided in the Figure 3, Figures S2 and S3 and Tables S2 and S3 in the Supplementary Material.
As shown in Figure 1, further PRS analysis also observed genetic associations between birth weight and brain volumes, such as gray matter (B = 0.027, p = 8.80 × 10−7), gray and white matter (B = 0.028, p = 2.41 × 10−7), and thalamus (B = 0.017, p = 0.004 for left and B = 0.020, p = 5.86 × 10−4 for right). More additional details were provided in Table S1 in the Supplementary Material. In men, we obtained the similar results that birth weight PRS was significant associated with each of the brain volumes, except the volume of gray matter in VIIIa cerebellum. And in women, birth weight PRS was only associated with the brain volume of gray matter, white matter, gray and white matter, left accumbens, and gray matter in frontal pole. Detailed information was provided in the Figure 3 Figures S2 and S3 and Tables S2 and S3 in the Supplementary Material.
Associations between depression and brain volumes
As shown in Figure 2, associations were observed between depression and brain volumes, such as gray matter (OR = 0.879, 95% CI = 0.805–0.960, p = 0.016), gray and white matter (OR = 0.702, 95% CI = 0.552–0.893, p = 0.016), left thalamus (OR = 0.891, 95% CI = 0.850–0.933, p = 4.46 × 10−5), and right thalamus (OR = 0.884, 95% CI = 0.841–0.928, p = 2.67 × 10−5). More additional details were provided in Table S4 in the Supplementary Material. In men, we only obtained that depression was significant associated with gray matter and gray matter in right VIIIa cerebellum. And in women, depression was associated with the brain volume of gray and white matter, thalamus, accumbens and gray matter in VIIIa cerebellum. Detailed information was provided in the Figure 3, Figures S4 and S5 and Tables S5 and S6 in the Supplementary Material.
Associations between birth weight and depression adjusting for brain volume as covariate
We observed no significant association between birth weight (both phenotype and PRS) and depression, when one of the brain volumes was used as a covariate (all p value >0.05). After gender classified, we observed the same results. Detailed information was provided in the Figure 3, Figures S6 and S7 and Tables S7–S9 in the Supplementary Material.
Summary of mediation effects of brain volumes between birth weight and depression
Through the causal steps approach, parts of the brain volumes were the complete mediators between birth weight and depression. For birth weight phenotype, the brain volumes included gray matter, gray and white matter, left thalamus, right thalamus, left accumbens, and gray matter in left VIIIa cerebellum. In men, gray matter and gray matter in left VIIIa cerebellum played the mediation effects between birth weight and depression. And in women, left and right thalamus, left and right accumbens and gray matter in left VIIIa cerebellum played the mediation effects between birth weight and depression.
For birth weight PRS, the brain volumes included the brain volumes included gray matter, gray and white matter, left and right thalamus, left accumbens, and gray matter in left VIIIa cerebellum. In men, only gray matter played the mediation effects between birth weight and depression. And in women, only left accumbens played the mediation effects between birth weight and depression. More additional details were provided in Figure 3 and Figures S6 and S7 in the Supplementary Material.
Sensitive analyses
We obtained the similar results from mediation analyses through R package. The indirect path from birth weight PRS that were based on a GWAS to depression via elevated levels of one of the brain volumes, including gray matter, gray and white matter, left thalamus, right thalamus, and left accumbens. More additional details were provided in Tables S10–S12 in the Supplementary Material.
Discussion
We conducted a large observational and genetic PRS study to evaluate the associations between birth weight, brain volume, and depression using UK Biobank population cohort. We observed phenotypic associations for birth weight versus brain volumes, birth weight versus depression, and depression versus brain volumes. Further PRS analysis also provided evidence for genetic predispositions between birth weight, depression and brain volumes. It is interesting that the correlation between birth weight and depression was not significant after adjusting the impact of brain volume, suggesting that brain volumes might be an important mediator between birth weight and depression.
Previous studies have demonstrated that birth weight was associated with brain development [Reference Gujral, Aizenstein, Reynolds, Butters and Erickson25–27], which is consistent with our results. Our results indicated the positive association between birth weight and multiple brain volumes. A meta-analysis of birth weight concluded that birth weight was associated with the volume of white matter, gray matter, gray and white matter, and thalamus [Reference Kesler, Ment, Vohr, Pajot, Schneider and Katz28]. In this study, we also detected several novel brain regions to test the association with birth weight, such as accumbens volume and frontal pole volume. Additionally, limited literature provided evidence for the genetic association between birth weight and brain volume. For instance, rs8756 and rs7968682 were reported to be associated with brain structure [Reference Warrington, Beaumont, Horikoshi, Day, Helgeland and Laurin4]. Through a GWAS of regional brain volume in 19,629 individuals, authors indicated a significant correlation between depressive symptoms and gray matter [Reference Zhao, Luo, Li, Li, Zhang and Shan13]. Our PRS analysis results provided novel evidence for the effects of birth weight related genetic factors on the variations of brain volumes. For the reason why low birth weight is risk factor for brain development, combining with previous studies [Reference Volpe29,Reference Bjuland, Rimol, Løhaugen and Skranes30], one possible explanation is that low birth weight may affect normal brain development and maturation after birth, eventually leading to abnormal brain volumes in adults.
Meanwhile, association was also observed between brain volume and depression in this study. It has been suggested that decreased volume of the hippocampus, anterior cingulate cortex, prefrontal cortex, striatum and amygdala, and impaired white matter integrity are often observed in depressed adults [Reference Gujral, Aizenstein, Reynolds, Butters and Erickson25]. Merz et al. conducted a significant correlation between higher levels of depressive symptoms and lower levels of hippocampal gray matter and larger volumes of pale gray matter in 328 children and adolescents between the ages of 7 and 21 [Reference Merz, He and Noble31]. In this study, we found that depression was associated with multiple brain volumes, such as gray matter, brain of gray and white matter, left and right thalamus, left and right accumbens, gray matter in left and right frontal pole, and VIIIa cerebellum. Interestingly, most of the brain volumes that were associated with birth weight were also associated with depression in this study.
As described in the introduction section, the association between birth weight and depression remains controversial [9–12]. Some of previous studies proved that low birth weight increased the risk of depressive disorder [Reference Lyall, Inskip, Mackay, Deary, McIntosh and Hotopf9,Reference Van Lieshout, Boyle and Favotto10], which was consistent with the present study. However, Levine et al. demonstrated that there is a lack of convincing evidence to support an association between very low birth weight and suspected major depression, taking into account well-documented confounding factors [Reference Levine11]. Our results shown that not only the increase of the birth weight phenotype, but also the PRS for birth weight can lead to the decreased prevalence of depression, which also meant the genetic liability between birth weight and depression. Combined with the insignificant results of other studies [Reference Levine11,Reference Wojcik, Lee, Colman, Hardy and Hotopf12], we may speculate that the effect of birth weight on depression was not direct but indirect.
Based on the results, higher birth weight was associated with higher brain volumes, higher brain volumes were associated with lower morbidity of depression, and the weak association between depression and birth weight, we inferred that the brain volumes might play an important role for the observed association between birth weight and depression. This hypothesis was supported by the fact that the association between birth weight and depression disappeared, after adjusting the one of brain volumes as a covariate in regression analysis and significant indirect path results from the mediation approach. The possible mechanism is that the low brain volume of low birth weight leads to neuropathologic changes in adulthood, resulting in the risk of depression. Also, birth weight may affect the risk of depression through other mediating variables, so that more researches are expected.
Using the individual-level phenotypic and genetic and data from the UK Biobank, the large sample size and a wide range of variables make our conclusions more rigorous and convincing. Using PRS method to explore the effect of birth weight on depression and brain volume is a powerful method to reduce the impact of nongenetic confounding factors on our study results.
There were some limitations to this study. First, depression was measured retrospectively by self-report through an online questionnaire and birth weight was collected from participants’ recall, which may increase the likelihood of measurement errors. The conclusions are mainly applicable to middle-aged and elderly European populations, lacking comparisons among the young and the interracial populations. Second, the study is a cross-sectional study. Therefore, a causal relationship between depression and brain volume cannot be established. Third, it could be at birth persons with low birth weight had different brain development at birth, through the brain dysplasia of low birth weight [Reference Volpe29,Reference Bjuland, Rimol, Løhaugen and Skranes30]. Although limited studies suggested that age had no effect on the difference in brain volume between low birth weight and normal weight [Reference Bjuland, Rimol, Løhaugen and Skranes32,Reference Kieviet, Zoetebier, van Elburg, Vermeulen and Oosterlaan33], further longitudinal studies are expected to prove the continuity of mediation. At last, there are some other confounders might influence depression. Although we chose Townsend deprivation index and alcohol use as covariates according to the previous studies [34–36], some another possibly confounders were not considered in this study and need further research.
Due to treatment for depression can only alleviate rather than cure the condition, and the burden of depression outcome continues to grow. Therefore, the identification and prevention of the high-risk group of depression is an effective and key measure. Our study showed that birth weight was positively correlated with depression, negatively correlated with brain volumes and brain volume might an important mediator between birth weight and depression. Therefore, this study suggested that paying more attention to the birth weight and brain development might be helpful for early prevention and treatment of depression.
Acknowledgments
This research has been conducted using the UK Biobank resource, which has been approved by the National Research Ethics Service Committee.
Conflict of Interest.
There are no conflicts of interest to declare.
Data Availability Statement
The UKB data are available through the UK Biobank Access Management System https://www.ukbiobank.ac.uk/. We will return the derived data fields following UKB policy; in due course, they will be available through the UK Biobank Access Management system.
Code Availability.
a<−lm (brain_volume~birth_weight+covariate).
b<−glm(depression~brain_volume+covariate,family=binomial).
c<−mediate(a,b,treat=‘birth_weight’, mediator=‘brain_volume’).
summary(c).
Supplementary Materials
To view supplementary material for this article, please visit http://dx.doi.org/10.1192/j.eurpsy.2020.74.
Comments
No Comments have been published for this article.