Characterization of pubertal development of girls in rural Bangladesh

Abstract

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing * E-mail: [email protected] Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America ORCID logo https://orcid.org/0000-0001-9458-362X Introduction Adolescence is a life stage during which individuals transition from childhood to adulthood, experiencing physical, physiological and psychological changes. In a birth cohort of >15,000 girls aged 9 to 15 years at enrollment living in a large, typical, population-dense, rural setting of Bangladesh, data collected between 2015 and 2017 were used to 1) describe the timing, sequence and progression of prominent features of female sexual maturation, namely menarche, breast development and pubic hair growth, and 2) relate observed patterns of development to those published from other adolescent populations around the world. Gaibandha is an area that reflects typical rural characteristics of Bangladesh [23], a region historically beset by low socioeconomic status [24], children being born low birth weight [25], high infant mortality [26] and chronic undernutrition among preschool children [27] and women in the reproductive years with underlying high levels of food insecurity [28, 29]. Pubertal benchmarks were assessed using questions that were adapted from the Pubertal Development Scale (PDS) [30] to the age range and levels of maturation expected of participants, and which included menarche plus secondary sex characteristics of breast development and pubic hair growth for which Tanner stages [31] exist.

About the Authors:

Jinhee Hur

Roles Formal analysis, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

ORCID logo https://orcid.org/0000-0003-0968-208X

Kerry J. Schulze

Roles Conceptualization, Writing – review & editing

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

Andrew L. Thorne-Lyman

Roles Writing – review & editing

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

Lee S.-F. Wu

Roles Data curation, Resources, Software

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

Saijuddin Shaikh

Roles Conceptualization, Investigation, Supervision

Affiliation: The JiVitA Project, Gaibandha, Bangladesh

Hasmot Ali

Roles Conceptualization, Investigation, Supervision

Affiliation: The JiVitA Project, Gaibandha, Bangladesh

ORCID logo https://orcid.org/0000-0003-0044-1992

Alain B. Labrique

Roles Conceptualization, Supervision

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

Keith P. West Jr.

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliation: Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

ORCID logo https://orcid.org/0000-0001-9458-362X

Introduction

Adolescence is a life stage during which individuals transition from childhood to adulthood, experiencing physical, physiological and psychological changes. The most salient milestone characterizing adolescence is attainment of reproductive capability through pubertal maturation. South Asia is a region undergoing rapid changes including growing wealth, shifting diets and decreasing physical activity [1]. These changes may influence the timing of sexual maturation of populations in the region [2], as indicated by an estimate of the onset of the primary sex characteristic (i.e. menarche) in rural Bangladesh occurring at age 13.0 years in 2017 [3], compared to age 15.8 years in 1976 [4]. A small number of studies have reported the timing of secondary sex characteristics (e.g. breast development and pubic hair growth) in the region since the late 1980s [5–9], with initial stages of breast development occurring at age ~10 years and pubic hair growth in the months that follow. Given the potential for regional and temporal variation in sexual maturation [10–12] against trends in recent years of generally improving nutritional conditions in historically, chronically undernourished rural Bangladesh [1], there is utility to periodically monitor the timing, sequence and progression of pubertal milestones as an indicator of adolescent maturity and early reproductive health [3, 4, 13, 14].

The timing of puberty in rural South Asian cultures can markedly affect trajectories of socialization, empowerment and agency [15], with menarche commonly viewed to indicate a young girl’s entry into adulthood and readiness for marriage [16, 17]. Coupled with infrequent contraception among newlyweds [18], early marriage can lead to mid-adolescent pregnancy [19], which can adversely affect nutritional status, decelerate growth [20], and thus potentially contribute to endemic short stature among rural South Asian women [20], in addition to often restricting educational achievement of young women [16]. These sequences suggest there is important information to gain by characterizing the timing and progression of pubertal maturation as one approach to monitor health and inform policies intending to delay early marriage and pregnancy, improve nutrition, growth and development, and foster social development among young women in high-risk populations.

In a birth cohort of >15,000 girls aged 9 to 15 years at enrollment living in a large, typical, population-dense, rural setting of Bangladesh, data collected between 2015 and 2017 were used to 1) describe the timing, sequence and progression of prominent features of female sexual maturation, namely menarche, breast development and pubic hair growth, and 2) relate observed patterns of development to those published from other adolescent populations around the world.

Methods

Study population

The JiVitA-1 Cohort Follow-Up Study was a prospective assessment of a birth cohort of 35,056 young adolescents, comprising 17,755 boys and 17,301 girls aged 9 to 15 at enrollment, a relevant age range during which a series of pubertal events take place [21] (Fig 1). The children were born to mothers who participated in a double-masked, cluster-randomized trial of weekly, antenatal vitamin A or β-carotene supplementation, JiVitA-1, conducted in Gaibandha District in northwestern Bangladesh from August 2001 through January 2007 [22]. Gaibandha is an area that reflects typical rural characteristics of Bangladesh [23], a region historically beset by low socioeconomic status [24], children being born low birth weight [25], high infant mortality [26] and chronic undernutrition among preschool children [27] and women in the reproductive years with underlying high levels of food insecurity [28, 29]. Maternal characteristics of the current study population have been described previously [22].

[Figure omitted. See PDF.]

