Joint models reveal genetic architecture of pubertal stage transitions and their association with BMI in admixed Chilean population.

Early or late pubertal onset can lead to disease in adulthood, including cancer, obesity, type 2 diabetes, metabolic disorders, bone fractures, and psychopathologies. Thus, knowing the age at which puberty is attained is crucial as it can serve as a risk factor for future diseases. Pubertal development is divided into five stages of sexual maturation in boys and girls according to the standardized Tanner scale. We performed genome-wide association studies (GWAS) on the "Growth and Obesity Chilean Cohort Study" cohort composed of admixed children with mainly European and Native American ancestry. Using joint models that integrate time-to-event data with longitudinal trajectories of body mass index (BMI), we identified genetic variants associated with phenotypic transitions between pairs of Tanner stages. We identified $42$ novel significant associations, most of them in boys. The GWAS on Tanner $3\rightarrow 4$ transition in boys captured an association peak around the growth-related genes LARS2 and LIMD1 genes, the former of which causes ovarian dysfunction when mutated. The associated variants are expression and splicing Quantitative Trait Loci regulating gene expression and alternative splicing in multiple tissues. Further, higher individual Native American genetic ancestry proportions predicted a significantly earlier puberty onset in boys but not in girls. Finally, the joint models identified a longitudinal BMI parameter significantly associated with several Tanner stages' transitions, confirming the association of BMI with pubertal timing.

Association between plasma leptin/adiponectin ratio and insulin resistance indexes in prepubertal children.

Objective: To assess the association between leptin/adiponectin ratio (LAR) and insulin resistance surrogates in prepubertal children. Materials and methods: Study based on data from the Growth and Obesity Chilean Cohort Study (GOCS) involving 968 Chilean prepubertal children. Plasma insulin, leptin, and adiponectin were determined by immunoassays. Several common insulin resistance surrogates were calculated, including the homeostasis model assessment of insulin resistance (HOMA-IR), triglyceride/HDL cholesterol index, triglyceride-glucose (TyG) index, and the TyG index corrected for body mass index (BMI; TyG-BMI) and waist circumference (WC; TyG-WC). Associations among variables were assessed using multiple linear and logistic regression analysis. Results: There was a significant direct association between plasma leptin and LAR with BMI z-score but no association between plasma adiponectin and adiposity. After adjustments for sex and age, LAR was significantly associated with all insulin resistance surrogates (which were categorized using the 75th percentile as the cutoff point), with the TyG-WC index emerging as the surrogate with the highest magnitude of association (odds ratio [OR] 2.44, 95% confidence interval [CI] 2.05-2.9). After additional adjustment for BMI z-score, only the association between LAR and TyG-WC remained significant (OR 1.64, 95% CI 1.27-2.12). Conclusion: Plasma leptin and LAR were strongly associated with several common insulin resistance surrogates in prepubertal children, most notably with the TyG-WC index. Associations between LAR and insulin resistance indexes were mainly driven by the effect of plasma leptin, which is also directly associated with increased adiposity.

Association between plasma leptin/adiponectin ratio and insulin resistance indexes in prepubertal children

ABSTRACT Objective: To assess the association between leptin/adiponectin ratio (LAR) and insulin resistance surrogates in prepubertal children. Subjects and methods: Study based on data from the Growth and Obesity Chilean Cohort Study (GOCS) involving 968 Chilean prepubertal children. Plasma insulin, leptin, and adiponectin were determined by immunoassays. Several common insulin resistance surrogates were calculated, including the homeostasis model assessment of insulin resistance (HOMA-IR), triglyceride/HDL cholesterol index, triglyceride-glucose (TyG) index, and the TyG index corrected for body mass index (BMI; TyG-BMI) and waist circumference (WC; TyG-WC). Associations among variables were assessed using multiple linear and logistic regression analysis. Results: There was a significant direct association between plasma leptin and LAR with BMI z-score but no association between plasma adiponectin and adiposity. After adjustments for sex and age, LAR was significantly associated with all insulin resistance surrogates (which were categorized using the 75th percentile as the cutoff point), with the TyG-WC index emerging as the surrogate with the highest magnitude of association (odds ratio [OR] 2.44, 95% confidence interval [CI] 2.05-2.9). After additional adjustment for BMI z-score, only the association between LAR and TyG-WC remained significant (OR 1.64, 95% CI 1.27-2.12). Conclusion: Plasma leptin and LAR were strongly associated with several common insulin resistance surrogates in prepubertal children, most notably with the TyG-WC index. Associations between LAR and insulin resistance indexes were mainly driven by the effect of plasma leptin, which is also directly associated with increased adiposity.

