Growth reference data presented here add to an expanding literature in this area, with the potential to benefit patient care and enhance research into qualities of growth among infants, children, and teens with achondroplasia. The large number of data points amassed enhance growth curves that our research group previously established based on a single data source [4–6]. We also explored the possibility of secular trends in weight in children by age group and decade of birth and utilized weight-for-age data to screen for extremes in weight that were excluded from our novel, prescriptive weight-for-age and weight-for-length/height curves. We expect these figures to have clinical utility for identifying children carrying excess weight for their stature. Additionally, we add to the literature by offering new head circumference curves for young children from the US. We expect the utility of this large achondroplasia anthropometric database in clinical care and research to be multi-fold.
As supported by the comparison of height trajectories to that of average stature children, the achondroplasia-specific FGFR3 mutation influence on linear growth is, of course, greater than any environmental or ethnic factor that could be derived from other studies of linear growth in achondroplasia. By comparing to other studies of linear growth in achondroplasia (e.g. from Europe, Argentina, Australia, Japan) [8, 10, 14, 16], more subtle environmental or ethnic influences on attained height for a given age across childhood could potentially be ascertained. Based on visual inspection of the most recently published height-for-age curves, median height was similar across populations at 5 and 10 years of age, ranging from ~ 85–87 cm in girls and 86–88 cm in boys at 5 years and 104–106 cm in girls and 106–108 cm in boys at 10 years. At 15 years, median height in Argentinian girls was at least 3 cm shorter than peers from the US, Australia, and Europe (116 vs 119–121 cm), as previously noted [14]. However, at 18 years, the oldest age for which data were available for all settings, both the US and Argentinian girls, at a median height of ~ 120 cm, were shorter than their European peers (124 cm). For boys, median heights were quite similar across settings at 15 years (123–125 cm), but variable at 18 years (127–134 cm; shortest in the US, tallest in Australia).
Greater variability in the older age groups could reflect underlying genetic differences in growth potential, environmental influences on growth, different timing of pubertal development (although any pubertal gains in height seem subtle), a relative paucity of data and selection bias at older ages, or different modelling approaches that contribute to variability in estimates of median height. All of these possibilities need to be studied in a more rigorous manner. For average stature children, a global growth standard is used to characterize length/height in children 0–5 years [22], in whom growth occurs quite comparably regardless of ethnic background when unrestricted by environmental constraints [23]. A reference based largely on US-derived data is also used globally to characterize heights of average stature children 5–20 years of age [24], although the assumption of comparability of growth trajectories through childhood and adolescence by background is not certain [25]. Ultimately, we may find that compiling height-for-age across populations could be used to generate a global standard for achondroplasia.
Additionally, whether puberty results in a linear growth spurt, and whether its timing, tempo, or magnitude varies by ethnicity or other environmental factors remains an area requiring resolution. The cross-sectional presentation of height-for-age data does not suggest a substantial pubertal growth spurt in achondroplasia, although others have speculated that one exists based on the application of height data collected serially in a small number of adolescents with achondroplasia to the Preece-Baines model [13]. That conclusion may be an artifact of the modeling approach, and if a pubertal spurt in linear growth does occur it is very modest. A harmonized effort across populations to collect longitudinal Tanner stage data in conjunction with height velocity would elucidate effects of puberty on height gains in achondroplasia, but this question is not resolvable in this analysis of retrospective growth data.
Accumulation of excess weight is a concern in the population with achondroplasia [5, 6]. The growth curves presented here, as shown previously [6], show more overlap in weight than height with the average stature population–findings consistent with a body morphology with a high ratio of trunk to extremity size. Additionally, however, excess adiposity is a concern in this population, with few tools available to guide what optimal weights should be. Given the known obesity epidemic in the US it is useful to compare weight-for-age across populations to ascertain whether the US population with achondroplasia is heavier than their global counterparts. In fact, also based on visual inspection of most recent growth curves [8, 10, 13, 14], median weights were similar across settings until 10 years of age in girls (14–16 kg at 5 years; 24–26 kg at 10 years) and boys (15–16 kg at 5 years; 24–26 kg at 10 years), while US girls from this study were heavier at 15 years (43 kg vs 36–40 kg) and 17 years (last age in common across settings; 47 kg in US vs 41–44 kg from other sites). Among boys, US and European boys were comparable at 15 years (42 kg vs 34 kg in Australia and 37 kg in Argentina) and 17 years (47 kg in the US, 48 kg in Europe vs 38 kg in Australia and 42 kg in Argentina). Although findings are subject to the same caveats described above that could account for variability in height-for-age across settings, they are also consistent with children with achondroplasia being heavier-set in the US, particularly for their stature, than elsewhere, with the possible exception of Europe. This is despite our efforts to screen out extremes in weight. We were surprised that secular trends in excess weight by decade of birth in our population were not found, but did note a large number of extreme weights for a given height. In order to avoid skewing reference data to higher than optimal weights, potentially giving patients the impression that a high weight is acceptable, we utilized the weight-for-age curves generated here to remove cases with extreme weights before generating novel weight-for-stature charts that we propose for clinical use in this population.
We propose new weight-for-stature charts to help clinicians integrate information on weight and height of patients to improve guidance for achieving appropriate weights in a patient population prone to obesity. We [6] and others [8, 10] had also established BMI-for-age curves, which integrate information on weight and height by age. However, BMI can be difficult to interpret, and it can be tempting to use inappropriate cutoffs that have been perpetuated in the literature for average stature individuals to misattribute overweight and obesity to patients with achondroplasia [26]. Additional studies are needed to quantify body composition in individuals with achondroplasia and correlate these values with BMI.
Sensitivity analysis of the head circumference data indicated there was no difference in the curves generated from all available data and those derived from data excluding those undergoing foramen magnum/C1/C2 decompression and/or ventricular shunting. Therefore, we opted to include all available head circumference data to derive these new curves. These curves may now be utilized in the clinic and research venues to ascertain deviation from normal cranial growth which should prompt further investigation.
Limitations of the study include the fact that data were not collected at pre-determined time intervals, due to the clinic-based and retrospective nature of the study. Regrettably, we do not have sufficient Tanner staging available in this cohort to compare landmarks of pubertal development against growth. Clearly this should be ascertained in future cohort studies.
Strengths of this study include the fact that in this large, well-phenotyped cohort of patients with achondroplasia, we could eliminate the potential confounding effect of surgical and medical interventions influencing linear growth by omitting the stature data collected after these procedures. In this regard, the curves represent a norm against which individuals and populations who undergo novel medical therapies and surgical interventions to ascertain effects of those treatments can be compared. Although these data were collected by multiple providers and therefore potentially subject to inconsistencies, these growth parameters were all collected in skeletal dysplasia clinics or affiliated clinical site where care providers are accustomed to performing anthropometry. Additionally, the large overall amount of anthropometric data and the fact that over 50% of the cohort contributed 5 or more data points for all 3 parameters means it is more longitudinal than cross-sectional in nature and therefore more reflective of longitudinal growth. Finally, the methodology employed to generate our curves by modeling separate splines over the age range accounted for the fewer data points available at older ages.
To the last point, we offer a note regarding the choice to use penalized cubic p-splines to estimate age (and height with the weight for height analyses) percentiles, means and standard deviation as compared to LMS (Box-Cox) employed by others. Our approach requires no assumptions about the age specific anthropometric measurement distributions. The LMS method uses a normalizing power transformation of age-specific anthropometric measures which requires the estimation of three age-specific parameters: L (lambda), the transforming power, M (mu), the age specific median, and S (sigma) the age specific standard deviation. The resulting percentile estimates based on this approach are sensitive to these estimates. The LMS approach is appropriate when there are a relatively large number of measurements in each 1-month age interval, allowing for estimates of these 3 parameters at each age. However, as noted in our data, the number of measures per age-interval is variable and decreases substantially with increasing age. Thus, our smoothing approach was appropriate given the characteristics of available data.