Fig 1. Flow diagram of the study participants.

https://doi.org/10.1371/journal.pone.0247762.g001

In March 2015, prior to the formal start of the adolescent follow-up study, an initial census was executed to update the vital and residential status of the entire birth cohort. Starting in September 2015, resident children were visited at home during their birth month by trained female interviewers to assess pubertal development, anthropometry, attained level of schooling and household socioeconomic status. With recorded birth dates, exact ages were known on the date of assessment. As the study was designed to assess children during their birth month, the distribution of ages was concentrated around each child’s known birth month across ages 9 to 15 years. The cohort study was designed to follow up children a year later, from ~10 to 16 years in order to assess progressions in growth, maturation, nutritional status and other factors.

In this analysis of female pubertal development, a total of 15,320 adolescent girls aged 9–15 years at the time of 1st interview served as the intended sample, of whom 1,263 were lost to follow-up by the time 2nd visit, leaving 14,057 girls for follow-up assessment by the same cadre of interviewers at a median (interquartile range, 10th and 90th percentiles) of 12.5 (12.0–13.3, 11.8 and 14.6) months later.

Pubertal development

At both baseline and follow-up visits, trained female staff interviewed girls in the presence of their mother or female guardian, providing an assuring interview setting and opportunity for parental guidance. Pubertal benchmarks were assessed using questions that were adapted from the Pubertal Development Scale (PDS) [30] to the age range and levels of maturation expected of participants, and which included menarche plus secondary sex characteristics of breast development and pubic hair growth for which Tanner stages [31] exist. Following an introduction during which interviewers explained the physical changes that take place during adolescence, girls were asked to report the extent of their breast development and growth of pubic hair. Specifically for pubic hair growth, a pictogram with a mark around the pubic region was shown to girls to facilitate ascertainment. Each girl was asked first whether or not she had experienced each pubertal event. If so, they were probed further about whether the event had recently begun or had been well underway, yielding a three-level ordinal response scale for each secondary sex characteristic: 1) prepubertal (not yet begun); 2) beginning; or 3) in an advanced stage of development. Although the original PDS had a response option for a postpubertal stage [30], it was not included in our adapted questionnaire due to the pubertal age range of participants.

Based on responses from the questionnaire, we categorized girls into maturational stages corresponding to those established by Tanner [31], where scales range successively from stage 1 indicating prematurity to stage 5 representing full maturity of each sex characteristic (Table 1). Girls who reported their breasts had not started to develop were categorized as stage 1 (B1). Stage 2 (B2) was assigned to girls who reported their breasts had just begun to develop. Stage 3–4 (B3-4) was assigned for those who answered their breast development had been well underway. We could not assign stage 5 as we did not ask if breast development was completed. The same classification was applied to gauge maturational progression of pubic hair growth (e.g. PH1, PH2 and PH3-4). We were not able to clinically verify girls’ responses primarily due to the sensitivity surrounding clinical ascertainment of puberty in a rural community setting.

[Figure omitted. See PDF.]

Table 1. Classification of maturational stages of breast development and pubic hair growth by Tanner and current study.

https://doi.org/10.1371/journal.pone.0247762.t001

For menarche, girls were asked to recall if they had experienced a first menstrual cycle and, if so, in what calendar year and month. To calculate menarcheal age, the default calendar day of 15 was assigned. If girls had reportedly attained menarche in the same month as when they were interviewed, a calendar day of 1 was assigned. For those who were premenarcheal at the baseline interview but reported to have attained menarche within the same month, as recalled at the follow-up assessment a year later, the last day of the month was assigned.

Statistical analysis

The distributions of stages of breast development and pubic hair growth by age are presented according to the 1st or 2nd visit and as total observations merged from both visits. Girls were grouped by age in years using mid-year as a cutoff; for example, girls aged 9.50–10.49 were considered to be 10-year-olds. Ages at stages of breast development and pubic hair growth were summarized as calculated means ± standard deviations (SDs). To maximize the age range presented and accommodate potential temporality in pubertal development, we pooled the data collected at baseline and follow-up assessments, with girls therefore contributing up to two maturational and age time points, yielding a total of 29,377 observations.

The prevalence of menarche and different maturational stages of breast development and pubic hair growth was calculated for each year of age at assessment. For breast development and pubic hair growth, proportions of girls at each age having achieved maturational stage 2 or greater (i.e. B2 or B3-4; PH2 or PH3-4) and stage 3 or greater (i.e. B3-4; PH3-4) were computed, respectively.

Probit regression was used to plot cumulative frequency curves for different stages of pubertal maturation according to age at assessment. In brief, a probit model fits the probability of the event, which follows an inverse standard normal distribution, as a form of linear regression with predictors of interest [32]. In our analysis, separate models were fit for different maturational stages of pubertal events with age in years as the independent variable of interest.

We next estimated ages at median (50th) and selective (3rd, 10th, 25th, 75th, 90th and 97th) percentiles of the onset of menarche and different maturational stages of breast development and pubic hair growth. A median age was obtained when the predicted probability from the probit model equaled 0.5, indicating the estimated age at which half of the study population reached a given maturational stage. As a supplementary analysis, we performed separate probit analysis for the data collected at baseline versus follow-up assessment. All statistical analyses and calculation were conducted using STATA 14 (StataCorp LP, College Station, Texas, USA).