Identifying patient-related predictors of permanent growth hormone deficiency.

Objective: Isolated childhood growth hormone deficiency (GHD) can persist into adulthood, and re-testing at the transition period is needed to determine whether continued growth hormone therapy is indicated. Here, our objective was to identify predictors of permanent GHD. Design: Retrospective single-centre study of patients with childhood-onset GHD who were re-tested after adult height attainment. Methods: Auxological, clinical, laboratory, and MRI data throughout follow-up were collected. Results: We included 101 patients. At GH treatment initiation, age was 8.1 ± 0.4 years, height -2.25 ± 0.8, and BMI -0.27 ± 0.1 SDS. The 29 (28.7%) patients with persistent GHD had lower height SDS (-2.57 ± 0.1 vs. -2.11 ± 0.1, p<0.001) and mean GH peaks (8.4 ± 1.0 vs.13.2 ± 0.5 mIU/L, p<0.001) at GHD diagnosis; at adult height, they had lower IGF1 (232 ± 19.9 vs. 331 ± 9.1 ng/mL, p<0.001) and higher BMI SDS (-0.15 ± 0.27 vs. -0.73 ± 0.13, p<0.005). By multivariate analysis, the best predictive model included height and BMI SDS, both GH peaks, and MRI findings at diagnosis. Patients with height at diagnosis <-3 SDS had a 7.7 (95% IC 1.4-43.1, p=0.02) fold higher risk of persistent GHD after adjustment on BMI SDS. An abnormal pituitary region by MRI was the strongest single predictor (7.2 times, 95% CI 2.7-19.8) and after multivariate analysis adjustment for GH peaks and height SDS at diagnosis, the risk increased to 10.6 (1.8 - 61.3) times. Conclusions: Height <-3 SDS at GHD diagnosis and pituitary MRI abnormalities should lead to a high index of suspicion for persistent GHD.

Association Between Markers of Adiposity During Childhood and Puberty Onset in Latino Girls.

CONTEXT: Prepubertal adiposity is associated with earlier puberty. It is unclear when this association starts, if all adiposity markers are similarly associated, and whether all pubertal milestones are similarly affected. OBJECTIVE: To evaluate the association between different adiposity markers during childhood and the timing of different pubertal milestones in Latino girls. DESIGN, SETTING, AND PARTICIPANTS: Longitudinal follow-up of 539 female participants of the Chilean Growth and Obesity Cohort recruited from childcare centers (mean age 3.5 years) from the southeast area of Santiago, Chile. Participants were singletons born between 2002 and 2003 within the normal birthweight range. Since 2006, a trained dietitian measured weight, height, waist circumference (WC) and skinfolds to estimate body mass index (BMI) Centers for Disease Control and Prevention percentiles, central obesity, percentage of fat mass (%FM), and fat mass index (FMI, fat mass/height2). MAIN OUTCOME: Since 2009, sexual maturation was assessed every 6 months to assess age at (1) thelarche, (2) pubarche, (3) menarche, and (4) peak height velocity (PHV). RESULTS: At thelarche, 12.5% were obese and 2% had central obesity. The median age of pubarche, menarche, and PHV were all associated with markers of adiposity at different time points during childhood whereas thelarche only with %FM and FMI. Adiposity clusters models showed that children with trajectories of high WC, %FM, and FMI during childhood were related with earlier thelarche, pubarche, menarche, and PHV but BMI trajectories only with menarche and PHV. CONCLUSIONS: Higher WC, %FM, and FMI were associated with earlier age at thelarche, pubarche, menarche, and PHV. The effect of BMI was less consistent.

Differential methylation pattern in pubertal girls associated with biochemical premature adrenarche.

