Introduction

The development of bone mass throughout childhood is important in determining the peak bone mass achieved in early adulthood [1], and simulation models have demonstrated the potential of small increases in peak bone mass to delay the onset of osteoporosis and therefore decrease the risk of fracture in the elderly [2]. Even during childhood, a reduction in bone mass is associated with an increased fracture risk from all levels of trauma severity [3, 4].

Maternal smoking during pregnancy has been shown to have a detrimental influence on the accrual of bone mass in utero. Two studies in the Southampton Women’s Survey reported associations between maternal smoking and decreased whole body bone mineral content (BMC) in neonatal offspring [5, 6]. The earlier of the two studies also found a similar relationship with bone mineral density (BMD) [5], but the more recent and larger study did not [6]. Little is known about longer term effects, although in a Tasmanian cohort of 330 participants, relationships were found between maternal smoking during pregnancy and reduced offspring femoral neck and lumbar spine BMC and BMD at age 8 years which remained after adjustment for current weight and height [7].

We assessed the associations of maternal smoking in pregnancy with the skeletal size and bone density at mean age 9.9 years of a large cohort of children: the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared the effects of maternal smoking with those of paternal smoking during pregnancy since the paternal exposure would not be expected to influence foetal development via an intrauterine mechanism. Hence, stronger maternal associations would provide evidence of a direct intrauterine effect on bone development, whilst similar-sized maternal and paternal associations would indicate relationships driven by shared familial, social, genetic and environmental factors. This method has been used effectively to study the influences of maternal smoking on other outcomes in the ALSPAC [810], and its validity is demonstrated by the much greater association of maternal compared with paternal smoking in pregnancy with offspring birth weight, which is known to be influenced by maternal smoking via an intrauterine mechanism [11].

Materials and methods

The ALSPAC

The ALSPAC is a prospective birth cohort study aiming to investigate environmental and inheritable influences on the health and development of children. It has been previously described in full elsewhere and on the web site www.alspac.bris.ac.uk. Pregnant women with expected delivery dates between 1 April 1991 and 31 December 1992 and living in a defined area of Avon including the city of Bristol were eligible for recruitment to the study. A total of 14,541 women were enrolled, and 13,678 of these had a singleton live birth. Ethical approval for the study was obtained from the ALSPAC Law and Ethics Committee and from local ethics committees. At age 9 years, all children with known addresses who were still participating were invited to a “Focus @ 9” clinic, and 7,121 of the singleton children attended. Of these, 6,868 underwent a full-body dual-energy X-ray absorptiometry (DXA) scan.

DXA measurements

Whole body DXA scans were carried out using a Lunar Prodigy scanner (GE Healthcare Bio-Sciences Corp., USA), and after exclusion of those containing artefacts, movement or skeletal irregularities, there remained 6,775 scans. The total body less head (TBLH) region was used to represent the child’s total bone mass, with the head excluded as bone development here is different from the rest of the skeleton and less likely to be influenced by environmental factors. In contrast to the overall skeleton which primarily comprises cortical bone, the spine subregion which was also analysed has a relatively high proportion of trabecular bone. One operator reanalysed the scans to check and adjust automated placement of body regions; in the case of the spinal region, the upper border comprises the cervicothoracic junction, the lower border the lumbosacral junction, and the lateral borders the bone/soft tissue interface. Since curvature in the image of the spine leads to contamination of the spinal region with the ribs, only images with minor or no curvature are included in the analysis of spinal outcomes. Measurements for TBLH and spine BMC, bone area (BA) and areal BMD were subsequently calculated. For both regions, area-adjusted BMC (ABMC) was also derived as a measure of volumetric BMD by using linear regression to adjust BMC for BA and adding the residuals to the mean BMC for the region. The coefficient of variation for TBLH BMD was 0.84% based on 122 pairs of scans repeated on the same day. At the same time as the DXA scan, the child’s standing height (without shoes) was measured using a Harpenden Stadiometer (Holtain Ltd., UK) and weight (unshod and in light clothing) was measured using a Tanita Body Fat Analyzer (model TBF 305, Tanita UK Limited, UK).

Maternal and paternal smoking