Ethical review and approval

Ahead of study enrollment, all participating adolescent girls gave verbal assent with parental or guardian’s written consent. Both the maternal trial [22] and current adolescent follow-up study were reviewed and approved by the institutional review board at the Johns Hopkins Bloomberg School of Public Health, Baltimore, MD and the Bangladesh Medical Research Council, Dhaka, Bangladesh.

Results

A total of 15,320 and 14,057 girls were eligible at baseline and follow-up assessments, respectively, contributing a total of 29,377 analyzed observations. At baseline, 38.3% of girls had a height-for-age z-score <-2 [33] and 25.5% had a body mass index-for-age z-score <-2 [33] (S1 Table). Nearly all girls (96.3%) were currently enrolled in school. Nutritional and socioeconomic status of 1,263 girls who were lost to follow-up was largely comparable to those who remained in the study.

The age distributions of participants at baseline and follow-up visits are presented in Table 2 according to 1-year age strata and maturational stages of breast development and pubic hair growth. The age range at baseline was 9–15 years with few (0.1%) 15-year-olds; whereas by the follow-up visit 11.3% of girls had aged to 15 years or older. Maturational stages for each secondary sex characteristic expectedly advanced with age at the follow-up visit. For example, initially 31.5% of girls were in the earliest stage of prepubertal breast development (B1), which declined to 12.6% at follow-up. Pubic hair was inapparent (PH1) in 71.0% of girls at baseline, which declined to 45.3%. Age distributions within each 1-year interval are older by ~0.1 year at the follow-up visit, an artifact of imbalance across age intervals, and likely explaining the older follow-up age distributions within each stage of breast development and public hair growth. Consequently, Table 2 also presents composite mean ± SD ages for the two visits combined within each age stratum and maturational stage in order to present average sexual maturation patterns for this population, including those for 9- and 16-year-old children, which are uniquely available from the baseline and follow-up assessments, respectively.

[Figure omitted. See PDF.]

Table 2. Distributions of age and maturational stages of breast development and pubic hair growth at baseline and follow-up visits of adolescent girls in rural Bangladesh.

https://doi.org/10.1371/journal.pone.0247762.t002

Fig 2 summarizes composite age-specific patterns of attained breast development and pubic hair growth (for stages 2+ and 3–4) and menarche. At age 9 (i.e. an interval encompassing 8.50–9.49 years), 10.6% reported onset of breast development (B2+), with the percent rising steeply with age such that by 14 (i.e. 13.50–14.49 years) virtually all children had begun breast development (Fig 2A), and by age 15 years 94.8% were in an advanced stage (B3-4) (Fig 2B). On the other hand, pubic hair growth lagged, evident by <2% of 9–10 (i.e. 8.50–10.49) year-olds reporting any growth (PH2+) and a slower rise in percent attaining this characteristic (Fig 2C). At 15 and 16 years, 95.9% and 100% of girls, respectively, reported having pubic hair, and 73.7% and 85.7%, respectively, in an advanced stage of hair growth (Fig 2D). Menarche lagged both secondary characteristics, with 2.7% of aged 11 (i.e. 10.50–11.49 years) and 15.4%, 45.8% and 78.8% of girls becoming postmenarcheal at average exact ages of 12, 13 and 14, respectively. Virtually all girls were postmenarcheal by age 16 (i.e. 15.50–16.49) years (Fig 2E).

[Figure omitted. See PDF.]

Fig 2. Prevalence of pubertal maturational stages by age at assessment of adolescent girls in rural Bangladesh.

A, Percent attainment of breast development stages B2+ and, B, B3-4. C, Percent attainment of pubic hair growth stages PH2+ and, D, PH3-4. E, Prevalence of menarche. All ages are centered on birth dates ± 6 months.

https://doi.org/10.1371/journal.pone.0247762.g002

Fig 3 demonstrates the cumulative, composite frequency curves plotted from a probit model revealing intervals between maturational stages of each pubertal milestone by age. We note that, on average, the onset of breast development initiates puberty, followed by pubic hair growth and menarche, with all except the advanced stage of pubic hair growth converging asymptotically by age 16 (i.e. 15.50–16.49) years.

[Figure omitted. See PDF.]

Fig 3. Cumulative frequency curves for pubertal development of adolescent girls in rural Bangladesh.

Abbreviations: B2, stage 2 for breast development; B3-4, stage 3–4 for breast development; PH2, stage 2 for pubic hair growth; PH3-4, stage 3–4 for pubic hair growth. Patterns represent composite of two assessments ~1 year apart. All ages are centered on birth dates ± 6 months.

https://doi.org/10.1371/journal.pone.0247762.g003

Complementing Fig 2, the age distributions in attained stages 2 and 3–4 of secondary sex characteristics and menarche, estimated from the probit model, are presented in Table 3. The median (95% confidence interval [CI]) age of onset of breast development (B2) was 11.02 (11.00–11.04) years. For the onset of pubic hair growth (PH2) and menarche, the median (95% CI) age was 12.93 (12.91–12.94) and 13.17 (13.15–13.19) years, respectively. The duration of overall pubertal progression, often designated as the interval between age at breast development and menarche, was 2.15 years. The median intervals from B2 to other stages of pubertal milestones were 1.80, 1.91 and 3.27 years for B3-4, PH2 and PH3-4, respectively. As the data presented in Fig 2 and Table 3 summarize the composite data, distributions of the timing, sequence and progression of attainment of each pubertal milestone are also presented in supplementary material, for baseline (S1 Fig and S2 Table) and follow-up (S2 Fig and S3 Table) assessments, showing that they mirror one another, differing only to minor degrees, and thus reinforce the representativeness of the more encompassing composite distributions of sexual maturation.