Biochemical premature adrenarche is defined by elevated serum DHEAS [≥40 µg/dL] before age 8 y in girls. This condition is receiving more attention due to its association with obesity, hyperinsulinemia, dyslipidemia, and polycystic ovary syndrome. Nevertheless, the link between early androgen excess and these risk factors remains unknown. Epigenetic modifications, and specifically DNA methylation, have been associated with the initiation and progression of numerous disorders, including obesity and insulin resistance. The aim of this study was to determine if prepubertal androgen exposure is associated with a different methylation profile in pubertal girls. Eighty-six healthy girls were studied. At age 7 y, anthropometric measurements were begun and DHEAS levels were determined. Girls were classified into Low DHEAS (LD) [<42 µg/dL] and High DHEAS (HD) [≥42 µg/dL] groups. At Tanner stages 2 and 4 a DNA methylation microarray was performed to identify differentially methylated CpG positions (DMPs) between HD and LD groups. We observed a differential methylation pattern between pubertal girls with and without biochemical PA. Moreover, a set of DNA methylation markers, selected by the LASSO method, successfully distinguished between HD and LD girls regardless of Tanner stage. Additionally, a subset of these markers were significantly associated with glucose-related measures such as insulin level, HOMA-IR, and glycaemia. This pilot study provides evidence consistent with the hypothesis that high DHEAS concentration, or its hormonally active metabolites, may induce a unique blood methylation signature in pubertal girls, and that this methylation pattern is associated with altered glucose metabolism.

New insights from GWAS on BMI-related growth traits in a longitudinal cohort of admixed children with Native American and European ancestry.

Body-mass index (BMI) is a hallmark of adiposity. In contrast with adulthood, the genetic architecture of BMI during childhood is poorly understood. The few genome-wide association studies (GWAS) on children have been performed almost exclusively in Europeans and at single ages. We performed cross-sectional and longitudinal GWAS for BMI-related traits on 904 admixed children with mostly Mapuche Native American and European ancestries. We found regulatory variants of the immune gene HLA-DQB3 strongly associated with BMI at 1.5 - 2.5 years old. A variant in the sex-determining gene DMRT1 was associated with the age at adiposity rebound (Age-AR) in girls (P = 9.8 × 10 - 9 ). BMI was significantly higher in Mapuche than in Europeans between 5.5 and 16.5 years old. Finally, Age-AR was significantly lower (P = 0.004 ) by 1.94 years and BMI at AR was significantly higher (P = 0.04 ) by 1.2 kg/ m 2 , in Mapuche children compared with Europeans.

Higher Prepubertal IGF-1 Concentrations Associate to Earlier Pubertal Tempo in Both Sexes.

INTRODUCTION: Pubertal onset is triggered by multiple neuroendocrine interactions. The role of prepubertal IGF-1 in this process has not been explored in both sexes. Our objective was to analyze the association of prepubertal IGF-1 concentration with age at thelarche (B2) and menarche (M) in girls and age at gonadarche (G2) in boys. METHODS: This is a longitudinal study (n = 1,196 boys and girls) within the Growth and Obesity Chilean Cohort Study (GOCS). At age ≈ 6.7 years, blood sample was taken for IGF-1. Subjects were divided into 4 groups according to the onset age of the pubertal event. RESULTS: Higher prepubertal IGF-1 levels were observed at earlier ages of B2 (p = 0.003) and M onset (p = 0.041). A taller prepubertal height was observed at younger ages of B2 and M (p=<0.001 and 0.002, respectively). The hazard proportional regression models (HR) showed that with an increase of 1 SD in IGF-1, the HR of presenting B2 at younger ages was 1.25, and this association was maintained when adjusted for confounding variables. Similarly, the HR of presenting M at earlier ages was 1.21. This association was maintained only when adjusting for body mass index but not using further confounders. In boys, prepubertal IGF-1 showed a tendency to be significantly higher in children with earlier G2 and taller height (both p < 0.001). The HR of presenting G2 at younger ages was 1.22, and this association was maintained after adjusting for confounders. CONCLUSIONS: Higher IGF-1 levels in mid-childhood are associated with earlier puberty onset. The role of IGF-1 in the onset of puberty requires further investigation.

Assessment of a consensus definition of obesity and metabolic health phenotypes in children at different pubertal stages.