At 18 weeks’ gestation, the mothers were sent a postal questionnaire which asked how many times per day they had smoked in the first trimester and in the last 2 weeks, representing smoking during the second trimester. At 32 weeks’ gestation, another postal questionnaire asked how many cigarettes per day the woman was currently smoking, representing smoking during the third trimester. Variables describing smoking in any trimester and smoking in all trimesters were derived from these responses, with one or more cigarettes smoked per day considered as smoking regularly. A questionnaire completed by the mother’s partner at 18 weeks’ gestation asked if he had smoked regularly at any time in the last 9 months. The mother was also asked if her partner smoked in her 18-week questionnaire, and a positive response from either the partner or the mother was assumed sufficient to indicate that the partner smoked regularly during the pregnancy.

Other variables

Maternal and paternal height, weight and highest educational qualifications, household social class, father’s age and the mother’s parity were obtained from questionnaires administered during pregnancy. Household social class was defined from the highest parental occupation, on a scale from I to V, with I indicating a professional/managerial role and V being unskilled manual. Maternal and paternal body mass index (BMI) were calculated as weight (kg)/height (m)2. The child’s sex was obtained at the time of birth, and the child’s birth weight, gestational age and the mother’s age at delivery were abstracted from obstetric records.

In the questionnaire administered at 18 weeks’ gestation, the mother was asked how many hours per week she spent engaging in strenuous physical activity. The questionnaire also asked the number of hours per week the mother spent in a number of specific types of leisure activity, each of which was assigned a MET score [12], and a weighted activity index was developed by multiplying the MET score by the number of hours of activity per week. Dietary information for the mothers was obtained from a food frequency questionnaire administered at 32 weeks’ gestation which asked how often they consumed each of the 43 food groups. Using nutrient information on standard-sized portions, the mother’s total weekly energy, carbohydrate, fat and protein intakes were derived [13]. Although the main analysis did not adjust for these variables, since the equivalent paternal information was not available, an additional analysis was performed in which the relationships of maternal smoking in pregnancy with offspring bone outcomes were adjusted for maternal physical activity (strenuous activity of 3 h or more per week and weighted activity index) and diet (weekly energy, carbohydrate, fat and protein intake) during pregnancy.

Pubertal stage data for the children were obtained from Tanner stage questionnaires administered to the parents at 116 months and were based on pubic hair development for boys and breast development for girls, or pubic hair development if this was unavailable. For girls, age at menarche was derived from a series of questionnaires administered between the ages of 8 and 17 years which asked if the daughter had started her menstrual periods and, if so, the age she was at her first menstrual period. Where there was disagreement between questionnaires, the age given on the earliest questionnaire was used. Most children (99% of boys and 96% of girls) with pubertal stage information were either pre- or early pubertal (Tanner stage 1 or 2). For this reason, and due to the high proportion of missing pubertal stage data, this has not been adjusted for in the main regression analysis, but an additional analysis was performed which adjusted for pubertal stage and, for girls, whether menarche occurred at age ≤10 years.

Paternity

If, when asked in a questionnaire administered in pregnancy, the mother had not confirmed her partner to be the child’s biological father, all paternal information (smoking status, BMI, age, height and education) was treated as missing.

Statistical analysis

We assessed maternal and paternal smoking associations with offspring bone outcomes separately and also in combined mutually adjusted regression models. We used sex-specific models due to evidence of interactions of maternal and paternal smoking with the child’s sex and adjusted first for the child’s age only, then additionally for the potential confounders of household social class, parity and maternal/paternal age, height, BMI and education. Maternal factors were included in maternal exposure models, paternal factors in paternal exposure models, and both maternal and paternal factors in combined models. To explore mediating relationships, we additionally adjusted for the child’s birth weight and gestational age and then finally included the child’s height and weight as potential mediators. Since there was little change in regression coefficients between the simple age-adjusted model and the model adjusting for all potential confounding factors (full results for all four models available from authors), only the confounder-adjusted model (age and all other potential confounders, model 1) and the two additional models exploring potential mediation by birth weight and gestational age (model 2) and by weight and height at age 9.9 (model 3) are presented. Sex-specific standard deviation (SD) scores of TBLH and spine BMC, BA, BMD and ABMC were used as outcomes.

We used multivariate multiple imputation of missing data to impute data for all children who attended the 9-year clinic and also analysed the complete cases with no missing data on any of the exposures, outcomes or covariates to compare findings from the fully observed data with those from partially imputed data. Multiple imputation was used to increase the efficiency of the model estimates and reduce selection bias, which can be present in complete case analysis when data are not missing completely at random. The multiple imputation method is valid provided that the reasons for missingness in the data can be explained by other observed variables [14]. Detailed methods for this procedure are described in the Electronic supplementary material (ESM). All analyses were carried out in Stata version 11.0 (StataCorp LP, USA).

