|Home | About | Journals | Submit | Contact Us | Français|
Elucidating the somatic and maturational influences on the biomechanical properties of bone in children is crucial for a proper understanding of bone strength and quality in childhood and later life, and has significant potential for predicting adult fracture and osteoporosis risks. The ability of a long bone to resist bending and torsion is primarily a function of its cross-sectional geometric properties, and is negatively impacted by smaller external bone diameter. In pubescent girls, elevated levels of estrogen impede subperiosteal bone growth and increase endosteal bone deposition, resulting in bones averaging a smaller external and internal diameter relative to boys. In addition, given a well-documented secular trend for an earlier menarche, the age at which the rate of subperiosteal bone deposition decreases may also be younger in more recent cohorts of girls.
In this study we examined the relationship between pubertal timing and subsequent bone strength in girls. Specifically, we investigated the effects of age at menarche on bone strength indicators (polar moment of inertia and section modulus) determined from cross-sectional geometry of the second metacarpal (MC2) using data derived from serial hand-wrist radiographs of female participants (N=223) in the Fels Longitudinal Study, with repeated measures of MC2 between the ages of 7 and 35 years. Using multivariate regression models, we evaluated the effects of age at menarche on associations between measures of bone strength in early adulthood and the same measures at a prepubertal age. Results indicate that later age at menarche is associated with stronger adult bone (in torsion and bending) when controlling for prepubertal bone strength (R2 ranged between 0.54 and 0.70, p<0.001). Since cross-sectional properties of bone in childhood may have long lasting implications, they should be considered along with pubertal timing in assessing risk for future fracture and in clinical recommendations.
Bone strength depends on its material properties, including bone mineral content (BMC), bone mineral density (BMD), and the spatial distribution of bone tissue (i.e., bone geometry) . Resistance to the principal loads that cause fracture (e.g., bending and torsion) is primarily a function of a long bone's cross-sectional geometric properties . In general, larger external bone diameters provide greater strength in bending and torsion than smaller external bone diameters, when all else is held equal (e.g., amount of cortex, BMD, BMC). Significant sex differences in bone diameter exist, with girls tending to have bones of smaller external diameter than boys, regardless of overall body size (i.e., height or weight) [3–6]. These sex differences arise during puberty, when exposure to estrogen promotes subperiosteal bone deposition in boys  but has the opposite effect in girls, where estrogen suppresses subperiosteal bone deposition and promotes endosteal bone deposition [3,5,6]. As a result, girls tend to have bones with both a smaller external diameter and a smaller medullary cavity, and possibly greater density , than boys [4,5,9,10]. From this it can be inferred that: 1) a later onset of exposure to estrogen in girls provides more time for periosteal bone deposition; 2) the result of this delayed exposure includes long bones with an increased external diameter; 3) increased external diameter ultimately leads to bones with higher strength under bending and torsion. Thus, it is plausible that among girls with equal prepubertal bone size, those with later age at menarche (i.e., later exposure to elevated estrogen) would have larger adult bone diameter than their peers with an earlier age at menarche.
Examination of the relationship between age at menarche and bone strength is particularly timely given the reported secular trend for younger age at menarche [11–16]. If menarche affects subperiosteal bone deposition as posited, the age at which subperiosteal bone deposition rate decreases may be occurring earlier as well. The resulting bone geometry could lead to lower bone strength, tracking into adulthood, and a greater proportion of women at risk for fracture. While many researchers have examined the relationship between age at menarche and BMD and BMC [17–34], the specific effects of age at menarche on cross-sectional geometric properties of bone remain largely unexplored.
Participants included in the current analysis are derived from the Fels Longitudinal Study, a large longitudinal study of human growth, development and body composition change over the life span, which has been running continuously since 1929 . The sample consisted of 223 female participants between the ages of 7 and 35 years who had at least one hand-wrist radiograph taken between age 7 years and the time of menarche, and at least one radiograph taken between the ages of 18 and 35 years. Fels Longitudinal Study participants were born primarily in Southwest Ohio or neighboring states, and are predominantly Caucasian (~98%). Participants in the current analysis were born between 1929 and 1988, and the year of birth was used for assessment of secular trends. Participants were not selected for any disease or bone-related trait, and therefore represent normal variation in such traits. For the age range of interest (7–35 years), the examination protocol was to see participants every 6 months until age 18 or until they reach skeletal maturity, and every 2 years thereafter . Each study exam included collecting measures of body habitus (e.g., stature, weight, and BMI), a hand-wrist radiograph, and a detailed health history (including menstrual history) via questionnaire. A maximum of 197 participants were included in any one analysis, though the number varied depending on the age range used in determining premenarcheal bone strength (Table 1). All data collection protocols were approved by the Wright State University Institutional Review Board.
Quantitative measures of bone geometry were obtained from the second metacarpal of the left hand in each individual from antero-posterior radiographs. The second metacarpal is the most cylindrical bone in the hand and can be assumed to have bone evenly distributed about its axis [3,37]. As it is free from weight-bearing stress that would influence bone morphology, it has been widely used in biomedical investigations of bone health [3,37–44]. The radiographs were digitized using a Vidar DosimetryPRO Advantage Radiographic Scanner, and analyzed with a dedicated image processing program written for this purpose in MatLab v.18.104.22.1681 . The program first determined diaphyseal length and maximum length (diaphysis + epiphysis) of MC2 along the long axis of the bone. Then, using pixel properties reflecting radiolucency (the differences between light and dark pixels representing cortical bone and medullary cavity, respectively), the program determined external bone diameter and medial and lateral cortical thicknesses perpendicular to the long axis at midshaft of the maximum bone length (Fig. 1A). All measurements were recorded in millimeters. Each radiograph was reviewed for accuracy of measurement placement and adjusted, if necessary, before the data were entered into a database. The radiographs were reviewed by two assessors, and a subset was reviewed by both assessors to assess inter-observer error. Inter-observer repeatability, expressed as Lin's concordance correlation coefficient , ranged from 0.93 to 0.99, with the mean measurement differences between 0.008 mm and 0.023 mm (N=256). Intra-observer repeatability ranged between 0.95 and 0.99, with the mean differences between 0.001 mm and 0.04 mm (N1=261, N2=248).
Two indicators of bone strength were calculated for the study presented here; polar moment of inertia (J) and section modulus (Z). J represents torsional rigidity of bone along an axis perpendicular to the plane of the cross-section . In case of a hollow cylinder, J is equal to twice the value of the second moment of area (I), or bending rigidity, along an axis in the plane of the cross-section. Z represents the ratio of the polar moment of area and the radius of the bone, and is also a measurement of bending rigidity. Assuming a cylindrical cross-section, J and Z were calculated from measurements described above using the following formulae : J = π*(T4 − t4)/32, and Z = J/(T/2), where T is the external bone diameter and t is the diameter of the medullary cavity (external diameter minus total cortical thickness) at the midshaft of maximum bone length (Figs. 1B and and22).
We used multivariate regression to test for associations between J and Z values in adulthood and age at menarche, adjusting for premenarcheal J and Z values, premenarcheal body mass index (BMI), and birth year (to account for any secular trend in bone strength). We also considered age-squared at menarche, as well as interactions between age at menarche terms and childhood bone strength (model formulae are given in Tables 2 and and3).3). Nominally significant (p<0.10) covariates were retained in the models. For each model, an overall model F-test was conducted and R2 was calculated. In all the regression models, bone strength (adult and childhood J and Z) was log-transformed in order to normalize the adult response variables and reduce the influence of outliers in the childhood predictor variables. Additionally, all covariates were centered at their mean value (or mean of log values for childhood J and Z). The analyses were carried out using PROC GLM in SAS software v9.2 .
We considered two approaches to modeling the relationship between adulthood and childhood bone strength. In both, an average of bone strength values between the ages of 18 and 35 years was used as the adult bone strength variable. The first approach we employed controls for childhood bone properties at a specific prepubertal age. This approach is prospective and easily translates to a clinical setting. To make our models applicable to a wide range of ages for timing of the premenarcheal bone strength measure, we considered five models. Each model controlled for childhood bone properties at a different premenarcheal age (the average of measurements in the age range 7 to <8 years, 8 to <9 years, 9 to <10 years, 10 to <11 years, or 11 to <12 years). Girls who had already experienced menarche were excluded from the models controlling for bone strength at subsequent ages.
The second approach controls for bone properties at a premenarcheal maturational stage, such as Tanner B2 (e.g., ), or at a certain time interval relative to menarche or another maturational marker, such as age at peak height velocity (PHV) . This approach is retrospective and can be less clinically applicable because it is limited to peripubertal girls, who will likely experience menarche in the very near future, thus reducing the window for potential intervention to improve bone strength. For comparison, and because the retrospective approach has been used in previous studies, we adapted this approach to our data to examine differences in outcome attributable to analytical method. Because Tanner staging was not available in this sample, we used an average of prepubertal bone strength values recorded between one and two years before menarche as the premenarcheal bone strength value. This interval was chosen as it is in the range of the Tanner B2 stage, and it includes the age at PHV .
Descriptive statistics of the prospective samples, including sample size, mean, standard deviation and range for observed age at menarche, BMI, and adult and childhood bone strength indicators, for each model are shown in Table 1. Table 1 also includes the mean, standard deviation and range for bone strength variables on the natural log scale.
Tables 2 and and33 present the estimated regression parameters for the significant covariates in the prospective models for adult polar moment of inertia (J) and section modulus (Z), respectively, each along with model R2 values and overall p-values. Graphic representations of these models, showing the relationship between adult bone strength and age at menarche at different levels of childhood bone strength (while controlling for all significant covariates), are shown in Fig. 3. In each panel of Fig. 3, the lines represent the regression of adult bone strength on age at menarche at childhood bone strength values at the mean, and at the mean ± one standard deviation.
Results of the multivariate regression modeling of adult J showed that all of our prospective models are statistically significant (F-test p <0.001), and account for 55% to 70% of the variability in adult bone strength. The models predicting adult Z are all statistically significant as well (F-test p <0.001), and account for 59% to 65% of the variability in adult bone strength. Adult J and Z are both, in general, positively associated with age at menarche. One of the most substantial components of bone strength measured by J and Z is the external diameter of bone (T; see Bone strength measurements). Like J and Z, adult external bone diameter exhibits a positive relationship with age at menarche when controlling for childhood bone diameter (r2 for different models between 0.55 and 0.67, p<0.001). Given equal childhood bone diameter, later age at menarche provides an additional augmentative effect resulting in a larger adult bone diameter.
The positive effect of later age at menarche on J and Z is clearly present for individuals with average to high childhood bone strength, but the effect is less clear for individuals with relatively low childhood bone strength. Statistically, the attenuation of the age of menarche effect for those with lower childhood bone strength is due to the presence of a significant interaction with childhood bone strength (see the 9 and 11-year old models in Table 2 and Fig. 3). We suspect that, in fact, there either is an interaction or not regardless of at what age one measures prepubertal bone strength, and the presence of an interaction in some models but not others is due to either Type I or Type II error. In any case, all five models consistently indicate a significant positive association between age at menarche and adult bone strength for girls with average to high childhood bone strength in our sample. This association may also be true for girls on the lower end of childhood bone strength in our sample, but our results are not conclusive for those girls.
In our dataset, participant birth years span almost six decades, so it is possible that there is an underlying secular trend (either in age at menarche, BMI, or a previously unidentified secular trend in bone strength) influencing our results. In this sample of girls, birth year and age at menarche are not significantly correlated (p-values between 0.15 and 0.91 for different models). Additionally, birth year in our sample is not significantly correlated with childhood bone strength (p-values between 0.17 and 0.98 for different models) but is weakly correlated with childhood BMI in 7-year-olds (r=0.16, p=0.03) and 9-year-olds (r=0.14, p=0.05). Childhood BMI was significantly and negatively associated with adult bone strength (both J and Z) after controlling for the other covariates. After controlling for all other covariates, birth year is a significant (p = 0.011) and negative predictor of adult J only when controlling for childhood J at age 8 years, whereas for adult Z it is a significant negative predictor in the 8-year old model (p=0.035), and marginally significant in the 10-year old (p = 0.087) and 11-year old models (0.057). However, in all instances, its absolute contribution to adult bone strength is small (Tables 2 and and33).
Descriptive statistics for the retrospective sample, providing the mean±standard deviation (range) for each variable, and the natural log values for bone strength variables are shown in Table 4. Regression coefficients for the retrospective models for adult J and Z are shown in Table 5. These models accounted for 63.1% and 60.8% of the variability in adult J and Z, respectively. When controlling for average childhood J in the period between one to two years before menarche, adult J has a curvilinear relationship with age at menarche, with later age at menarche associated with, in general, lower adult J (see Fig. 4A). Birth year and BMI were significant negative predictors of adult J (see Table 5). Similar results were obtained for Z (see Table 5 and Fig. 4B).
Menarcheal timing has lasting effects on the female skeleton. Our principal finding using a prospective analysis shows that, for a majority of girls, a later age at menarche is associated with an increase in external bone diameter and a concomitant increase in bone strength under bending and torsion. A later age at menarche appears to be particularly beneficial for adult bone strength in girls of average and high childhood bone strength in our sample, regardless of age at which childhood bone strength is measured. For girls on the lower end of the childhood bone strength distribution in our sample, the effect of age at menarche is less clear (Fig. 3).
Previous work examining the effects of age at menarche on bone strength has generally focused on bone mineral content (BMC) and density (BMD), rather than bone geometry, with differing results [17–34]. Most studies report that late menarche is associated with lower BMD [17–22,24,28,29,32–34] and higher fracture risks for different skeletal sites [23,25,26,30,31] (but see [27,51]). While this may seem at odds with our results for bone strength, we believe the two outcomes are mutually compatible. If early exposure to estrogen at menarche halts periosteal expansion (and limits cross-sectional strength) in early puberty, but mineralization continues, the pattern of bone deposition may be altered [52,53]. Increased endosteal deposition after menarche compensates in part for the reduced periosteal expansion in girls. This configuration, however, does not benefit bending and torsional strength of the bone as much as does periosteal expansion. Thickening the cortex through endosteal expansion (if the external diameter cannot be increased) improves bone strength, but is less biomechanically beneficial than increasing the external diameter of bone (e.g., by delaying exposure to estrogen), even when cortical thickness is held constant (Fig. 2). The exponentiation of bone diameter values (both T and t) to the fourth power in calculations of polar moment of inertia and section modulus means that even small increases in the external diameter of bone due to a later age at menarche exponentially increase both measures of bone strength.
Endosteal expansion of bone after menarche is anatomically and physiologically limited by the necessary functions of the medullary cavity. It is also possible that new bone mineralization in postmenarcheal girls with early age at menarche could increase the BMD and BMC of bone. This may, at least in part, explain why some studies find skeletal benefits to a late menarche (increased BMD/BMC), and some studies show skeletal benefits to early menarche (bending/torsional strength).
Another explanation of divergent conclusions in the literature regarding the effect of age at menarche on bone strength is the differing approaches to statistical modeling used in each study. The primary statistical model presented here is prospective and, in controlling for premenarcheal values of the outcome, it utilizes bone strength measures at a specific prepubertal chronological age, regardless of age at menarche. Other studies have used a more retrospective approach, controlling for bone measures at a certain time point relative to a maturational marker, such as breast development or age at PHV.
To examine the differences between modeling approaches in our own data, we conducted analyses using the retrospective criteria. Using this retrospective approach, we found that age at menarche and bone strength were negatively associated (i.e., later age at menarche was associated with lower adult bone strength). This result is the opposite of what we found using the prospective approach. This discrepancy, however, stems from the fact that these two models statistically condition on childhood bone strength in two very different ways. The prospective model controls for childhood bone strength at a specific chronological age, whereas the retrospective model controls for childhood bone strength at a specific time point relative to a sexual maturation event, but does not include information on chronological age. Fig. 5 presents these two scenarios in two hypothetical individuals and serves to illustrate a possible reason why the prospective and retrospective approaches/models result in different conclusions regarding the relationship between age at menarche and adult bone strength.
In the prospective approach, the model is conditioned on childhood bone strength at a specific age, but subjects are not equally distant from menarche, as some girls will experience it earlier, and some later. Fig. 5A illustrates a hypothetical girl with earlier menarche (point of flexion in the figure) who started with the same premenarcheal bone strength at age 8 years as another hypothetical girl with later menarche. Assuming the rate of increase in bone strength diminishes at menarche (i.e., greater estrogen exposure), the girl with later menarche ends up with higher adult bone strength, ostensibly because she had more time before menarche for increasing bone diameter and bone strength.
In the retrospective approach, our model is conditioned on childhood bone strength at a specific number of years prior to menarche, but subjects are of different chronological ages. Fig. 5B presents one of the possible scenarios where early age at menarche could be seen to have a positive effect on cross-sectional bone strength in two hypothetical girls. Here, bone strength is measured (and is equal) at the same number of years before menarche, but each girl is measured at a different chronological age since age at menarche differs between them. In both girls, the event of menarche (point of flexion in the figure) slows the rate of increase in bone strength by the same degree. The girl with earlier menarche has higher adult bone strength, and the girl with later menarche has lower adult bone strength. However, in this model, the difference in adult bone strength is heavily based on a higher premenarcheal bone strength at any given chronological age in the girl with earlier menarche, and not necessarily because of the timing of menarche.
The prospective and retrospective modeling approaches are equally valid mathematically but, biologically, they seem to contradict one another. Age at menarche is either positively or negatively associated with adult bone strength, given childhood bone strength. The contradiction is resolved with a better understanding of what the models are conditioned on, how the pre-menarcheal age for bone strength measures is determined, and how that may influence the results. In the prospective model, when evaluating the effects of menarcheal timing on bone strength, girls are compared at the same age with similar bone strength and we prospectively determine which girl ends up with higher adult bone strength. This approach accounts for chronological age, bone strength at that age, and age at menarche. The retrospective model, however, controls for bone strength at a certain time point relative to the age at menarche, but does not account for chronological age. This leads to a comparison of girls of different chronological ages whose developmental timing is already known to be different, which makes them less comparable to each other than are the girls in the prospective model. Since childhood growth in normal children is largely a function of age, and to a lesser extent of sexual maturation, omitting chronological age when evaluating the effects of subsequent maturational events on childhood bone growth may introduce bias into the interpretation. Moreover, the retrospective approach may be less clinically applicable, as it applies only to peripubertal girls who would experience menarche in the near future, reducing the possible intervention window regarding geometric properties of bone. The prospective approach, on the other hand, can be applied to girls of any age regardless of their pubertal status, and provides a greater window of opportunity for potential intervention to increase cross-sectional bone strength given the predicted outcome.
Secular trends have been noted in a variety of maturational milestones in the U.S. and worldwide [11,12,35,54]. It is, therefore, possible that there is an underlying secular trend influencing our results, such as changes in age at menarche over the course of the 20th century [11–16], or a previously unrecognized secular trend in childhood or adult bone strength. Using a large sample of girls participating in the Fels Longitudinal Study, Demerath et al.  found a 3–6 month decrease in the mean age at menarche starting in the 1980s. In our smaller Fels dataset, birth year is not significantly correlated with age at menarche. Given our need for a matched measure of bone strength in adulthood, our sample did not include many girls who entered puberty in or after the 1980s, the timeframe in which Demerath et al. detected the trend. We nonetheless felt it was prudent to test for secular trend in our sample.
In our sample, in only a few models is birth year significantly associated with adult bone strength, and even then its contribution to the adult bone strength outcome is very small. Likewise, birth year is not significantly correlated with childhood bone strength, and is weakly correlated with childhood BMI. Birth year, when used as a predictor, also does not have a large impact on adult bone strength. Thus, there does not appear to be a consistent secular trend in bone strength during the 59-year birth year span captured in this sample.
In addition to the potential factors affecting the relationship between menarche and bone strength already discussed, the effects of BMI on bone strength and quality are another important consideration. In children, higher BMI generally has a positive effect on BMC and different measures of BMD [55–57], but the relationship between higher BMI and cortical thickness and density is less straightforward . Furthermore, higher adiposity generally has a negative impact on bone strength, including bone geometry [55,56,58,59]. In the current study, higher childhood BMI is associated with lower adult torsional and bending bone strength when controlling for prepubertal bone strength and age at menarche. Given the increasing proportion of overweight children in the United States, further studies focusing on the effects of BMI on various aspects of bone quality and strength are highly warranted.
Pubertal timing and bone strength are each under genetic influence [60–62]. Studies of genetic underpinnings of skeletal traits frequently show heritabilities of 50–65% [37,63,64], with the other proportion of trait variance due to the environment. In this sense, the term “environment” typically encompasses any number of non-genetic unmeasured covariates, such as diet or nutrition, physical activity, or hormonal influences. Because menarche directly contributes to one of these unmeasured environmental influences, namely hormones (e.g., estrogen), it is difficult to identify independent genetic influence on the two, versus shared genetic underpinnings (pleiotropy). Nevertheless, covariation between these two traits, whether it is genetic in origin or not, is an important consideration. This is especially true given the secular trend in menarche (discussed above) and the implications for lifelong skeletal health.
Our analyses demonstrate that later age at menarche may be beneficial for bone strength under bending and torsion. A likely mechanism for this is the delayed exposure to estrogen and longer period available for subperiosteal bone expansion in girls with later menarche. Therefore, although late menarche may result in lower BMC or BMD, at the same time it leads to bones of greater diameter, and consequently greater bone strength. Furthermore, based on our results, the strength of bone in childhood is a key predictor of adult bone strength, as measured by its cross-sectional geometric properties. Thus, any factors that negatively affect subperiosteal bone deposition in childhood (including early menarche), leading to bones of a smaller external diameter, have consequences into adulthood.
We thank the editor and two reviewers for their feedback. We also thank Drs. Siervogel and Czerwinski for their input on study design in the early stages of this project. Kimberly Lever and Sharon Lawrence collected the cortical data, which were compiled by Joseph Wagner, under the supervision of Dr. Duren. Anthropometric data were collected by a number of data collection staff of the Fels Longitudinal Study since 1929. We also thank the IBMS for the Alice L. Jee Young Investigator Award and the ASBMR for the Young Investigator Travel Award to Dr. Šešelj. Finally, we thank all the Fels Longitudinal Study participants. This work was supported by the National Institute of Health grants NIH R01HD056247, R01HD056247-S1, and R01HD012252.
Role of funding source: The funding source (National Institutes of Health) did not have any involvement in study design, data collection, analysis or interpretation, in writing of the report or the decision to submit the article for publication.
Conflict of interest: All authors state that they have no conflicts of interest.