Not all individuals with obesity develop metabolic complications, which has brought about the concepts of metabolically healthy and unhealthy obesity (MHO/MUO). However, inconsistent definitions of these conditions have limited their understanding. We assessed whether a recently-proposed consensus definition for MHO/MUO correlates with adiposity and reflects metabolic risk parameters during puberty. Low-middle income children from the Growth and Obesity Cohort Study (Santiago, Chile) were included (n = 949; 1692 visits at Tanner (T)2, T4 and/or one-year post menarche (1YPM)). Anthropometry, body composition and metabolic parameters were compared between MUO and MHO, and also in children without obesity. The risk for presenting MUO phenotype was significantly elevated with higher waist-height ratio (T2), zBMI (T2, T4), trunk fat, and C-reactive protein (T4). Elevated cardiometabolic indices were important predictors of the "unhealthy" phenotype allocation in children with or without obesity. Our observations suggest that the consensus definition in children at T2, T4 and 1YPM reflects metabolic risk and central obesity. Metabolic health phenotype allocation by this equation enables easy detection of risk factors that call for action to prevent long-term metabolic derangements in children with obesity and, importantly, also those without obesity.

Habitual Phytoestrogen Intake Is Associated with Breast Composition in Girls at 2 Years after Menarche Onset.

BACKGROUND: High phytoestrogen intake during adolescence is associated with a reduced risk of breast cancer. Breast density (BD) is a strong predictor of breast cancer and can be considered an early marker. We aim to assess the association between the mean habitual intake of isoflavones, lignans, and total phytoestrogens intake during puberty until 2 years after menarche onset and absolute fibroglandular volume (AFGV) and percentage of fibroglandular volume (%FGV) in Hispanic girls at the end of puberty. METHODS: Longitudinal study set up in the Growth and Obesity Chilean Cohort Study (GOCS). We included 329 girls with dietary data (multiple 24-hours recalls) from puberty until 2 years after menarche onset (81% had 2-4 recalls). Two international datasets were used to estimate isoflavones, lignans, and total phytoestrogens in the diet. Breast composition was measured by dual energy X-ray absorptiometry at 2 years after menarche. Multiple linear regression models were used to assess the association between isoflavones, lignans, and total phytoestrogens intake and AFGV and %FGV. RESULTS: The average total phytoestrogen intake was 1 mg/day and %FGV was 50.7% (SD = 15.2) and AFGV 218.8 cm3 (SD = 79.3). An inverse association was found between consumption of isoflavones and AFGV, as well as, with total phytoestrogens [Q4 vs. Q1 adjusted model ß = -49.2 cm3; 95% CI (-85.5 to -13.0)]. CONCLUSIONS: Girls with a higher intake of total phytoestrogens and isoflavones during puberty until 2 years after menarche onset had significantly lower AFGV. IMPACT: Although the intake of phytoestrogens is low in Western populations, higher consumption of them during a critical period of life like puberty could be beneficial to reduce breast cancer during adulthood.

Novel OTX2 loss of function variant associated with congenital hypopituitarism without eye abnormalities.

OBJECTIVES: The normal development of the pituitary gland requires multiple induction signals and transcription factors encoded by more than 30 genes, including OTX2. OTX2 mutations have been described with eye abnormalities and variable congenital hypopituitarism, but rarely with hypopituitarism without ocular manifestations. CASE PRESENTATION: We report a girl with hypopituitarism associated with pituitary hypoplasia and pituitary stalk atrophy, without ocular manifestations. NGS revealed a novel heterozygous mutation in OTX2 c.426dupC:p.(Ser143Leufs*2). CONCLUSIONS: Mutations in the transcription factor OTX2 have been associated with ocular, craniofacial, and pituitary development anomalies. Here we describe a novel mutation in OTX2 associated with hypopituitarism without an ocular phenotype.

Biallelic POC1A variants cause syndromic severe insulin resistance with muscle cramps.