[Figure omitted. See PDF.]

Table 3. Ages at which different stages of secondary sex characteristics and menarche were attained1.

https://doi.org/10.1371/journal.pone.0247762.t003

Discussion

In a large, contemporary birth cohort of young adolescent girls living in a typical rural setting in northwestern Bangladesh, we have described the sequence and tempo of pubertal development, reflected by breast development, pubic hair growth and menarche, for which international references exist for comparison. We observed the onset of these three events to occur, on average, at ages 11.02, 12.93 and 13.17 years.

The timing of the onset of breast development and menarche was found to be comparable to what has been reported from a well-published British population studied in the 1960s in which these pubertal milestones occurred at ages 11.15 and 13.47 years, respectively [21]. We have also placed our findings in context of relevant studies in the literature since the late 1970s, prioritized by their multiple pubertal milestones estimated with the use of probit analysis (Table 4). Also relevant, are findings from a recent meta-analysis of 27 studies in adolescent populations born in or after 1998 in low and middle income countries (LMICs) that reported pooled mean ages of 10.4 and 12.3 years for breast development and menarche, respectively [2]. From these comparisons, rural Bangladeshi girls appear to be maturing later than other adolescent populations in LMICs. Of note, 19 out of 27 studies included in the meta-analysis came from urban settings, which could affect developmental timing. Also, a later menarche was observed in studies from African than Asian countries, at ages 13.8 versus 12.3 years, respectively [2]. Together, these studies reveal variation in pubertal timing by country over a period of five decades within LMICs and that pubertal timing in our study population plausibly reflects later pubertal development in rural contexts. Beyond potential differences associated with study design, inter-population and decadal variation may be explained by suboptimal nutritional and socioeconomic status that often exist in rural Bangladesh [3, 34].

[Figure omitted. See PDF.]

Table 4. Mean or median age at the onset of secondary sex characteristics reported from different populations.

https://doi.org/10.1371/journal.pone.0247762.t004

Our findings also allow trends in age at menarche (AAM) to be examined over time within rural Bangladesh. In 1976, five years following independence and two years after a major flood and subsequent famine, menarche was reported to be occurring at age 15.8 in Matlab Thana in southeast Bangladesh [4]. By 1996, a study in nearby Narayanganj District reported a mean AAM of 13.0 years [13], which is comparable to our current estimate of 13.2 years in Gaibandha District in the northwest region. Data from our field site collected a decade before the current study reported a mean AAM of 12.8 years [14]. This was derived from an adolescent population in which 88% of girls had attained menarche; thus, by excluding premenarcheal adolescents, the attained estimate could have been biased towards an earlier age. A study conducted between 2016 and 2017 in Matlab Thana reported a median AAM of 13.0 years [3]. Taken together, our findings suggest that menarche has stabilized at ~13 years of age in rural Bangladesh over the last 2 decades. As a downward secular trend in AAM is considered to have plateaued at ~12 years of age as of the 1970s in industrialized countries, a further lowering in AAM may still occur in LMICs with improved living standards and nutritional conditions [10, 51], consistent with an earlier age of breast development being observed around the world [12]. Little data are available to be able to compare trends in pubic hair development within populations over time.

Girls in rural Bangladesh progressed through puberty evident by a predictable sequence of breast development, followed by pubic hair growth on average 1.91 years later, and after another 0.24 years by menarche [31]. A much shorter interval of 0.30–0.55 years between breast development and pubic hair growth, calculated from probit analysis, has been reported from the US and Denmark (Table 4) [36–38]. On the other hand, studies in Thailand [44] and China [39, 46] that used the same statistical approach as ours have reported intervals between 1.70–1.96 years, similar to this study. Whether this variation can be explained by differences in biological phenotype, such as a lack of pubic hair follicles, or production and circulating levels of adrenal androgens [39], remains to be elucidated.

Our observations of the timing of pubertal onset and its distribution, indicated by stage 2 of breast development (B2), were similar to what Marshall and Tanner observed in a longitudinal study of British girls reported in the 1960s [21]. Their pioneering data has evolved to become the standard for assessing extent and progression of pubertal development, including providing a definition of precocious puberty as starting at age 8 or younger, usually assessed by onset of breast development [21]. We estimate that <3% of girls in our study population experienced breast development before age 8.

The interval between ages at breast development and menarche is considered a reliable indicator of the duration of puberty [21]. In our study population, despite a relatively late pubertal onset at 11.02 years, girls progressed from breast development stage B2 to menarche over a period of 2.15 years, comparable to that reported in other populations [7, 35]. In the UK, Marshall and Tanner reported a similar pubertal duration between early versus late maturing girls [21], and more recent studies from Europe have reported that later onset of puberty was followed by a shorter duration of maturation [52, 53]. Of note, the interval between consecutive pubertal events derived from cross-sectional data is likely to appear longer than that from longitudinal data, due to higher variance in estimates between individuals [46].

