- Research article
- Open access
- Published:
Both indirect maternal and direct fetal genetic effects reflect the observational relationship between higher birth weight and lower adult bone mass
BMC Medicine volume 20, Article number: 361 (2022)
Abstract
Background
Birth weight is considered not only to undermine future growth, but also to induce lifelong diseases; the aim of this study is to explore the relationship between birth weight and adult bone mass.
Methods
We performed multivariable regression analyses to assess the association of birth weight with bone parameters measured by dual-energy X-ray absorptiometry (DXA) and by quantitative ultrasound (QUS), independently. We also implemented a systemic Mendelian randomization (MR) analysis to explore the causal association between them with both fetal-specific and maternal-specific instrumental variables.
Results
In the observational analyses, we found that higher birth weight could increase the adult bone area (lumbar spine, β-coefficient= 0.17, P < 2.00 × 10−16; lateral spine, β-coefficient = 0.02, P = 0.04), decrease bone mineral content-adjusted bone area (BMCadjArea) (lumbar spine, β-coefficient= − 0.01, P = 2.27 × 10−14; lateral spine, β-coefficient = − 0.05, P = 0.001), and decrease adult bone mineral density (BMD) (lumbar spine, β-coefficient = − 0.04, P = 0.007; lateral spine; β-coefficient = − 0.03, P = 0.02; heel, β-coefficient = − 0.06, P < 2.00 × 10−16), and we observed that the effect of birth weight on bone size was larger than that on BMC. In MR analyses, the higher fetal-specific genetically determined birth weight was identified to be associated with higher bone area (lumbar spine; β-coefficient = 0.15, P = 1.26 × 10−6, total hip, β-coefficient = 0.15, P = 0.005; intertrochanteric area, β-coefficient = 0.13, P = 0.0009; trochanter area, β-coefficient = 0.11, P = 0.03) but lower BMD (lumbar spine, β-coefficient = − 0.10, P = 0.01; lateral spine, β-coefficient = − 0.12, P = 0.0003, and heel β-coefficient = − 0.11, P = 3.33 × 10−13). In addition, we found that the higher maternal-specific genetically determined offspring birth weight was associated with lower offspring adult heel BMD (β-coefficient = − 0.001, P = 0.04).
Conclusions
The observational analyses suggested that higher birth weight was associated with the increased adult bone area but decreased BMD. By leveraging the genetic instrumental variables with maternal- and fetal-specific effects on birth weight, the observed relationship could be reflected by both the direct fetal and indirect maternal genetic effects.
Background
Osteoporosis is a common skeletal disease characterized by the loss of bone mass and the risk of brittle fracture. It has become a major global public health problem, affecting about 200 million people around the world [1]. Osteoporosis could be clinically diagnosed by measuring the bone mineral density (BMD), which is a highly heritable trait [2]. Some modifiable environmental factors such as body weight [3], physical activity [4], sleep behavior [5], and inflammatory disease [6] could have an impact on bone mass gain and the development of osteoporosis. Other non-modifiable factors included sex, age, and genetic factors. The genome-wide association studies (GWASs) have discovered hundreds of genetic loci that are associated with BMD, osteoporosis, and fracture in the past decade [7]. However, there is yet a considerable proportion of the variance in bone mass that cannot be explained by known genetic and environmental determinants, e.g., the role of the growth during the prenatal period on adult osteoporosis was often underestimated.
The Developmental Origins of Health and Disease (DOHaD) hypothesis was initially proposed in the 1990s [8], which proposed that the growth during the prenatal period might play a critical role in an individual’s short- and long-term health [9]. Although the relationship between birth weight and risk of osteoporosis in later life has been investigated in observational studies, the results were controversial. Some studies detected a positive association between birth weight and bone mass [10,11,12,13], but some other studies found an inverse association [14,15,16]; these conflicting results might be in part due to the limited sample size and other unmeasured confounding factors. With the development of genomic medicine, the Mendelian randomization (MR) approach somehow provides a possibility to reduce the confounding bias or reverse causation by leveraging genetic variants to instrument a potential exposure and therefore leads to more robust results than the conventional observational approach [17]. Recently, by conducting genome-wide association analyses of own birth weight versus own genetic factors and offspring birth weight versus mother genetic factors, Warrington et al. have successfully partitioned the genetic effects on birth weight into direct fetal genetic component and indirect maternal genetic component [18]. Consequently, by using the maternal and fetal genetic effects separately, Moen et al. suggested that the maternal intrauterine environment, as proxied by maternal genetic variants that influence offspring birthweight, was unlikely to be a major determinant of adverse cardiometabolic outcomes, but offspring birth weight determined by own genetic factors was found to be associated with offspring cardiometabolic risk factors [19].
Therefore, the aim of this study is to investigate the effect of birth weight on adult BMD by integrating the observational data and both fetal and maternal genetic evidence. We first performed multivariable regression analyses to assess the association of birth weight with BMD-related parameters measured by dual-energy X-ray absorptiometry (DXA) and by quantitative ultrasound (QUS) in the UK Biobank. We also applied the MR approach to test whether the observational relationship between birth weight and adult BMD could be reflected by either indirect maternal or direct fetal genetic effects or both.
Methods
Data sources and phenotypes
An overview of the study design is illustrated in Fig. 1. In the observational study and weighted genetic risk score analysis, we used individual-level data from the UK Biobank (Application 41376) as we used before [6, 20]. In the present study, the exposure feature of interest was birth weight (Field ID 20022), and we obtained 280,255 participants with birth weight measurements in the UK Biobank dataset. Next, we applied the exclusion criterion step by step: (1) 2 time points of birth weight difference > 1 kg (N = 84), (2) birth weight < 2.5 kg, (3) birth weight > 4.5 kg (N = 42,613), and (4) non-Europeans (N = 29,767). We calculated the mean value of the birth weight if there were more than one visit. The outcomes were the spine DXA BMD and ultrasound-derived heel BMD. The spine DXA data in the UK Biobank were provided as lumbar spine (LS) and lateral spine (LaS), the scanning of the lumbar spine (from L1 to L4) requested the lower legs of the participant be placed on a polystyrene block, bringing the hips and knees to 80° flexion. The scanning of the lateral spine (from L4 to T4) requested the participant to lie on their side and acquire the spine in the lateral plane. Here, we extracted the corresponding bone parameters including the LS BMD (Field ID 23204), LS bone mineral content (BMC) (Field ID 23203), LS bone area (BA) (Field ID 23200), LaS BMD (Field ID 23234), LaS BMC (Field ID 23312), and LaS BA (Field ID 23311). Finally, we got 20,367 participants with both LS data and birth weight measurements, and 21,364 participants with both LaS data and birth weight measurements. Furthermore, we calculated the area-adjusted bone mineral content (BMCadjArea) as an additional outcome; to generate the measurement of BMCadjArea, we first performed a linear regression for BMC with the independent variable bone area, and the residual of BMC for each individual was predicted, then, the BMCadjArea of each individual was calculated as the sum of the mean of BMC and the residual. The heel BMD (Field ID 3084, 3148, and 4105) was estimated from the quantitative ultrasound (QUS) measurement with the following formula: heel estimated BMD = 0.0025926 × (bone ultrasound attenuation + speed of sound) − 3.687. We excluded the unusually large or small heel BMD value (the heel BMD value > mean + 3SD or mean − 3SD), and we also excluded the individuals which were measured by DXA scan, leaving 177,675 individuals with both QUS heel BMD and birth weight measurements.
We extracted the potential confounding factors from the questionnaires such as sex (Field ID 31), age (Field ID 21003), weight (Field ID 21002), smoking status (Field ID 20116), alcohol intake frequency (Field ID 1558), physical activity (Field ID 884, 894, 904, 914), and assessment center (Field ID 54). The physical activity was assessed whether the individual accorded with at least one of the 4 criteria: (1) at least moderate physical activity 5 days a week, (2) vigorous activity once a week, (3) more than 150 min of moderate activity per week, and (4) more than 75 min of vigorous activity per week.
The selection of instrumental variables
The fetal single nucleotide polymorphisms (SNPs) associated with own birth weight (sample size: 321,223) and the maternal SNPs associated with offspring birth weight (sample size: 230,069 mothers) were extracted from the large-scale GWAS meta-analyses from the Early Growth Genetics (EGG) Consortium and the UK Biobank dataset [18]. Specifically, 205 autosomal SNPs were independently (r2 < 0.1) associated with birth weight at genome-wide significance P < 6.6 × 10−9. Of which, 63 autosomal SNPs that have a fetal-only effect on own birth weight represent the effect of own genotype on own birth weight after adjusting the maternal effect on birth weight using the weighted linear regression (WLM method; an approximation of the structural equation model), and 31 autosomal SNPs that had a maternal-only effect on offspring birth weight represent the effect of maternal genotype on offspring birth weight after adjusting the fetal effect on birth weight; 41 autosomal SNPs have both fetal and maternal effects (Additional file 1: Table S1). Here, we selected 63 autosomal SNPs as the fetal instrument variables to investigate the causal association between the genetically determined birth weight and BMD. As a complement, MR analyses with 104 (= 63 + 41) and 205 SNPs were also performed. Similarly, we selected 31 maternal-only effect SNPs as the maternal instrumental variables to represent the intrauterine environment in the parent-offspring pair MR analysis. We also performed MR analyses with 72 (= 31 + 41) and 205 SNPs as additional evidence. Furthermore, we evaluated the strength of the abovementioned instrument variables using two parameters: the proportion of variance explained (R2) [21], which was calculated from the formula ∑β2 × 2 × MAF × (1 − MAF), where the β was the estimated effect, MAF was the minor allele frequency, and the F statistic, which was calculated from the formula F = [(N – K − 1)/K] × R2/(1 − R2), where the N was the sample size, K was the number of included SNPs, and the R2 was the proportion of variance explained by the genetic variants [21]. Typically, the SNPs with an F-statistic > 10 were considered as the strong and validated instrumental variables [22]. Additionally, the GWAS summary-level data for femoral neck BMD (FN BMD, n = 49,988) [23] was extracted from the GEnetic Factors for Osteoporosis (GEFOS) Consortium (http://www.gefos.org/). We also obtained the GWAS summary-level data for bone areas including the femoral neck area, intertrochanteric area, trochanter area, and total hip area [24] from the deCODE genetics (https://www.decode.com/summarydata). The abovementioned datasets were shown in Additional file 1: Table S2.
Observational analyses
Firstly, we attempted to test the potential non-linear relationship between the birth weight and each of the bone parameters (BMC, bone area, BMCadjArea, and BMD) using the multivariable restricted cubic spline (RCS) model with four knots (using the “ols” function from the “rms” library in R) adjusting for age, gender, weight, alcohol status, smoking status, and physical status. Then, we performed a multivariable linear regression to explore the association of birth weight with the spine DXA scan results including the lumbar spine and the lateral spine sites in 20,367 and 21,364 participants, respectively; here, the outcomes were LS/LaS BMC, LS/LaS BA, LS/LaS BMCadjArea, and LS/LaS BMD. The confounding factors were sex, age, weight, smoking, alcohol, and physical activity. The baseline characteristics of the DXA measured datasets were shown in Additional file 1: Table S3. Besides, we also implemented the multivariable linear regressions with the same confounders to explore the relationship between birth weight and QUS measured heel BMD in 177,675 individuals. Of note, these individuals did not overlap with the DXA-measured participants, which could be served as an independent replication. The baseline characteristics of the QUS-measured heel BMD dataset were shown in Additional file 1: Table S4. To be noted, we applied a standardized approach (z-scores) for the outcomes to give a mean = 0 and standard deviation = 1 and then conducted the above analyses. All observational analyses were performed using the R software (https://www.r-project.org/), and P < 0.05 was set as the significant threshold.
Mendelian randomization
Fetal SNPs and own birth weight
We performed a weighted genetic risk score (wGRS) analysis in the UK Biobank using the individual-level data to test the relationship between genetically determined birth weight and BMC, bone area, BMCadjArea, and BMD. Here, we constructed the wGRS with the selected instrumental variables (63, 104, 205 SNPs described in “The selection of instrumental variables” section) [18] using the PLINK software (http://zzz.bwh.harvard.edu/plink/) in the UK Biobank individuals; briefly, the wGRS was calculated by the summation of the individual’s effect alleles, weighted by corresponding effect sizes [25], and we used the option --score sum to obtain the sum of valid per-allele scores [6]. Then, we applied the multivariable linear regression adjusting for sex, age, weight, smoking, alcohol, and physical activity to assess the genetically determined birth weight on BMD-related parameters (LS/LaS BMC, LS/LaS BA, LS/LaS BMCadjArea, and LS/LaS BMD) in 19,851 participants with LS BMD measurement, 20,874 participants with LaS BMD measurement, and 173,874 participants with heel BMD measurement. In addition, we implemented a two-sample MR analysis to detect the association of genetically determined birth weight with hip bone areas and femoral neck BMD. The information on the instrumental variables was shown in Additional file 1: Table S5-S7. The inverse variance weighted (IVW), MR-Egger, MR-Egger intercept, and MR-PRESSO methods were used for the two-sample MR analyses [26,27,28]. All MR analyses were conducted using the “TwoSampleMR” and “MR-PRESSO” packages in R [27, 29]. In this analysis, P < 0.05 was considered as the significant level. Additionally, we also calculated the statistical power for MR analyses using the “online sample size and power calculator for Mendelian randomization” (https://sb452.shinyapps.io/power/), with the parameters of the sample size, the proportion of variance explained (R2), the causal effect, and the significance level (0.05).
Genetic pleiotropy assessment
We conducted a further sensitivity analysis to address the horizontal pleiotropy assumption of the MR analysis. According to the existing data in the UK Biobank, we selected the bone metabolic markers including alkaline phosphatase (Field ID 30610), vitamin D (Field ID 30890), oestradiol (Field ID 30800), calcium (Field ID 30680), and inflammatory markers including C-reactive protein (Field ID 30710), rheumatoid factor (Field ID 30820), and urea (Field ID 30670) as the confounding risk factors. We performed a linear regression to test the associations between the genetically determined birth weight (wGRS) and the potential risk factors of bone parameters in a one-sample MR analysis. Additionally, we also evaluated the pleiotropic association of instrumental variables (SNPs associated with birth weight) with potential confounders (P < 5 × 10−8). Previous studies have suggested some risk factors for bone mineral density, including anorexia nervosa [30], alcohol [31, 32], body mass index (BMI) [33], Crohn’s disease [34], education years [35], inflammatory bowel disease [36, 37], rheumatoid arthritis [38], smoking [39], ulcerative colitis [36, 37], and 25-hydroxyvitamin D [40]. We checked the association of the SNPs in the GWAS catalog (https://www.ebi.ac.uk/gwas/) to exclude any of the pleiotropic instrumental variables (Additional file 1: Table S8) [41,42,43,44,45,46,47,48,49].
Maternal SNPs and offspring birth weight
With the UK Biobank individual-level dataset, we also examined whether the intrauterine environment (proxied by the maternal SNPs that influence offspring birth weight) was associated with the offspring adult BMD. Firstly, we confirmed the parent-offspring pairs in the UK Biobank with the released kinship information file “ukbA_rel_sP.txt” (https://biobank.ndph.ox.ac.uk/ukb/ukb/docs/ukbgeneinstruct.html). The kinship information file contained kinship coefficients and estimates of the proportion of SNPs with zero identical-by-state (IBS0). The parent-offspring pairs were defined using the formula \(\frac{1}{2^{5/2}}<\mathrm{kinship}\ \mathrm{coefficient}<\frac{1}{2^{3/2}}\) and IBS0 < 0.1 as recommended by Manichaikul et al. [50]. We then excluded the parent-offspring pairs who had less than 15 years of age difference between the mother’s/father’s age and the son’s/daughter’s age. We kept all the offsprings when one mother/father had one or more children. Finally, we obtained 6219 parent-offspring pairs of European ancestry, that is, 4284 mother-offspring pairs and 1935 father-offspring pairs. The detailed process of quality control was provided in Fig. 1. Eventually, we calculated the unweighted genetic risk scores (GRS) with 31 maternal-only effect SNPs [18] for the mothers/fathers in the parent-offspring pairs. Then, we examined the association of the maternal SNP-determined birth weight (exposure) and offspring heel BMD (outcome), either controlling for the offspring GRS or not, in the linear mixed-effects model (with a random effects term to account for within sibling variance) adjusting for the offspring age, sex, weight, and assessment centers in the analysis. Of note, some parent-offspring pairs might miss the covariates data. We also performed MR analyses with 72 and 205 SNPs (see the “The selection of instrumental variables” section) as additional evidence. All statistical analyses were performed with the “lme4” package in R.
Results
The association of birth weight with BMC, bone area, BMCadjArea, and BMD
The results from the restricted cubic spline analysis demonstrated that there was no evidence of a non-linear relationship between birth weight and BMD-related parameters (BMC, bone area, BMCadjArea, and BMD) (P > 0.05, Additional file 1: Fig. S1-S3). In the multivariable linear regression analyses at the site of the lumbar spine (LS), we found that higher birth weight was associated with higher adult LS BMC (β-coefficient = 0.06, P = 2.87 × 10−7) and higher LS area (β-coefficient = 0.17, P < 2.00 × 10−16) (Table 1). However, when BMC was adjusted for the bone area (BMCadjArea), the direction of association between birth weight and BMCadjArea turned negative with significance (β-coefficient = − 0.01, P = 2.27 × 10−14) (Table 1). Similarly, if we looked at the BMD directly, we found a higher birth weight was associated with lower adult LS BMD (β-coefficient = − 0.04, P = 0.007) (Table 1). At the site of the lateral spine (LaS), the association between birth weight and LaS BMC was not observed (P = 0.83); however, when the area was adjusted, we found a strong negative association between birth weight and LaS BMCadjArea (β-coefficient = − 0.05, P = 0.001) (Table 1). In addition, we found that higher birth weight was associated with lower adult LaS BMD (β-coefficient = − 0.03, P = 0.02) (Table 1). At the site of the heel, the higher birth weight was associated with QUS-estimated BMD (β-coefficient = − 0.06, P < 2.00 × 10−16), and the direction of the effect was consistent with the findings from DXA-measured spine BMD (Table 1).
Fetal SNP-determined birth weight, BMC, bone area, BMCadjArea, and BMD
A total of 63 fetal-only effect SNPs explained 1.7% of the variance (R2) of the own birth weight, and the corresponding F-statistic was 88. We calculated the weighted genetic risk score (wGRS) for each individual with these 63 fetal-only effect SNPs within the UK Biobank dataset, and the weighted effect size (beta value) of each SNP on birth weight was extracted from Warrington et al. [18]. We conducted a multivariable linear regression analysis to test the relationship between fetal-only effect SNP-determined birth weight and spine DXA scan parameters and QUS-estimated heel BMD. We observed that higher genetically determined birth weight was significantly associated with higher LS area (β-coefficient = 0.15, P = 1.26 × 10−6) but associated with lower LS BMCadjArea (β-coefficient = − 0.19, P = 6.03 × 10−5) and lower LS BMD (β-coefficient = − 0.10, P = 0.01) (Table 2). At the site of the lateral spine (LaS), we also found that higher fetal SNP-determined birth weight was associated with lower LaS BMCadjArea (β-coefficient = − 0.17, P = 0.0001) and lower LaS BMD (β-coefficient = − 0.12, P = 0.0003) (Table 2). Similarly, the higher fetal-only effect SNP-determined birth weight was found to be associated with lower QUS-estimated BMD at the site of the heel (β-coefficient = − 0.11, P = 3.33 × 10−13) (Table 2). Additionally, we also performed the multivariable linear regression analysis of the wGRS calculated with 104 and 205 SNPs (see the “Methods” section) and obtained similar trends of associations (Table 2). The 104 and 205 SNPs explained 2.4% and 3.6% of the variance (R2) of the own birth weight, respectively, and the corresponding F-statistics were 58 and 75. Furthermore, we conducted a sensitivity analysis to address the horizontal pleiotropy assumption of the MR analysis. Within the selected confounding risk factors such as the bone metabolic markers (alkaline phosphatase, vitamin D, oestradiol, calcium) and inflammatory markers (C-reactive protein, rheumatoid factor, urea), we only found a significant relationship between the fetal-only effect SNP-determined birth weight (wGRS) and serum calcium (Additional file 1: Table S9). Then, we included serum calcium in the one-sample MR model (adjusting the covariates sex, age, weight, smoking, alcohol, physical activity, and serum calcium), and we obtained similar results as before (Additional file 1: Table S10). Therefore, by leveraging the genetic instrumental variables with fetal-SNPs, the wGRS analysis provided additional evidence to support the observational findings, suggesting that the observational findings that higher birth weight was associated with lower BMD might not be biased by potential confounding factors or reverse causality.
Because bone area data at the hip cannot be downloaded from the UK Biobank, we tested the association of 63 fetal-only effect SNP-determined birth weight with bone area and BMD at the hip using a two-sample MR analysis. We found that the 63 fetal-only effect SNP-determined higher birth weight was associated with the higher total hip area (β-coefficient = 0.15, P = 0.005), the intertrochanteric area (β-coefficient = 0.13, P = 0.0009), and the trochanter area (β-coefficient = 0.11, P = 0.03) in two-sample MR-PRESSO analysis (Fig. 2 and Additional file 1: Table S11). Additionally, a higher 63 fetal-only effect SNP-determined birth weight was associated with lower femoral neck BMD at a nominal significance (β-coefficient = − 0.09, P = 0.06). Furthermore, we also performed a two-sample MR analysis with 104 and 205 SNPs (see the “Methods” section) as the instrumental variables and obtained similar trends of association (Fig. 2). The MR-Egger intercept also showed no evidence of directional horizontal pleiotropy for the above three types of fetal instruments (MR-Egger intercept; P > 0.05), and the statistical power was 94–100% as evaluated by the online tool (Additional file 1: Table S11). Additionally, we evaluated the pleiotropic association of instrumental variables (SNPs associated with birth weight) with potential confounders (such as anorexia nervosa, alcohol, BMI, Crohn’s disease, education year, inflammatory bowel disease, rheumatoid arthritis, smoking, ulcerative colitis, and 25-hydroxyvitamin D). We found that rs1547669, rs6911024, and rs9267812 were associated with rheumatoid arthritis, rs516246 was associated with alcohol and Crohn’s disease, and rs7903146 was associated with BMI. Likewise, excluding these SNPs did not significantly change the results of the MR analysis (Additional file 1: Table S12).
Maternal SNP-determined birth weight and BMD
In order to test whether maternal-only effect SNP-determined offspring birth weight was associated with the offspring BMD, we calculated the unweighted genetic risk scores (GRS) with 31 maternal-only effect SNPs in 4284 mother-offspring pairs and 1935 father-offspring pairs. We first checked whether the maternal GRS was associated with offspring birth weight and found that the maternal GRS was indeed strongly associated with increased offspring birth weight before (β-coefficient = 0.01, P = 5.07 × 10−7) and after conditioning on offspring GRS (β-coefficient = 0.01, P = 0.0003) in 4284 mother-offspring pairs (Table 3). In the linear mixed model (exposure: maternal GRS, outcome: offspring heel BMD) adjusting for age, sex, weight, and assessment centers, we found that higher maternal-only effect SNP-determined offspring birth weight was associated with lower adult BMD in both before (β-coefficient = − 0.001, P = 0.03) and after (β-coefficient = − 0.001, P = 0.04) adjusting for offspring GRS in 4284 mother-offspring pairs (Table 3). Additionally, we also performed a linear mixed model analysis of the GRS calculated with 72 and 205 SNPs (see the “Methods” section) and obtained similar trends of association (Table 3). Here, 31, 72, and 205 SNPs explained 1.0%, 1.8%, and 3.0% of the variance (R2) of the offspring birth weight, respectively, and the corresponding F-statistics were 74, 58, and 34. We did not detect any significant relationship in 1935 father-offspring pairs with these types of maternal SNPs (all P > 0.05) (Table 3).
Discussion
In this study, the observational analyses suggested that higher birth weight was associated with increased adult spine bone area but decreased spine BMCadjArea, spine BMD, and heel BMD. By leveraging fetal-only effect SNPs as the genetic instrumental variables, the birth weight determined by own genetic factors was found to be positively associated with offspring bone area but negatively associated with BMD. Furthermore, by using the parent-offspring pairs, our results suggested that the maternal intrauterine environment, as proxied by maternal genetic variants that influence offspring birthweight, was also a determinant of adult BMD. In other words, the observational relationship between birth weight and BMD could be reflected by both the direct fetal and indirect maternal genetic effects.
Lower birth weight often occurs along with inferior intrauterine environments, for example, malnutrition or stress during pregnancy. The Developmental Origin of Health and Disease (DOHaD) hypothesis proposes that exposure to a hostile uterine environment would change the fetus’ structure, function, and metabolism in ways that lead to significantly adverse consequences on an individual’s short- and long-term health, followed by the reduced fetal growth rate, low body weight at birth, and the high risk of future related diseases [51, 52]. The hypothesis was applied well in the studies of cardiovascular diseases [9, 53]. Other epidemiologic observations also showed that lower birth weight was associated with an increased risk of stroke, type 2 diabetes mellitus (T2DM), obesity, and hypertension [54,55,56,57].
A previous study has tried to investigate the relationship between birth weight and osteoporosis risk, but as they listed as a limitation, they failed to analyze the DXA data, which is the gold standard for the diagnosis of osteoporosis [58]. In our study design, the observation that higher birth weight was associated with lower adult BMD was consistently shown in two independent study samples: the DXA scan data (20,367 and 21,364 participants) and QUS measurement data (177,675 samples). With DXA scan data, we could also look at the association between birth weight and BMC/bone area. It is better to pay close attention to BMC/bone area because Dennison et al. demonstrated that birth weight had greater contributions to bone size and mineral content than to bone mineral density [10], and birth weight was previously detected to be associated with the forearm and tibial bone size in the old population aged 65–73 years [59]. In our study, we found that higher birth weight was associated with higher adult LS BMC and higher LS area. Interestingly, the effect of birth weight on bone size was much larger than the effect on BMC after standardization. At the site of the lateral spine (LaS), the association between birth weight and LaS BMC was not observed, but birth weight was positively associated with the LaS area. These results suggested that higher birth weight would lead to lower adult BMD; this might, to some extent, be because of the larger effect of birth weight on bone size than on BMC. Coincidentally, Steer et al. reported similar results in teenagers that birth weight was positively related to bone size but inversely related to cortical BMD, but they suggested that the relationship was mediated by the effects on late onset of puberty and bone resorption [14]. In addition, we included an alternative parameter BMCadjArea (BMC adjusted for bone area) in our analysis, and as expected, the results for BMCadjArea were consistent with the result for BMD, because BMCadjArea is conceptually similar to BMD.
Mendelian randomization (MR) is an approach to test the causality of an observed relationship between an exposure and an outcome [17]. Recently, some studies have tried to leverage the MR method to detect the relationship between birth weight and non-communicable diseases (i.e., atrial fibrillation, type 2 diabetes, and cardiovascular diseases) to address the DOHaD hypothesis [60,61,62]. However, these MR analyses only took the offspring genetic variants to instrument the birth weight. An optimal way is to also take the maternal genotypes that effect the offspring’s birth weight to proxy the intrauterine environment to test aspects of DOHaD [63]. Benefiting from the study by Warrington et al., the genetic effects on birth weight were partitioned into direct fetal genetic component and indirect maternal genetic component [18]. They identified a total of 205 autosomal SNPs associated with birth weight. Of which, 63 autosomal SNPs that have a fetal-only effect on own birth weight represent the effect of own genotype on own birth weight after adjusting the maternal effect on birth weight, 31 autosomal SNPs that had a maternal-only effect on offspring birth weight represent the intrauterine environment, and 41 autosomal SNPs have both fetal and maternal effects [18]. By leveraging these SNPs with different effects, we found that higher fetal-only effect SNP-determined birth weight was associated with lower LS/FN BMD and heel BMD. We also found that higher maternal-only effect SNP-determined offspring birth weight was associated with lower adult BMD after adjusting for offspring genetic effects but did not find any association evidence in father-offspring pairs, which is suggestive of DOHaD mechanisms (related to lower birthweight) through the effect on bone size. Here, we calculated the unweighted GRS in mothers because we were ignorant of the effect size of each allele to avoid the inflation in type 1 error [19]. Previously, Moen et al. suggested that the maternal intrauterine environment, as proxied by maternal genetic variants that influence offspring birth weight, was unlikely to be a major determinant of adverse cardiometabolic outcomes, but offspring birth weight determined by own genetic factors was found to be associated with offspring cardiometabolic risk factors [19]. However, in our study, we suggested that both indirect maternal and direct fetal genetic effects reflected the observational relationship between higher birth weight and lower adult bone mass.
Acknowledged as a limitation of this study, we excluded the individuals with extreme birth weight (very low birth weight < 1.5 kg, low birth weight > 1.5 kg and < 2.5 kg, and birth weight > 4.5 kg) in our study to exclude likely pre-term births and fetal macrosomia. Hovi et al. found that adults born with very low birth weight (< 1.5 kg) had lower LS BMD compared with the participants born at normal birth weight at ages 18.5 to 27.1 [64]. When we checked back in the UK Biobank dataset, we found that the participants born with very low birth weight had lower lumbar spine BMD than the participants born with normal birth weight (1.166 vs 1.179 g/cm2). Therefore, the results should be interpreted with caution when facing the extreme value of birth weight.
Conclusions
In summary, we observed that higher birth weight would lead to lower bone mass in later adulthood; this might be because of the larger effect of birth weight on bone size than on BMC. By leveraging the genetic instrumental variables with maternal- and fetal-specific SNP information, we detected the fetal genetic contributions to birth weight had a casually positive effect on the bone area but a negative effect on BMD. Furthermore, the maternal intrauterine environment, as proxied by maternal genetic variants that influence offspring birth weight, was also a determinant of adult BMD.
Availability of data and materials
The genetic summary statistics for BMD and bone size were downloaded from the GEnetic Factors for OSteoporosis (GEFOS) consortium (http://www.gefos.org/) and deCODE Genetics (https://www.decode.com/summarydata/). The individual-level genetic and phenotype data require permission from the UK Biobank (Application 41376) (https://www.ukbiobank.ac.uk/).
Abbreviations
- BMC:
-
Bone mineral contents
- BMCadjArea:
-
BMC-adjusted bone area
- BMD:
-
Bone mineral density
- DXA:
-
Dual-energy X-ray absorptiometry
- GRS:
-
Genetic risk score
- LaS:
-
Lateral spine
- LS:
-
Lumbar spine
- MR:
-
Mendelian randomization
- QUS:
-
Quantitative ultrasound
References
Sözen T, Özışık L, Başaran NÇ. An overview and management of osteoporosis. Eur J Rheumatol. 2017;4:46.
Zhu X, Zheng H. Factors influencing peak bone mass gain. Front Med. 2021;15:53–69.
Zhu XW, Liu KQ, Yuan CD, Xia JW, Qian Y, et al. General and abdominal obesity operate differently as influencing factors of fracture risk in old adults. iScience. 2022;25:104466.
Hagen KB, Dagfinrud H, Moe RH, Østerås N, Kjeken I, et al. Exercise therapy for bone and muscle health: an overview of systematic reviews. BMC Med. 2012;10:167.
Qian Y, Xia J, Liu KQ, Xu L, Xie SY, et al. Observational and genetic evidence highlight the association of human sleep behaviors with the incidence of fracture. Commun Biol. 2021;4:1339.
Xia J, Xie SY, Liu KQ, Xu L, Zhao PP, et al. Systemic evaluation of the relationship between psoriasis, psoriatic arthritis and osteoporosis: observational and Mendelian randomisation study. Ann Rheum Dis. 2020;79:1460–7.
Zhu X, Bai W, Zheng H. Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications. Bone Res. 2021;9:23.
Hales CN, Barker DJ, Clark PM, Cox LJ, Fall C, et al. Fetal and infant growth and impaired glucose tolerance at age 64. BMJ (Clinical research ed). 1991;303:1019–22.
Barker DJ. The origins of the developmental origins theory. J Intern Med. 2007;261:412–7.
Dennison EM, Syddall HE, Sayer AA, Gilbody HJ, Cooper C. Birth weight and weight at 1 year are independent determinants of bone mass in the seventh decade: the Hertfordshire cohort study. Pediatr Res. 2005;57:582–6.
Jensen RB, Vielwerth S, Frystyk J, Veldhuis J, Larsen T, et al. Fetal growth velocity, size in early life and adolescence, and prediction of bone mass: association to the GH-IGF axis. J Bone Miner Res. 2008;23:439–46.
Jones G, Dwyer T. Birth weight, birth length, and bone density in prepubertal children: evidence for an association that may be mediated by genetic factors. Calcif Tissue Int. 2000;67:304–8.
Saito T, Nakamura K, Okuda Y, Nashimoto M, Yamamoto N, et al. Weight gain in childhood and bone mass in female college students. J Bone Miner Metab. 2005;23:69–75.
Steer CD, Sayers A, Kemp J, Fraser WD, Tobias JH. Birth weight is positively related to bone size in adolescents but inversely related to cortical bone mineral density: findings from a large prospective cohort study. Bone. 2014;65:77–82.
Leunissen RW, Stijnen T, Boot AM, Hokken-Koelega AC. Influence of birth size and body composition on bone mineral density in early adulthood: the PROGRAM study. Clin Endocrinol. 2008;69:386–92.
Steer CD, Tobias JH. Insights into the programming of bone development from the Avon longitudinal study of parents and children (ALSPAC). Am J Clin Nutr. 2011;94:1861s–4s.
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89–98.
Warrington NM, Beaumont RN, Horikoshi M, Day FR, Helgeland Ø, et al. Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors. Nat Genet. 2019;51:804–14.
Moen GH, Brumpton B, Willer C, Åsvold BO, Birkeland KI, et al. Mendelian randomization study of maternal influences on birthweight and future cardiometabolic risk in the HUNT cohort. Nat Commun. 2020;11:5404.
Bai WY, Wang L, Ying ZM, Hu B, Xu L, et al. Identification of PIEZO1 polymorphisms for human bone mineral density. Bone. 2020;133:115247.
Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21:223–42.
Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40:755–64.
Zheng HF, Forgetta V, Hsu YH, Estrada K, Rosello-Diez A, et al. Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture. Nature. 2015;526:112–7.
Styrkarsdottir U, Stefansson OA, Gunnarsdottir K, Thorleifsson G, Lund SH, et al. GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures. Nat Commun. 2019;10:2054.
Choi SW, Mak TS-H, O’Reilly PF. Tutorial: a guide to performing polygenic risk score analyses. Nat Protoc. 2020;15:2759–72.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol. 2015;44:512–25.
Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–8.
Burgess S, Thompson SG. Bias in causal estimates from Mendelian randomization studies with weak instruments. Stat Med. 2011;30:1312–23.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, et al. The MR-base platform supports systematic causal inference across the human phenome. eLife. 2018;7:e34408.
Legroux I, Cortet B. Factors influencing bone loss in anorexia nervosa: assessment and therapeutic options. RMD Open. 2019;5:e001009.
Jang HD, Hong JY, Han K, Lee JC, Shin BJ, et al. Relationship between bone mineral density and alcohol intake: a nationwide health survey analysis of postmenopausal women. PLoS One. 2017;12:e0180132.
McLernon DJ, Powell JJ, Jugdaohsingh R, Macdonald HM. Do lifestyle choices explain the effect of alcohol on bone mineral density in women around menopause? Am J Clin Nutr. 2012;95:1261–9.
Shanbhogue VV, Mitchell DM, Rosen CJ, Bouxsein ML. Type 2 diabetes and the skeleton: new insights into sweet bones. Lancet Diabetes Endocrinol. 2016;4:159–73.
Hoffmann P, Krisam J, Kasperk C, Gauss A. Prevalence, risk factors and course of osteoporosis in patients with Crohn’s disease at a tertiary referral center. J Clin Med. 2019;8:2178.
Ho SC, Chen YM, Woo JL. Educational level and osteoporosis risk in postmenopausal Chinese women. Am J Epidemiol. 2005;161:680–90.
Scott EM, Gaywood I, Scott BB. Guidelines for osteoporosis in coeliac disease and inflammatory bowel disease. Brit Soc Gastroenterol Gut. 2000;46(Suppl 1):i1–8.
Bjarnason I, Macpherson A, Mackintosh C, Buxton-Thomas M, Forgacs I, et al. Reduced bone density in patients with inflammatory bowel disease. Gut. 1997;40:228–33.
Lodder MC, de Jong Z, Kostense PJ, Molenaar ET, Staal K, et al. Bone mineral density in patients with rheumatoid arthritis: relation between disease severity and low bone mineral density. Ann Rheum Dis. 2004;63:1576–80.
Kanis JA, Johnell O, Odén A, Johansson H, De Laet C, et al. Smoking and fracture risk: a meta-analysis. Osteoporos Int. 2005;16:155–62.
Bischoff-Ferrari HA, Kiel DP, Dawson-Hughes B, Orav JE, Li R, et al. Dietary calcium and serum 25-hydroxyvitamin D status in relation to BMD among US adults. J Bone Miner Res. 2009;24:935–42.
Walters RK, Polimanti R, Johnson EC, McClintick JN, Adams MJ, et al. Transancestral GWAS of alcohol dependence reveals common genetic underpinnings with psychiatric disorders. Nat Neurosci. 2018;21:1656–69.
Watson HJ, Yilmaz Z, Thornton LM, Hübel C, Coleman JR, et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nat Genet. 2019;51:1207–14.
Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518:197–206.
Liu M, Jiang Y, Wedow R, Li Y, Brazel DM, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. 2019;51:237–44.
Rietveld CA, Medland SE, Derringer J, Yang J, Esko T, et al. GWAS of 126,559 individuals identifies genetic variants associated with educational attainment. Science. 2013;340:1467–71.
Okbay A, Beauchamp JP, Fontana MA, Lee JJ, Pers TH, et al. Genome-wide association study identifies 74 loci associated with educational attainment. Nature. 2016;533:539–42.
Goyette P, Boucher G, Mallon D, Ellinghaus E, Jostins L, et al. High-density mapping of the MHC identifies a shared role for HLA-DRB1* 01: 03 in inflammatory bowel diseases and heterozygous advantage in ulcerative colitis. Nat Genet. 2015;47:172–9.
Manousaki D, Mitchell R, Dudding T, Haworth S, Harroud A, et al. Genome-wide association study for vitamin D levels reveals 69 independent loci. Am J Hum Genet. 2020;106:327–37.
Eyre S, Bowes J, Diogo D, Lee A, Barton A, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. 2012;44:1336–40.
Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, et al. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010;26:2867–73.
Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of in utero and early-life conditions on adult health and disease. N Engl J Med. 2008;359:61–73.
Mandy M, Nyirenda M. Developmental origins of health and disease: the relevance to developing nations. Int Health. 2018;10:66–70.
Arima Y, Fukuoka H. Developmental origins of health and disease theory in cardiology. J Cardiol. 2020;76:14–7.
Hales CN, Barker DJ. Type 2 (non-insulin-dependent) diabetes mellitus: the thrifty phenotype hypothesis. Diabetologia. 1992;35:595–601.
Goldstein LB, Adams R, Alberts MJ, Appel LJ, Brass LM, et al. Primary prevention of ischemic stroke: a guideline from the American Heart Association/American Stroke Association stroke council: cosponsored by the atherosclerotic peripheral vascular disease interdisciplinary working group; cardiovascular nursing council; clinical cardiology council; nutrition, physical activity, and metabolism council; and the quality of care and outcomes research interdisciplinary working group: the American Academy of Neurology affirms the value of this guideline. Stroke. 2006;37:1583–633.
Kensara OA, Wootton SA, Phillips DI, Patel M, Jackson AA, et al. Fetal programming of body composition: relation between birth weight and body composition measured with dual-energy X-ray absorptiometry and anthropometric methods in older Englishmen. Am J Clin Nutr. 2005;82:980–7.
Curhan GC, Willett WC, Rimm EB, Spiegelman D, Ascherio AL, et al. Birth weight and adult hypertension, diabetes mellitus, and obesity in US men. Circulation. 1996;94:3246–50.
Yu XH, Wei YY, Zeng P, Lei SF. Birth weight is positively associated with adult osteoporosis risk: observational and Mendelian randomization studies. J Bone Miner Res. 2021;36:1469–80.
Oliver H, Jameson KA, Sayer AA, Cooper C, Dennison EM. Growth in early life predicts bone strength in late adulthood: the Hertfordshire cohort study. Bone. 2007;41:400–5.
Zanetti D, Tikkanen E, Gustafsson S, Priest JR, Burgess S, et al. Birthweight, type 2 diabetes mellitus, and cardiovascular disease: addressing the Barker hypothesis with Mendelian randomization. Circul: Genom Precis Med. 2018;11:e002054.
Kember RL, Levin MG, Cousminer DL, Tsao N, Judy R, et al. Genetically determined birthweight associates with atrial fibrillation: a Mendelian randomization study. Circul: Genom Precis Med. 2020;13:e002553.
Wang T, Huang T, Li Y, Zheng Y, Manson JE, et al. Low birthweight and risk of type 2 diabetes: a Mendelian randomisation study. Diabetologia. 2016;59:1920–7.
D’Urso S, Wang G, Hwang LD, Moen GH, Warrington NM, et al. A cautionary note on using Mendelian randomization to examine the Barker hypothesis and developmental origins of health and disease (DOHaD). J Dev Orig Health Dis. 2021;12:688–93.
Hovi P, Andersson S, Järvenpää AL, Eriksson JG, Strang-Karlsson S, et al. Decreased bone mineral density in adults born with very low birth weight: a cohort study. PLoS Med. 2009;6:e1000135.
Acknowledgements
We thank the UK Biobank database, the GEnetic Factors for OSteoporosis (GEFOS) Consortium, and deCODE Consortium. We also thank the High-Performance Computing Center at Westlake University for the facility support and technical assistance.
Funding
This work was supported by the National Natural Science Foundation of China (81871831) and supported by the Chinese National Key Technology R&D Program, Ministry of Science and Technology (2021YFC2501702). This work was also supported by the joint program of the National Natural Science Foundation of China (NSFC) and the Israel Science Foundation (ISF) (NSFC#32061143019 and ISF#3408/20).
Author information
Authors and Affiliations
Contributions
H.-F.Z. conceptualized and designed the study. J-W.X. conducted the analyses and drafted the paper. L.Z, J.L., C-D.Y., X.-W.Z., Y.Q., S.K., J.-X.G., L.X., S-Y.X., J-H.G., K-Q.L., G-B.C., and D.K. helped in the data acquisition, result interpretation, and manuscript editing. The authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All individuals provided written informed consent. The North West Multi-Centre Research Ethics Committee approved the UK Biobank ethical application (reference number: 16/NW/0274).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Fig. S1.
Restricted cubic spline analysis to show the relationship between birth weight and Lumbar spine bone parameters. Fig. S2. Restricted cubic spline analysis to show the relationship between birth weight and Lateral spine bone parameters. Fig. S3. Restricted cubic spline analysis to show the relationship between birth weight and heel BMD. Table S1. Characteristics of the single nucleotide polymorphisms associated with the birth weight. Table S2. Detailed information on the GWAS summary-statistic data for birth weight and bone parameters. Table S3. Characteristics of participants from the UK Biobank dataset in the observational analysis for the Spine DXA data. Table S4. Characteristics of participants from the UK Biobank dataset in the model 3 observational analysis for heel BMD measured by QUS. Table S5. Summary statistics for the 63 SNPs used in instrumental variable analysis to assess the effect of fetal birth weight on the hip bone parameters. Table S6. Summary statistics for the 104 SNPs used in instrumental variable analysis to assess the effect of fetal birth weight on the hip bone parameters. Table S7. Summary statistics for the 205 SNPs used in instrumental variable analysis to assess the effect of fetal birth weight on the hip bone parameters. Table S8. P values for the associations of 205 birth weight variants with potential risk factors. Table S9. The association between fetal SNPs determined birth weight and potential risk factors in UK biobank dataset. Table S10. The association between fetal SNPs determined birth weight and adult BMC, Bone area, BMCadjArea and BMD in UK biobank dataset (confounders including calcium level). Table S11. Causal associations of the fetal genotype effects on birth weight (63, 104, 205 instrument variables) with hip bone parameters in two-sample MR analyses. Table S12. Causal associations of the fetal genotype effects on birth weight (62, 102, 200 instrumental variables which excluded the potential pleiotropic variants) with hip bone parameters in two-sample MR analyses.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Xia, JW., Zhang, L., Li, J. et al. Both indirect maternal and direct fetal genetic effects reflect the observational relationship between higher birth weight and lower adult bone mass. BMC Med 20, 361 (2022). https://doi.org/10.1186/s12916-022-02531-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12916-022-02531-w