Objective: To describe clinical, laboratory, and genetic characteristics of three unrelated cases from Chile, Portugal, and Saudi Arabia with severe insulin resistance, SOFT syndrome, and biallelic pathogenic POC1A variants. Design: Observational study. Methods: Probands' phenotypes, including short stature, dysmorphism, and insulin resistance, were compared with previous reports. Results: Cases 1 (female) and 3 (male) were homozygous for known pathogenic POC1A variants: c.649C>T, p.(Arg217Trp) and c.241C>T, p.(Arg81*), respectively. Case 2 (male) was compound heterozygous for p.(Arg217Trp) variant and the rare missense variant c.370G>A, p.(Asp124Asn). All three cases exhibited severe insulin resistance, acanthosis nigricans, elevated serum triglycerides and decreased HDL, and fatty liver, resembling three previously reported cases. All three also reported severe muscle cramps. Aggregate analysis of the six known cases with biallelic POC1A variants and insulin resistance showed decreased birth weight and length mean (s.d.): -2.8 (0.9) and -3.7 (0.9) SDS, respectively), severe short stature mean (s.d.) height: -4.9 (1.7) SDS) and moderate microcephaly (mean occipitofrontal circumference -3.0 (range: -4.7 to -1.2)). These findings were similar to those reported for patients with SOFT syndrome without insulin resistance. Muscle biopsy in Case 3 showed features of muscle involvement secondary to a neuropathic process. Conclusions: Patients with SOFT syndrome can develop severe dyslipidaemic insulin resistance, independent of the exonic position of the POC1A variant. They also can develop severe muscle cramps. After diagnosis, patients should be regularly screened for insulin resistance and muscle complaints.

Genome-Wide Association Study and Polygenic Risk Scores of Serum DHEAS Levels in a Chilean Children Cohort.

CONTEXT: Adrenarche reflects the developmental growth of the adrenal zona reticularis, which produces increasing adrenal androgen secretion (eg, dehydroepiandrosterone [DHEA]/dehydroepiandrosterone sulfate [DHEAS]) from approximately age 5 to 15 years. OBJECTIVE: We hypothesized that the study of the genetic determinants associated with variations in serum DHEAS during adrenarche might detect genetic variants influencing the rate or timing of this process. METHODS: Genome-wide genotyping was performed in participants of the Chilean pediatric Growth and Obesity Chilean Cohort Study (GOCS) cohort (n = 788). We evaluated the genetic determinants of DHEAS levels at the genome-wide level and in targeted genes associated with steroidogenesis. To corroborate our findings, we evaluated a polygenic risk score (PRS) for age at pubarche, based on the discovered variants, in children from the same cohort. RESULTS: We identified one significant variant at the genome-wide level in the full cohort, close to the GALR1 gene (P = 3.81 × 10-8). In addition, variants suggestive of association (P < 1 × 10-5) were observed in PRLR, PITX1, PTPRD, NR1H4, and BCL11B. Stratifying by sex, we found variants suggestive of association in SERBP1 and CAMTA1/VAMP3 for boys and near ZNF98, TRPC6, and SULT2A1 for girls. We also found significant reductions in age at pubarche in those children with higher PRS for greater DHEAS based on these newly identified variants. CONCLUSION: Our results disclose one variant associated with DHEAS concentrations at the level of genome-wide association study significance, and several variants with a suggestive association that may be involved in the genetic regulation of adrenarche.

High DHEAS in girls and metabolic features throughout pubertal maturation.

CONTEXT: An association between premature adrenarche and metabolic syndrome at presentation has been described. Our aim was to assess whether the presence of high dehydroepiandrosterone sulphate (DHEAS [HD]) at the adrenarche determines the risk of metabolic syndrome during puberty, taking into account body mass index (BMI) and birth weight. DESIGN: Prospective observational. PATIENTS: Five hundred four girls from the Growth and Obesity Chilean Cohort Study were followed from birth through puberty. At age ~7, subjects were classified by DHEAS concentrations into the HD (>75th percentile) or normal DHEAS (ND, ≤75th percentile) subgroups. MEASUREMENTS: Anthropometrics, semiannual clinical pubertal staging and hormonal and metabolic levels. The relationships among DHEAS at age ~7, metabolic syndrome, and each of its components independently, were analyzed by linear and logistic regression models during puberty and 1-year postmenarche, adjusted by confounders. RESULTS: Girls with HD at 7 years exhibited higher BMI, more central fat and higher serum androgen and insulin like growth factor (IGF)-I levels throughout puberty. Also, girls with HD had a greater prevalence of hyperglycemia at B2 and B4 breast stages, and of low HDL at B4. At 1 year after menarche, HD girls had a higher prevalence of metabolic syndrome, and those with BMI > 1 SD score had a higher metabolic score and insulin levels than ND girls with similar BMI. CONCLUSIONS: Our observations suggest that girls with HD at the age of adrenarche may be at greater risk for metabolic syndrome at adolescence, especially in those who are overweight or obese. Our results emphasize the importance of lifestyle interventions for childhood overweight and obesity among girls with HD.