To our knowledge, this is the first study concurrently characterizing multiple pubertal milestones of girls in rural Bangladesh and comparing findings with those in studied adolescent populations across different regions of the world. With the use of probit analysis, we were able to estimate pubertal timing of the entire, large sample of young adolescents, without biasing our estimates that can result from restricting analysis to those who have attained a pubertal milestone, and improve comparability with other studies that employed probit analysis [4, 5, 7, 8, 13, 36–39, 42–49]. Given the present study comprises a birth cohort with dates of birth known, girls’ ages were calculated accurately to the year, month and day.

As a limitation, rather than being ascertained via clinical examination, pubertal maturation was assessed by lay interviewers administering a questionnaire adapted from the Pubertal Development Scale (PDS) [30], with drawn anatomic site on pubic region to help children answer the questionnaire without showing explicit maturational images. This approach allowed pubertal maturation to be assessed in a very large number of children in a culturally sensitive and acceptable way. Several studies have noted a moderate to high level of agreement (e.g. weighted kappa or Kendall tau of 0.57 to 0.73) between the PDS and clinically ascertained Tanner stages, with greater agreement when the PDS was mapped onto the three-level pubertal scales [54–57], as used in the current study. Of relevance, a recent study of similarly aged and nourished girls in rural Pakistan [57] reported excellent agreement (weighted kappa of 0.73) between the PDS and physician-assessed Tanner stages, supporting the clinical validity of this self-reported method to gauge pubertal maturation in large-scale epidemiologic studies. Our observations of increasing maturation of secondary sex characteristics with advancing age and sequential timing of pubertal milestones were consistent with a well-established pattern of pubertal development [21]. Only ~2% of children reportedly “regressed” in recording breast development and pubic hair growth stages from baseline to the follow-up assessment, while 44% and 40% advanced, and 54% and 57% remained at the same reported stage, respectively, offering evidence of reliability (data not shown). In turn, menarche is a salient and thus particularly memorable event across a series of pubertal milestones; however, a short interval between an interview and past menarche is likely to yield more accurate data on recalled AAM [58]. Among menarcheal girls in our study population, the mean interval between the date of menarche and interview was 6.4 months, which we believe lies within the valid time period to support the accuracy of self-reported AAM [3, 58].

By characterizing multiple pubertal milestones including both primary and secondary sex characteristics in a large, contemporary birth cohort of adolescent girls in rural Bangladesh, the current study contributes knowledge to the literature on pubertal maturation of girls in South Asia [3–9, 13, 43, 44, 50]. Additional evidence on great detail of pubertal milestones of adolescent populations in LMICs, particularly in rural settings, will help establish the generalizability of our findings to other populations. Future studies elucidating nutritional, socioeconomic or environmental determinants of puberty in this context can be expected to contribute to a better understanding of maturational patterns of adolescent girls in South Asia and help inform the designs of interventions intended to benefit their growth, health and development.

In conclusion, adolescent girls in rural Bangladesh progressed through puberty following a well-documented sequence of breast development, pubic hair growth and menarche. The timing and pace of reaching these pubertal milestones appear to be similar to that reported in the western population in the 1960s. Observed onsets of phenotypical pubertal milestones were among the latest in a limited set of studies globally; whereas the tempo for the overall pubertal progression, that is from breast development to menarche, was in a comparable range.

Supporting information

S1 Checklist. STROBE statement—Checklist of items that should be included in reports of cohort studies.

https://doi.org/10.1371/journal.pone.0247762.s001

(DOCX)

S1 Fig. Cumulative frequency curves for pubertal development of adolescent girls in rural Bangladesh using data collected at baseline visit.

Abbreviations: B2, stage 2 for breast development; B3-4, stage 3–4 for breast development; PH2, stage 2 for pubic hair growth; PH3-4, stage 3–4 for pubic hair growth. All ages are centered on birth dates ± 6 months.

https://doi.org/10.1371/journal.pone.0247762.s002

(TIF)

S2 Fig. Cumulative frequency curves for pubertal development of adolescent girls in rural Bangladesh using data collected at follow-up visit.

Abbreviations: B2, stage 2 for breast development; B3-4, stage 3–4 for breast development; PH2, stage 2 for pubic hair growth; PH3-4, stage 3–4 for pubic hair growth. All ages are centered on birth dates ± 6 months.

https://doi.org/10.1371/journal.pone.0247762.s003

(TIF)

S1 Table. Baseline characteristics of adolescent Bangladeshi girls (9–15 years) who were included in the analysis in comparison to those lost to follow-up.

Abbreviations: BMI, body mass index calculated as weight (kg) / height (m)2; BMIZ, body mass index-for-age z-score [33]; HAZ, height-for-age z-score [33]; SD, standard deviation. 1 Eligibility defined as girls who were alive, met, gave consent and agreed to interview.

https://doi.org/10.1371/journal.pone.0247762.s004

(DOCX)

S2 Table. Ages at which different stages of secondary sex characteristics and menarche were attained using data collected at baseline visit.