Results

Table 1 shows the characteristics of the 7,121 children who attended the 9-year clinic. There were 6,101 sets of parents for whom both maternal and paternal smoking information was available; for 3,576 (58.6%) of these neither parent smoked, for 369 (6.0%) only the mother smoked, for 1,313 (21.5%) only the father smoked, and for 843 (13.8%) both parents smoked. Mothers who smoked at any time during pregnancy were younger and shorter on average, more likely to be of a manual social class and less likely to have an A-level or higher qualification than mothers who did not smoke (ESM Web Table 2). Pre-pregnancy BMI did not differ between mothers who smoked and those who did not. Children of mothers who smoked were lighter at birth and older, heavier and had higher fat mass at the time of the DXA scan on average.

Table 1 Characteristics of children who attended the clinic at age 9 years and their parents (N = 7,121)

Pairwise correlations of total body and spinal bone measures are given in ESM Web Table 3, and correlations of these measures with child and parental characteristics are shown in ESM Web Table 4. The child’s height and weight were strongly positively correlated with TBLH and spine BMC and BA and moderately with TBLH and spine BMD. Higher birth weight, longer gestation and greater age at DXA scan were all associated with increased TBLH BMC, BA and BMD.

Multiple imputation analysis of maternal and paternal smoking in relation to TBLH bone outcomes is shown in Table 2 and the analysis of spinal outcomes shown in Table 3. In girls, maternal smoking in any trimester was associated with increased TBLH and spine BMC, BA and BMD in confounder-adjusted models (model 1). Mean differences in TBLH BMC and BA were slightly higher for mothers who smoked in all trimesters of pregnancy, but other associations were similar. In boys, maternal smoking in any trimester was not robustly associated with any TBLH or spinal bone outcomes. P values for sex differences were 0.007, 0.003 and 0.085 for TBLH BMC, BA and BMD and 0.036, 0.035 and 0.119 for spine BMC, BA and BMD, respectively, in models with maternal smoking in any trimester as the exposure. In combined confounder-adjusted models (model 1) for girls, there were greater paternal smoking associations with TBLH BMC, BA and BMD and spine BMD compared with those for maternal smoking, whilst maternal associations were larger than paternal associations with spine BMC and BA. However, all P values for differences between maternal and paternal effects were >0.15. There were no strong associations of paternal smoking with bone outcomes in boys. On additional adjustment for the child’s birth weight and gestational age (model 2), there were increases in maternal associations, whilst paternal associations did not change. In boys, maternal smoking in all trimesters was positively associated with TBLH BMC, BA and BMD after adjustment for birth weight and gestational age. In fully adjusted models including offspring height and weight at age 9.9 years (model 3), all maternal relationships attenuated to the null, although a weak association remained with spine BA in girls. Paternal associations were similarly attenuated, and although evidence remained of an association with TBLH BA, this weakened in combined models. There were no associations between parental smoking during pregnancy and TBLH or spine ABMC, except for a weak positive association between paternal smoking and spine ABMC in girls. These models are not included in the tables (full data available from authors on request). Variance inflation factors for maternal smoking and paternal smoking in mutually adjusted models were all <2, indicating that there was low colinearity between the two parental smoking variables in these models.

Table 2 Sex-specific associations of maternal and paternal smoking with total body less head bone outcomes at age 9.9 years in multiple imputation analysis (boys N = 3,530; girls N = 3,591)
Table 3 Sex-specific associations of maternal and paternal smoking with spinal bone outcomes at age 9.9 years in multiple imputation analysis (boys N = 2,772; girls N = 2,715)

To examine whether the associations in girls were mediated mainly by the child’s height or weight, we adjusted for each individually (plus confounders). On adjustment for height without weight, mean differences in TBLH BMC, BA and BMD associated with maternal smoking in any trimester were 0.13 SD, 0.12 SD and 0.12 SD, respectively (all P < 0.01). However, on adjustment for weight without height, mean differences were −0.02 SD, −0.03 SD and 0.00 SD (all P > 0.2), suggesting that the positive associations of maternal smoking with offspring bone mass are driven by the child’s weight at age 9.9 years. Mean differences in TBLH BMC, BA and BMD associated with paternal smoking on adjustment for height without weight were 0.10 SD, 0.10 SD and 0.10 SD (all P < 0.01), and adjusting for weight without height were 0.01 SD, 0.01 SD and 0.03 SD, respectively (all P > 0.2). A similar pattern occurred in spine BMC, BA and BMD.

