Prospective association of dietary soy and fibre intake with puberty timing: a cohort study among Chinese children
BMC Medicine volume 20, Article number: 145 (2022)
Dietary phytoestrogens have been suggested to influence puberty timing, a critical stage for well-being in adulthood. We hypothesized that childhood soy intake might prospectively influence puberty timing and that dietary fibre and the key isoflavone metabolite equol might play roles.
Cox proportional hazard regression models were performed in 4781 children (2152 girls and 2629 boys) aged 6–8 years old from the Chinese Adolescent Cohort Study for whom a food frequency questionnaire at baseline and information about potential confounders were available. Anthropometry and pubertal status including age at Tanner stage 2 for breast development (B2) or age at the initiation of gonadal growth (G2), and age at menarche (M) or voice break (VB) were assessed annually. Equol excretion was determined by urine samples from 1311 participants.
Among girls and boys, higher soy intake was associated with later puberty timing (hazard ratio (HR)-B2: 0.88 (95% CI, 0.80–0.96), p=0.02; HR-M, 0.87 (0.77–0.94), p=0.01; HR-G2, 0.91 (0.82–0.98), p=0.013; HR-VB, 0.90 (0.82–0.9), p=0.02), independent of prepubertal body fatness and fibre intake. These associations were more pronounced among children with a high urinary equol level (pfor-interaction ≤ 0.04) or with a high cereal fibre intake (pfor-interaction ≤ 0.06). Intake of dietary fibre or its subtype was not prospectively associated with puberty onset after adjusting for dietary soy intake (p≥0.06).
Higher childhood soy intake is prospectively associated with later puberty timing in both Chinese girls and boys, independent of prepubertal body fatness, and the association is particularly pronounced among individuals with a higher urinary equol level.
Puberty is a critical period of time encompassing sequential dramatic developmental changes to reach mature reproductive functions. Early puberty onset is recognized as an established risk factor for all-cause mortality, hormone-related cancers, insulin resistance and obesity later in life [1, 2]. Thus, factors influencing puberty timing have been increasingly acknowledged [3, 4]. Over the last 10 years, the impact of childhood dietary phytoestrogen intake on puberty timing has been discussed. Isoflavones, the major class of dietary phytoestrogen , is structurally and functionally similar to endogenous oestrogen . Whether isoflavone intake in childhood is implicated in the timing of puberty is currently a controversially debated issue. Support for its relevance has come from prospective studies, primarily conducted in girls, which showed that higher dietary isoflavone intake was associated with later breast development  or menarche  among German  and US  girls. In contrast, Wolff et al. [9, 10] showed that urinary isoflavone levels were not associated with breast or pubic hair development, and dietary isoflavone intake was not associated with voice break . To date, the above studies were all conducted in Western countries. Since soy and its products, which are the major food sources of dietary isoflavones, are substantially consumed in developing countries, prospective investigation and understanding of the role of dietary isoflavones in puberty timing in these developing populations have important implications.
Moreover, food rich in isoflavones often contains large amounts of dietary fibre. Findings on the association of dietary fibre intake with pubertal development have been inconsistent among prospective cohorts, and evidence from non-Western populations has been lacking: high fibre consumption was associated with later menarche among Dutch  and Canadian  girls, while a null association was reported in US girls  and German children . Therefore, the relevance of dietary isoflavones and fibre intake with puberty timing merits simultaneous evaluation. Interestingly, dietary fibre has been suggested to be associated with isoflavone bioavailability . Isoflavones are transformed into aglycones with the help of gut microbes, leading to easier absorbance and higher biological activities . Equol is an isoflavone bacterial metabolite and a key player in soy-related health benefits , and its production relies on equol-producing bacteria and is supported by fibre intake. More equol producers were reported in populations with higher dietary fibre intake, e.g., approximately 30% of omnivorous Caucasians , 59% of vegetarians  and 60.4% of Asians , presumably due to microbiota heterogeneity and dietary preferences for fibre subtypes favouring equol-producing bacterial growth. Nonetheless, existing studies have not considered the impacts of equol excretion or fibre subtype on pubertal development.
Over the last 40 years, a noteworthy secular trend in early puberty onset has been observed in both Chinese boys (from 16.1 years old in 1979  to 14.3 years old in 2010 ) and girls (from 13.5 years old in 1979  to 12.3 years old in 2014 ). Given that children with early puberty onset are associated with a spectrum of diseases leading to increasing disease burden in China [24,25,26], the influences of dietary isoflavones and fibre intake on puberty timing, like those investigated in our study, which include a large representative sample of Chinese children, will have important public health implications. Because the Chinese food component database is under construction, we conducted this analysis with dietary soy intake rather than calculating dietary isoflavones. Using prospective data from the Chinese Adolescent Cohort Study (CAC), we investigated the hypothesis that higher childhood dietary soy and fibre intake was associated with later puberty timing (as indicated by the early pubertal markers of age at Tanner stage 2 for breast development in girls (B2) and age at the initiation of gonadal growth in boys (G2) and the late pubertal markers of age at menarche (M) and age at voice break (VB)) and that this association could be modified both by urinary equol excretion and cereal fibre intake.
Study design and participants
The CAC study details have been described elsewhere . Cooperative and voluntary children aged 6–8 years old in 23 selected schools have been recruited yearly since 2013. Baseline information included sociodemographic issues, dietary intake and eating behaviours, physical activity and sedentary behaviours, anthropometry and pubertal development. Follow-up data on nutrition, growth, metabolism, and health status were collected at regular intervals until the children were 15 years old: anthropometry and puberty assessments were conducted annually, and dietary intake and physical activity data were collected biennially. This study was approved by the Ethics Committee of Sichuan University, and all of the parents of the participants provided their written confirmed consent before enrolment. All examinations and questionnaires were administered with parental consent.
Between January 2013 and December 2018, 6967 children aged 6–8 years old were included at baseline. Of these children, 5439 had completed at least 2 follow-up assessments by the end of 2020. Since we were interested in the prospective relevance of diet to puberty timing, 389 children who had already reached B2/G2 at baseline were excluded from our current analysis. Among them, 141 participants with implausible energy intakes (less than or greater than age- and sex-specific cut-offs)  and 128 children with incomplete information on potential confounders were further excluded. In total, 4781 children (2152 girls, 2629 boys) were eligible (Additional file 1: Fig. S1), and of them, 1311 children provided first morning voided midstream urine samples.
Nutrition data were collected via a validated food frequency questionnaire (FFQ) by trained investigators . This FFQ included 17 categories of the 53 most representative local foods or food groups among children: whole grains, refined grains, tubers, vegetables, fruits, nuts, meat, fish and shrimp, animal viscera, eggs, dairy and products, total soy (soybean and products), fried foods, sugary snacks, sugar-sweetened beverages, fruit juices and dietary supplements. The participants reported their frequency (never, daily, weekly, monthly or annually) for each item and estimated portion sizes using food models and picture aids. During the interviews, the investigators checked FFQs for potentially incorrect responses and made clarifications when necessary. Dietary intake data were converted into energy and nutrient intake data using the continuously updated in-house nutrient database based on NCCW software (version 11.0, 2014), which reflects food composition in China.
This study investigated individual mean daily intakes of total soy (soy and soy products), dietary fibre and major fibre subtypes (cereal fibre: cereals, noodle, rice, tubers, cookies and cakes; fruit fibre: fruits and its products; vegetable fibre: vegetables and its products).
Detailed instructions on collecting first morning voided midstream urine samples were carefully given to parents and children. All urine samples were stored immediately at − 20 °C before transportation and then at − 80 °C until analysis. Equol levels were determined using a previously validated gas chromatography-mass spectrometry method . The detection limit was 3.8 ng/ml. All laboratory equipment was calibrated, and blinded duplicate samples were used. All of the data were double entered into the database.
According to Tanner stage standardized criteria , B2 and pubic hair (girls and boys) were assessed at each examination by investigators. G2 was assessed by comparative palpation with a Prader orchidometer. If the volumes of the two testes varied, the larger volume was recorded. Testicular volume less than 1 mL was recorded as 1. Moreover, children were asked whether M or VB occurred during the annual physical examination; if so, respective months and years were recorded.
An ultrasonic weight and height metre (DHM-30, Dingheng Ltd, Zhengzhou, China) was used to assess standing height to the nearest 0.1 cm and weight to the nearest 0.1 kg with the subject lightly dressed and barefoot. Triceps skinfold thicknesses and subscapular angle sites were measured on the right side to the nearest 0.1 mm using Holtain callipers (Holtain Ltd, Crymych, UK). All measurements were performed twice to calculate averages. Body mass index (BMI) sex- and age-independent BMI standard deviation scores (SDS) were calculated using Chinese reference curves . Overweight was defined according to the International Obesity Task Force (IOTF) BMI cut-offs for children, which corresponds to an adult BMI of 25 kg/m2 . The percent body fat (%BF) was calculated using Slaughter equations .
Information on the frequency, duration and type of physical activity in various settings among children was collected by a validated physical activity questionnaire with 38 items (e.g., walking, running, climbing stairs, ball games, dancing) . The participants reported typical time spent on sedentary behaviours associated with television, computers, smartphones and homework.
Furthermore, parents provided information about pregnancy and infancy (i.e., children's birth weight, exclusive breastfeeding duration, timing of complementary feeding) and domestic characteristics (i.e., residency, income, family size, smoking status, parental age, occupations and education level).
SAS® procedures (version 9.4, SAS Inc., Cary, NC, USA) and Stata 14 (Stata Corp., College Station, TX, USA) were used for data analyses. All analyses were performed with a significance level of p<0.05. Although there was no statistical interaction between dietary soy intake and sex, in theory, dietary oestrogen, similar to endogenous oestrogen, might differentially impact the course of puberty in girls and boys . Data from girls and boys were thus analysed separately.
Since energy intake has been suggested to influence pubertal development, and energy intake is dependent on age , intake of total soy was expressed as age-specific residuals from the regression of soybean and its product intake on energy intake. Similarly, dietary fibre intake was expressed as age-specific residuals from the regression of fibre intake on energy intake. To examine the potential associations of total soy intake or dietary fibre intake with puberty timing, their distributions were grouped into tertiles (T1–T3).
The Kolmogorov-Smirnov and Shapiro-Wilk tests were conducted to examine the data for normality. Baseline birth weight and %BF were nonnormally distributed and presented as medians together with the interquartile rages; other continuous variables were normally distributed and presented as the means with their standard deviations (SD). Differences in anthropometric, sociodemographic and nutritional data between tertiles were analysed using an ANOVA for normally distributed continuous variables, the Kruskal-Wallis test for nonnormally distributed continuous variables, and the chi-square test for categorical variables. Statistical models and descriptive tables were stratified by sex.
Cox proportional hazard regression models were used to investigate the prospective relevance of total soy or dietary fibre (and its types) intake at baseline with age at B2/G2 or M/VB. Censoring occurred at the age of reaching B2/G2 and M/VB or age at the last follow-up if puberty events had not been reported.
In the basic models, the tertiles of total soy intake (residuals) or dietary fibre (and its types) intake (residuals) at baseline were the main independent fixed effects. These following potential confounders were considered for the Cox regression models: birth weight, age at baseline, school location, physical activity, body composition (Z scores of BMI, overweight (Y/N), %BF), parental/paternal/maternal educational level, family income, mother’s age at menarche, smoking status in the household, and total energy intake at baseline, as well as dietary fibre intake (residual) at baseline (in the total soy intake model) and total soy intake (residual) at baseline (in the dietary fibre intake model). In addition, we conceptualized confounders using the DAG platform  to validate and justify the potential confounders. Each potential confounder was initially considered separately and was included if it was associated with both the dietary index and indicators of puberty timing and if it substantially altered the estimate by more than 10% . As high levels of isoflavones and dietary fibre often coexist in food, Model 2 was adjusted for parental education level, energy intake at baseline, mother’s age at menarche, and fibre intake (residuals) at baseline (in the total soy intake model) or total soy intake (residuals) at baseline (in the dietary fibre intake model). In the final model, we controlled for confounding and/or mediation by percent body fat at baseline (Model 3), because it has been proposed that body composition in childhood might be relevant to the timing of puberty . Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated by comparing the 2nd and 3rd tertiles to the 1st tertile in these models. We assessed the linear trends by entering the value of dietary fibre/soy intake as a continuous variable in the above models.
To explore potential nonlinear relationships, we examined the associations (based on Model 3) of dietary soy intake and fibre intake with pubertal markers using restricted cubic spline models (four knots, according to Harrell’s recommendation ) among all of the participants. Four knots offer an adequate fit of the model and constitute a good compromise between flexibility and loss of precision caused by overfitting.
Moreover, we tested the potential interactions of urinary equol level (or fibre and its subtype intake) on the relationship between dietary soy intake and puberty timing. Further stratified analyses were conducted if the p for interaction was < 0.05.
To test the robustness of our results, we re-run our analyses using mixed model (PROC MIXED in SAS) with school clustering as a random effect, to investigate the associations of total soy intakes or fibre intakes in childhood with puberty timing.
Characteristics are presented according to tertiles of soy intake (Table 1). The mean baseline age was 7.2 (0.7) years old for girls and 7.3 (0.6) years old for boys. Among girls, 1748 (81.2%) and 1162 (54.0%) reached B2 and M, respectively, soy intake varied from 0 to 69 g/day, and the highest tertile had a lower baseline %BF and higher fibre intake. In boys, 1233 (46.9%) and 829 (31.5%) reached G2 and VB, respectively, soy intake ranged from 0 to 82.6 g/day and the highest tertile had a lower baseline %BF and a higher parental education level. According to the Chinese food component database , the estimated isoflavone intake was approximately 65 mg/day in girls and 87 mg/day in boys, based on the mean soy intake of our sample. There were no significant differences in age, BMI SDS or %BF between eligible and excluded participants.
The relevance of dietary soy intake with puberty timing
In both girls and boys, higher total soy intake was correlated with later puberty timing, independent of prepubertal %BF and fibre intake (Table 2, p for trend ≤ 0.03): girls with high dietary total soy intake had a 12% lower risk of reaching B2 or a 13% lower risk of experiencing M than girls with low dietary total soy intake, and boys with high dietary total soy intake had a 9% lower risk of reaching G2 or a 10% lower risk of experiencing VB than boys with low dietary total soy intake. Similar results were observed in the sensitivity analysis using mixed model with school clustering as a random effect (Additional file 2: Tables S1 and S2).
Interestingly, these associations were modified by the urinary equol level (p for interaction ≤ 0.04): girls with a high urinary equol level and a high dietary total soy intake had an approximately 13% lower risk of reaching B2 or a 16% lower risk of experiencing M than those in the opposite group (Fig. 1A, B), and boys with a high urinary equol level and a high dietary total soy intake had a 12% lower risk of reaching G2 or a 13% lower risk of experiencing VB than those in the opposite group (Fig. 1C, D). Similar results were observed in the participants with urinary equol data (Additional file 2: Table S3). In addition, we found evidence of nonlinear associations (p value for nonlinearity < 0.005) for dietary soy intake, with strong inverse associations with pubertal markers at low intake levels, but weaker associations at moderate to higher intake levels (Additional file 1: Fig. S2: A, B, C and D).
Furthermore, the association of total soy intake with puberty timing was modified by the cereal fibre intake (p for interaction ≤ 0.06): girls with a high cereal fibre intake and a high dietary total soy intake had approximately 13% lower risk of reaching B2 or an 18% lower risk of experiencing M than those in the opposite group (Fig. 2A, B). Boys with a high cereal fibre intake and a high dietary total soy intake had an 11% lower risk of reaching G2 or a 12% lower risk of experiencing VB than those in the opposite group (Fig. 2C, D). Similar effects were not observed for vegetable and fruit fibre intake (data not shown).
The relevance of dietary fibre intake with puberty timing
Total fibre intake was not significantly associated with puberty timing after adjusting for soy intake (Table 3, p for trend ≥ 0.06). Results remained unchanged in the sensitivity analysis using mixed model with school clustering as a random effect (Additional file 2: Table S4). Similar results were observed with subtype fibre intake (Additional file 2: Tables S5 and S6). Moreover, the associations between dietary fibre intake and pubertal markers appeared to be largely linear (Additional file 1: Fig. S2: E, F, G and H).
In the present analysis, girls and boys with higher total soy intake reached puberty later than children with lower soy intake, independent of prepubertal body fat and fibre consumption. These prospective associations were particularly pronounced among those with a higher urinary equol level or those with a high cereal fibre intake.
Puberty timing is a developmental milestone for children. Prospective cohort studies of soy or isoflavone intake and puberty timing have been limited. Our findings are in line with those from 230 US girls with higher phytoestrogen intake and delayed age at menarche  and from 119 German girls from the Dortmund Nutritional and Anthropometric Longitudinally Designed (DONALD) study with higher isoflavone intake and later breast development . However, the associations between soy or isoflavone intake and age at menarche were not consistent between girls from our study and those in the DONALD study. This inconsistency might partly be due to the small sample size, population bias toward high socioeconomic status, and lack of control over potential maternal genetic influences in the DONALD study. In addition, we speculate that the average level of individual isoflavone intake might also play a role. In our study, the estimated median value of soy isoflavone intake in girls was higher than that in Western populations (approximately 65 mg/day vs. 0.04–19.2 mg/day [7, 42]). Therefore, we are led to believe that the association between isoflavone intake and age at menarche in German girls might be obscured by a low isoflavone intake level. Similarly, the relevance between isoflavone intake and puberty timing was absent in 108 boys from the DONALD Study . However, we have evidently observed that higher soy intake was correlated with later age at gonadal growth and voice break in Chinese boys, possibly because the estimated isoflavone intake level in boys from our study was markedly higher than the reported isoflavone intake level in boys from the DONALD study (approximately 87 mg/day vs. 0.8–62.7 mg/day [7, 43]).
Although the mechanism underlying the association between soy intake (isoflavone) and puberty timing remains unknown, we speculate that possible explanations are: isoflavone is structurally and functionally similar to endogenous oestrogen  and could inhibit the activity of aromatase, the rate-limiting enzyme of oestrogen biosynthesis , therefore possibly affecting endocrine homeostasis and further pubertal development. Conversely, isoflavones have been reported to directly bind to and influence the expression of oestrogen receptors (ERs)  and impact the hypothalamus-pituitary-gonadal (HPG) axis . Menarche is a late pubertal stage primarily governed by the HPG axis, whereas breast development represents an early pubertal stage influenced by both the HPG axis and peripheral factors, including oestrogen-related enzymes and ERs on gonadal cells . Therefore, we speculate that low levels of isoflavone intake is able to influence breast development, which was observed in girls from the DONALD study , while the relevance between isoflavone intake and menarche was absent because a higher isoflavone concentration might be needed to influence the HPG axis. However, further investigations of the HPG axis and peripheral pubertal development are needed to elucidate possible underlying mechanisms.
Food rich in isoflavones typically contains appreciable amounts of dietary fibre. The present analysis demonstrated that dietary fibre and its subtypes were not independently associated with puberty timing. Our data further suggested that the associations between total soy intake and pubertal markers were more pronounced in children with higher cereal fibre intake. In our analysis, cereal fibre was defined as fibre from cereals, noodles, rice, tubers, cookies and cakes, representing a major group of influential nutrients known as microbiota-accessible carbohydrates (MACs). MACs vary distinctly in Western and traditional fibre-rich diets, leading to changes in gut microbiota, microbial functionality and bacteria-host interactions . Cereal fibre intake was not independently related to puberty timing in our participants; however, it was intriguing to observe a modified effect of dietary cereal fibre on soy-puberty associations. One possible explanation is that high cereal fibre intake, hence elevated MAC levels, facilitates the establishment of a “healthy” gut microbiota favouring an isoflavone-friendly microbial ecosystem in which soy is efficiently digested, absorbed and utilized .
Gut microbes play a critical role in transforming glucose-conjugated soy isoflavones into aglycones, which are easier to absorb and retain higher biological activities than the glycosylated parent compounds . Additionally, equol production occurs in the intestine via bacterial reductases from a number of gut microbes, and the list of equol-producing bacteria is constantly expanding with newly characterized members . More equol producers have been reported in populations with higher fibre intake [17,18,19]. Since equol is generally considered to be the main conductor for soy-related health benefits , the relationship between equol and puberty onset is of particular interest. However, observational data on this matter are scarce. In this study, associations of total soy intake with puberty timing were more prominent in children with higher urinary equol levels, which might indicate that a pro-equol-production gut microbiota might reinforce the impact of soy consumption on puberty timing. Combined with the findings on cereal fibre intake, we believe that a higher intake of foods rich in MACs might build an optimal intestinal microflora, likely indicated by urinary equol levels, to efficiently deploy the puberty-influential effects of soy isoflavones.
Our study has several strengths. In contrast to most studies which focused only on girls, we observed both girls and boys. Our participants and their parents/family were representative of the general population in age, economic and educational status according to regional statistical books . The prospective nature and repeated detailed measurements of anthropometric, pubertal and dietary data in participants, in conjunction with the ability to adjust for a number of major potential confounders both in children and in their parents, represent substantial strengths. Notably, prepubertal body fatness was considered in our analysis because childhood body composition might potentially influence puberty onset .
One of our study’s limitations is that the observational design of our study could not establish causal relationships, but possible factors have been provided for further causal investigation. Although we considered possible confounders in our statistical adjustment, there might be residual confounders that we could not capture. Third, given that our study interest was dietary fibre/soy intake at baseline, we were unable to detect eating behaviour changes during the study follow-ups, which might confound these associations. In addition, because soy and its products are widely and substantially consumed in Asian diets, isoflavone intake among Asians is considerably higher than that in Western populations , so that generalization of our findings might be limited in low-soy-consumption populations. Furthermore, urinary equol level might merely represent the quantity of excreted equol in the urine, rather than one’s ability to produce or utilize equol. The soy challenging test could provide valuable information about equol-production status in future cohorts.
In conclusion, higher childhood soy intake is prospectively associated with later puberty timing in both Chinese girls and boys, independent of prepubertal body fatness, and is particularly pronounced among individuals with a higher urinary equol level.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its supplementary files.
The Chinese Adolescent Cohort Study
Age at Tanner stage 2 for breast development or age at the initiation of gonadal growth
Age at menarche or voice break
Food frequency questionnaire
Percent body fat
Body mass index
Day F, Thompson D, Helgason H, Chasman D, Finucane H, Sulem P, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49(6):834–41. https://doi.org/10.1038/ng.3841.
Cancer CGoHFiB. Menarche, menopause, and breast cancer risk: individual participant meta-analysis, including 118 964 women with breast cancer from 117 epidemiological studies. Lancet Oncol. 2012;13(11):1141–51.
Beszterda M, Frański R. Endocrine disruptor compounds in environment: as a danger for children health. Pediatr Endocrinol Diabetes Metab. 2018;24(2):88–95 https://doi.org/10.18544/pedm-24.02.0107.
Cheng G, Buyken A, Shi L, Karaolis-Danckert N, Kroke A, Wudy S, et al. Beyond overweight: nutrition as an important lifestyle factor influencing timing of puberty. Nutr Rev. 2012;70(3):133–52. https://doi.org/10.1111/j.1753-4887.2011.00461.x.
Adlercreutz H. Phytoestrogens and breast cancer. J Steroid Biochem Mol Biol. 2002;83:113–8. https://doi.org/10.1016/s0960-0760(02)00273-x.
Knight DC, Eden JA. Phytoestrogens--a short review. Maturitas. 1995;22(3):167–75. https://doi.org/10.1016/0378-5122(95)00937-g.
Cheng G, Remer T, Prinz-Langenohl R, Blaszkewicz M, Degen G, Buyken A. Relation of isoflavones and fiber intake in childhood to the timing of puberty. Am J Clin Nutr. 2010;92(3):556–64. https://doi.org/10.3945/ajcn.2010.29394.
Kissinger DG, Sanchez A. The association of dietary factors with the age of menarche. Nutr Res. 1987;7(5):471–9.
Wolff M, Teitelbaum S, McGovern K, Pinney S, Windham G, Galvez M, et al. Environmental phenols and pubertal development in girls. Environ Int. 2015;84:174–80. https://doi.org/10.1016/j.envint.2015.08.008.
Wolff M, Teitelbaum S, Pinney S, Windham G, Liao L, Biro F, et al. Investigation of relationships between urinary biomarkers of phytoestrogens, phthalates, and phenols and pubertal stages in girls. Environ Health Perspect. 2010;118(7):1039–46. https://doi.org/10.1289/ehp.0901690.
de Ridder CM, Thijssen JH, Van’t Veer P, van Duuren R, Bruning PF, Zonderland ML, et al. Dietary habits, sexual maturation, and plasma hormones in pubertal girls: a longitudinal study. Am J Clin Nutr. 1991;54(5):805–13. https://doi.org/10.1093/ajcn/54.5.805.
Koo MM, Rohan TE, Jain M, McLaughlin JR, Corey PN. A cohort study of dietary fibre intake and menarche. Public Health Nutr. 2002;5(2):353–60. https://doi.org/10.1079/PHN2002261.
Maclure M, Travis LB, Willett W, MacMahon B. A prospective cohort study of nutrient intake and age at menarche. Am J Clin Nutr. 1991;54(4):649–56. https://doi.org/10.1093/ajcn/54.4.649.
Nielsen I, Williamson G. Review of the factors affecting bioavailability of soy isoflavones in humans. Nutr Cancer. 2007;57(1):1–10. https://doi.org/10.1080/01635580701267677.
Zaheer K, Humayoun Akhtar M. An updated review of dietary isoflavones: nutrition, processing, bioavailability and impacts on human health. Crit Rev Food Sci Nutr. 2017;57(6):1280–93. https://doi.org/10.1080/10408398.2014.989958.
Mayo B, Vázquez L, Flórez A. Equol: a bacterial metabolite from the daidzein isoflavone and its presumed beneficial health effects. Nutrients. 2019;11(9). https://doi.org/10.3390/nu11092231.
Franke A, Lai J, Halm B. Absorption, distribution, metabolism, and excretion of isoflavonoids after soy intake. Arch Biochem Biophys. 2014;559:24–8. https://doi.org/10.1016/j.abb.2014.06.007.
Setchell K, Cole S: Method of defining equol-producer status and its frequency among vegetarians. J Nutr 2006, 136(8):2188-2193. http://doi.org/https://doi.org/10.1093/jn/136.8.2188
Liu B, Qin L, Liu A, Uchiyama S, Ueno T, Li X, et al. Prevalence of the equol-producer phenotype and its relationship with dietary isoflavone and serum lipids in healthy Chinese adults. J Epidemiol. 2010;20(5):377–84. https://doi.org/10.2188/jea.je20090185.
Ma HM, Chen SK, Chen RM, Zhu C, Xiong F, Li T, et al. Pubertal development timing in urban Chinese boys. Int J Androl. 2011;34(5 Pt 2):e435–45. https://doi.org/10.1111/j.1365-2605.2011.01173.x.
Song Y, Ma J, Wang HJ, Wang Z, Lau PW, Agardh A. Age at spermarche: 15-year trend and its association with body mass index in Chinese school-aged boys. Pediatr Obes. 2016;11(5):369–74. https://doi.org/10.1111/ijpo.12073.
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. https://doi.org/10.1542/peds.2008-2638.
Lei Y, Luo D, Yan X, Zhang J, Hu P, Ma J, et al. The mean age of menarche among Chinese schoolgirls declined by 6 months from 2005 to 2014. Acta Paediatr. 2021;110(2):549–55. https://doi.org/10.1111/apa.15441.
Liu M, Liu SW, Wang LJ, Bai YM, Zeng XY, Guo HB, et al. Burden of diabetes, hyperglycaemia in China from to 2016: Findings from the 1990 to 2016, global burden of disease study. Diabetes Metab. 2019;45(3):286–93. https://doi.org/10.1016/j.diabet.2018.08.008.
Liu S, Li Y, Zeng X, Wang H, Yin P, Wang L, et al. Burden of cardiovascular diseases in China, 1990-2016: findings from the 2016 Global Burden of Disease Study. JAMA Cardiol. 2019;4(4):342–52. https://doi.org/10.1001/jamacardio.2019.0295.
Sun D, Li H, Cao M, He S, Lei L, Peng J, et al. Cancer burden in China: trends, risk factors and prevention. Cancer Biol Med. 2020;17(4):879–95 https://doi.org/10.20892/j.issn.2095-3941.2020.0387.
Duan R, Wang X, Shan S, Zhao L, Xiong J, Libuda L, et al. The Chinese Adolescent Cohort Study: design, implementation, and major findings. Front Nutr. 2021;8(904). https://doi.org/10.3389/fnut.2021.747088.
Duan R, Qiao T, Chen Y, Chen M, Xue H, Zhou X, et al. The overall diet quality in childhood is prospectively associated with the timing of puberty. Eur J Nutr. 2021;60(5):2423–34. https://doi.org/10.1007/s00394-020-02425-8.
Zhang X, Chen M, Duan R, Xue H, Luo J, Lv X, et al. The Nutrition and Health in Southwest China (NHSC) study: design, implementation, and major findings. Eur J Clin Nutr. 2021;75(2):299–306. https://doi.org/10.1038/s41430-020-00703-6.
Li H, Ji CY, Zong XN, Zhang YQ. Body mass index growth curves for Chinese children and adolescents aged 0 to 18 years. Zhonghua Er Ke Za Zhi. 2009;47(7):493.
Cole T, Bellizzi M, Flegal K, Dietz W. Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ (Clinical research ed). 2000;320(7244):1240–3. https://doi.org/10.1136/bmj.320.7244.1240.
Slaughter M, Lohman T, Boileau R, Horswill C, Stillman R, Van Loan M, et al. Skinfold equations for estimation of body fatness in children and youth. Hum Biol. 1988;60(5):709–23.
Matthews CE, Shu XO, Yang G, Jin F, Ainsworth BE, Liu D, et al. Reproducibility and validity of the Shanghai Women's Health Study physical activity questionnaire. Am J Epidemiol. 2003;158(11):1114–22. https://doi.org/10.1093/aje/kwg255.
Brook C, Grumbach MM. Clinical paediatric endocrinology. London: Blackwell Publishing Ltd; 1981.
Meyer F, Moisan J, Marcoux D, Bouchard C. Dietary and physical determinants of menarche. Epidemiology. 1990;1(5):377–81. https://doi.org/10.1097/00001648-199009000-00007.
Textor J, van der Zander B, Gilthorpe MS, Liskiewicz M, Ellison GT. Robust causal inference using directed acyclic graphs: the R package 'dagitty'. Int J Epidemiol. 2016;45(6):1887–94. https://doi.org/10.1093/ije/dyw341.
Maldonado G, Greenland S. Simulation study of confounder-selection strategies. Am J Epidemiol. 1993;138(11):923–36. https://doi.org/10.1093/oxfordjournals.aje.a116813.
Kaplowitz PB. Link between body fat and the timing of puberty. Pediatrics. 2008;121(Suppl 3):S208–17. https://doi.org/10.1542/peds.2007-1813F.
Harrell FE Jr. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis: springer; 2015.
Yuexin Y, Peking NFS, Guangya W, Xingchang P. China food composition. 1st ed; 2002.
China NBoSo. China Statistical Yearbook. Beijing: China Statistics Press; 2020.
Wolff M, Britton J, Boguski L, Hochman S, Maloney N, Serra N, et al. Environmental exposures and puberty in inner-city girls. Environ Res. 2008;107(3):393–400. https://doi.org/10.1016/j.envres.2008.03.006.
Segovia-Siapco G, Pribis P, Oda K, Sabaté J. Soy isoflavone consumption and age at pubarche in adolescent males. Eur J Nutr. 2018;57(6):2287–94. https://doi.org/10.1007/s00394-017-1504-1.
Nielsen AJ, McNulty J. Polyphenolic natural products and natural product-inspired steroidal mimics as aromatase inhibitors. Med Res Rev. 2019;39(4):1274–93. https://doi.org/10.1002/med.21536.
Barnes S, Kim H, Darley-Usmar V, Patel R, Xu J, Boersma B, et al. Beyond ERalpha and ERbeta: estrogen receptor binding is only part of the isoflavone story. J Nutr. 2000;130(3):656S–7S. https://doi.org/10.1093/jn/130.3.656S.
Arispe S, Adams B, Adams T. Effect of phytoestrogens on basal and GnRH-induced gonadotropin secretion. J Endocrinol. 2013;219(3):243–50. https://doi.org/10.1530/joe-13-0158.
Sultan C, Gaspari L, Maimoun L, Kalfa N, Paris F. Disorders of puberty. Best Pract Res Clin Obstet Gynaecol. 2018;48:62–89. https://doi.org/10.1016/j.bpobgyn.2017.11.004.
Sonnenburg E, Sonnenburg J. Starving our microbial self: the deleterious consequences of a diet deficient in microbiota-accessible carbohydrates. Cell Metab. 2014;20(5):779–86. https://doi.org/10.1016/j.cmet.2014.07.003.
O'Keeffe L, Frysz M, Bell J, Howe L, Fraser A. Puberty timing and adiposity change across childhood and adolescence: disentangling cause and consequence. Hum Reprod (Oxford, England). 2020;35(12):2784–92. https://doi.org/10.1093/humrep/deaa213.
The participation of all of the children and their parents is gratefully acknowledged.
The study was supported by a research grant from the National Key R&D Program of China (2020YFC2006300). The funders of this study played no role in the study design, collection, analysis, or interpretation of data, in the writing of the report, or in the decision to submit the paper for publication.
Ethics approval and consent to participate
The CAC cohort was approved by the Ethics Committee of West China School of Public Health of Sichuan University. All of the parents of the participants provided their written confirmed consent before enrolment.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flowchart for the study sample. Figure S2. Shape of the associations between total soy intakes and fibre intakes in childhood with puberty timing. In the analyses, 2152 girls and 2629 boys were included. B2, Tanner stage 2 for breast development; G2, the initiation of gonadal growth.
Associations of total soy intakes in childhood with puberty timing using mixed model. Table S2. Associations of total soy intakes in childhood with puberty timing among participants with urinary equol data using mixed model. Table S3. Associations of total soy intakes in childhood with puberty timing among participants with urinary equol data. Table S4. Associations of total dietary fibre intakes in childhood with puberty timing using mixed model. Table S5. Associations of cereal fibre intakes in childhood with puberty timing. Table S6. Associations of vegetable and fruit fibre intakes in childhood with puberty timing.
About this article
Cite this article
Xiong, J., Xu, Y., Liu, X. et al. Prospective association of dietary soy and fibre intake with puberty timing: a cohort study among Chinese children. BMC Med 20, 145 (2022). https://doi.org/10.1186/s12916-022-02320-5