Abbreviations: B2, stage 2 for breast development; B3-4, stage 3–4 for breast development; CI, confidence interval; P3-97, 3rd-97th percentile, respectively; PH2, stage 2 for pubic hair growth; PH3-4, stage 3–4 for pubic hair growth. 1 Ages at attainment were estimated from probit analysis using baseline assessment at which adolescent girls were aged 9–15 years.

https://doi.org/10.1371/journal.pone.0247762.s005

(DOCX)

S3 Table. Ages at which different stages of secondary sex characteristics and menarche were attained using data collected at follow-up visit.

Abbreviations: B2, stage 2 for breast development; B3-4, stage 3–4 for breast development; CI, confidence interval; P3-97, 3rd-97th percentile, respectively; PH2, stage 2 for pubic hair growth; PH3-4, stage 3–4 for pubic hair growth. 1 Ages at attainment were estimated from probit analysis using follow-up assessment at which adolescent girls were aged 10–16 years.

https://doi.org/10.1371/journal.pone.0247762.s006

(DOCX)

Acknowledgments

We gratefully acknowledge participating investigators whose contributions were critical to the study: Kelsey Alland and Sucheta Mehra coordinated field activities; Parul Christian contributed to conceptualization and protocol development; Maithilee Mitra and Rezwanul Haque directed and oversaw data management and constructed the database. We also thank the study participants, their parents and field interviewers involved in data collection.

Citation: Hur J, Schulze KJ, Thorne-Lyman AL, Wu LS-F, Shaikh S, Ali H, et al. (2021) Characterization of pubertal development of girls in rural Bangladesh. PLoS ONE 16(4): e0247762. https://doi.org/10.1371/journal.pone.0247762

References

1. Bishwajit G. Nutrition transition in South Asia: the emergence of non-communicable chronic diseases. F1000Res. 2015;4:8. pmid:26834976

2. Moodie JL, Campisi SC, Salena K, Wheatley M, Vandermorris A, Bhutta ZA. Timing of pubertal milestones in low- and middle-income countries: a systematic review and meta-analysis. Adv Nutr. 2020;11(4):951–9. pmid:32027344

3. Svefors P, Pervin J, Islam Khan A, Rahman A, Ekstrom EC, El Arifeen S, et al. Stunting, recovery from stunting and puberty development in the MINIMat cohort, Bangladesh. Acta Paediatr. 2020;109(1):122–33. pmid:31283046

4. Chowdhury AK, Huffman SL, Curlin GT. Malnutrition, menarche, and marriage in rural Bangladesh. Soc Biol. 1977;24(4):316–25. pmid:607404

5. Agarwal DK, Agarwal KN, Upadhyay SK, Mittal R, Prakash R, Rai S. Physical and sexual growth pattern of affluent Indian children from 5 to 18 years of age. Indian Pediatr. 1992;29(10):1203–82. pmid:1286886

6. Gupta N, Singh MP, Dhillon BS, Saxena NC. Preparing for adulthood—patterns of physical growth, sexual maturity and menarche of adolescent girls in selected urban slums and rural areas. J Indian Med Assoc. 2007;105(3):119–22, 26. pmid:17824461

7. Qamra SR, Mehta S, Deodhar SD. A mixed-longitudinal study on the pattern of pubertal growth: relationship to socioeconomic status and caloric-intake—IV. Indian Pediatr. 1991;28(2):147–56. pmid:2055629

8. Khadgawat R, Marwaha RK, Mehan N, Surana V, Dabas A, Sreenivas V, et al. Age of onset of puberty in apparently healthy school girls from northern India. Indian Pediatr. 2016;53(5):383–7. pmid:27254044

9. Campisi SC, Humayun KN, Rizvi A, Lou W, Soder O, Bhutta ZA. Later puberty onset among chronically undernourished adolescents living in a Karachi slum, Pakistan. Acta Paediatr. 2020;109(5):1019–25. pmid:31606907

10. Parent AS, Teilmann G, Juul A, Skakkebaek NE, Toppari J, Bourguignon JP. The timing of normal puberty and the age limits of sexual precocity: variations around the world, secular trends, and changes after migration. Endocr Rev. 2003;24(5):668–93. pmid:14570750

11. Parent AS, Franssen D, Fudvoye J, Gerard A, Bourguignon JP. Developmental variations in environmental influences including endocrine disruptors on pubertal timing and neuroendocrine control: revision of human observations and mechanistic insight from rodents. Front Neuroendocrinol. 2015;38:12–36. pmid:25592640

12. Eckert-Lind C, Busch AS, Petersen JH, Biro FM, Butler G, Brauner EV, et al. Worldwide secular trends in age at pubertal onset assessed by breast development among girls: a systematic review and meta-analysis. JAMA Pediatr. 2020;174(4):e195881. pmid:32040143

13. Chowdhury S, Shahabuddin AK, Seal AJ, Talukder KK, Hassan Q, Begum RA, et al. Nutritional status and age at menarche in a rural area of Bangladesh. Ann Hum Biol. 2000;27(3):249–56. pmid:10834289

14. Rah JH, Shamim AA, Arju UT, Labrique AB, Rashid M, Christian P. Age of onset, nutritional determinants, and seasonal variations in menarche in rural Bangladesh. J Health Popul Nutr. 2009;27(6):802–7. pmid:20099764