Primary Ovarian Failure in Addition to Classical Clinical Features of Coats Plus Syndrome in a Female Carrying 2 Truncating Variants of CTC1.

Coats plus syndrome is an autosomal recessive multisystemic and pleiotropic disorder affecting the eyes, brain, bone, and gastrointestinal tract, usually caused by compound heterozygous variants of the conserved telomere maintenance component 1 gene (CTC1), involved in telomere homeostasis and replication. So far, most reported patients are compound heterozygous for a truncating mutation and a missense variant. The phenotype is believed to result from telomere dysfunction, with accumulation of DNA damage, cellular senescence, and stem cell depletion. Here, we report a 23-year-old female with prenatal and postnatal growth retardation, microcephaly, osteopenia, recurrent fractures, intracranial calcification, leukodystrophy, parenchymal brain cysts, bicuspid aortic valve, and primary ovarian failure. She carries a previously reported maternally inherited pathogenic variant in exon 5 (c.724_727del, p.(Lys242Leufs*41)) and a novel, paternally inherited splice site variant (c.1617+5G>T; p.(Lys480Asnfs*17)) in intron 9. CTC1 transcript analysis showed that the latter resulted in skipping of exon 9. A trace of transcripts was normally spliced resulting in the presence of a low level of wild-type CTC1 transcripts. We speculate that ovarian failure is caused by telomere shortening or chromosome cohesion failure in oocytes and granulosa cells, with early decrease in follicular reserve. This is the first patient carrying 2 truncating CTC1 variants and the first presenting primary ovarian failure.

Novel loci and Mapuche genetic ancestry are associated with pubertal growth traits in Chilean boys.

Puberty is a complex developmental process that varies considerably among individuals and populations. Genetic factors explain a large proportion of the variability of several pubertal traits. Recent genome-wide association studies (GWAS) have identified hundreds of variants involved in traits that result from body growth, like adult height. However, they do not capture many genetic loci involved in growth changes over distinct growth phases. Further, such GWAS have been mostly performed in Europeans, but it is unknown how these findings relate to other continental populations. In this study, we analyzed the genetic basis of three pubertal traits; namely, peak height velocity (PV), age at PV (APV) and height at APV (HAPV). We analyzed a cohort of 904 admixed Chilean children and adolescents with European and Mapuche Native American ancestries. Height was measured on roughly a [Formula: see text]month basis from childhood to adolescence between 2006 and 2019. We predict that, in average, HAPV is 4.3 cm higher in European than in Mapuche adolescents (P = 0.042), and APV is 0.73 years later in European compared with Mapuche adolescents (P = 0.023). Further, by performing a GWAS on 774, 433 single-nucleotide polymorphisms, we identified a genetic signal harboring 3 linked variants significantly associated with PV in boys (P [Formula: see text]). This signal has never been associated with growth-related traits.

Total and Central Adiposity Are Associated With Age at Gonadarche and Incidence of Precocious Gonadarche in Boys.

CONTEXT: A close link between body mass index (BMI) and female pubertal onset is well established. However, observations in boys remain inconclusive. OBJECTIVE: We aim to determine whether BMI as well as total and central adiposity in prepubertal Chilean boys is associated with pubertal timing. METHODS: We performed a longitudinal study in which 494 boys from the Growth and Obesity Chilean Cohort Study were followed starting from birth and throughout puberty, including 5 prepubertal visits. The main outcome measures included anthropometric data and semi-annual clinical pubertal staging. The association between BMI, obesity (BMI standard deviation score [SDS] ≥ 2) and central adiposity (waist circumference ≥ 90th centile) with precocious puberty and age at gonadarche was analyzed using survival- and logistic regression models. RESULTS: BMI, prevalence of total obesity, and central obesity increased throughout childhood. Among the study population, 45 boys entered puberty before the age of 9 years (9.1%). Obesity at 4 to 7 years and childhood mean BMI SDS were significantly associated with precocious gonadarche. Mean age at testicular enlargement (≥4 mL), was 11.0 years (95% CI, 10.9-11.1) and was inversely associated with BMI SDS, waist circumference, and percentage fat mass in almost all prepubertal visits. Age at testicular enlargement in normal weight, overweight, and obese boys was 11.2 (11.0-11.3), 10.9 (10.6-11.1) and 10.7 (10.4-11.1) years, respectively. CONCLUSION: Our observation supports the association of BMI SDS and obesity with pubertal timing and precocious gonadarche in boys, respectively. Early intervention controlling the obesity epidemic could be useful in decreasing detrimental impact on later health.

