Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Association Between Age at Puberty and Bone Accrual From 10 to 25 Years of Age

Association Between Age at Puberty and Bone Accrual From 10 to 25 Years of Age Key Points Question Is puberty timing associated IMPORTANCE Bone health in early life is thought to influence the risk of osteoporosis in later life. with growth-related bone accrual up to adulthood? OBJECTIVE To examine whether puberty timing is associated with bone mineral density accrual up Findings In this cohort study of 6389 to adulthood. participants who underwent repeated bone density scans from ages 10 to 25 DESIGN, SETTING, AND PARTICIPANTS This cohort study used data from the Avon Longitudinal years, later puberty was associated with Study of Parents and Children, a prospective population-based birth cohort initiated in 1991 to 1992 persistently lower bone mineral density, in southwest England. The participants were 6389 healthy British people who underwent regular despite some catch-up during puberty. follow-up, including up to 6 repeated bone density scans from ages 10 to 25 years. Data analysis was performed from June 2018 to June 2019. Meaning People with older pubertal age should be advised on how to EXPOSURES Age at puberty from estimated age at peak height velocity (years). maximize bone density and minimize its decrease in later life to help prevent MAIN OUTCOMES AND MEASURES Gains per year in whole-body bone mineral density (grams per fracture and osteoporosis. square centimeter), assessed by dual-energy x-ray absorptiometry at ages 10, 12, 14, 16, 18, and 25 years and modeled using linear splines. Supplemental content RESULTS A total of 6389 participants (3196 [50.0%] female) were included. The mean (SD) age at Author affiliations and article information are listed at the end of this article. peak height velocity was 13.5 (0.9) years for male participants and 11.6 (0.8) years for female participants. Male participants gained bone mineral density at faster rates than did female participants, with the greatest gains in both male participants (0.139 g/cm /y; 95% CI, 0.127-0.151 2 2 2 g/cm /y) and female participants (0.106 g/cm /y; 95% CI, 0.098-0.114 g/cm /y) observed between the year before and 2 years after peak height velocity. When aligned by chronological age, per 1-year older age at puberty was associated with faster subsequent gains in bone mineral density; the magnitudes of faster gains were greatest between ages 14 and 16 years in both male participants 2 2 2 (0.013 g/cm /y; 95% CI, 0.011-0.015 g/cm /y) and female participants (0.014 g/cm /y; 95% CI, 2 2 0.014-0.015 g/cm /y), were greater in male participants (0.011 g/cm /y; 95% CI, 0.010-0.013 2 2 2 g/cm /y) than in female participants (0.003 g/cm /y; 95% CI, 0.003-0.004 g/cm /y) between ages 16 and 18 years, and were least in both male participants (0.002 g/cm /y; 95% CI, 0.001-0.003 2 2 2 g/cm /y) and female participants (0.000 g/cm /y; 95% CI, −0.001 to 0.000 g/cm /y) between ages 18 and 25 years. Despite faster gains, older age at puberty was associated with persistently lower 2 2 bone mineral density, changing from 0.050 g/cm (95% CI, −0.056 to −0.045 g/cm )loweratage14 2 2 years to 0.047 g/cm (95% CI, −0.051 to −0.043 g/cm ) lower at age 25 years in male participants 2 2 2 and from 0.044 g/cm (95% CI, −0.046 to −0.041 g/cm ) to 0.034 g/cm (95% CI, −0.036 to −0.032 g/cm ) lower at the same ages in female participants. (continued) Open Access. This is an open access article distributed under the terms of the CC-BY License. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 1/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Abstract (continued) CONCLUSIONS AND RELEVANCE People with older pubertal age should be advised on how to maximize bone mineral density and minimize its decrease in later life to help prevent fracture and osteoporosis. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 Introduction Peak bone mass at the end of growth is thought to be an important determinant of later-life risk of 1-5 6,7 fracture and osteoporosis, a bone loss disorder with substantial and increasing health costs. For example, bone remodeling simulations showed that a 10% increase in peak bone mineral density (BMD) would delay osteoporosis by 13 years. Puberty is a key early life milestone that is characterized by endocrine-initiated reproductive maturation and a dramatic skeletal growth spurt 8,9 10-12 in height. Although bone mass potential and puberty timing are both strongly heritable, 13-16 evidence indicates that later puberty may lead to lower BMD in adolescence and adulthood and, thus, an increased risk of osteoporosis later in life. However, the association between puberty timing and long-term bone accrual from early life up to adulthood, including the extent and duration of any catch-up bone accrual by later maturing adolescents, has not been described, to our knowledge. Age at peak height velocity (APHV) is an accurate and precise marker of puberty timing that allows direct 13,17 comparisons of this association between male and female individuals. Yet, to our knowledge, only 2 studies, which included 230 participants who underwent repeated bone scans and 500 male participants assessed once at baseline and once at follow-up, have examined the association between APHV and bone accrual. Therefore, the aim of this study was to examine the association between age at puberty (measured by APHV) and bone accrual from ages 10 to 25 years in a large birth cohort of male and female participants. We assessed the rates of BMD and bone mineral content (BMC) accrual from before puberty up to adulthood relative to pubertal age and examined the association between age at puberty and subsequent BMD and BMC accrual. Our secondary aim was to examine the associations between age at puberty and site-specific bone accrual. Methods Study Sample 20,21 The Avon Longitudinal Study of Parents and Children (ALSPAC) is a prospective birth cohort study that recruited all pregnant women residing within the catchment area of 3 National Health Service authorities in southwest England with an expected date of delivery between April 1991 and December 1992. In total, 15 247 eligible pregnancies were enrolled in ALSPAC (75% response), resulting in 14 973 live births, of whom 14 899 were alive at 1 year of age. Detailed information has been collected from offspring and parents using questionnaires, data extraction from medical records, linkage to health records, and dedicated clinic assessments up to the last completed contact in 2018. Details of all available data can be found in the ALSPAC study website, which includes a fully searchable data dictionary and variable search tool. Ethics approval was obtained from the ALSPAC law and ethics committee and the local National Health Service research ethics committee. Written informed consent was obtained from all participants. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline. Height Measurements Numerous height measures have been obtained from various sources from birth to age 25 years, including from routine data collected from midwives, health visitors, linkage to child health records, JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 2/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual and ALSPAC research clinic visits. Height from research clinic visits (annually or more frequently up to age 14 years and then at mean ages of 16, 18, and 25 years) was measured by accredited fieldworkers to the nearest 0.1 cm using a Harpenden stadiometer (Holtain Ltd). Measured heights were supplemented by extensive maternal and self-reported height records collected throughout the study. For these analyses, we restricted APHV estimation to individuals with height measurements taken between ages 5 and 20 years and at least one height measurement after age 9 years because these were relevant to assessment of APHV. Dual-Energy X-ray Absorptiometry Scans All participants were invited to undergo whole-body dual-energy x-ray absorptiometry (DXA) scans as part of clinic assessments at mean ages 10, 12, 14, 16, 18, and 25 years. Scans were performed using a Lunar Prodigy scanner (Lunar Radiation Corp) and analyzed according to the manufacturer’s standard scanning software and positioning protocols. Scans were reanalyzed as necessary to ensure optimal placement of borders between adjacent subregions, and scans with anomalies were 23,24 excluded. For our primary outcomes, we extracted whole-body (except for the head) BMD (grams per square centimeters) and BMC (grams) at each age. Our secondary outcomes were site- specific BMD and BMC (arms, legs, trunk, ribs, spine, and pelvis) from whole-body DXA, in addition to total hip and femoral neck BMD from up to 3 repeated hip DXA scans performed at ages 14, 18, and 25 years. Confounding Variables Birth weight, ethnicity, socioeconomic position, body mass index (calculated as the weight in kilograms divided by height in meters squared), and diet in early life were hypothesized to potentially confound associations between puberty timing and bone outcomes and were included as model adjustments. Birth weights were recorded to the nearest gram and were extracted from hospital records. Ethnicity was based on the mother’s ethnic background and was reported during a general prenatal questionnaire. Early life socioeconomic position was based on the mother’s highest educational qualifications reported in pregnancy. Body mass index was derived from height and weight measurements taken at the age 7 years research clinic visit. Dietary intake in early life was based on the child’s daily energy intake (kilojoules per day), which was derived from food frequency questionnaires completed by the parents when the child was aged 7 years. Statistical Analysis Data analysis was performed from June 2018, to June 2019. Statistical analyses were performed in 2 stages. The first stage involved measuring age at puberty by estimating APHV, which was done using 13,26 Superimposition by Translation and Rotation (SITAR) growth curve analysis. The SITAR models reduce complex growth data into clinically relevant parameters that represent the timing, intensity, and duration of the pubertal growth spurt, which, in turn, simplifies comparison between individuals and between male and female participants. We estimated APHV through transformation of the random age intercept that reflects individual differences in the timing of the growth spurt. The second stage involved investigating the association between age at puberty and whole- body and site-specific BMD and BMC accrual up to adulthood using 2 sets of mixed-effects linear 27,28 spline regression models. These models summarize nonlinear change through a series of linear splines joined at knot points by estimating person-specific rates of accrual during each growth 27,28 period. We included APHV-by-age splines interaction terms to investigate whether bone accrual differs by age at puberty, and models with sex-by-age splines interaction terms were used to test differences in bone accrual between male and female participants. In the first set of these analyses, we examined gains in BMD and BMC from childhood (prepuberty) to adulthood relative to pubertal age by centering chronological age at APHV. Results from these models can be interpreted as rates of BMD and BMC accrual during 4 periods: 8 years up to 1 year before APHV, 1 year before APHV to 2 years after APHV, 2 to 4 years after APHV, and 4 to JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 3/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 16 years after APHV. The estimated trajectories were plotted for male and female participants. In the second set of analyses, we excluded bone scans taken before APHV to investigate the associations between age at puberty and subsequent bone accrual up to adulthood. Results from these models can be interpreted as rate of BMD and BMC accrual per older age at puberty during 4 periods; up to age 14 years, 14 to 16 years, 16 to 18 years, and 18 to 25 years. The estimated trajectories from these models were plotted for individuals in the 10th, 50th, and 90th sex-specific APHV percentiles. We assessed whether adolescent differences in BMD and BMC according to pubertal age persisted into adulthood by regressing the estimated bone values at age 25 years on APHV. We fitted initial unadjusted linear spline models followed by models that included adjustment for birth weight, ethnicity, maternal education, body mass index, and diet; unadjusted and adjusted estimates were similar, and only the adjusted estimates are presented. Participants were excluded from the linear spline analyses if they did not have a valid DXA scan from at least 1 age or complete data on pubertal age and confounders. In sensitivity analyses, the linear spline models were refitted after (1) reestimating APHV from measured heights only (ie, excluding any parental- or self-report of height) and (2) using age at menarche (reported prospectively by the mother in years and months) instead of APHV as a marker of pubertal age in female participants. All analyses were performed in R statistical software version 3.5.1 (R Project for Statistical Computing). We fitted SITAR models with the sitar package and linear spline models with the nlme package. A more detailed description of the statistical methods is provided in eMethods, eTable 1, and eTable 2 in the Supplement. Results Participant Characteristics A total of 6389 participants (3196 [50.0%] female) had bone measures from at least 1 age (26 202 BMD or BMC measures in total; median [interquartile range] measurements per individual, 4 [2-6]), and data on APHV, birth weight, ethnicity, maternal education, and childhood body mass index and diet (Figure 1). Of these, 5477 participants (2975 female [54.3%]) had bone measurements from at least 1 age after peak height velocity (Figure 1). The numbers of participants with data from each age are presented in Figure 1 and in a footnote to the Table. Participants with missing data on covariates (28% of those potentially eligible) were excluded from the analyses. The APHV occurred around 2 years earlier in female participants (mean [SD], 11.6 [0.8] years) than male participants (mean [SD], 13.5 [0.9] years) (Table and Figure 2). The BMD and BMC increased over follow-up; their values were similar between male and female participants at assessment ages 10, 12, and 14 years and were higher in male participants at ages 16, 18, and 25 years (Table and eFigure 1 in the Supplement). BMD and BMC Accrual From Childhood (Prepuberty) to Adulthood Figure 3A and B presents the rates of whole-body BMD and BMC accrual relative to pubertal age (APHV) during each growth period from childhood (prepuberty) to adulthood. Both male and female participants had gains in BMD during all 4 growth periods (ie, from 8 years before and up to 16 years after age of peak height growth) and gains in BMC during the first 3 periods (ie, up to 4 years after APHV). Male participants gained BMD and BMC at faster rates than did female participants from childhood and up to 4 years after APHV. The fastest gains in BMD and BMC in both male and female participants were between the year before APHV and up to 2 years after APHV (0.139 g/cm /y [95% 2 2 2 CI, 0.127-0.151 g/cm /y] for male participants vs 0.106 g/cm /y [95% CI, 0.098-0.114 g/cm /y] for female participants). No further gains in BMC in either male or female participants were observed from 4 years after APHV, with some evidence of modest loss in BMC in male participants during this period (eTable 3 in the Supplement). Figure 4A and B shows the mean estimated trajectories from childhood to adulthood for male and female participants, which illustrate these rapid pubertal gains in BMD and BMC. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 4/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Age at Puberty and Subsequent BMD and BMC Accrual Figure 3C and D shows the rates of whole-body BMD and BMC accrual per 1-year older age at puberty (ie, APHV) during each growth period from the time of peak height velocity up to adulthood. Older age at puberty was associated with faster gains in BMD and BMC in both male and female participants. Gains in BMD continued at a similar pace up to age 18 years in male participants (up to 2 2 age 14 years, 0.012 g/cm /y [95% CI, 0.003-0.020 g/cm /y]; between ages 14 and 16 years, 0.013 2 2 2 g/cm /y [95% CI, 0.011-0.015 g/cm /y]; between ages 16 and 18 years, 0.011 g/cm /y [95% CI, 0.010- 0.013 g/cm /y]) but had slowed considerably between ages 16 and 18 years in female participants 2 2 (up to age 14 years, 0.010 g/cm /y [95% CI, 0.009-0.011 g/cm /y]; between ages 14 and 16 years, 2 2 2 0.014 g/cm /y [95% CI, 0.014-0.015 g/cm /y]; between ages 16 and 18 years, 0.003 g/cm /y [95% CI, 0.003-0.004 g/cm /y]). Between ages 18 and 25 years, BMD gains had slowed substantially in the 2 2 male participants (0.002 g/cm /y; 95% CI, 0.001-0.003 g/cm /y) and had ceased completely in the 2 2 female participants (0.000 g/cm /y; 95% CI, −0.001 to 0.000 g/cm /y) with older pubertal age. Gains in BMC in male participants continued at similar peak levels up to age 16 years, somewhat slowing between ages 16 and 18 years before a more substantial decrease in pace between 18 and 25 years. In female participants, the gains in BMC peaked between 14 and 16 years and then continued slowing considerably between ages 16 to 18 and 18 to 25 years. Examining associations between age at puberty and estimated bone values showed evidence of persisting differences in BMD and of catch-up in BMC by age 25 years. For example, in male and female participants, respectively, per 1-year older APHV was associated with 0.050 g/cm (95% CI, 2 2 2 −0.056 to −0.045 g/cm ) and 0.044 g/cm (95% CI, −0.046 to −0.041 g/cm ) lower BMD at age 14 2 2 2 years and with 0.047 g/cm (95% CI, −0.051 to −0.043 g/cm ) and 0.034 g/cm (95% CI, −0.036 to −0.032 g/cm ) lower BMD at age 25 years. In contrast, in the male and female participants combined, per 1-year older APHV was associated with 185.5 g (95% CI, −196.4 to −174.6 g) lower BMC at age 14 years and with 116.7 g (95% CI, 111.7-121.8 g) higher BMC at age 25 years. Figure 1. Study Flowchart for the Avon Longitudinal Study of Parents and Children (ALSPAC), 1991 to 2018 14 062 Live births recruited to the ALSPAC core sample 13 973 Alive at age 1 y 7711 Invited to 7129 Invited to 6148 Invited to 10 545 With APHV 5509 Invited to 5217 Invited to 4021 Invited to DXA scan at DXA scan at DXA scan at data DXA scan at DXA scan at DXA scan at age 10 y age 12 y age 14 y age 16 y age 18 y age 25 y 7376 Had scan 7051 Had scan 6075 Had scan 5151 Had scan 4904 Had scan 3896 Had scan 7333 With valid 7000 With valid 5982 With valid 4684 With valid 4851 With valid 3525 With valid scan scan scan scan scan scan 8939 With at least 1 valid DXA scan and data on APHV 2550 Missing birth weight, ethnicity, maternal education, and BMI and diet at age 7 y 6389 With at least 1 valid DXA scan and complete data on APHV, birth weight, ethnicity, maternal education, and BMI and diet at age 7 y 5477 With at least 1 DXA scan recorded after APHV APHV indicates age at peak height velocity; BMI, body mass index; and DXA, dual-energy x-ray absorptiometry. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 5/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Figure 4C and D shows the mean estimated BMD and BMC trajectories for individuals in the 10th, 50th, and 90th APHV percentiles. The oldest pubertal age group in both sexes (ie, the 90th percentile of APHV) gained BMD and BMC at faster rates over follow-up, which partially reduced the differences in BMD and more substantially reduced differences in BMC between groups by age 25 years. Assessing associations with estimated bone values at age 25 years showed that male and female participants in the oldest pubertal age group had 0.054 g/cm (95% CI, −0.60 to −0.49 g/cm ) lower BMD (equivalent to 0.5-SD difference in BMD) and 90.9 g (95% CI, 79.6-102.1 g) higher BMC when compared with those in the youngest pubertal age group. eFigure 2 in the Supplement shows that the trajectories in Figure 4C and D were largely unchanged when APHV was estimated from measured heights only. eFigure 3 in the Supplement shows the estimated trajectories for age at menarche percentile groups in female participants, which were very similar to the trajectories for APHV percentile groups presented in Figure 4C and D. Table. Characteristics of Participants From the Avon Longitudinal Study of Parents and Children With Relevant Data, 1991 to 2018 Mean (SD) Male Participants Female Participants Characteristic (n = 3193) (n = 3196) Age at peak height velocity, y 13.5 (0.9) 11.6 (0.8) Age at dual-energy x-ray absorptiometry scan, y 10 9.8 (0.3) 9.8 (0.3) 12 11.7 (0.2) 11.7 (0.2) 14 13.8 (0.2) 13.8 (0.2) 16 15.4 (0.3) 15.4 (0.3) 18 17.8 (0.4) 17.8 (0.4) 25 24.5 (0.8) 24.4 (0.8) Whole-body bone mineral density, g/cm , by age at assessment, y 10 0.78 (0.1) 0.77 (0.1) 12 0.85 (0.1) 0.85 (0.1) 14 0.95 (0.1) 0.96 (0.1) 16 1.05 (0.1) 1.00 (0.1) 18 1.14 (0.1) 1.04 (0.1) 25 1.31 (0.1) 1.19 (0.1) Whole-body bone mineral content, g, by age at assessment, y 10 905.0 (174.6) 880.0 (191.1) 12 1186.1 (250.6) 1239.5 (298.4) 14 1721.3 (405.1) 1722.2 (343.1) 16 2226.5 (449.7) 1921.3 (344.8) 18 2562.8 (467.3) 2043.8 (373.2) 25 3047.7 (407.0) 2341.3 (274.7) Abbreviation: BMI, body mass index (calculated as the Birth weight, g 3474.2 (574.0) 3375.8 (506.3) weight in kilograms divided by height in Race, No. (%) meters squared). White 3137 (98.3) 3136 (98.1) a Numbers of participants at each age: 10 years, 2839 Nonwhite 56 (1.8) 60 (1.9) male and 2864 female; 12 years, 2691 male and 2730 female; 14 years, 2360 male and 2440 female; 16 Maternal education, No. (%) years, 1791 male and 2000 female; 18 years, 1674 Certificate of secondary education 362 (11.3) 386 (12.1) male and 2072 female; and 25 years, 1056 male and Vocational 282 (8.8) 256 (8.0) 1685 female. Numbers of participants at each age General certificate of education, ordinary level 1160 (36.3) 1107 (34.6) after removing bone measures recorded before age General certificate of education, advanced level 862 (27.0) 891 (27.8) at peak height velocity: 10 years, 0 male and 57 female; 12 years, 55 male and 1518 female; 14 years, Degree or higher 527 (16.5) 556 (17.4) 1528 male and 2474 female; 16 years, 1741 male and BMI at age 7 y 16.1 (1.9) 16.3 (2.1) 2036 female; 18 years, 1674 male and 2116 female; Daily energy intake at age 7 y, kJ/d 7785.5 (1878.8) 7477.7 (1759.8) and 25 years, 1056 male and 1725 female. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 6/14 Height Velocity, cm/y Height Velocity, cm/y JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Age at Puberty and Site-Specific Bone Trajectories eFigures 4 to 9 in the Supplement show the mean estimated BMD and BMC trajectories for arms, legs, trunk, spine, ribs, and pelvis, respectively, in male and female participants in the 10th, 50th, and 90th APHV percentiles. Overall, these were similar to the main whole-body trajectories in that those in the oldest APHV groups started lower and exhibited some catch-up over follow-up. eFigure 10 in the Supplement shows similar mean estimated hip BMD trajectories. Male and female participants in the oldest APHV group had lower total hip and femoral neck BMD at age 14 years, and these differences were more substantially reduced in male than female participants at ages 18 and 25 years (eFigure 10 in the Supplement). However, the lack of hip measures from earlier age makes it difficult to discern whether these patterns are different from trajectories of whole-body parameters. Discussion We investigated the association between age at puberty (measured by APHV) and bone accrual from ages 10 to 25 years in male and female participants from a large prospective British birth cohort. The BMD and BMC increased over follow-up with sex differences emerging during puberty. Male participants accrued BMD and BMC at faster rates than did female participants from before puberty and up to 4 years after APHV. The fastest gains in BMD and BMC in both male and female participants were during the year before APHV to 2 years after APHV. Older pubertal age was associated with an Figure 2. Height Growth and Height Velocity Curves A Male participants 200 12 Mean –1.96 SD +1.96 SD 10 100 0 5 10 15 20 Age, y Female participants 200 12 Mean height growth curves (solid lines), mean height velocity curves (dashed lines), and mean age at peak height velocity (plus estimated curves and age at peak height velocity at 1.96 SD and −1.96 SD) (dotted 100 0 vertical lines) are shown for male (A) and female (B) 5 10 15 20 participants included in the Superimposition by Age, y Translation and Rotation models. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 7/14 Height, cm Height, cm JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual initially accelerating and subsequently decelerating faster BMD and BMC accrual over follow-up. This catch-up was greater in male participants than in female participants between ages 16 and 18 years for BMD and between ages 14 and 25 years for BMC. Despite these faster gains, older pubertal age remained associated with lower BMD throughout follow-up to age 25 years and with a BMC that was lower at younger ages but higher by age 25 years. Findings were overall similar for site-specific BMD and BMC accrual. The findings that the greatest gains in BMD and BMC occurred during the year before and 2 years after APHV and that gains in BMC ceased 4 years after APHV agree with those of a recent study of the temporal association between peak height and peak bone accrual; however, that study did not test the association between APHV and bone accrual. That study, which included 2000 participants aged 5 to 19-years who underwent annual DXA scans for up to 7 years, showed that up to 36% of BMC is acquired in the 2 years before and 2 years after APHV. Our study expands on this by showing that, in contrast to BMC, gains in BMD continued into adulthood. We also showed that male participants had faster gains in both BMD and BMC from before puberty than did female participants but that sex differences in BMD and BMC emerged during puberty, which is consistent with the suggestion that sex hormones are driving sexual dimorphism in body composition. Our study contributes to previous research by documenting a process of transient pubertal catch-up in BMD and BMC among those with older pubertal age that did not fully eliminate differences in BMD. This agrees with findings that differences in BMD and BMC by APHV group were Figure 3. Gains in Bone Mineral Density (BMD) and Bone Mineral Content (BMC) A Gains in BMD from childhood to adulthood B Gains in BMC from childhood to adulthood 0.15 700 Male participants 0.12 Female participants 0.09 0.06 0.03 0 0 –8 to –1 –1 to 2 2 to 4 4 to 16 –8 to –1 –1 to 2 2 to 4 4 to 16 Time Before or After Puberty, y Time Before or After Puberty, y C Gains in BMD from puberty to adulthood D Gains in BMC from puberty to adulthood 0.025 140 0.020 0.015 0.010 0.005 0 0 ≤14 14-16 16-18 18-25 ≤14 14-16 16-18 18-25 Age, y Age, y Gains in BMD and BMC are shown from childhood (prepuberty) to adulthood, with age denote means, and whiskers denote 95% confidence intervals. Estimates were adjusted centered at pubertal age (ie, age at peak height velocity [APHV]) (A, BMD; B, BMC) and for birth weight, ethnicity, maternal education, and childhood body mass index and from puberty to adulthood per 1-year older age at puberty (C, BMD; D, BMC). Circles dietary intake. These estimates are presented in eTable 3 in the Supplement. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 8/14 Whole-body BMD, g/cm /y Whole-body BMD, g/cm /y Whole-body BMC, g/y Whole-body BMC, g/y JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual attenuated at follow-up in a mixed sample of 230 participants aged 8 to 14 years with long-term follow-up and in 500 male participants aged 18 years who were assessed at baseline and 5 years later. Our work adds to these studies by having considerably larger numbers and repeated measures of BMD from ages 10 to 25 years and by including female and male participants. This allowed us to demonstrate transient pubertal catch-up gains in BMD among those with older pubertal age, which were greater for male than for female participants. In contrast to BMD, we found that later puberty was associated with more prolonged catch-up in BMC, leading to higher total levels in adulthood. This finding may be explained by those reaching puberty later ultimately ending up 17,32 taller and, thus, having bigger (albeit less dense) bones, but further research is required to disentangle the complex associations between pubertal growth and bone mineralization up to adulthood. Given the persisting associations between later puberty and lower BMD reported here and evidence that peak bone mineralization lags behind peak height accrual, our findings suggest that 15,19,30 adolescents who mature later may be at higher risk of fractures throughout adolescence. Because differences in BMD persisted up to age 25 years, including at the spine and hip sites, which are preferentially affected by osteoporotic fracture in later life, those with older pubertal age could also be at increased risk of osteoporosis in later life, although continued follow-up of this cohort is required to identify how associations might vary through adult life. Figure 4. Estimated Bone Mineral Density (BMD) and Bone Mineral Content (BMC) Trajectories A Estimated BMD trajectory from childhood to adulthood B Estimated BMC trajectory from childhood to adulthood 1.4 3500 1.3 1.2 1.1 1.0 0.9 Male participants 1000 Male participants 0.8 Female participants Female participants 0.7 0.6 0 –8 –6 –4 –2 0 2 4 6 8 10 12 14 16 18 –8 –6 –4 –2 0 2 4 6 8 10 12 14 16 18 Time Since Peak Height Velocity, y Time Since Peak Height Velocity, y C Estimated BMD trajectory from puberty to adulthood D Estimated BMC trajectory from puberty to adulthood Male participants Female participants Male participants Female participants 1.4 3500 1.3 APHV percentile APHV percentile th th 10 10 1.2 th th 50 50 th 2500 th 90 90 1.1 1.0 0.9 0.8 1000 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 Age, y Age, y Mean estimated BMD and BMC trajectories are shown from childhood (prepuberty) to the 10th (<12.4 years in male and <10.6 years in female participants), 50th (13.3-13.5 years adulthood (A, BMD; B, BMC) against time since puberty (ie, age at peak height velocity in male and 11.4-11.6 years in female participants), and 90th (>14.5 years in male and >12.7 [APHV], represented by vertical line), and from puberty to adulthood for individuals in years in female participants) pubertal age (APHV) percentiles (C, BMD; D, BMC). JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 9/14 2 2 Whole-body BMD, g/cm Whole-body BMD, g/cm Whole-body BMC, g Whole-body BMC, g JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Limitations To our knowledge, this is the largest study on BMD in adolescents to date and has a longer follow-up than previous studies, but some limitations should be noted. Dual-energy x-ray absorptiometry is routinely used in clinical practice to assess BMD and osteoporosis; however, measurement error and changing soft-tissue distribution might influence findings. Newer imaging approaches could be used to investigate structural changes in bone microarchitecture, including growth-associated cortical porosity. Using APHV as measure of age at puberty is a key strength because it allows the same analyses in male and female participants and is less prone to differential measurement error than pubertal staging using reports of sexual maturity, such as Tanner stages. Measurement error, particularly from the self or parental reports of height, may have influenced our results. However, because parents would be unaware of their child’s subsequent bone measurements, that is unlikely to have been systematic. Furthermore, restricting analyses to research-clinic measures of height only produced similar results. In addition, similar results were obtained when using age at menarche as an alternative marker of pubertal age in female participants. The linear spline models reduced bias resulting from missing bone data by including all 27,28 participants with at least 1 bone measure under the missing-at-random assumption (ie, that the probability of missing bone data is assumed to be associated with model covariates and not associated with unmeasured factors). Although it is not possible to fully test that the assumption of missing at random holds, the probability of missing bone data among those included in the linear spline analysis was associated with the variables included in the models (eTable 1 in the Supplement). Those excluded from the analyses because of missing all 6 DXA scans were socioeconomically different from the analytic sample (eTable 2 in the Supplement), which might limit the generalizability of our findings. Participants with missing data on covariates (28% of those potentially eligible) were excluded from the analyses, which might introduce a bias if they had systematically different bone measurements; however, this seems unlikely because participants would not have 35,36 known their bone measurements. Our study sample were mostly of white British ethnicity, which means that the findings may not be generalizable to participants of other ethnicities, with potentially different bone accrual rates. Also, our analyses were performed in 2 stages rather than 37,38 as a single joint model, which might introduce a bias ; however, given the complex nonlinear SITAR models, joint modeling is unlikely to be feasible for this study. Conclusions This large and long-running prospective follow-up cohort study showed that later puberty was associated with persistently lower BMD from ages 10 to 25 years in male and female participants, despite faster gains in BMD during puberty in those with older pubertal age. Future studies should seek more robust evidence of associations between puberty timing and bone accrual from large emerging collaborations such as the European Union’s Child Cohort Network, which can provide more repeated measures over longer follow-up in populations with differing confounding structures and with larger sample sizes that can support methods that are statistically inefficient, such as within 40,41 sibship and mendelian randomization analyses. Our findings suggest that advice on how to 42,43 increase and maintain BMD such as through physical activity should be offered to people with older pubertal age. ARTICLE INFORMATION Accepted for Publication: June 17, 2019. Published: August 9, 2019. doi:10.1001/jamanetworkopen.2019.8918 Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2019 Elhakeem A et al. JAMA Network Open. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 10/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Corresponding Author: Ahmed Elhakeem, PhD, Medical Research Council Integrative Epidemiology Unit at University of Bristol, Population Health Sciences, Bristol Medical School, Oakfield House, Oakfield Grove, Bristol BS8 2BN, United Kingdom (a.elhakeem@bristol.ac.uk). Author Affiliations: Medical Research Council Integrative Epidemiology Unit at University of Bristol, Population Health Sciences, Bristol Medical School, Bristol, United Kingdom (Elhakeem, Tilling, Tobias, Lawlor); Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom (Elhakeem, Tilling, Lawlor); Musculoskeletal Research Unit, Translational Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom (Frysz, Tobias). Author Contributions: Dr Elhakeem had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: Elhakeem, Tobias, Lawlor. Acquisition, analysis, or interpretation of data: Elhakeem, Frysz, Tilling, Lawlor. Drafting of the manuscript: Elhakeem, Lawlor. Critical revision of the manuscript for important intellectual content: Elhakeem, Frysz, Tilling, Tobias. Statistical analysis: Elhakeem, Tilling. Obtained funding: Lawlor. Administrative, technical, or material support: Frysz. Supervision: Lawlor. Conflict of Interest Disclosures: Dr Elhakeem reported grants from European Union’s Horizon 2020 research and innovation program outside the submitted work. Dr Tilling reported grants from the Economic and Social Research Council and grants from the UK Medical Research Council (MRC) during the conduct of the study. Dr Lawlor reported grants from the MRC, Wellcome, and European Union’s Horizon 2020 program during the conduct of the study, and grants from several national and international government and charitable research funders and Roche Diagnostics and Medtronic Ltd outside the submitted work. No other disclosures were reported. Funding/Support: This work was supported by the European Union’s Horizon 2020 research and innovation program under grant agreement 733206 (LifeCycle), which pays the salary of Dr Elhakeem. Drs Elhakeem and Lawlor work in a unit that receives funds from the University of Bristol and UK MRC (grant MC_UU_00011/6). Dr Lawlor is a National Institute of Health Research Senior Investigator (grant NF-SI-0611-10196). The UK MRC and Wellcome (grant r102215/2/13/2) and the University of Bristol provide core support for the Avon Longitudinal Study of Parents and Children (ALSPAC). A comprehensive list of grant funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. Disclaimer: Views expressed in this article are those of the authors and not necessarily any funder. Additional Contributions: We thank all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, and nurses. Additional Information: Details of all available data and the processes and procedures involved in accessing the ALSPAC resource can be found in the ALSPAC study website, which includes a fully searchable data dictionary and variable search tool (http://www.bristol.ac.uk/alspac/researchers/our-data/). REFERENCES 1. Hansen MA, Overgaard K, Riis BJ, Christiansen C. Role of peak bone mass and bone loss in postmenopausal osteoporosis: 12 year study. BMJ. 1991;303(6808):961-964. doi:10.1136/bmj.303.6808.961 2. Harvey N, Dennison E, Cooper C. Osteoporosis: a lifecourse approach. J Bone Miner Res. 2014;29(9):1917-1925. doi:10.1002/jbmr.2286 3. Weaver CM, Gordon CM, Janz KF, et al. The National Osteoporosis Foundation’s position statement on peak bone mass development and lifestyle factors: a systematic review and implementation recommendations. Osteoporos Int. 2016;27(4):1281-1386. doi:10.1007/s00198-015-3440-3 4. Hernandez CJ, Beaupré GS, Carter DR. A theoretical analysis of the relative influences of peak BMD, age-related bone loss and menopause on the development of osteoporosis. Osteoporos Int. 2003;14(10):843-847. doi:10. 1007/s00198-003-1454-8 5. Heaney RP, Abrams S, Dawson-Hughes B, et al. Peak bone mass. Osteoporos Int. 2000;11(12):985-1009. doi:10. 1007/s001980070020 JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 11/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 6. Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005-2025. J Bone Miner Res. 2007;22(3):465-475. doi:10. 1359/jbmr.061113 7. Hopkins RB, Burke N, Von Keyserlingk C, et al. The current economic burden of illness of osteoporosis in Canada. Osteoporos Int. 2016;27(10):3023-3032. doi:10.1007/s00198-016-3631-6 8. Rogol AD, Roemmich JN, Clark PA. Growth at puberty. J Adolesc Health. 2002;31(6)(suppl):192-200. doi:10. 1016/S1054-139X(02)00485-8 9. Patton GC, Viner R. Pubertal transitions in health. Lancet. 2007;369(9567):1130-1139. doi:10.1016/S0140-6736 (07)60366-3 10. Cousminer DL, Widén E, Palmert MR. The genetics of pubertal timing in the general population: recent advances and evidence for sex-specificity. Curr Opin Endocrinol Diabetes Obes. 2016;23(1):57-65. doi:10.1097/ MED.0000000000000213 11. Day FR, Bulik-Sullivan B, Hinds DA, et al. Shared genetic aetiology of puberty timing between sexes and with health-related outcomes. Nat Commun. 2015;6:8842. doi:10.1038/ncomms9842 12. Kemp JP, Morris JA, Medina-Gomez C, et al. Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis. Nat Genet. 2017;49(10):1468-1475. doi:10.1038/ ng.3949 13. Cole TJ, Kuh D, Johnson W, et al. Using Super-Imposition by Translation And Rotation (SITAR) to relate pubertal growth to bone health in later life: the Medical Research Council (MRC) National Survey of Health and Development. Int J Epidemiol. 2016;45(4):1125-1134. 14. Cousminer DL, Mitchell JA, Chesi A, et al. Genetically determined later puberty impacts lowered bone mineral density in childhood and adulthood. J Bone Miner Res. 2018;33(3):430-436. doi:10.1002/jbmr.3320 15. Bonjour JP, Chevalley T. Pubertal timing, bone acquisition, and risk of fracture throughout life. Endocr Rev. 2014;35(5):820-847. doi:10.1210/er.2014-1007 16. Day FR, Elks CE, Murray A, Ong KK, Perry JR. Puberty timing associated with diabetes, cardiovascular disease and also diverse health outcomes in men and women: the UK Biobank study. Sci Rep. 2015;5:11208. doi:10.1038/ srep11208 17. Cousminer DL, Berry DJ, Timpson NJ, et al; ReproGen Consortium; Early Growth Genetics (EGG) Consortium. Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity. Hum Mol Genet. 2013;22(13):2735-2747. doi:10.1093/hmg/ddt104 18. Jackowski SA, Erlandson MC, Mirwald RL, et al. Effect of maturational timing on bone mineral content accrual from childhood to adulthood: evidence from 15 years of longitudinal data. Bone. 2011;48(5):1178-1185. doi:10. 1016/j.bone.2011.02.010 19. Darelid A, Ohlsson C, Nilsson M, Kindblom JM, Mellström D, Lorentzon M. Catch up in bone acquisition in young adult men with late normal puberty. J Bone Miner Res. 2012;27(10):2198-2207. doi:10.1002/jbmr.1675 20. Fraser A, Macdonald-Wallis C, Tilling K, et al. Cohort profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int J Epidemiol. 2013;42(1):97-110. doi:10.1093/ije/dys066 21. Boyd A, Golding J, Macleod J, et al. Cohort profile: the ‘children of the 90s’—the index offspring of the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2013;42(1):111-127. doi:10.1093/ije/dys064 22. Avon Longitudinal Study of Parents and Children. Explore data and samples. http://www.bristol.ac.uk/alspac/ researchers/our-data/. Accessed May 1, 2019. 23. Lawlor DA, Wills AK, Fraser A, Sayers A, Fraser WD, Tobias JH. Association of maternal vitamin D status during pregnancy with bone-mineral content in offspring: a prospective cohort study. Lancet. 2013;381(9884): 2176-2183. doi:10.1016/S0140-6736(12)62203-X 24. Warrington NM, Kemp JP, Tilling K, Tobias JH, Evans DM. Genetic variants in adult bone mineral density and fracture risk genes are associated with the rate of bone mineral density acquisition in adolescence. Hum Mol Genet. 2015;24(14):4158-4166. doi:10.1093/hmg/ddv143 25. Emmett PM, Jones LR, Northstone K. Dietary patterns in the Avon Longitudinal Study of Parents and Children. Nutr Rev. 2015;73(suppl 3):207-230. doi:10.1093/nutrit/nuv055 26. Cole TJ, Donaldson MD, Ben-Shlomo Y. SITAR—a useful instrument for growth curve analysis. Int J Epidemiol. 2010;39(6):1558-1566. doi:10.1093/ije/dyq115 27. Howe LD, Tilling K, Matijasevich A, et al. Linear spline multilevel models for summarising childhood growth trajectories: a guide to their application using examples from five birth cohorts. Stat Methods Med Res. 2016;25 (5):1854-1874. doi:10.1177/0962280213503925 JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 12/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 28. Tilling K, Macdonald-Wallis C, Lawlor DA, Hughes RA, Howe LD. Modelling childhood growth using fractional polynomials and linear splines. Ann Nutr Metab. 2014;65(2-3):129-138. doi:10.1159/000362695 29. Frysz M, Howe LD, Tobias JH, Paternoster L. Using SITAR (SuperImposition by Translation and Rotation) to estimate age at peak height velocity in Avon Longitudinal Study of Parents and Children. Wellcome Open Res. 2018;3:90. doi:10.12688/wellcomeopenres.14708.1 30. McCormack SE, Cousminer DL, Chesi A, et al. Association between linear growth and bone accrual in a diverse cohort of children and adolescents. JAMA Pediatr. 2017;171(9):e171769. doi:10.1001/jamapediatrics.2017.1769 31. Wells JC. Sexual dimorphism of body composition. Best Pract Res Clin Endocrinol Metab. 2007;21(3):415-430. doi:10.1016/j.beem.2007.04.007 32. Carel JC, Lahlou N, Roger M, Chaussain JL. Precocious puberty and statural growth. Hum Reprod Update. 2004;10(2):135-147. doi:10.1093/humupd/dmh012 33. Kanis JA. Diagnosis of osteoporosis and assessment of fracture risk. Lancet. 2002;359(9321):1929-1936. doi: 10.1016/S0140-6736(02)08761-5 34. Farr JN, Khosla S. Skeletal changes through the lifespan—from growth to senescence. Nat Rev Endocrinol. 2015;11(9):513-521. doi:10.1038/nrendo.2015.89 35. Munafò MR, Tilling K, Taylor AE, Evans DM, Davey Smith G. Collider scope: when selection bias can substantially influence observed associations. Int J Epidemiol. 2018;47(1):226-235. doi:10.1093/ije/dyx206 36. Hughes RA, Heron J, Sterne JAC, Tilling K. Accounting for missing data in statistical analyses: multiple imputation is not always the answer. Int J Epidemiol. 2019;dyz032. doi:10.1093/ije/dyz032 37. Sayers A, Heron J, Smith A, et al. Joint modelling compared with two stage methods for analysing longitudinal data and prospective outcomes: a simulation study of childhood growth and BP. Stat Methods Med Res. 2017;26 (1):437-452. doi:10.1177/0962280214548822 38. Macdonald-Wallis C, Lawlor DA, Palmer T, Tilling K. Multivariate multilevel spline models for parallel growth processes: application to weight and mean arterial pressure in pregnancy. Stat Med. 2012;31(26):3147-3164. doi: 10.1002/sim.5385 39. LifeCycle. European Union Child Cohort Network. https://lifecycle-project.eu. Accessed May 1, 2019. 40. Lawlor DA, Tilling K, Davey Smith G. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016;45(6): 1866-1886. 41. Munafò MR, Davey Smith G. Robust research needs many lines of evidence. Nature. 2018;553(7689): 399-401. doi:10.1038/d41586-018-01023-3 42. Mitchell JA, Chesi A, McCormack SE, et al. Physical activity and bone accretion: isotemporal modeling and genetic interactions. Med Sci Sports Exerc. 2018;50(5):977-986. doi:10.1249/MSS.0000000000001520 43. Muthuri SG, Ward KA, Kuh D, Elhakeem A, Adams JE, Cooper R. Physical activity across adulthood and bone health in later life: the 1946 British birth cohort. J Bone Miner Res. 2019;34(2):252-261. SUPPLEMENT. eMethods. Supplemental Methods. eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV eReferences eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty eFigure 1. Observed Whole-body BMD and BMC at Each Follow-up Assessment From 10 to 25 Years eFigure 2. Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 13/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 14/14 Supplementary Online Content Elhakeem A, Frysz M, Tilling K, Tobias JH, Lawlor DA. Association between age at puberty and bone accrual from 10 to 25 years of age. JAMA Netw Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 eMethods. Supplemental Methods eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV eReferences eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty eFigure 1. Observed Whole-body BMD and BMC at Each Follow-up Assessment From 10 to 25 Years eFigure 2. Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for Iindivudals in the 10th, 50th, and 90th APHV Percentiles © 2019 Elhakeem A et al. JAMA Network Open. eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles This supplementary material has been provided by the authors to give readers additional information about their work. © 2019 Elhakeem A et al. JAMA Network Open. eMethods. Supplemental Methods. APHV was estimated using Superimposition by Translation and Rotation (SITAR) growth curve analysis (1, 2). SITAR is a shape invariant model consisting of a natural cubic spline mean curve of height vs age. SITAR fits this curve with three person-specific random-effects that reflect individual differences from the mean curve in size (random height intercept for person-specific shift up or down in the spline curve), velocity (random age scaling that reflects differences in the duration of the growth spurt in individuals) and tempo (random age intercept reflecting individual differences in the timing of the growth spurt). SITAR models were fitted separately for males and females with 5 degrees of freedom. APHV estimation was restricted to individuals with height measurements between ages 5 and 20 years and at least one height measure after age 9 years. Data cleaning similar to previous ALSPAC analyses was performed to remove data points with velocity exceeding 4 SDs and standardized residuals exceeding 4 in absolute value (2, 3, 4), leaving a total of 5223 males (with 55 752 height measurements) and 5322 females (with 59 318 height measurements) with an estimate of APHV. The final SITAR models explained 95.7% and 95.6% of variance in height in males and females respectively. APHV in years was estimated by transforming the random age intercept that reflects individual differences in the tempo (timing) of the growth spurt. We subsequently investigated the relationship between APHV and whole-body BMD and BMC accrual using two sets of mixed-effects linear spline regression analyses. Linear spline models facilitate the summarising of nonlinear bone accrual through a series of linear splines joined at knot points (5, 6). They estimate mean and person-specific bone values at the first assessment and mean and person-specific rates of accrual during each of the different growth periods. Smoothing functions, prior subject knowledge and fit statistics were used to identify the best fitting models (5, 6) and subsequently knot point selection was based on the mean age at DXA scans. Models were fitted using maximum likelihood and included random intercepts and random age slopes. A st 1 order autoregressive correlation structure was used to account for temporal autocorrelation, and a fixed variance structure was used to model heteroskedasticity. In both sets of linear spline models, we examined the relationship between APHV and rates of bone accrual by including APHV-by-age splines interaction terms and included sex-by-age splines interactions terms to test for statistical evidence of different bone trajectories for males and females. All models were adjusted for birth weight, ethnicity, maternal education and childhood BMI and energy intake. the observed and predicted BMD and BMC values (without random effects) were 0.84 and 0.83. © 2019 Elhakeem A et al. JAMA Network Open. In the first analyses, age was centred at APHV in order to estimate rates of bone accrual relative to APHV; i.e. from before puberty up to adulthood. Knot points were placed at the mean ages of the 12, 14, and 16 -year scans which equated to summarising rate of accrual from 7.7 up to 0.9 years prior to APHV, 0.9 year prior up to 2.1 years following APHV, 2.1 to 3.7 years following APHV, and 3.7 to 16 years after APHV. The mean predicted trajectories relative to APHV were plotted for males and females. For the second analyses, we excluded individuals with accrual from before APHV in order to investigate associations between APHV and subsequent rate of bone accrual up to adulthood. Here, knot points were placed at the mean ages of the 14 (centred at age 14), 16, and 18-year assessments, which then estimated associations with rate of accrual up to age 14, between 14 16, 16 18 and 18 25 years. To assess if difference in bone outcomes persisted into adulthood, predicted bone values at age 25 were regressed on APHV (after adjusting for predicted values at younger ages) and compared to associations with predicted values at age 14. To visualise associations between APHV and bone th th th accrual, predicted mean trajectories were plotted for individuals in the 10 , 50 , and 90 sex-specific APHV percentiles. In sensitivity analyses, we refitted models after re-estimating APHV from measured heights only (3) and using age at menarche (reported prospectively by the mother in years and months) instead of APHV as an alternative marker of age at puberty in females. For our secondary outcomes, we investigated associations between APHV and subsequent site-specific BMD and BMC accrual from whole-body DXA, and with total hip and femoral neck BMD accrual from hip DXA. Given that hips were only scanned at three different ages (14, 18 and 25), we fitted hip models using a single knot point placed at the mean ages of the 18-year assessment. The predicted mean trajectories were plotted for all secondary outcomes. All analyses were performed in R (R Foundation for Statistical Computing, Vienna); SITAR models were fitted wit The linear spline models reduced bias due to missing bone data by including all measures under the missing at random assumption (i.e. the probability of missing bone data assumed to be related to measured factors but unrelated to unmeasured factors). eTable 1 shows differences between those included in the analysis with complete and partial DXA data, and suggests that the probability of missing at least one bone measure was related to the model covariates, and thus possibly missing at random. Those missing all DXA measures were excluded from the analysis. eTable 2 shows differences between those included in the linear spline models and excluded due to missing all DXA data, and shows that compared with the participants in this study, those © 2019 Elhakeem A et al. JAMA Network Open. excluded due to missing bone data had lower socioeconomic position and had slightly lower birthweight and slightly higher BMI (eTable 2). Those missing data on model covariates were also excluded from the analysis. © 2019 Elhakeem A et al. JAMA Network Open. eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure Males Females Included in Included in Included in Included in analyses: all six analyses: missing analyses: all six analyses: missing DXA measures DXA measures available (n=585) (n=2608) available (n=993) (n=2203) Birth weight g [mean (SD)] 3484.1 (548.2) 3472.0 (597.8) 3375.3 (499.7) 3376.1 (509.3) Ethnicity [No. (%)] White 575 (98.3) 2562 (98.2) 972 (97.9) 2164 (98.2) Non-white 10 (1.7) 46 (1.8) 21 (2.1) 39 (1.8) Maternal education [No. (%)] Certificate of Secondary Education 37 (6.3) 325 (12.5) 58 (5.8) 328 (14.9) Vocational 34 (5.8) 248 (9.5) 54 (5.4) 202 (9.2) General Certificate of Education: Ordinary Level 175 (29.9) 985 (37.8) 343 (34.5) 764 (34.7) General Certificate of Education: Advanced Level 179 (30.6) 683 (26.2) 310 (31.2) 581 (26.4) Degree or higher 160 (27.4) 367 (14.1) 228 (23.0) 328 (14.9) BMI at age 7y kg/m [mean (SD)] 15.8 (1.7) 16.1 (1.9) 16.0 (1.8) 16.4 (2.2) Daily energy intake at age 7y kJ/day [mean (SD)] 7712.9 (1603.2) 7801.7 (1935.2) 7431.4 (1612.6) 7498.5 (1821.7) APHV years [mean (SD)] 13.4 (0.9) 13.5 (0.9) 11.7 (0.8) 11.6 (0.8) © 2019 Elhakeem A et al. JAMA Network Open. eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV Males Females Included in linear Excluded due to Included in linear Excluded due to spline models missing all DXA spline models missing all (n=3193) (n=851) (n=3196) DXA (n=755) Birth weight g [mean (SD)] 3474.2 (574.0) 3449.4 (578.1) 3375.8 (506.3) 3372.7 (506.3) Ethnicity [No. (%)] White 3137 (98.3) 727 (96.7) 3136 (98.1) 621 (97.8) Non-white 56 (1.8) 25 (3.3) 60 (1.9) 14 (2.2) Maternal education [No. (%)] Certificate of Secondary Education 362 (11.3) 161 (21.3) 386 (12.1) 172 (26.8) Vocational 282 (8.8) 74 (9.8) 256 (8.0) 75 (11.7) General Certificate of Education: Ordinary Level 1160 (36.3) 266 (35.2) 1107 (34.6) 217 (33.8) General Certificate of Education: Advanced Level 862 (27.0) 157 (20.8) 891 (27.8) 113 (17.6) Degree or higher 527 (16.5) 97 (12.9) 556 (17.4) 66 (10.3) BMI at age 7y kg/m [mean (SD)] 16.1 (1.9) 16.2 (1.9) 16.3 (2.1) 16.8 (2.9) Daily energy intake at age 7y kJ/day [mean (SD)] 7785.5 (1878.8) 7731.9 (2068.8) 7477.7 (1759.8) 7459.9 (2029.1) APHV years [mean (SD)] 13.5 (0.9) 13.5 (0.7) 11.6 (0.8) 11.7 (0.8) © 2019 Elhakeem A et al. JAMA Network Open. eReferences 1. Cole TJ, Donaldson MD, Ben-Shlomo Y. SITAR--a useful instrument for growth curve analysis. Int J Epidemiol. 2010;39(6):1558-66. 2. Cole TJ, Kuh D, Johnson W, Ward KA, Howe LD, Adams JE, et al. Using Super-Imposition by Translation And Rotation (SITAR) to relate pubertal growth to bone health in later life: the Medical Research Council (MRC) National Survey of Health and Development. Int J Epidemiol. 2016. 3. Frysz M, Howe L, Tobias J, Paternoster L. Using SITAR (SuperImposition by Translation and Rotation) to estimate age at peak height velocity in Avon Longitudinal Study of Parents and Children. Wellcome Open Res; 2018. 4. Mahmoud O, Granell R, Tilling K, Minelli C, Garcia-Aymerich J. Association of Height Growth in Puberty with Lung Function: A Longitudinal Study. Am J Respir Crit Care Med. 2018;doi: 10.1164/rccm.201802 - 0274OC 5. Tilling K, Macdonald-Wallis C, Lawlor DA, Hughes RA, Howe LD. Modelling childhood growth using fractional polynomials and linear splines. Annals of nutrition & metabolism. 2014;65(2-3):129-38. 6. Howe LD, Tilling K, Matijasevich A, Petherick ES, Santos AC, Fairley L, et al. Linear spline multilevel models for summarising childhood growth trajectories: A guide to their application using examples from five birth cohorts. Statistical methods in medical research. 2016;25(5):1854-74. © 2019 Elhakeem A et al. JAMA Network Open. eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty BMD g/cm /yr. (95% CI) BMC g/yr. (95% CI) Males Females Males Females A. Rates relative to APHV (n= 6389)* -7.7 to -0.9 yrs. 0.054 (0.047 to 0.062) 0.014 (0.000 to 0.028) 289.109 (256.887 to 321.331) 74.065 (5.229 to 139.537) -0.9 to 2.1 yrs 0.137 (0.125 to 0.150) 0.106 (0.098 to 0.115) 644.480 (592.822 to 696.138) 428.786 (390.064 to 463.701) 2.1 to 3.7 yrs. 0.078 (0.060 to 0.096) 0.047 (0.037 to 0.058) 456.124 (378.401 to 533.847) 312.841 (263.411 to 306.351) 3.7 to 16.0 yrs. 0.013 (0.005 to 0.021) 0.005 (0.001 to 0.009) -35.423 (-69.913 to -0.934) 16.057 (-2.495 to 34.609) B. Rates per 1-year older APHV (n= 5477)** up to 14 yrs 0.012 (0.003 to 0.020) 0.010 (0.009 to 0.011) 92.259 (49.562 to 135.955) 40.902 (35.538 to 46.265) 14 to 16 yrs 0.013 (0.011 to 0.015) 0.014 (0.014 to 0.015) 85.125 (76.005 to 94.244) 61.378 (56.995 to 65.761) 16 to 18 yrs 0.011 (0.010 to 0.013) 0.003 (0.003 to 0.004) 64.252 (58.846 to 69.658) 19.374 (15.398 to 23.350) 18 to 25 yrs 0.002 (0.001 to 0.003) 0.000 (-0.001 to 0.000) 13.368 (10.823 to 15.913) 5.793 (3.902 to 7.683) * Knots at the 12, 14, 16 year assessments. Age centred at APHV (i.e. 0=APHV): models includes APHV:age splines interaction terms, and adjsutment for birth weight, maternal education, childhood BMI and diet. ** Knots at the 14, 16 and 18 year assessments. Models includes APHV:age splines interaction terms, and adjsutment for birth weight, ethnicity, maternal education and childhood BMI and diet. These models excluded bone measures recorded before APHV. © 2019 Elhakeem A et al. JAMA Network Open. eFigure 1 Observed Whole-body BMD and BMC at Each Follow-up Assessment from 10 to 25 Years © 2019 Elhakeem A et al. JAMA Network Open. eFigure 2 Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for IIndividuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Individuals in the 10 th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png JAMA Network Open American Medical Association

Association Between Age at Puberty and Bone Accrual From 10 to 25 Years of Age

Loading next page...
 
/lp/american-medical-association/association-between-age-at-puberty-and-bone-accrual-from-10-to-25-jOBbnNjUWz

References (67)

Publisher
American Medical Association
Copyright
Copyright 2019 Elhakeem A et al. JAMA Network Open.
eISSN
2574-3805
DOI
10.1001/jamanetworkopen.2019.8918
Publisher site
See Article on Publisher Site

Abstract

Key Points Question Is puberty timing associated IMPORTANCE Bone health in early life is thought to influence the risk of osteoporosis in later life. with growth-related bone accrual up to adulthood? OBJECTIVE To examine whether puberty timing is associated with bone mineral density accrual up Findings In this cohort study of 6389 to adulthood. participants who underwent repeated bone density scans from ages 10 to 25 DESIGN, SETTING, AND PARTICIPANTS This cohort study used data from the Avon Longitudinal years, later puberty was associated with Study of Parents and Children, a prospective population-based birth cohort initiated in 1991 to 1992 persistently lower bone mineral density, in southwest England. The participants were 6389 healthy British people who underwent regular despite some catch-up during puberty. follow-up, including up to 6 repeated bone density scans from ages 10 to 25 years. Data analysis was performed from June 2018 to June 2019. Meaning People with older pubertal age should be advised on how to EXPOSURES Age at puberty from estimated age at peak height velocity (years). maximize bone density and minimize its decrease in later life to help prevent MAIN OUTCOMES AND MEASURES Gains per year in whole-body bone mineral density (grams per fracture and osteoporosis. square centimeter), assessed by dual-energy x-ray absorptiometry at ages 10, 12, 14, 16, 18, and 25 years and modeled using linear splines. Supplemental content RESULTS A total of 6389 participants (3196 [50.0%] female) were included. The mean (SD) age at Author affiliations and article information are listed at the end of this article. peak height velocity was 13.5 (0.9) years for male participants and 11.6 (0.8) years for female participants. Male participants gained bone mineral density at faster rates than did female participants, with the greatest gains in both male participants (0.139 g/cm /y; 95% CI, 0.127-0.151 2 2 2 g/cm /y) and female participants (0.106 g/cm /y; 95% CI, 0.098-0.114 g/cm /y) observed between the year before and 2 years after peak height velocity. When aligned by chronological age, per 1-year older age at puberty was associated with faster subsequent gains in bone mineral density; the magnitudes of faster gains were greatest between ages 14 and 16 years in both male participants 2 2 2 (0.013 g/cm /y; 95% CI, 0.011-0.015 g/cm /y) and female participants (0.014 g/cm /y; 95% CI, 2 2 0.014-0.015 g/cm /y), were greater in male participants (0.011 g/cm /y; 95% CI, 0.010-0.013 2 2 2 g/cm /y) than in female participants (0.003 g/cm /y; 95% CI, 0.003-0.004 g/cm /y) between ages 16 and 18 years, and were least in both male participants (0.002 g/cm /y; 95% CI, 0.001-0.003 2 2 2 g/cm /y) and female participants (0.000 g/cm /y; 95% CI, −0.001 to 0.000 g/cm /y) between ages 18 and 25 years. Despite faster gains, older age at puberty was associated with persistently lower 2 2 bone mineral density, changing from 0.050 g/cm (95% CI, −0.056 to −0.045 g/cm )loweratage14 2 2 years to 0.047 g/cm (95% CI, −0.051 to −0.043 g/cm ) lower at age 25 years in male participants 2 2 2 and from 0.044 g/cm (95% CI, −0.046 to −0.041 g/cm ) to 0.034 g/cm (95% CI, −0.036 to −0.032 g/cm ) lower at the same ages in female participants. (continued) Open Access. This is an open access article distributed under the terms of the CC-BY License. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 1/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Abstract (continued) CONCLUSIONS AND RELEVANCE People with older pubertal age should be advised on how to maximize bone mineral density and minimize its decrease in later life to help prevent fracture and osteoporosis. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 Introduction Peak bone mass at the end of growth is thought to be an important determinant of later-life risk of 1-5 6,7 fracture and osteoporosis, a bone loss disorder with substantial and increasing health costs. For example, bone remodeling simulations showed that a 10% increase in peak bone mineral density (BMD) would delay osteoporosis by 13 years. Puberty is a key early life milestone that is characterized by endocrine-initiated reproductive maturation and a dramatic skeletal growth spurt 8,9 10-12 in height. Although bone mass potential and puberty timing are both strongly heritable, 13-16 evidence indicates that later puberty may lead to lower BMD in adolescence and adulthood and, thus, an increased risk of osteoporosis later in life. However, the association between puberty timing and long-term bone accrual from early life up to adulthood, including the extent and duration of any catch-up bone accrual by later maturing adolescents, has not been described, to our knowledge. Age at peak height velocity (APHV) is an accurate and precise marker of puberty timing that allows direct 13,17 comparisons of this association between male and female individuals. Yet, to our knowledge, only 2 studies, which included 230 participants who underwent repeated bone scans and 500 male participants assessed once at baseline and once at follow-up, have examined the association between APHV and bone accrual. Therefore, the aim of this study was to examine the association between age at puberty (measured by APHV) and bone accrual from ages 10 to 25 years in a large birth cohort of male and female participants. We assessed the rates of BMD and bone mineral content (BMC) accrual from before puberty up to adulthood relative to pubertal age and examined the association between age at puberty and subsequent BMD and BMC accrual. Our secondary aim was to examine the associations between age at puberty and site-specific bone accrual. Methods Study Sample 20,21 The Avon Longitudinal Study of Parents and Children (ALSPAC) is a prospective birth cohort study that recruited all pregnant women residing within the catchment area of 3 National Health Service authorities in southwest England with an expected date of delivery between April 1991 and December 1992. In total, 15 247 eligible pregnancies were enrolled in ALSPAC (75% response), resulting in 14 973 live births, of whom 14 899 were alive at 1 year of age. Detailed information has been collected from offspring and parents using questionnaires, data extraction from medical records, linkage to health records, and dedicated clinic assessments up to the last completed contact in 2018. Details of all available data can be found in the ALSPAC study website, which includes a fully searchable data dictionary and variable search tool. Ethics approval was obtained from the ALSPAC law and ethics committee and the local National Health Service research ethics committee. Written informed consent was obtained from all participants. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline. Height Measurements Numerous height measures have been obtained from various sources from birth to age 25 years, including from routine data collected from midwives, health visitors, linkage to child health records, JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 2/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual and ALSPAC research clinic visits. Height from research clinic visits (annually or more frequently up to age 14 years and then at mean ages of 16, 18, and 25 years) was measured by accredited fieldworkers to the nearest 0.1 cm using a Harpenden stadiometer (Holtain Ltd). Measured heights were supplemented by extensive maternal and self-reported height records collected throughout the study. For these analyses, we restricted APHV estimation to individuals with height measurements taken between ages 5 and 20 years and at least one height measurement after age 9 years because these were relevant to assessment of APHV. Dual-Energy X-ray Absorptiometry Scans All participants were invited to undergo whole-body dual-energy x-ray absorptiometry (DXA) scans as part of clinic assessments at mean ages 10, 12, 14, 16, 18, and 25 years. Scans were performed using a Lunar Prodigy scanner (Lunar Radiation Corp) and analyzed according to the manufacturer’s standard scanning software and positioning protocols. Scans were reanalyzed as necessary to ensure optimal placement of borders between adjacent subregions, and scans with anomalies were 23,24 excluded. For our primary outcomes, we extracted whole-body (except for the head) BMD (grams per square centimeters) and BMC (grams) at each age. Our secondary outcomes were site- specific BMD and BMC (arms, legs, trunk, ribs, spine, and pelvis) from whole-body DXA, in addition to total hip and femoral neck BMD from up to 3 repeated hip DXA scans performed at ages 14, 18, and 25 years. Confounding Variables Birth weight, ethnicity, socioeconomic position, body mass index (calculated as the weight in kilograms divided by height in meters squared), and diet in early life were hypothesized to potentially confound associations between puberty timing and bone outcomes and were included as model adjustments. Birth weights were recorded to the nearest gram and were extracted from hospital records. Ethnicity was based on the mother’s ethnic background and was reported during a general prenatal questionnaire. Early life socioeconomic position was based on the mother’s highest educational qualifications reported in pregnancy. Body mass index was derived from height and weight measurements taken at the age 7 years research clinic visit. Dietary intake in early life was based on the child’s daily energy intake (kilojoules per day), which was derived from food frequency questionnaires completed by the parents when the child was aged 7 years. Statistical Analysis Data analysis was performed from June 2018, to June 2019. Statistical analyses were performed in 2 stages. The first stage involved measuring age at puberty by estimating APHV, which was done using 13,26 Superimposition by Translation and Rotation (SITAR) growth curve analysis. The SITAR models reduce complex growth data into clinically relevant parameters that represent the timing, intensity, and duration of the pubertal growth spurt, which, in turn, simplifies comparison between individuals and between male and female participants. We estimated APHV through transformation of the random age intercept that reflects individual differences in the timing of the growth spurt. The second stage involved investigating the association between age at puberty and whole- body and site-specific BMD and BMC accrual up to adulthood using 2 sets of mixed-effects linear 27,28 spline regression models. These models summarize nonlinear change through a series of linear splines joined at knot points by estimating person-specific rates of accrual during each growth 27,28 period. We included APHV-by-age splines interaction terms to investigate whether bone accrual differs by age at puberty, and models with sex-by-age splines interaction terms were used to test differences in bone accrual between male and female participants. In the first set of these analyses, we examined gains in BMD and BMC from childhood (prepuberty) to adulthood relative to pubertal age by centering chronological age at APHV. Results from these models can be interpreted as rates of BMD and BMC accrual during 4 periods: 8 years up to 1 year before APHV, 1 year before APHV to 2 years after APHV, 2 to 4 years after APHV, and 4 to JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 3/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 16 years after APHV. The estimated trajectories were plotted for male and female participants. In the second set of analyses, we excluded bone scans taken before APHV to investigate the associations between age at puberty and subsequent bone accrual up to adulthood. Results from these models can be interpreted as rate of BMD and BMC accrual per older age at puberty during 4 periods; up to age 14 years, 14 to 16 years, 16 to 18 years, and 18 to 25 years. The estimated trajectories from these models were plotted for individuals in the 10th, 50th, and 90th sex-specific APHV percentiles. We assessed whether adolescent differences in BMD and BMC according to pubertal age persisted into adulthood by regressing the estimated bone values at age 25 years on APHV. We fitted initial unadjusted linear spline models followed by models that included adjustment for birth weight, ethnicity, maternal education, body mass index, and diet; unadjusted and adjusted estimates were similar, and only the adjusted estimates are presented. Participants were excluded from the linear spline analyses if they did not have a valid DXA scan from at least 1 age or complete data on pubertal age and confounders. In sensitivity analyses, the linear spline models were refitted after (1) reestimating APHV from measured heights only (ie, excluding any parental- or self-report of height) and (2) using age at menarche (reported prospectively by the mother in years and months) instead of APHV as a marker of pubertal age in female participants. All analyses were performed in R statistical software version 3.5.1 (R Project for Statistical Computing). We fitted SITAR models with the sitar package and linear spline models with the nlme package. A more detailed description of the statistical methods is provided in eMethods, eTable 1, and eTable 2 in the Supplement. Results Participant Characteristics A total of 6389 participants (3196 [50.0%] female) had bone measures from at least 1 age (26 202 BMD or BMC measures in total; median [interquartile range] measurements per individual, 4 [2-6]), and data on APHV, birth weight, ethnicity, maternal education, and childhood body mass index and diet (Figure 1). Of these, 5477 participants (2975 female [54.3%]) had bone measurements from at least 1 age after peak height velocity (Figure 1). The numbers of participants with data from each age are presented in Figure 1 and in a footnote to the Table. Participants with missing data on covariates (28% of those potentially eligible) were excluded from the analyses. The APHV occurred around 2 years earlier in female participants (mean [SD], 11.6 [0.8] years) than male participants (mean [SD], 13.5 [0.9] years) (Table and Figure 2). The BMD and BMC increased over follow-up; their values were similar between male and female participants at assessment ages 10, 12, and 14 years and were higher in male participants at ages 16, 18, and 25 years (Table and eFigure 1 in the Supplement). BMD and BMC Accrual From Childhood (Prepuberty) to Adulthood Figure 3A and B presents the rates of whole-body BMD and BMC accrual relative to pubertal age (APHV) during each growth period from childhood (prepuberty) to adulthood. Both male and female participants had gains in BMD during all 4 growth periods (ie, from 8 years before and up to 16 years after age of peak height growth) and gains in BMC during the first 3 periods (ie, up to 4 years after APHV). Male participants gained BMD and BMC at faster rates than did female participants from childhood and up to 4 years after APHV. The fastest gains in BMD and BMC in both male and female participants were between the year before APHV and up to 2 years after APHV (0.139 g/cm /y [95% 2 2 2 CI, 0.127-0.151 g/cm /y] for male participants vs 0.106 g/cm /y [95% CI, 0.098-0.114 g/cm /y] for female participants). No further gains in BMC in either male or female participants were observed from 4 years after APHV, with some evidence of modest loss in BMC in male participants during this period (eTable 3 in the Supplement). Figure 4A and B shows the mean estimated trajectories from childhood to adulthood for male and female participants, which illustrate these rapid pubertal gains in BMD and BMC. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 4/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Age at Puberty and Subsequent BMD and BMC Accrual Figure 3C and D shows the rates of whole-body BMD and BMC accrual per 1-year older age at puberty (ie, APHV) during each growth period from the time of peak height velocity up to adulthood. Older age at puberty was associated with faster gains in BMD and BMC in both male and female participants. Gains in BMD continued at a similar pace up to age 18 years in male participants (up to 2 2 age 14 years, 0.012 g/cm /y [95% CI, 0.003-0.020 g/cm /y]; between ages 14 and 16 years, 0.013 2 2 2 g/cm /y [95% CI, 0.011-0.015 g/cm /y]; between ages 16 and 18 years, 0.011 g/cm /y [95% CI, 0.010- 0.013 g/cm /y]) but had slowed considerably between ages 16 and 18 years in female participants 2 2 (up to age 14 years, 0.010 g/cm /y [95% CI, 0.009-0.011 g/cm /y]; between ages 14 and 16 years, 2 2 2 0.014 g/cm /y [95% CI, 0.014-0.015 g/cm /y]; between ages 16 and 18 years, 0.003 g/cm /y [95% CI, 0.003-0.004 g/cm /y]). Between ages 18 and 25 years, BMD gains had slowed substantially in the 2 2 male participants (0.002 g/cm /y; 95% CI, 0.001-0.003 g/cm /y) and had ceased completely in the 2 2 female participants (0.000 g/cm /y; 95% CI, −0.001 to 0.000 g/cm /y) with older pubertal age. Gains in BMC in male participants continued at similar peak levels up to age 16 years, somewhat slowing between ages 16 and 18 years before a more substantial decrease in pace between 18 and 25 years. In female participants, the gains in BMC peaked between 14 and 16 years and then continued slowing considerably between ages 16 to 18 and 18 to 25 years. Examining associations between age at puberty and estimated bone values showed evidence of persisting differences in BMD and of catch-up in BMC by age 25 years. For example, in male and female participants, respectively, per 1-year older APHV was associated with 0.050 g/cm (95% CI, 2 2 2 −0.056 to −0.045 g/cm ) and 0.044 g/cm (95% CI, −0.046 to −0.041 g/cm ) lower BMD at age 14 2 2 2 years and with 0.047 g/cm (95% CI, −0.051 to −0.043 g/cm ) and 0.034 g/cm (95% CI, −0.036 to −0.032 g/cm ) lower BMD at age 25 years. In contrast, in the male and female participants combined, per 1-year older APHV was associated with 185.5 g (95% CI, −196.4 to −174.6 g) lower BMC at age 14 years and with 116.7 g (95% CI, 111.7-121.8 g) higher BMC at age 25 years. Figure 1. Study Flowchart for the Avon Longitudinal Study of Parents and Children (ALSPAC), 1991 to 2018 14 062 Live births recruited to the ALSPAC core sample 13 973 Alive at age 1 y 7711 Invited to 7129 Invited to 6148 Invited to 10 545 With APHV 5509 Invited to 5217 Invited to 4021 Invited to DXA scan at DXA scan at DXA scan at data DXA scan at DXA scan at DXA scan at age 10 y age 12 y age 14 y age 16 y age 18 y age 25 y 7376 Had scan 7051 Had scan 6075 Had scan 5151 Had scan 4904 Had scan 3896 Had scan 7333 With valid 7000 With valid 5982 With valid 4684 With valid 4851 With valid 3525 With valid scan scan scan scan scan scan 8939 With at least 1 valid DXA scan and data on APHV 2550 Missing birth weight, ethnicity, maternal education, and BMI and diet at age 7 y 6389 With at least 1 valid DXA scan and complete data on APHV, birth weight, ethnicity, maternal education, and BMI and diet at age 7 y 5477 With at least 1 DXA scan recorded after APHV APHV indicates age at peak height velocity; BMI, body mass index; and DXA, dual-energy x-ray absorptiometry. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 5/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Figure 4C and D shows the mean estimated BMD and BMC trajectories for individuals in the 10th, 50th, and 90th APHV percentiles. The oldest pubertal age group in both sexes (ie, the 90th percentile of APHV) gained BMD and BMC at faster rates over follow-up, which partially reduced the differences in BMD and more substantially reduced differences in BMC between groups by age 25 years. Assessing associations with estimated bone values at age 25 years showed that male and female participants in the oldest pubertal age group had 0.054 g/cm (95% CI, −0.60 to −0.49 g/cm ) lower BMD (equivalent to 0.5-SD difference in BMD) and 90.9 g (95% CI, 79.6-102.1 g) higher BMC when compared with those in the youngest pubertal age group. eFigure 2 in the Supplement shows that the trajectories in Figure 4C and D were largely unchanged when APHV was estimated from measured heights only. eFigure 3 in the Supplement shows the estimated trajectories for age at menarche percentile groups in female participants, which were very similar to the trajectories for APHV percentile groups presented in Figure 4C and D. Table. Characteristics of Participants From the Avon Longitudinal Study of Parents and Children With Relevant Data, 1991 to 2018 Mean (SD) Male Participants Female Participants Characteristic (n = 3193) (n = 3196) Age at peak height velocity, y 13.5 (0.9) 11.6 (0.8) Age at dual-energy x-ray absorptiometry scan, y 10 9.8 (0.3) 9.8 (0.3) 12 11.7 (0.2) 11.7 (0.2) 14 13.8 (0.2) 13.8 (0.2) 16 15.4 (0.3) 15.4 (0.3) 18 17.8 (0.4) 17.8 (0.4) 25 24.5 (0.8) 24.4 (0.8) Whole-body bone mineral density, g/cm , by age at assessment, y 10 0.78 (0.1) 0.77 (0.1) 12 0.85 (0.1) 0.85 (0.1) 14 0.95 (0.1) 0.96 (0.1) 16 1.05 (0.1) 1.00 (0.1) 18 1.14 (0.1) 1.04 (0.1) 25 1.31 (0.1) 1.19 (0.1) Whole-body bone mineral content, g, by age at assessment, y 10 905.0 (174.6) 880.0 (191.1) 12 1186.1 (250.6) 1239.5 (298.4) 14 1721.3 (405.1) 1722.2 (343.1) 16 2226.5 (449.7) 1921.3 (344.8) 18 2562.8 (467.3) 2043.8 (373.2) 25 3047.7 (407.0) 2341.3 (274.7) Abbreviation: BMI, body mass index (calculated as the Birth weight, g 3474.2 (574.0) 3375.8 (506.3) weight in kilograms divided by height in Race, No. (%) meters squared). White 3137 (98.3) 3136 (98.1) a Numbers of participants at each age: 10 years, 2839 Nonwhite 56 (1.8) 60 (1.9) male and 2864 female; 12 years, 2691 male and 2730 female; 14 years, 2360 male and 2440 female; 16 Maternal education, No. (%) years, 1791 male and 2000 female; 18 years, 1674 Certificate of secondary education 362 (11.3) 386 (12.1) male and 2072 female; and 25 years, 1056 male and Vocational 282 (8.8) 256 (8.0) 1685 female. Numbers of participants at each age General certificate of education, ordinary level 1160 (36.3) 1107 (34.6) after removing bone measures recorded before age General certificate of education, advanced level 862 (27.0) 891 (27.8) at peak height velocity: 10 years, 0 male and 57 female; 12 years, 55 male and 1518 female; 14 years, Degree or higher 527 (16.5) 556 (17.4) 1528 male and 2474 female; 16 years, 1741 male and BMI at age 7 y 16.1 (1.9) 16.3 (2.1) 2036 female; 18 years, 1674 male and 2116 female; Daily energy intake at age 7 y, kJ/d 7785.5 (1878.8) 7477.7 (1759.8) and 25 years, 1056 male and 1725 female. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 6/14 Height Velocity, cm/y Height Velocity, cm/y JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Age at Puberty and Site-Specific Bone Trajectories eFigures 4 to 9 in the Supplement show the mean estimated BMD and BMC trajectories for arms, legs, trunk, spine, ribs, and pelvis, respectively, in male and female participants in the 10th, 50th, and 90th APHV percentiles. Overall, these were similar to the main whole-body trajectories in that those in the oldest APHV groups started lower and exhibited some catch-up over follow-up. eFigure 10 in the Supplement shows similar mean estimated hip BMD trajectories. Male and female participants in the oldest APHV group had lower total hip and femoral neck BMD at age 14 years, and these differences were more substantially reduced in male than female participants at ages 18 and 25 years (eFigure 10 in the Supplement). However, the lack of hip measures from earlier age makes it difficult to discern whether these patterns are different from trajectories of whole-body parameters. Discussion We investigated the association between age at puberty (measured by APHV) and bone accrual from ages 10 to 25 years in male and female participants from a large prospective British birth cohort. The BMD and BMC increased over follow-up with sex differences emerging during puberty. Male participants accrued BMD and BMC at faster rates than did female participants from before puberty and up to 4 years after APHV. The fastest gains in BMD and BMC in both male and female participants were during the year before APHV to 2 years after APHV. Older pubertal age was associated with an Figure 2. Height Growth and Height Velocity Curves A Male participants 200 12 Mean –1.96 SD +1.96 SD 10 100 0 5 10 15 20 Age, y Female participants 200 12 Mean height growth curves (solid lines), mean height velocity curves (dashed lines), and mean age at peak height velocity (plus estimated curves and age at peak height velocity at 1.96 SD and −1.96 SD) (dotted 100 0 vertical lines) are shown for male (A) and female (B) 5 10 15 20 participants included in the Superimposition by Age, y Translation and Rotation models. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 7/14 Height, cm Height, cm JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual initially accelerating and subsequently decelerating faster BMD and BMC accrual over follow-up. This catch-up was greater in male participants than in female participants between ages 16 and 18 years for BMD and between ages 14 and 25 years for BMC. Despite these faster gains, older pubertal age remained associated with lower BMD throughout follow-up to age 25 years and with a BMC that was lower at younger ages but higher by age 25 years. Findings were overall similar for site-specific BMD and BMC accrual. The findings that the greatest gains in BMD and BMC occurred during the year before and 2 years after APHV and that gains in BMC ceased 4 years after APHV agree with those of a recent study of the temporal association between peak height and peak bone accrual; however, that study did not test the association between APHV and bone accrual. That study, which included 2000 participants aged 5 to 19-years who underwent annual DXA scans for up to 7 years, showed that up to 36% of BMC is acquired in the 2 years before and 2 years after APHV. Our study expands on this by showing that, in contrast to BMC, gains in BMD continued into adulthood. We also showed that male participants had faster gains in both BMD and BMC from before puberty than did female participants but that sex differences in BMD and BMC emerged during puberty, which is consistent with the suggestion that sex hormones are driving sexual dimorphism in body composition. Our study contributes to previous research by documenting a process of transient pubertal catch-up in BMD and BMC among those with older pubertal age that did not fully eliminate differences in BMD. This agrees with findings that differences in BMD and BMC by APHV group were Figure 3. Gains in Bone Mineral Density (BMD) and Bone Mineral Content (BMC) A Gains in BMD from childhood to adulthood B Gains in BMC from childhood to adulthood 0.15 700 Male participants 0.12 Female participants 0.09 0.06 0.03 0 0 –8 to –1 –1 to 2 2 to 4 4 to 16 –8 to –1 –1 to 2 2 to 4 4 to 16 Time Before or After Puberty, y Time Before or After Puberty, y C Gains in BMD from puberty to adulthood D Gains in BMC from puberty to adulthood 0.025 140 0.020 0.015 0.010 0.005 0 0 ≤14 14-16 16-18 18-25 ≤14 14-16 16-18 18-25 Age, y Age, y Gains in BMD and BMC are shown from childhood (prepuberty) to adulthood, with age denote means, and whiskers denote 95% confidence intervals. Estimates were adjusted centered at pubertal age (ie, age at peak height velocity [APHV]) (A, BMD; B, BMC) and for birth weight, ethnicity, maternal education, and childhood body mass index and from puberty to adulthood per 1-year older age at puberty (C, BMD; D, BMC). Circles dietary intake. These estimates are presented in eTable 3 in the Supplement. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 8/14 Whole-body BMD, g/cm /y Whole-body BMD, g/cm /y Whole-body BMC, g/y Whole-body BMC, g/y JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual attenuated at follow-up in a mixed sample of 230 participants aged 8 to 14 years with long-term follow-up and in 500 male participants aged 18 years who were assessed at baseline and 5 years later. Our work adds to these studies by having considerably larger numbers and repeated measures of BMD from ages 10 to 25 years and by including female and male participants. This allowed us to demonstrate transient pubertal catch-up gains in BMD among those with older pubertal age, which were greater for male than for female participants. In contrast to BMD, we found that later puberty was associated with more prolonged catch-up in BMC, leading to higher total levels in adulthood. This finding may be explained by those reaching puberty later ultimately ending up 17,32 taller and, thus, having bigger (albeit less dense) bones, but further research is required to disentangle the complex associations between pubertal growth and bone mineralization up to adulthood. Given the persisting associations between later puberty and lower BMD reported here and evidence that peak bone mineralization lags behind peak height accrual, our findings suggest that 15,19,30 adolescents who mature later may be at higher risk of fractures throughout adolescence. Because differences in BMD persisted up to age 25 years, including at the spine and hip sites, which are preferentially affected by osteoporotic fracture in later life, those with older pubertal age could also be at increased risk of osteoporosis in later life, although continued follow-up of this cohort is required to identify how associations might vary through adult life. Figure 4. Estimated Bone Mineral Density (BMD) and Bone Mineral Content (BMC) Trajectories A Estimated BMD trajectory from childhood to adulthood B Estimated BMC trajectory from childhood to adulthood 1.4 3500 1.3 1.2 1.1 1.0 0.9 Male participants 1000 Male participants 0.8 Female participants Female participants 0.7 0.6 0 –8 –6 –4 –2 0 2 4 6 8 10 12 14 16 18 –8 –6 –4 –2 0 2 4 6 8 10 12 14 16 18 Time Since Peak Height Velocity, y Time Since Peak Height Velocity, y C Estimated BMD trajectory from puberty to adulthood D Estimated BMC trajectory from puberty to adulthood Male participants Female participants Male participants Female participants 1.4 3500 1.3 APHV percentile APHV percentile th th 10 10 1.2 th th 50 50 th 2500 th 90 90 1.1 1.0 0.9 0.8 1000 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 10 12 14 16 18 20 22 24 26 Age, y Age, y Mean estimated BMD and BMC trajectories are shown from childhood (prepuberty) to the 10th (<12.4 years in male and <10.6 years in female participants), 50th (13.3-13.5 years adulthood (A, BMD; B, BMC) against time since puberty (ie, age at peak height velocity in male and 11.4-11.6 years in female participants), and 90th (>14.5 years in male and >12.7 [APHV], represented by vertical line), and from puberty to adulthood for individuals in years in female participants) pubertal age (APHV) percentiles (C, BMD; D, BMC). JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 9/14 2 2 Whole-body BMD, g/cm Whole-body BMD, g/cm Whole-body BMC, g Whole-body BMC, g JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Limitations To our knowledge, this is the largest study on BMD in adolescents to date and has a longer follow-up than previous studies, but some limitations should be noted. Dual-energy x-ray absorptiometry is routinely used in clinical practice to assess BMD and osteoporosis; however, measurement error and changing soft-tissue distribution might influence findings. Newer imaging approaches could be used to investigate structural changes in bone microarchitecture, including growth-associated cortical porosity. Using APHV as measure of age at puberty is a key strength because it allows the same analyses in male and female participants and is less prone to differential measurement error than pubertal staging using reports of sexual maturity, such as Tanner stages. Measurement error, particularly from the self or parental reports of height, may have influenced our results. However, because parents would be unaware of their child’s subsequent bone measurements, that is unlikely to have been systematic. Furthermore, restricting analyses to research-clinic measures of height only produced similar results. In addition, similar results were obtained when using age at menarche as an alternative marker of pubertal age in female participants. The linear spline models reduced bias resulting from missing bone data by including all 27,28 participants with at least 1 bone measure under the missing-at-random assumption (ie, that the probability of missing bone data is assumed to be associated with model covariates and not associated with unmeasured factors). Although it is not possible to fully test that the assumption of missing at random holds, the probability of missing bone data among those included in the linear spline analysis was associated with the variables included in the models (eTable 1 in the Supplement). Those excluded from the analyses because of missing all 6 DXA scans were socioeconomically different from the analytic sample (eTable 2 in the Supplement), which might limit the generalizability of our findings. Participants with missing data on covariates (28% of those potentially eligible) were excluded from the analyses, which might introduce a bias if they had systematically different bone measurements; however, this seems unlikely because participants would not have 35,36 known their bone measurements. Our study sample were mostly of white British ethnicity, which means that the findings may not be generalizable to participants of other ethnicities, with potentially different bone accrual rates. Also, our analyses were performed in 2 stages rather than 37,38 as a single joint model, which might introduce a bias ; however, given the complex nonlinear SITAR models, joint modeling is unlikely to be feasible for this study. Conclusions This large and long-running prospective follow-up cohort study showed that later puberty was associated with persistently lower BMD from ages 10 to 25 years in male and female participants, despite faster gains in BMD during puberty in those with older pubertal age. Future studies should seek more robust evidence of associations between puberty timing and bone accrual from large emerging collaborations such as the European Union’s Child Cohort Network, which can provide more repeated measures over longer follow-up in populations with differing confounding structures and with larger sample sizes that can support methods that are statistically inefficient, such as within 40,41 sibship and mendelian randomization analyses. Our findings suggest that advice on how to 42,43 increase and maintain BMD such as through physical activity should be offered to people with older pubertal age. ARTICLE INFORMATION Accepted for Publication: June 17, 2019. Published: August 9, 2019. doi:10.1001/jamanetworkopen.2019.8918 Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2019 Elhakeem A et al. JAMA Network Open. JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 10/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual Corresponding Author: Ahmed Elhakeem, PhD, Medical Research Council Integrative Epidemiology Unit at University of Bristol, Population Health Sciences, Bristol Medical School, Oakfield House, Oakfield Grove, Bristol BS8 2BN, United Kingdom (a.elhakeem@bristol.ac.uk). Author Affiliations: Medical Research Council Integrative Epidemiology Unit at University of Bristol, Population Health Sciences, Bristol Medical School, Bristol, United Kingdom (Elhakeem, Tilling, Tobias, Lawlor); Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom (Elhakeem, Tilling, Lawlor); Musculoskeletal Research Unit, Translational Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom (Frysz, Tobias). Author Contributions: Dr Elhakeem had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: Elhakeem, Tobias, Lawlor. Acquisition, analysis, or interpretation of data: Elhakeem, Frysz, Tilling, Lawlor. Drafting of the manuscript: Elhakeem, Lawlor. Critical revision of the manuscript for important intellectual content: Elhakeem, Frysz, Tilling, Tobias. Statistical analysis: Elhakeem, Tilling. Obtained funding: Lawlor. Administrative, technical, or material support: Frysz. Supervision: Lawlor. Conflict of Interest Disclosures: Dr Elhakeem reported grants from European Union’s Horizon 2020 research and innovation program outside the submitted work. Dr Tilling reported grants from the Economic and Social Research Council and grants from the UK Medical Research Council (MRC) during the conduct of the study. Dr Lawlor reported grants from the MRC, Wellcome, and European Union’s Horizon 2020 program during the conduct of the study, and grants from several national and international government and charitable research funders and Roche Diagnostics and Medtronic Ltd outside the submitted work. No other disclosures were reported. Funding/Support: This work was supported by the European Union’s Horizon 2020 research and innovation program under grant agreement 733206 (LifeCycle), which pays the salary of Dr Elhakeem. Drs Elhakeem and Lawlor work in a unit that receives funds from the University of Bristol and UK MRC (grant MC_UU_00011/6). Dr Lawlor is a National Institute of Health Research Senior Investigator (grant NF-SI-0611-10196). The UK MRC and Wellcome (grant r102215/2/13/2) and the University of Bristol provide core support for the Avon Longitudinal Study of Parents and Children (ALSPAC). A comprehensive list of grant funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. Disclaimer: Views expressed in this article are those of the authors and not necessarily any funder. Additional Contributions: We thank all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, and nurses. Additional Information: Details of all available data and the processes and procedures involved in accessing the ALSPAC resource can be found in the ALSPAC study website, which includes a fully searchable data dictionary and variable search tool (http://www.bristol.ac.uk/alspac/researchers/our-data/). REFERENCES 1. Hansen MA, Overgaard K, Riis BJ, Christiansen C. Role of peak bone mass and bone loss in postmenopausal osteoporosis: 12 year study. BMJ. 1991;303(6808):961-964. doi:10.1136/bmj.303.6808.961 2. Harvey N, Dennison E, Cooper C. Osteoporosis: a lifecourse approach. J Bone Miner Res. 2014;29(9):1917-1925. doi:10.1002/jbmr.2286 3. Weaver CM, Gordon CM, Janz KF, et al. The National Osteoporosis Foundation’s position statement on peak bone mass development and lifestyle factors: a systematic review and implementation recommendations. Osteoporos Int. 2016;27(4):1281-1386. doi:10.1007/s00198-015-3440-3 4. Hernandez CJ, Beaupré GS, Carter DR. A theoretical analysis of the relative influences of peak BMD, age-related bone loss and menopause on the development of osteoporosis. Osteoporos Int. 2003;14(10):843-847. doi:10. 1007/s00198-003-1454-8 5. Heaney RP, Abrams S, Dawson-Hughes B, et al. Peak bone mass. Osteoporos Int. 2000;11(12):985-1009. doi:10. 1007/s001980070020 JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 11/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 6. Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005-2025. J Bone Miner Res. 2007;22(3):465-475. doi:10. 1359/jbmr.061113 7. Hopkins RB, Burke N, Von Keyserlingk C, et al. The current economic burden of illness of osteoporosis in Canada. Osteoporos Int. 2016;27(10):3023-3032. doi:10.1007/s00198-016-3631-6 8. Rogol AD, Roemmich JN, Clark PA. Growth at puberty. J Adolesc Health. 2002;31(6)(suppl):192-200. doi:10. 1016/S1054-139X(02)00485-8 9. Patton GC, Viner R. Pubertal transitions in health. Lancet. 2007;369(9567):1130-1139. doi:10.1016/S0140-6736 (07)60366-3 10. Cousminer DL, Widén E, Palmert MR. The genetics of pubertal timing in the general population: recent advances and evidence for sex-specificity. Curr Opin Endocrinol Diabetes Obes. 2016;23(1):57-65. doi:10.1097/ MED.0000000000000213 11. Day FR, Bulik-Sullivan B, Hinds DA, et al. Shared genetic aetiology of puberty timing between sexes and with health-related outcomes. Nat Commun. 2015;6:8842. doi:10.1038/ncomms9842 12. Kemp JP, Morris JA, Medina-Gomez C, et al. Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis. Nat Genet. 2017;49(10):1468-1475. doi:10.1038/ ng.3949 13. Cole TJ, Kuh D, Johnson W, et al. Using Super-Imposition by Translation And Rotation (SITAR) to relate pubertal growth to bone health in later life: the Medical Research Council (MRC) National Survey of Health and Development. Int J Epidemiol. 2016;45(4):1125-1134. 14. Cousminer DL, Mitchell JA, Chesi A, et al. Genetically determined later puberty impacts lowered bone mineral density in childhood and adulthood. J Bone Miner Res. 2018;33(3):430-436. doi:10.1002/jbmr.3320 15. Bonjour JP, Chevalley T. Pubertal timing, bone acquisition, and risk of fracture throughout life. Endocr Rev. 2014;35(5):820-847. doi:10.1210/er.2014-1007 16. Day FR, Elks CE, Murray A, Ong KK, Perry JR. Puberty timing associated with diabetes, cardiovascular disease and also diverse health outcomes in men and women: the UK Biobank study. Sci Rep. 2015;5:11208. doi:10.1038/ srep11208 17. Cousminer DL, Berry DJ, Timpson NJ, et al; ReproGen Consortium; Early Growth Genetics (EGG) Consortium. Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity. Hum Mol Genet. 2013;22(13):2735-2747. doi:10.1093/hmg/ddt104 18. Jackowski SA, Erlandson MC, Mirwald RL, et al. Effect of maturational timing on bone mineral content accrual from childhood to adulthood: evidence from 15 years of longitudinal data. Bone. 2011;48(5):1178-1185. doi:10. 1016/j.bone.2011.02.010 19. Darelid A, Ohlsson C, Nilsson M, Kindblom JM, Mellström D, Lorentzon M. Catch up in bone acquisition in young adult men with late normal puberty. J Bone Miner Res. 2012;27(10):2198-2207. doi:10.1002/jbmr.1675 20. Fraser A, Macdonald-Wallis C, Tilling K, et al. Cohort profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int J Epidemiol. 2013;42(1):97-110. doi:10.1093/ije/dys066 21. Boyd A, Golding J, Macleod J, et al. Cohort profile: the ‘children of the 90s’—the index offspring of the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2013;42(1):111-127. doi:10.1093/ije/dys064 22. Avon Longitudinal Study of Parents and Children. Explore data and samples. http://www.bristol.ac.uk/alspac/ researchers/our-data/. Accessed May 1, 2019. 23. Lawlor DA, Wills AK, Fraser A, Sayers A, Fraser WD, Tobias JH. Association of maternal vitamin D status during pregnancy with bone-mineral content in offspring: a prospective cohort study. Lancet. 2013;381(9884): 2176-2183. doi:10.1016/S0140-6736(12)62203-X 24. Warrington NM, Kemp JP, Tilling K, Tobias JH, Evans DM. Genetic variants in adult bone mineral density and fracture risk genes are associated with the rate of bone mineral density acquisition in adolescence. Hum Mol Genet. 2015;24(14):4158-4166. doi:10.1093/hmg/ddv143 25. Emmett PM, Jones LR, Northstone K. Dietary patterns in the Avon Longitudinal Study of Parents and Children. Nutr Rev. 2015;73(suppl 3):207-230. doi:10.1093/nutrit/nuv055 26. Cole TJ, Donaldson MD, Ben-Shlomo Y. SITAR—a useful instrument for growth curve analysis. Int J Epidemiol. 2010;39(6):1558-1566. doi:10.1093/ije/dyq115 27. Howe LD, Tilling K, Matijasevich A, et al. Linear spline multilevel models for summarising childhood growth trajectories: a guide to their application using examples from five birth cohorts. Stat Methods Med Res. 2016;25 (5):1854-1874. doi:10.1177/0962280213503925 JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 12/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual 28. Tilling K, Macdonald-Wallis C, Lawlor DA, Hughes RA, Howe LD. Modelling childhood growth using fractional polynomials and linear splines. Ann Nutr Metab. 2014;65(2-3):129-138. doi:10.1159/000362695 29. Frysz M, Howe LD, Tobias JH, Paternoster L. Using SITAR (SuperImposition by Translation and Rotation) to estimate age at peak height velocity in Avon Longitudinal Study of Parents and Children. Wellcome Open Res. 2018;3:90. doi:10.12688/wellcomeopenres.14708.1 30. McCormack SE, Cousminer DL, Chesi A, et al. Association between linear growth and bone accrual in a diverse cohort of children and adolescents. JAMA Pediatr. 2017;171(9):e171769. doi:10.1001/jamapediatrics.2017.1769 31. Wells JC. Sexual dimorphism of body composition. Best Pract Res Clin Endocrinol Metab. 2007;21(3):415-430. doi:10.1016/j.beem.2007.04.007 32. Carel JC, Lahlou N, Roger M, Chaussain JL. Precocious puberty and statural growth. Hum Reprod Update. 2004;10(2):135-147. doi:10.1093/humupd/dmh012 33. Kanis JA. Diagnosis of osteoporosis and assessment of fracture risk. Lancet. 2002;359(9321):1929-1936. doi: 10.1016/S0140-6736(02)08761-5 34. Farr JN, Khosla S. Skeletal changes through the lifespan—from growth to senescence. Nat Rev Endocrinol. 2015;11(9):513-521. doi:10.1038/nrendo.2015.89 35. Munafò MR, Tilling K, Taylor AE, Evans DM, Davey Smith G. Collider scope: when selection bias can substantially influence observed associations. Int J Epidemiol. 2018;47(1):226-235. doi:10.1093/ije/dyx206 36. Hughes RA, Heron J, Sterne JAC, Tilling K. Accounting for missing data in statistical analyses: multiple imputation is not always the answer. Int J Epidemiol. 2019;dyz032. doi:10.1093/ije/dyz032 37. Sayers A, Heron J, Smith A, et al. Joint modelling compared with two stage methods for analysing longitudinal data and prospective outcomes: a simulation study of childhood growth and BP. Stat Methods Med Res. 2017;26 (1):437-452. doi:10.1177/0962280214548822 38. Macdonald-Wallis C, Lawlor DA, Palmer T, Tilling K. Multivariate multilevel spline models for parallel growth processes: application to weight and mean arterial pressure in pregnancy. Stat Med. 2012;31(26):3147-3164. doi: 10.1002/sim.5385 39. LifeCycle. European Union Child Cohort Network. https://lifecycle-project.eu. Accessed May 1, 2019. 40. Lawlor DA, Tilling K, Davey Smith G. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016;45(6): 1866-1886. 41. Munafò MR, Davey Smith G. Robust research needs many lines of evidence. Nature. 2018;553(7689): 399-401. doi:10.1038/d41586-018-01023-3 42. Mitchell JA, Chesi A, McCormack SE, et al. Physical activity and bone accretion: isotemporal modeling and genetic interactions. Med Sci Sports Exerc. 2018;50(5):977-986. doi:10.1249/MSS.0000000000001520 43. Muthuri SG, Ward KA, Kuh D, Elhakeem A, Adams JE, Cooper R. Physical activity across adulthood and bone health in later life: the 1946 British birth cohort. J Bone Miner Res. 2019;34(2):252-261. SUPPLEMENT. eMethods. Supplemental Methods. eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV eReferences eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty eFigure 1. Observed Whole-body BMD and BMC at Each Follow-up Assessment From 10 to 25 Years eFigure 2. Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 13/14 JAMA Network Open | Pediatrics Association Between Age at Puberty and Bone Accrual eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles JAMA Network Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 (Reprinted) August 9, 2019 14/14 Supplementary Online Content Elhakeem A, Frysz M, Tilling K, Tobias JH, Lawlor DA. Association between age at puberty and bone accrual from 10 to 25 years of age. JAMA Netw Open. 2019;2(8):e198918. doi:10.1001/jamanetworkopen.2019.8918 eMethods. Supplemental Methods eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV eReferences eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty eFigure 1. Observed Whole-body BMD and BMC at Each Follow-up Assessment From 10 to 25 Years eFigure 2. Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for Iindivudals in the 10th, 50th, and 90th APHV Percentiles © 2019 Elhakeem A et al. JAMA Network Open. eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Indivudals in the 10th, 50th, and 90th APHV Percentiles This supplementary material has been provided by the authors to give readers additional information about their work. © 2019 Elhakeem A et al. JAMA Network Open. eMethods. Supplemental Methods. APHV was estimated using Superimposition by Translation and Rotation (SITAR) growth curve analysis (1, 2). SITAR is a shape invariant model consisting of a natural cubic spline mean curve of height vs age. SITAR fits this curve with three person-specific random-effects that reflect individual differences from the mean curve in size (random height intercept for person-specific shift up or down in the spline curve), velocity (random age scaling that reflects differences in the duration of the growth spurt in individuals) and tempo (random age intercept reflecting individual differences in the timing of the growth spurt). SITAR models were fitted separately for males and females with 5 degrees of freedom. APHV estimation was restricted to individuals with height measurements between ages 5 and 20 years and at least one height measure after age 9 years. Data cleaning similar to previous ALSPAC analyses was performed to remove data points with velocity exceeding 4 SDs and standardized residuals exceeding 4 in absolute value (2, 3, 4), leaving a total of 5223 males (with 55 752 height measurements) and 5322 females (with 59 318 height measurements) with an estimate of APHV. The final SITAR models explained 95.7% and 95.6% of variance in height in males and females respectively. APHV in years was estimated by transforming the random age intercept that reflects individual differences in the tempo (timing) of the growth spurt. We subsequently investigated the relationship between APHV and whole-body BMD and BMC accrual using two sets of mixed-effects linear spline regression analyses. Linear spline models facilitate the summarising of nonlinear bone accrual through a series of linear splines joined at knot points (5, 6). They estimate mean and person-specific bone values at the first assessment and mean and person-specific rates of accrual during each of the different growth periods. Smoothing functions, prior subject knowledge and fit statistics were used to identify the best fitting models (5, 6) and subsequently knot point selection was based on the mean age at DXA scans. Models were fitted using maximum likelihood and included random intercepts and random age slopes. A st 1 order autoregressive correlation structure was used to account for temporal autocorrelation, and a fixed variance structure was used to model heteroskedasticity. In both sets of linear spline models, we examined the relationship between APHV and rates of bone accrual by including APHV-by-age splines interaction terms and included sex-by-age splines interactions terms to test for statistical evidence of different bone trajectories for males and females. All models were adjusted for birth weight, ethnicity, maternal education and childhood BMI and energy intake. the observed and predicted BMD and BMC values (without random effects) were 0.84 and 0.83. © 2019 Elhakeem A et al. JAMA Network Open. In the first analyses, age was centred at APHV in order to estimate rates of bone accrual relative to APHV; i.e. from before puberty up to adulthood. Knot points were placed at the mean ages of the 12, 14, and 16 -year scans which equated to summarising rate of accrual from 7.7 up to 0.9 years prior to APHV, 0.9 year prior up to 2.1 years following APHV, 2.1 to 3.7 years following APHV, and 3.7 to 16 years after APHV. The mean predicted trajectories relative to APHV were plotted for males and females. For the second analyses, we excluded individuals with accrual from before APHV in order to investigate associations between APHV and subsequent rate of bone accrual up to adulthood. Here, knot points were placed at the mean ages of the 14 (centred at age 14), 16, and 18-year assessments, which then estimated associations with rate of accrual up to age 14, between 14 16, 16 18 and 18 25 years. To assess if difference in bone outcomes persisted into adulthood, predicted bone values at age 25 were regressed on APHV (after adjusting for predicted values at younger ages) and compared to associations with predicted values at age 14. To visualise associations between APHV and bone th th th accrual, predicted mean trajectories were plotted for individuals in the 10 , 50 , and 90 sex-specific APHV percentiles. In sensitivity analyses, we refitted models after re-estimating APHV from measured heights only (3) and using age at menarche (reported prospectively by the mother in years and months) instead of APHV as an alternative marker of age at puberty in females. For our secondary outcomes, we investigated associations between APHV and subsequent site-specific BMD and BMC accrual from whole-body DXA, and with total hip and femoral neck BMD accrual from hip DXA. Given that hips were only scanned at three different ages (14, 18 and 25), we fitted hip models using a single knot point placed at the mean ages of the 18-year assessment. The predicted mean trajectories were plotted for all secondary outcomes. All analyses were performed in R (R Foundation for Statistical Computing, Vienna); SITAR models were fitted wit The linear spline models reduced bias due to missing bone data by including all measures under the missing at random assumption (i.e. the probability of missing bone data assumed to be related to measured factors but unrelated to unmeasured factors). eTable 1 shows differences between those included in the analysis with complete and partial DXA data, and suggests that the probability of missing at least one bone measure was related to the model covariates, and thus possibly missing at random. Those missing all DXA measures were excluded from the analysis. eTable 2 shows differences between those included in the linear spline models and excluded due to missing all DXA data, and shows that compared with the participants in this study, those © 2019 Elhakeem A et al. JAMA Network Open. excluded due to missing bone data had lower socioeconomic position and had slightly lower birthweight and slightly higher BMI (eTable 2). Those missing data on model covariates were also excluded from the analysis. © 2019 Elhakeem A et al. JAMA Network Open. eTable 1. Comparison on Model Covariates Between Participants Included in the Linear Spline Models That Had All Six DXA Measures With Those Also Included in the Analysis But Were Missing at Least One DXA Measure Males Females Included in Included in Included in Included in analyses: all six analyses: missing analyses: all six analyses: missing DXA measures DXA measures available (n=585) (n=2608) available (n=993) (n=2203) Birth weight g [mean (SD)] 3484.1 (548.2) 3472.0 (597.8) 3375.3 (499.7) 3376.1 (509.3) Ethnicity [No. (%)] White 575 (98.3) 2562 (98.2) 972 (97.9) 2164 (98.2) Non-white 10 (1.7) 46 (1.8) 21 (2.1) 39 (1.8) Maternal education [No. (%)] Certificate of Secondary Education 37 (6.3) 325 (12.5) 58 (5.8) 328 (14.9) Vocational 34 (5.8) 248 (9.5) 54 (5.4) 202 (9.2) General Certificate of Education: Ordinary Level 175 (29.9) 985 (37.8) 343 (34.5) 764 (34.7) General Certificate of Education: Advanced Level 179 (30.6) 683 (26.2) 310 (31.2) 581 (26.4) Degree or higher 160 (27.4) 367 (14.1) 228 (23.0) 328 (14.9) BMI at age 7y kg/m [mean (SD)] 15.8 (1.7) 16.1 (1.9) 16.0 (1.8) 16.4 (2.2) Daily energy intake at age 7y kJ/day [mean (SD)] 7712.9 (1603.2) 7801.7 (1935.2) 7431.4 (1612.6) 7498.5 (1821.7) APHV years [mean (SD)] 13.4 (0.9) 13.5 (0.9) 11.7 (0.8) 11.6 (0.8) © 2019 Elhakeem A et al. JAMA Network Open. eTable 2. Comparison on Model Covariates Between Those Included in the Linear Spline Models With Those That Were Excluded From the Analysis Due to Missing All Six DXA Measures, Despite Having Data on APHV Males Females Included in linear Excluded due to Included in linear Excluded due to spline models missing all DXA spline models missing all (n=3193) (n=851) (n=3196) DXA (n=755) Birth weight g [mean (SD)] 3474.2 (574.0) 3449.4 (578.1) 3375.8 (506.3) 3372.7 (506.3) Ethnicity [No. (%)] White 3137 (98.3) 727 (96.7) 3136 (98.1) 621 (97.8) Non-white 56 (1.8) 25 (3.3) 60 (1.9) 14 (2.2) Maternal education [No. (%)] Certificate of Secondary Education 362 (11.3) 161 (21.3) 386 (12.1) 172 (26.8) Vocational 282 (8.8) 74 (9.8) 256 (8.0) 75 (11.7) General Certificate of Education: Ordinary Level 1160 (36.3) 266 (35.2) 1107 (34.6) 217 (33.8) General Certificate of Education: Advanced Level 862 (27.0) 157 (20.8) 891 (27.8) 113 (17.6) Degree or higher 527 (16.5) 97 (12.9) 556 (17.4) 66 (10.3) BMI at age 7y kg/m [mean (SD)] 16.1 (1.9) 16.2 (1.9) 16.3 (2.1) 16.8 (2.9) Daily energy intake at age 7y kJ/day [mean (SD)] 7785.5 (1878.8) 7731.9 (2068.8) 7477.7 (1759.8) 7459.9 (2029.1) APHV years [mean (SD)] 13.5 (0.9) 13.5 (0.7) 11.6 (0.8) 11.7 (0.8) © 2019 Elhakeem A et al. JAMA Network Open. eReferences 1. Cole TJ, Donaldson MD, Ben-Shlomo Y. SITAR--a useful instrument for growth curve analysis. Int J Epidemiol. 2010;39(6):1558-66. 2. Cole TJ, Kuh D, Johnson W, Ward KA, Howe LD, Adams JE, et al. Using Super-Imposition by Translation And Rotation (SITAR) to relate pubertal growth to bone health in later life: the Medical Research Council (MRC) National Survey of Health and Development. Int J Epidemiol. 2016. 3. Frysz M, Howe L, Tobias J, Paternoster L. Using SITAR (SuperImposition by Translation and Rotation) to estimate age at peak height velocity in Avon Longitudinal Study of Parents and Children. Wellcome Open Res; 2018. 4. Mahmoud O, Granell R, Tilling K, Minelli C, Garcia-Aymerich J. Association of Height Growth in Puberty with Lung Function: A Longitudinal Study. Am J Respir Crit Care Med. 2018;doi: 10.1164/rccm.201802 - 0274OC 5. Tilling K, Macdonald-Wallis C, Lawlor DA, Hughes RA, Howe LD. Modelling childhood growth using fractional polynomials and linear splines. Annals of nutrition & metabolism. 2014;65(2-3):129-38. 6. Howe LD, Tilling K, Matijasevich A, Petherick ES, Santos AC, Fairley L, et al. Linear spline multilevel models for summarising childhood growth trajectories: A guide to their application using examples from five birth cohorts. Statistical methods in medical research. 2016;25(5):1854-74. © 2019 Elhakeem A et al. JAMA Network Open. eTable 3. Gains in BMD and BMC During Each Period (A) From Childhood (Prepuberty) to Adulthood, With Age Centered at Pubertal Age (ie, APHV) and (B) From Puberty to Adulthood Per 1-Year Older Age at Puberty BMD g/cm /yr. (95% CI) BMC g/yr. (95% CI) Males Females Males Females A. Rates relative to APHV (n= 6389)* -7.7 to -0.9 yrs. 0.054 (0.047 to 0.062) 0.014 (0.000 to 0.028) 289.109 (256.887 to 321.331) 74.065 (5.229 to 139.537) -0.9 to 2.1 yrs 0.137 (0.125 to 0.150) 0.106 (0.098 to 0.115) 644.480 (592.822 to 696.138) 428.786 (390.064 to 463.701) 2.1 to 3.7 yrs. 0.078 (0.060 to 0.096) 0.047 (0.037 to 0.058) 456.124 (378.401 to 533.847) 312.841 (263.411 to 306.351) 3.7 to 16.0 yrs. 0.013 (0.005 to 0.021) 0.005 (0.001 to 0.009) -35.423 (-69.913 to -0.934) 16.057 (-2.495 to 34.609) B. Rates per 1-year older APHV (n= 5477)** up to 14 yrs 0.012 (0.003 to 0.020) 0.010 (0.009 to 0.011) 92.259 (49.562 to 135.955) 40.902 (35.538 to 46.265) 14 to 16 yrs 0.013 (0.011 to 0.015) 0.014 (0.014 to 0.015) 85.125 (76.005 to 94.244) 61.378 (56.995 to 65.761) 16 to 18 yrs 0.011 (0.010 to 0.013) 0.003 (0.003 to 0.004) 64.252 (58.846 to 69.658) 19.374 (15.398 to 23.350) 18 to 25 yrs 0.002 (0.001 to 0.003) 0.000 (-0.001 to 0.000) 13.368 (10.823 to 15.913) 5.793 (3.902 to 7.683) * Knots at the 12, 14, 16 year assessments. Age centred at APHV (i.e. 0=APHV): models includes APHV:age splines interaction terms, and adjsutment for birth weight, maternal education, childhood BMI and diet. ** Knots at the 14, 16 and 18 year assessments. Models includes APHV:age splines interaction terms, and adjsutment for birth weight, ethnicity, maternal education and childhood BMI and diet. These models excluded bone measures recorded before APHV. © 2019 Elhakeem A et al. JAMA Network Open. eFigure 1 Observed Whole-body BMD and BMC at Each Follow-up Assessment from 10 to 25 Years © 2019 Elhakeem A et al. JAMA Network Open. eFigure 2 Predicted Mean Trajectories of Whole-body BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles, Using APHV From Measured Heights (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 3. Predicted Mean Trajectories of Whole-body BMD and BMC From Menarche to Adulthood for Females in the 10th, 50th, and 90th Age at Menarche Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 4. Predicted Mean Trajectories of Arms BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 5. Predicted Mean Trajectories of Legs BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 6. Predicted Mean Trajectories of Trunk BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 7. Predicted Mean Trajectories of Spine BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 8. Predicted Mean Trajectories of Ribs BMD and BMC From APHV to Adulthood for IIndividuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 9. Predicted Mean Trajectories of Pelvis BMD and BMC From APHV to Adulthood for Individuals in the 10th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open. eFigure 10. Predicted Mean Trajectories of Total Hip and Femur Neck BMD From APHV to Adulthood for Individuals in the 10 th, 50th, and 90th APHV Percentiles (Dotted points represent predicted values.) © 2019 Elhakeem A et al. JAMA Network Open.

Journal

JAMA Network OpenAmerican Medical Association

Published: Aug 9, 2019

There are no references for this article.