15. United Nations Fund for Population Activities. Marrying too young: end child marriage. New York: UNFPA; 2012.

16. Field E, Ambrus A. Early marriage, age of menarche, and female schooling attainment in Bangladesh. Journal of Political Economy. 2008;116(5):881–930.

17. Ibitoye M, Choi C, Tai H, Lee G, Sommer M. Early menarche: a systematic review of its effect on sexual and reproductive health in low- and middle-income countries. PLoS One. 2017;12(6):e0178884. pmid:28591132

18. Shahabuddin AS, Nostlinger C, Delvaux T, Sarker M, Bardaji A, Brouwere VD, et al. What influences adolescent girls’ decision-making regarding contraceptive methods use and childbearing? a qualitative exploratory study in Rangpur District, Bangladesh. PLoS One. 2016;11(6):e0157664. pmid:27336673

19. Sayem AM, Nury AT. Factors associated with teenage marital pregnancy among Bangladeshi women. Reprod Health. 2011;8:16. pmid:21599904

20. Rah JH, Christian P, Shamim AA, Arju UT, Labrique AB, Rashid M. Pregnancy and lactation hinder growth and nutritional status of adolescent girls in rural Bangladesh. J Nutr. 2008;138(8):1505–11. pmid:18641198

21. Marshall WA, Tanner JM. Variations in pattern of pubertal changes in girls. Arch Dis Child. 1969;44(235):291–303. pmid:5785179

22. West KP Jr., Christian P, Labrique AB, Rashid M, Shamim AA, Klemm RD, et al. Effects of vitamin A or beta carotene supplementation on pregnancy-related mortality and infant mortality in rural Bangladesh: a cluster randomized trial. JAMA. 2011;305(19):1986–95. pmid:21586714

23. Labrique AB, Christian P, Klemm RD, Rashid M, Shamim AA, Massie A, et al. A cluster-randomized, placebo-controlled, maternal vitamin A or beta-carotene supplementation trial in Bangladesh: design and methods. Trials. 2011;12:102. pmid:21510905

24. Gunnsteinsson S, Labrique AB, West KP Jr., Christian P, Mehra S, Shamim AA, et al. Constructing indices of rural living standards in Northwestern Bangladesh. J Health Popul Nutr. 2010;28(5):509–19. pmid:20941903

25. Klemm RD, Merrill RD, Wu L, Shamim AA, Ali H, Labrique A, et al. Low-birthweight rates higher among Bangladeshi neonates measured during active birth surveillance compared to national survey data. Matern Child Nutr. 2015;11(4):583–94. pmid:23647669

26. West KP Jr., Shamim AA, Mehra S, Labrique AB, Ali H, Shaikh S, et al. Effect of maternal multiple micronutrient vs iron-folic acid supplementation on infant mortality and adverse birth outcomes in rural Bangladesh: the JiVitA-3 randomized trial. JAMA. 2014;312(24):2649–58. pmid:25536256

27. Christian P, Shaikh S, Shamim AA, Mehra S, Wu L, Mitra M, et al. Effect of fortified complementary food supplementation on child growth in rural Bangladesh: a cluster-randomized trial. Int J Epidemiol. 2015;44(6):1862–76. pmid:26275453

28. Na M, Mehra S, Christian P, Ali H, Shaikh S, Shamim AA, et al. Maternal dietary diversity decreases with household food insecurity in rural Bangladesh: a longitudinal analysis. J Nutr. 2016;146(10):2109–16. pmid:27581578

29. Na M, Shamim AA, Mehra S, Labrique A, Ali H, Wu LS, et al. Maternal nutritional status mediates the linkage between household food insecurity and mid-infancy size in rural Bangladesh. Br J Nutr. 2020;123(12):1415–25. pmid:32102702

30. Petersen AC, Crockett L, Richards M, Boxer A. A self-report measure of pubertal status: reliability, validity, and initial norms. J Youth Adolesc. 1988;17(2):117–33. pmid:24277579

31. Tanner JM. Growth at adolescence: with a general consideration of the effects of hereditary and environmental factors upon growth and maturation from birth to maturity. 2nd ed. Oxford: Blackwell Scientific Publications; 1962.

32. Hosmer DW, Lemeshow S, Sturdivant RX. Applied logistic regression. 3rd ed. New Jersey: John Wiley & Sons; 2013.

33. de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, Siekmann J. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007;85(9):660–7. pmid:18026621

34. Thorne-Lyman AL, Shaikh S, Mehra S, Wu LSF, Ali H, Alland K, et al. Dietary patterns of >30,000 adolescents 9–15 years of age in rural Bangladesh. Ann N Y Acad Sci. 2020;1468(1):3–15. pmid:31403718

35. Biro FM, Huang B, Crawford PB, Lucky AW, Striegel-Moore R, Barton BA, et al. Pubertal correlates in black and white girls. J Pediatr. 2006;148(2):234–40. pmid:16492435

36. Wu T, Mendola P, Buck GM. Ethnic differences in the presence of secondary sex characteristics and menarche among US girls: the third National Health and Nutrition Examination Survey, 1988–1994. Pediatrics. 2002;110(4):752–7. pmid:12359790

