Abstract
Background: Higher socioeconomic position (SEP) has been associated with increased risk of breast cancer. Its relationship with earlier age of pubertal onset, a risk factor for breast cancer, is less clear.
Methods: We studied the relationship of SEP to pubertal onset in a multiethnic cohort of 1,237 girls ages 6 to 8 years at baseline. Girls in three U.S. cities were followed for 5 to 8 years with annual clinical examinations from 2004 to 2012. SEP measures were examined for associations with pubertal onset, assessed by breast budding (thelarche) and pubic hair development (adrenarche). Analyses were conducted with accelerated failure time models using a Weibull distribution, with left, right, and interval censoring.
Results: Higher body mass index percentage at entry to the study and black or Hispanic race/ethnicity were the strongest predictors of age at pubertal onset. An SEP index comprising household family income, mother's education, and home ownership was an independent predictor of thelarche in adjusted models for all girls together and for white and Latina, separately, but not black girls, and the relationship varied by study site. The SEP index was not related to adrenarche in adjusted models. Overall, girls from the lowest quintile of SEP entered puberty on average 6% earlier than girls from the highest quintile (time ratio = 0.94; 95% confidence interval 0.91–0.97) in adjusted models.
Conclusions: Our results suggest that early-life SEP may influence the timing of pubertal development.
Impact: Factors related to lower SEP in childhood can adversely affect early development in ways that may increase the risk of breast cancer. Cancer Epidemiol Biomarkers Prev; 26(12); 1714–21. ©2017 AACR.
Introduction
Socioeconomic position (SEP) at both an individual and population level has been repeatedly associated with increased breast cancer incidence in high-income countries (1, 2). The strongest explanation of this association is that women with higher SEP tend to have reproductive patterns and practices that are directly related to breast cancer incidence such as younger age at menarche (3, 4), older age at menopause (3, 5), older age at first full-term pregnancy (4, 5), lower parity (4, 6), and shorter duration of breast feeding (4, 6). Other breast cancer risk factors, such as use of hormone therapy and higher average alcohol intake, have also been linked to higher SEP (5, 7, 8). Among the reproductive risk factors, a younger age at menarche has been associated with a higher risk of breast cancer in multiple types of epidemiologic studies in many countries around the world (5). Earlier menarche is associated with early onset of ovulatory cycles with increased hormone exposures over a lifetime and has a long-lasting influence on breast cancer risk (9, 10). A meta-analysis of pooled data from studies of breast cancer and reproductive risk factors in Europe and North American estimated that breast cancer risk increases by 5% for each younger year of age at menarche (3).
Long-term trends in the average age of menarche from European countries over the past 100 years and from Korea in the past 50 years have shown a progressive drop in age of up to 5 years (11, 12). Causes for this striking downward trend are thought to be related to improved social circumstances, including better nutrition (13). Other possible influences on the downward trend of the age of menarche include environmental chemicals that disrupt hormone pathways (e.g., endocrine disrupting chemicals), lower physical activity, and psychosocial stressors (14–16).
Studies that examined SEP and age at menarche are inconsistent, with some finding lower SEP directly related to earlier menarche (17–20) and others the inverse (21–23). This inconsistency may be at least in part due to the changing relationship of SEP with body size over time (24, 25) and the association of obesity with race and ethnicity (26). The age of menarche is consistently related to increased subcutaneous fat and body mass in the prepubertal period (27), but the trends toward increased overweight and obesity in children (28–30) have been primarily among those of lower, rather than higher, SEP in more recent years (26, 28, 31), emphasizing the importance of considering the effects of SEP, race/ethnicity, and obesity together (32).
Although the age of menarche dropped substantially in the past century, the decrease has slowed. The age of pubertal onset, on the other hand, seems still to be decreasing (29, 33, 34). The reasons for the apparent drop in the age of pubertal onset have been a topic of intense concern and investigation (32, 33, 35), given earlier onset of female reproductive maturity is associated not only with increased breast cancer rates in adulthood, but also with more immediate negative consequences in adolescence, including sexual abuse, sexually transmitted diseases, and mental health issues such as depression (34, 36). Breast budding, or thelarche, is driven by activation of the hypothalamic–pituitary–ovarian axis and mediated by estrogen and is best measured by both observation and palpation, whereas the appearance of pubic hair and other secondary sexual characteristics, or adrenarche, is driven by androgen secretion from the adrenal cortex and assessed visually (29). In the current prospective cohort study, we sought to understand the relationship of SEP to pubertal onset in girls, which tends to precede menarche by about 2 years. The duration of time between pubertal onset and menarche is variable (33). The determinants of pubertal onset may be different from those of menarche and a more sensitive indicator of environmental and hormonal influences on pubertal maturation (33, 34). Menarche has been used in most epidemiologic studies, especially case–control studies, as it is an event more easily recalled by adult women. A subsequent report in our ongoing analyses will cover the relationship of SEP to menarche.
Materials and Methods
The purpose and study design of the Breast Cancer and the Environment Research Program (BCERP) have been described previously (14), and details of results from the BCERP Puberty Study have been published (34, 37–39). Briefly, the BCERP Puberty Study is a consortium of three collaborative prospective studies examining predictors of the onset of puberty in girls. Between 2004 and 2007, 1,239 socioeconomically and race/ethnically diverse girls, aged 6 to 8 years, were enrolled from three locations: the Greater San Francisco Bay Area, the Greater Cincinnati Area, and East Harlem in New York City, hereafter referred to as California, Ohio, and New York. The parent, legal guardian, or primary caregiver gave informed consent, and assent was obtained from the girl. Institutional review boards at each participating institution [Kaiser Permanente (Oakland, CA), Cincinnati Children's Hospital (Cincinnati, OH), Mount Sinai School of Medicine (New York, NY), and UCSF (San Francisco, CA)] approved the study protocols and procedures. Inclusion criteria were age (6–8 years), female sex, no underlying endocrine medical conditions, and in New York, black, or Hispanic race/ethnicity. The parents of 2 girls requested that data from their daughters not be included in the study. We followed 1,237 girls to the time of pubertal onset and included those with observed outcomes in this analysis.
Measurements
Data were obtained from questionnaires completed annually by a parent, legal guardian, or caregiver, either by in-person interviews (California and New York) or by self-administration for the first 5 years and then by interview (Ohio) for the 5- to 8-year duration of follow-up for this analysis, through March 2012. Interviews were conducted in English or Spanish according to the parent/guardian's preference. Anthropometry and assessment of pubertal signs were performed at annual visits in California and New York and semiannual visits in Ohio by clinical research assistants, nurse practitioners, or physicians trained with a standard protocol developed by expert pediatricians across the three sites (37).
Determination of SEP
We examined the following measures of SEP as well as parent/primary caregiver–reported race/ethnicity of the girl from the baseline questionnaire: education of the mother, household income, occupation of the primary financial provider, home ownership, and female as the head of household (defined as financial support only by one or more adult females). Household income combined the incomes of all wage earners and was recorded in categories as <$25,000, $25–$50,000, $50–$100,000, and >$100,000/year. Occupation of the primary financial provider was coded according to status categories as professional, nonmanual, and manual occupations for analyses (40). Mother's education was categorized as high school or less, some college or vocational school, bachelor's degree, and master's degree or higher. We constructed an SEP index that included household income, mother's education, and home ownership by standardizing each variable to mean zero and SD one and summing the standardized variables. Occupation of the primary financial provider had a substantial number of missing values and did not improve reliability of the SEP index (as measured by Cronbach alpha) and was not included.
Other independent variables
The girl's body mass index [weight (kg)/height (m)2] percentile (BMI%) was based on age- and sex-specific growth charts from the Centers for Disease Control and Prevention for 2000 (41) at baseline (categorized as <50, 50–85, and ≥85th percentiles). BMI percentile was based on the baseline height and weight measurement for each girl using calibrated scales and stadiometers by research staff that had been trained and certified uniformly across all three sites. Race/ethnicity was categorized as black, white, Hispanic or Latino, and Asian American in hierarchical order following an algorithm that made each race/ethnic group mutually exclusive. Other variables examined included mother's place of birth (United States/Canada, Mexico, other Latin America/Puerto Rico/Virgin Islands, rest of the world), the study site, and mother's age at menarche (<12, 12–13, ≥14), which may reflect a genetic influence on pubertal development.
Pubertal onset
Girls' pubertal development was assessed using standard methods of Tanner staging conducted by trained staff at each in-person clinic visit using an established 5-stage classification scheme for describing the onset and progression of breast and pubic hair changes by inspection (29) and palpation (42). Details of training, certification, and assessment procedures are reported elsewhere (37). The outcomes were onset of signs of puberty as assessed both by observation and palpation of breast budding for stages B2 or higher (breast) and by observation of stages PH2 or higher (pubic hair).
Statistical analysis
Our a priori hypothesis about SEP and pubertal onset was that girls from lower SEP families would go through puberty earlier, and that relationship would be stronger among overweight or obese (BMI% > 85) girls. Race/ethnicity, which is known to be strongly related to pubertal onset independent of BMI (37, 43), was seen as a potential confounder in this relationship being related both to SEP and to pubertal onset. To explore this hypothesis, we first compared the characteristics of the girls in terms of the independent variables by site using χ2 tests. We then computed a correlation matrix of multiple available SEP measures and assessed reliability using Cronbach alpha. The relationship of the SEP index with BMI was evaluated using polytomous logistic regression to model BMI percentile above the 85th percentile and between the 50th and 85th percentiles versus <50th percentile as a function of SEP quintile, first adjusted only for age at BMI measurement followed by adjustment for race/ethnicity. The relationship of the SEP index and other variables with age at onset of B2 and PH2 was then evaluated in unadjusted accelerated failure time models using a Weibull distribution (SAS PROC LIFEREG, SAS v. 9.3). Left and right censoring was used to account for pubertal transitions taking place outside the period of observation, and interval censoring accounted for pubertal transitions between examination visits. The number in each race/ethnic group was sufficient to allow estimation of interactions between race/ethnicity and other variables, and to produce stable estimates for black, white, and Hispanic girls. Subsequent models of the association between the SEP index and pubertal onset adjusted for BMI%, race/ethnicity, and interactions between BMI% and race/ethnicity, which were included to allow the association between pubertal onset and race/ethnicity to vary by BMI%. The effects of BMI%, race/ethnicity, and their interaction were modeled as the effect of BMI% for whites and the effect of race/ethnicity by BMI% category. We also estimated race/ethnicity–specific SEP effects in a model of the age at B2 to the SEP index, including a race/ethnicity–SEP interaction, and finally a site-specific model of age at B2 that included site–SEP, site–BMI%, and site–race/ethnicity interactions to estimate site-specific associations. Additional models were created to estimate trends across SEP levels.
For each model, time ratios (TRs) and estimated median age at onset for girls in the reference category of all variables in the model were computed, along with their 95% confidence intervals (CIs). The TR indicates how much earlier or later the estimated onset of puberty occurred relative to the reference category of a particular variable. For instance, if the median age at onset was 10 years in the reference category for a particular model, a TR of 0.95 indicates a 5% earlier onset, which corresponds to 6 months (i.e., 0.5 year). In a multivariable model, the number of months difference indicated by a TR depends on the levels of all independent variables in the model.
Results
All but 2 girls (N = 1,235) contributed to the assessment of the interval for the pubertal transition for breast development; 14% were left censored (B2 at baseline), 72% interval censored, and 13% right censored (still B1 when last observed). A pubertal transition interval for pubic hair development could be determined for 1,230 girls, with 12% left censored, 71% interval censored, and 17% right censored.
The study sample had a high degree of variability in SEP measures and covariates across sites (Table 1). All measures of SEP were highly skewed by site with lower SEP girls in New York and higher SEP level girls in the other two sites. Only black and Hispanic participants recruited in New York by design. Mothers of girls in New York were more likely to be first-generation immigrants from Mexico and Latin America, and more mothers of girls in New York had gone through menarche at ages less than 12 years than the other two sites. New York girls were also more likely to be obese (BMI% ≥ 85th percentile) at baseline compared with girls in other sites. Missing data ranged from 0.2% (BMI%) to 11.6% (occupation) of participants. All measures in the SEP indices were strongly correlated with each other except for occupational status.
In models adjusted only for age, BMI% above the 85th percentile was associated with lower SEP [quintile 1 (Q1) vs. Q5: OR = 2.6; 95% CI, 1.6–4.3; Q2 vs. Q5: OR = 2.6; 95% CI, 1.6–4.3; Q3 vs. Q5: OR = 1.8; 95% CI, 1.1–3.0; Q4 vs. Q5: OR = 1.8; 95% CI, 1.1–3.0]. BMI% above the 85th percentile was more common among black (OR = 2.1; 95% CI, 1.3–3.3) and Hispanic girls (OR = 2.9; 95% CI, 1.8–4.8) and less common among Asian (OR = 0.3; 95% CI, 0.1–0.9) girls compared with whites.
In unadjusted models, lower SEP index or any measure, except occupation, predicted earlier onset of B2 and PH2 by a substantial amount (Table 2). For example, a TR of 0.95 for girls with a household income of <$25K is equivalent to 5.8 months earlier onset of B2 compared with girls with household incomes >$100,000. For the SEP index, a TR of 0.94 for the lowest quintile of the SEP index is equivalent to 7.0 months earlier onset of B2 compared with girls in the highest quintile of SEP index. The unadjusted association of the SEP index with age at B2 and at PH2 was statistically significant (Ptrend < 0.0001 for both). In adjusted models, the SEP index association was diminished only slightly after adjustment for race/ethnicity with race/ethnic–specific BMI% (interaction) effects; however, the association of SEP index with PH2 was not longer present (Table 3).
To further examine the impact of SEP on B2 by race/ethnicity, we estimated SEP effects separately among white, Hispanic, black, and Asian American girls (Table 4). After adjustment for BMI% and the interaction BMI% with race/ethnicity and with SEP, pubertal onset in Hispanic girls had the strongest association with SEP, although the trend for whites was also significant. No relationship was seen for black girls, and there were too few Asian girls to produce stable results. Finally, in stratified models, we found variation among the sites in the association of SEP with age at B2 adjusted for BMI% and race/ethnicity, with SEP significantly associated with age at B2 only in San Francisco (Table 5).
Discussion
To our knowledge, this is the first study to prospectively examine the influence of SEP on the onset of puberty determined by physical examinations in a multiethnic population of U.S. girls. In previously published analyses of this cohort of girls, we reported the median age of onset of breast budding (Tanner breast stage 2) to be at age 8.8, 9.3, 9.7, and 9.7 years for black, Latina, white, and Asian girls, respectively (34), compared with 8.9 and 10.0 years for black and white girls in the Pediatric Research in Office Settings (PROS; ref. 43), and BMI% predicted earlier puberty more strongly in the current study in the PROS performed 10 to 20 years earlier (43). In this analysis, we were interested in examining the interrelationship of social position with the age of pubertal onset, which is viewed as a window of susceptibility in the life course perspective of breast cancer etiology (44, 45). We found that a prepubertal girl's SEP as measured by multiple variables, but especially by a lower SEP index, which was comprised of household income, mother's education, and home ownership, was associated with earlier age at the onset of pubertal signs of breast but not pubic hair development in models adjusted for BMI%, race/ethnicity, and their interaction. Specifically, girls in the lowest SEP index quintile developed pubertal signs of breast budding a full 7 months earlier than girls in the highest SEP quintile. Adjustment for BMI%, with or without race/ethnicity and their interaction, somewhat attenuated but did not substantially change the relationship of SEP with the onset of breast development compared with the unadjusted model. We did not find a statistically significant relationship of the SEP index to adrenarche as assessed by pubic hair development. This relationship deserves further exploration, but it may be that influences on pubertal development associated with social disadvantage act more strongly on the hypothalamic–pituitary–ovarian axis than that on the development of the adrenal cortex.
Our a priori concept was that the effect of SEP on age at puberty worked through, or was mediated by, its association with obesity. Whereas the decreasing age of menarche internationally has long been associated with higher standards of living and better nutrition, many countries are now experiencing an epidemic of pediatric obesity, which is associated with earlier age at puberty and is more prevalent in lower, not higher, SEP groups (28, 31). In our study population, although lower SEP was associated with earlier onset of B2, it was relatively independent of obesity (i.e., BMI%) and of race/ethnicity as well. This is consistent with the findings from the National Longitudinal Study of Youth (20) and suggests that there may be other pathways through which SEP influences early pubertal development. For example, one such pathway relates to psychosocial stress in early childhood (46). There is evidence that the absence of fathers is associated with earlier pubertal onset, at least in higher SEP girls (47, 48). Other factors may relate to exposures to environmental chemicals (38, 39) and to the built environment (49) associated with higher levels of pollution and to obesity, food insecurity, fewer playgrounds, and less opportunities for physical activity that were not examined in this analysis.
Consistent with the hypothesis, in addition to BMI, race/ethnicity was a strong predictor of the age of pubertal onset, with black girls entering puberty substantially before girls in other groups. In a model with race/ethnicity–specific SEP effects, a lower SEP index was related to earlier onset of puberty measured by breast budding in Hispanic and white girls, but not in black girls. For black girls, BMI% proved a strong predictor of pubertal onset, whereas SEP did not. The absence of a relationship between SEP and pubertal onset among black girls has also been observed by Braithwaite and colleagues (18) in the National Growth and Health Study and by Krieger and colleagues examining long terms trends in the National Health and Nutritional Examination Survey (32), although it was not seen in the National Longitudinal Study of Youth (20).
Our results are consistent with one other recent prospective study of SEP and puberty in a large Australian cohort of both girls and boys that relied on parental self-report of pubertal onset (50). In 1,770 girls assessed at age 10 to 11 years, the rate of early puberty was increased 2-fold (OR = 1.96; 95% CI, 1.08–3.56) for girls with low household SEP. In other studies of pubertal onset as measured by menarche, our results are compatible with Windham and colleagues (17), Braithwaite and colleagues (18), James-Todd and colleagues (19), and Deardorff and colleagues (20), where lower SEP was associated with an earlier age at menarche. For example, in the National Growth and Health Study, Braithwaite and colleagues found that higher SEP (measured by household income) white girls were more likely to go through menarche later than lower SEP white girls. In contrast, black girls of higher SEP went through menarche earlier than lower SEP black girls.
Our study has the advantage of a prospective design and repeated direct examination of girls as they entered the pubertal transition. Most of the literature relates to SEP and age of menarche obtained retrospectively (3, 4), as menarche is more easily recalled by adult women in epidemiologic studies than pubertal onset. Inferences in this study were strengthened by its longitudinal design and annual (or semiannual in Ohio) data collection, including physical examination with breast palpation. Breast palpation is preferred as it reduces errors in assessment from observation alone where there may appear to be breast development in overweight and obese girls (51). The study participants were also race/ethnically and socioeconomically diverse and more than 80% of the girls showed signs of pubertal development before the end of the data collection period or loss to follow-up.
Limitations of the study include fairly long (usually annual) intervals between observations and left censoring of a proportion of pubertal transitions. Cincinnati girls were examined every 6 months, and thus, the dating of pubertal onset was more precise at this site. That the other two sites examined girls annually could be considered a limitation, but it is not clear how serious an effect this was as mean pubertal onset in Cincinnati was estimated at a time intermediate between the other two sites (Table 3). It should be noted that we assessed the onset of signs of puberty in this study, but are aware that these signs, breast budding and pubic hair growth, may not be reflections of true pubertal onset. Rather, they may reflect changes in body phenotype due to external environmental exposures, such as those potentially brought on by endocrine-disrupting chemicals that may also vary by SEP (52, 53). Also, these results may not be generalizable to the entire U.S. population, even though we used data from an integrated cohort in three national geographic locations; participants were selected primarily from urban areas and were a selected sample in the sense that they agreed to participate for multiple visits over many years.
The relevance of these findings to breast cancer incidence in adulthood is a matter of speculation, as previously conducted longitudinal studies on the relationship are inconsistent. Higher childhood SEP, as measured by either father's occupation or education, was not related to breast cancer incidence in a large Dutch study (54), but, as measured by higher early family income, it was associated with greater risk of breast cancer in the Wisconsin Longitudinal Study (55). In the future, the increasing prevalence of pediatric obesity in the United States (and other industrialized countries; ref. 56) may play a role in the relationship of SEP and breast cancer (57). Higher obesity may increase breast cancer risk by reducing the age of puberty, and thus lifelong exposure to estrogen. This hypothesis is not supported by data from the Nurse's Health Study where recalled childhood obesity was associated with lower, not higher, risk of breast cancer (58), but that was a cohort of women who grew up in the early part of the last century, and the assessment of obesity is subject to recall bias. If current shifts in the association of SEP with pubertal onset and menarche reverse the traditional relationship such that higher childhood SEP is associated with lower rates of breast cancer, the currently understood pattern of reproductive risk factors for breast cancer could be altered in future decades.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Disclaimer
The contents of this article are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS or NCI, the NIH, the Centers for Disease Control and Prevention, or the California Department of Public Health.
Authors' Contributions
Conception and design: R.A. Hiatt, S.L. Stewart, L.H. Kushi, G.C. Windham, F.M. Biro, S.M. Pinney, M.S. Wolff, S.L. Teitelbaum, D. Braithwaite
Development of methodology: R.A. Hiatt, L.H. Kushi, G.C. Windham, F.M. Biro, S.M. Pinney, M.S. Wolff, S.L. Teitelbaum, D. Braithwaite
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K.S. Hoeft, L.H. Kushi, G.C. Windham, F.M. Biro, S.M. Pinney, M.S. Wolff, S.L. Teitelbaum
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): R.A. Hiatt, S.L. Stewart, L.H. Kushi, S.M. Pinney, M.S. Wolff, D. Braithwaite
Writing, review, and/or revision of the manuscript: R.A. Hiatt, S.L. Stewart, K.S. Hoeft, L.H. Kushi, G.C. Windham, F.M. Biro, S.M. Pinney, M.S. Wolff, S.L. Teitelbaum, D. Braithwaite
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): K.S. Hoeft
Study supervision: R.A. Hiatt, F.M. Biro, S.L. Teitelbaum
Other (obtained funding for the project): D. Braithwaite
Acknowledgments
This work was supported by a NCIRO3 CA143936-01 awarded to D. Braithwaite and made possible by the Breast Cancer and the Environment Research Program award numbers U01ES012770 (to F.M. Biro and S.M. Pinney), U01ES012771 (to M.S. Wolff and S.L. Teitelbaum), U01ES012800 and U01ES012801 (to R.A. Hiatt, S.L. Stewart, L.H. Kushi, and G.C. Windham), U01ES019435 (to L.H. Kushi and G.C. Windham), U01ES019453 (to F.M. Biro and S.M. Pinney), U01ES019454 (to M.S. Wolff and S.L. Teitelbaum), and U01ES019457 (to R.A. Hiatt and S.L. Stewart) from the National Institute of Environmental Health Sciences (NIEHS) and the NCI; P01ES009584 and P30ES006096 from NIEHS; UL1RR024131, UL1RR029887, and UL1RR026314 from the National Center for Research Resources; and F31 DE023282 from NIDCR (to K. S. Hoeft).
We gratefully acknowledge support of the Avon Foundation for community outreach and education.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.