In complete case analysis (ESM Web Tables 5 and 6), associations of maternal smoking with TBLH and spinal BMC, BA and BMD were equivalent to those using multiple imputation, but associations of paternal smoking were generally smaller in girls (by up to 0.07 SD). No strong associations of maternal or paternal smoking in pregnancy with bone outcomes were found in boys in the complete case in confounder-adjusted models. In combined confounder-adjusted models for TBLH bone outcomes in girls in the complete case maternal and paternal smoking associations were of a similar size, with little evidence for a difference between parental effects, as in multiple imputation models. However, in models for spinal bone outcomes, there were greater maternal compared with paternal associations, and there was statistical evidence for a difference between parental smoking associations with spinal BA. ESM Web Tables 7 and 8 compare the characteristics of multiply imputed and complete case datasets for TBLH and spinal bone outcomes, respectively, and show that parental educational qualifications tended to be higher in the complete case. We thus investigated the relationships between maternal and paternal smoking and TBLH and spinal BMC, BA and BMD in girls in the complete case and stratified the analysis into two subgroups: families where neither parent had an A-level or higher qualification and families where one or both parents was qualified to A level or above (data not shown). In TBLH models, paternal associations were greater than maternal associations in the stratum with lower parental qualifications, whilst maternal associations were greater in the stratum with higher parental qualifications. In the stratum with less educated parents, there were similar-sized parental smoking associations with spinal bone outcomes, but greater maternal associations in the higher educated stratum. This suggests that the dissimilarities between the multiple imputation and complete case analyses are likely to be due to the differing parental education distributions in the datasets.

Of the 6,741 children whose ethnicity was known, 6,470 (96.0%) were white. Restricting the analysis to children of known white ethnicity did not meaningfully change the model coefficients. Including maternal diet and physical activity during pregnancy in the multiple imputation process and additionally adjusting for these variables in models with maternal smoking as the exposure did not alter the findings. When we repeated the multiple imputation process with pubertal stage (for both boys and girls) and age of menarche (for girls only) included and additionally adjusted for these variables, model coefficients were similar for boys. In models with maternal smoking as the exposure for girls, associations were attenuated by up to 0.07 SD compared with the original multiple imputation analysis, whilst associations of paternal smoking were unchanged.

Discussion

We compared the relationships of maternal and paternal smoking during pregnancy with offspring bone mass at mean age 9.9 years in a large birth cohort and found similar-sized associations of smoking in both parents with increased total body and spinal BMC, BA and areal BMD in girls, but little evidence for any associations in boys.

Maternal smoking during pregnancy was associated with 0.10–0.13 SD increases in TBLH and spinal BMC, BA and BMD in daughters. These relationships were masked by the negative association of maternal smoking with the child’s birth weight and gestational age and increased on adjustment for these factors, whilst effect sizes associated with paternal smoking did not change. This may be due to the negative intrauterine effect on the accrual of bone mass by the foetus [5, 6], which is unique to the maternal smoking exposure. Maternal smoking during pregnancy is known to lead to a smaller child at birth, both through an increased risk of preterm birth and through intrauterine growth retardation [15, 16], and a positive relationship has been reported between birth weight and BMD at the femoral neck and lumbar spine in 8-year-old children [17].