37. Juul A, Teilmann G, Scheike T, Hertel NT, Holm K, Laursen EM, et al. Pubertal development in Danish children: comparison of recent European and US data. Int J Androl. 2006;29(1):247–55. pmid:16466546

38. Herman-Giddens ME, Slora EJ, Wasserman RC, Bourdony CJ, Bhapkar MV, Koch GG, et al. Secondary sexual characteristics and menses in young girls seen in office practice: a study from the Pediatric Research in Office Settings network. Pediatrics. 1997;99(4):505–12. pmid:9093289

39. Huen KF, Leung SS, Lau JT, Cheung AY, Leung NK, Chiu MC. Secular trend in the sexual maturation of southern Chinese girls. Acta Paediatr. 1997;86(10):1121–4. pmid:9350897

40. Miller EM. Iron status and reproduction in US women: National Health and Nutrition Examination Survey, 1999–2006. PLoS One. 2014;9(11):e112216. pmid:25375360

41. Biro FM, Greenspan LC, Galvez MP, Pinney SM, Teitelbaum S, Windham GC, et al. Onset of breast development in a longitudinal cohort. Pediatrics. 2013;132(6):1019–27. pmid:24190685

42. Aksglaede L, Sorensen K, Petersen JH, Skakkebaek NE, Juul A. Recent decline in age at breast development: the Copenhagen Puberty Study. Pediatrics. 2009;123(5):e932–9. pmid:19403485

43. Rao S, Joshi S, Kanade A. Height velocity, body fat and menarcheal age of Indian girls. Indian Pediatr. 1998;35(7):619–28. pmid:10216671

44. Mahachoklertwattana P, Suthutvoravut U, Charoenkiatkul S, Chongviriyaphan N, Rojroongwasinkul N, Thakkinstian A, et al. Earlier onset of pubertal maturation in Thai girls. J Med Assoc Thai. 2002;85 Suppl 4:S1127–34. pmid:12549786

45. Rabbani A, Khodai S, Mohammad K, Sotoudeh A, Karbakhsh M, Nouri K, et al. Pubertal development in a random sample of 4,020 urban Iranian girls. J Pediatr Endocrinol Metab. 2008;21(7):681–7. pmid:18780603

46. Ma HM, Du ML, Luo XP, Chen SK, Liu L, Chen RM, et al. Onset of breast and pubic hair development and menses in urban Chinese girls. Pediatrics. 2009;124(2):e269–77. pmid:19651567

47. Kashani HH, Kavosh MS, Keshteli AH, Montazer M, Rostampour N, Kelishadi R, et al. Age of puberty in a representative sample of Iranian girls. World J Pediatr. 2009;5(2):132–5. pmid:19718536

48. Dasgupta P. Physical growth, body composition and nutritional status of Bengali school aged children, adolescents and young adults of Calcutta, India: effects of socioeconomic factors on secular trends. Neys-van Hoogstraten Foundation, The Netherlands and Indian Statistical Institute Kolkata, India, 2015.

49. Opare-Addo PM, Stowe M, Ankobea-Kokroe F, Zheng T. Menarcheal and pubertal development and determining factors among schoolgirls in Kumasi, Ghana. J Obstet Gynaecol. 2012;32(2):159–65. pmid:22296429

50. Rahmawati NT, Hastuti J, Suriyanto RA. Relationship between somatotype and age at menarche among adolescent girls in Yogyakarta province, Indonesia. Pakistan Journal of Nutrition. 2017;16(5):351–8.

51. Gluckman PD, Hanson MA. Evolution, development and timing of puberty. Trends Endocrinol Metab. 2006;17(1):7–12. pmid:16311040

52. Marti-Henneberg C, Vizmanos B. The duration of puberty in girls is related to the timing of its onset. J Pediatr. 1997;131(4):618–21. pmid:9386670

53. Llop-Vinolas D, Vizmanos B, Closa Monasterolo R, Escribano Subias J, Fernandez-Ballart JD, Marti-Henneberg C. Onset of puberty at eight years of age in girls determines a specific tempo of puberty but does not affect adult height. Acta Paediatr. 2004;93(7):874–9. pmid:15303800

54. Chan NP, Sung RY, Nelson EA, So HK, Tse YK, Kong AP. Measurement of pubertal status with a Chinese self-report Pubertal Development Scale. Matern Child Health J. 2010;14(3):466–73. pmid:19517073

55. Pompeia S, Zanini GAV, Freitas RS, Inacio LMC, Silva FCD, Souza GR, et al. Adapted version of the Pubertal Development Scale for use in Brazil. Rev Saude Publica. 2019;53:56. pmid:31432913

56. Koopman-Verhoeff ME, Gredvig-Ardito C, Barker DH, Saletin JM, Carskadon MA. Classifying pubertal development using child and parent report: comparing the pubertal development scales to Tanner staging. J Adolesc Health. 2020;66(5):597–602. pmid:31974011

57. Campisi SC, Humayun KN, Wasan Y, Soofi S, Islam M, Lou W, et al. Self-assessed puberty is reliable in a low-income setting in rural Pakistan. J Pediatr Endocrinol Metab. 2020;33(9):1191–6. pmid:32817583

58. Koo MM, Rohan TE. Accuracy of short-term recall of age at menarche. Ann Hum Biol. 1997;24(1):61–4. pmid:9022907