Reproductive hormones during pubertal transition in girls with transient Thelarche.

CONTEXT: Transient thelarche (TT), that is, the appearance, regression and subsequent reappearance of breast buds, is a frequent phenomenon, but little is known about pubertal transition in these girls. OBJECTIVE: To describe pubertal progression, growth, genotypes, reproductive hormones and growth factors in girls with TT compared to those who do not present TT (non-TT). DESIGN: Retrospective analysis of a longitudinal population-based study. PATIENTS OR OTHER PARTICIPANTS: Girls (n = 508) of the Chilean Growth and Obesity cohort. MEASUREMENTS: Pubertal progression, reproductive hormones, follicle stimulating hormone (FSH) beta subunit/FSH receptor gene single nucleotide polymorphisms and growth. RESULTS: Thirty-seven girls (7.3%) were presented TT. These girls entered puberty by pubarche more frequently (51%) than girls with normal progression (non-TT; n = 471; 23%, P = .005). Girls with TT who were under 8 years old had lower androgens, anti-Müllerian hormone (AMH), luteinizing hormone (LH) and oestradiol (all P < .05) than older girls with TT. At the time of Tanner breast stage 2 (B2), girls with TT had higher androgens, LH, FSH, IGF1, LH, insulin and oestradiol (P < .01) than at the time of TT. TT girls were older at B2 (10.3 ± 1.1 vs. 9.2 ± 1.2 years, P < .001) and menarche (12.3 ± 0.8 vs. 12.0 ± 1.0 years, P = .040) than their counterparts (non-TT). No differences in anthropometric variables or FSHB/FSHR genotypes were detected. CONCLUSION: Transient thelarche is a frequent phenomenon that does not appear to be mediated by hypothalamic-pituitary-gonadal axis activation or by adiposity. Hormonal differences between earlier TT and later TT suggest that their mechanisms are different.

A Polygenic Risk Score Suggests Shared Genetic Architecture of Voice Break With Early Markers of Pubertal Onset in Boys.

CONTEXT: Voice break, as a landmark of advanced male puberty in genome-wide association studies (GWAS), has revealed that pubertal timing is a highly polygenic trait. Although voice break is easily recorded in large cohorts, it holds quite low precision as a marker of puberty. In contrast, gonadarche and pubarche are early and clinically well-defined measures of puberty onset. OBJECTIVE: To determine whether a polygenic risk score (PRS) of alleles that confer risk for voice break associates with age at gonadarche (AAG) and age at pubarche (AAP) in Chilean boys. EXPERIMENTAL DESIGN: Longitudinal study. SUBJECTS AND METHODS: 401 boys from the Growth and Obesity Chilean Cohort Study (n = 1194; 49.2% boys). MAIN OUTCOME MEASURES: Biannual clinical pubertal staging including orchidometry. AAG and AAP were estimated by censoring methods. Genotyping was performed using the Multi-Ethnic Global Array (Illumina). Using GWAS summary statistics from the UK-Biobank, 29 significant and independent single nucleotide polymorphisms associated with age at voice break were extracted. Individual PRS were computed as the sum of risk alleles weighted by the effect size. RESULTS: The PRS was associated with AAG (ß=0.01, P = 0.04) and AAP (ß=0.185, P = 0.0004). In addition, boys within the 20% highest PRS experienced gonadarche and pubarche 0.55 and 0.67 years later than those in the lowest 20%, respectively (P = 0.013 and P = 0.007). CONCLUSIONS: Genetic variants identified in large GWAS on age at VB significantly associate with age at testicular growth and pubic hair development, suggesting that these events share a genetic architecture across ethnically distinct populations.