Conversely, relationships of maternal and paternal smoking with offspring bone mass attenuated to the null when the child’s height and weight were included in regression models. BMC, BA and BMD are all related to bone size (as BMD is incompletely adjusted for bone area) and therefore correlate strongly with height and weight. Since no relationships were found between maternal smoking and ABMC, which reflects ‘volumetric’ BMC, it appears that the associations are working through skeletal size rather than density. The relationships were driven mainly by offspring weight, concurring with studies which have demonstrated an association between maternal smoking in pregnancy and increased BMI and risk of overweight in childhood [15, 1825], whilst the child’s height deficit at birth has been shown to track to age 8 years [22]. Studies comparing the associations of maternal and paternal smoking in pregnancy with offspring BMI also report similar effects for both parental exposures [9, 24]. Furthermore, a previous study in ALSPAC found an inverse relationship of parental social position with offspring BMC and BA at age 9.9 years, also acting via the pathway of offspring weight [26]. It therefore seems most plausible that our associations are not explained by intrauterine effects, but rather that unmeasured aspects of the shared family environment which are associated with parental smoking, such as diet or level of physical activity, influence increased weight gain and greater bone mass in the children. Studies have shown that overweight children and adolescents have higher whole body and spinal bone mass [2729] and that BMC is positively related to both lean and fat mass in childhood [30, 31]. Fat mass has been demonstrated to stimulate bone growth in prepubertal children previously in the ALSPAC [32, 33].

There has been a greater association reported between fat mass and bone mineral accrual in girls than in boys during puberty [34, 35], which may in part explain why we found no associations in boys, although one study suggests that this sex difference is not present in prepubertal children [35]. In our cohort, there was also a weaker univariate relationship between maternal smoking and offspring weight in sons than in daughters, so it is also possible that the social characteristics in families where parents smoke have a lesser influence on adiposity in boys than girls. In analysis adjusted for pubertal stage (both genders) and age at menarche (in girls), the associations between maternal smoking and bone outcomes in girls were attenuated, whereas the paternal associations remained similar. This suggests that these positive maternal associations may partly be explained by the association between maternal smoking in pregnancy and earlier age at menarche, which has been shown previously in ALSPAC [36]. Adjustment for pubertal stage in boys did not affect the associations between parental smoking and bone outcomes, and parental smoking was not related to pubertal stage at age 10 years in boys.

Our findings conflict with the study by Jones et al. [7] which indicated negative relationships between maternal smoking in pregnancy and bone mass in 8-year-olds for the total body, femoral neck and lumbar spine, with relationships at the femoral neck and lumbar spine remaining after adjustment for the child’s height and weight. However, they studied a Tasmanian cohort identified at birth as at increased risk of sudden infant death syndrome which contained 65% male offspring and a higher prevalence of maternal smoking during pregnancy (49%) compared with ours (21%). Children of mothers who smoked were lighter at age 8 years in Jones’ study, whereas we found a strong positive relationship between maternal smoking and offspring weight. Jones et al. do not make comparison with paternal smoking or give sex-specific findings.

The differing associations of parental smoking with spinal bone outcomes which we found in the complete case and multiple imputation analyses are probably due to the different populations represented by the two datasets since the complete case includes only around half of the children who attended the 9-year clinic. The rationale for comparing maternal and paternal smoking associations with offspring bone mass was that there is likely to be residual confounding in these relationships from unmeasured factors. Differing distributions of unmeasured confounders in the complete case and multiply imputed datasets could explain the difference between associations seen. Since there were differing educational distributions between the complete case and multiply imputed datasets and we found that parental smoking associations in the complete case differed between strata of parental education levels despite adjusting for all observed confounders, it seems that residual confounding is a possible explanation. Another possible reason for the difference is violation of the multiple imputation assumption that the missing data mechanisms can be explained by other observed variables. However, we verified that missingness in each of the variables with missing data was strongly associated with other observed variables and included a number of predictors of missingness in prediction equations to impute missing data. We therefore expect the multiply imputed datasets to be more representative of the study population and analyses based on these data more accurate.

A limitation to our study was the self-report of smoking by the mothers and fathers. Maternal smoking could be affected by reporting bias since mothers may be aware of the harmful effects of smoking and less likely to respond affirmatively. Nevertheless, where both the mother and father provided information about the father’s smoking status, there was agreement in 94.5% of couples. The study benefitted from its large size, the ability to control for a number of potential confounders and the ability to compare associations of bone outcomes with both maternal and paternal exposures to assess the level of residual confounding.

Conclusions

Our study has found positive associations of maternal smoking during pregnancy with offspring total body and spinal bone mass in girls, with minimal evidence for any associations in boys, and our multivariable analyses and parental comparisons suggest that these associations are largely driven by familial characteristics related to childhood adiposity and unlikely to be due to intrauterine mechanisms. Although our findings do not demonstrate negative effects of maternal smoking in pregnancy on offspring bone mass, its known adverse effects for mothers and offspring health mean than women should be encouraged not to smoke.