banner
News center
All-embracing business

Metabolic characteristics and pathogenesis of precocious puberty in girls: the role of perfluorinated compounds

Nov 06, 2023

BMC Medicine volume 21, Article number: 323 (2023) Cite this article

68 Accesses

Metrics details

Precocious puberty (PP) in girls is traditionally defined as the onset of breast development before the age of 8 years. The specific biomarkers of premature thelarche (PT) and central precocious puberty (CPP) girls are uncertain, and little is known about their metabolic characteristics driven by perfluorinated compounds (PFCs) and clinical phenotype. This study aimed to screen specific biomarkers of PT and CPP and elucidate their underlying pathogenesis. The relationships of clinical phenotype-serum PFCs-metabolic characteristics were also explored to reveal the relationship between PFCs and the occurrence and development of PT and CPP.

Nuclear magnetic resonance (NMR)-based cross-metabolomics strategy was performed on serum from 146 PP (including 30 CPP, 40 PT, and 76 unspecified PP) girls and 64 healthy girls (including 36 prepubertal and 28 adolescent). Specific biomarkers were screened by the uni- and multivariate statistical analyses. The relationships between serum PFCs and clinical phenotype were performed by correlation analysis and weighted gene co-expression network analysis to explore the link of clinical phenotype-PFCs-metabolic characteristics in PT and CPP.

The disordered trend of pyruvate and butyrate metabolisms (metabolites mapped as formate, ethanol, and 3-hydroxybutyrate) were shared and kept almost consistent in PT and CPP. Eight and eleven specific biomarkers were screened for PT and CPP, respectively. The area under curve of specific biomarker combination was 0.721 in CPP vs. prepubertal, 0.972 in PT vs. prepubertal, 0.646 in CPP vs. prepubertal integrated adolescent, and 0.822 in PT vs. prepubertal integrated adolescent, respectively. Perfluoro-n-heptanoic acid and perfluoro-n-hexanoic acid were statistically different between PT and CPP. Estradiol and prolactin were significantly correlated with PFCs in CPP and PT. Clinical phenotypes and PFCs drive the metabolic characteristics and cause metabolic disturbances in CPP and PT.

The elevation of formate, ethanol, and 3-hydroxybutyrate may serve as the early diagnostic indicator for PP in girls. But the stratification of PP still needs to be further determined based on the specific biomarkers. Specific biomarkers of CPP and PT exhibited good sensitivity and can facilitate the classification diagnosis of CPP and PT. PFC exposure is associated with endocrine homeostasis imbalance. PFC exposure and/or endocrine disturbance directly or indirectly drive metabolic changes and form overall metabolic network perturbations in CPP and PT.

Peer Review reports

Pubertal timing is usually regulated by complex interplay of genetic, environmental, nutritional, and epigenetic factors. Therefore, the criteria for normal pubertal timing and thus the definition of precocious puberty are hard to determine. Precocious puberty (PP) in girls is traditionally defined as the onset of breast development before the age of 8 years [1]. Its underlying pathophysiology may be gonadotropin-releasing hormone (GnRH)-dependent for central precocious puberty (CPP) girls or GnRH-independent for premature thelarche (PT) girls. CPP is mainly induced by the continuous pulse secretion of GnRH to prematurely activate the hypothalamic-pituitary-gonadal (HPG) axis; however, the exact mechanisms remain unclear. The main clinical manifestation of the PT girls is simple breast development due to exposure to the peripheral estrogen environment. When PT is accompanied by the significant advance growth of bone age, it is more likely to evolve into secondary CPP. CPP can lead to short- and long-term complications in girls, including increased risk of psychosocial distress, short stature, obesity, cardiovascular disease, and type 2 diabetes in adulthood [2]. Therefore, it is vital to understand the etiology of PT and CPP for accurate diagnosis and prompt intervention.

Some researchers tried to quantify PP with the help of clinical phenotype such as luteinizing hormone (LH), follicular-stimulating hormone (FSH), and estradiol to determine the index threshold for CPP diagnosis, but it is still confronted with great controversy and challenge at this moment [3, 4]. Some evidences have indicated the changes of metabolic profile during puberty. Qi et al. found that catecholamine metabolic pathway, tryptophan metabolic pathway, and TCA cycle were disturbed in CPP girls by GC/LC-MS-based urinary metabolomics analysis [5]. Yang et al. used LC-MS technology to characterize the urinary metabolomes of CPP girls and found that amino acids, especially aromatic amino acids, were closely related to the pathogenesis of CPP by activating the HPG axis and inhibiting the hypothalamic-pituitary-adrenal axis [6]. However, the clinical differential diagnosis of CPP and PT is still in a vague interface, and the lack of powerful molecular biomarkers is a long-term bottleneck in the clinical diagnosis and evaluation of PP.

Recently, ubiquitous exposures to polyfluorinated compounds (PFCs) have attracted concerns regarding their possible harmful effects during critical periods of development in early-life and long-term consequences on health in consideration of their persistence and bioaccumulation potential. Massive researches have shown that PFCs can interfere with estrogen homeostasis and pose a risk of endocrine-disrupting effects [7,8,9], and they are association with dyslipidemia, renal function, and age at menarche [7, 10,11,12], but there is still inconsistency in the research results [13] as well as certain gender differences [14,15,16,17]. In addition, evidences have shown that PFCs can affect the HPG axis [18, 19] or directly affect the gonad axis through their weak estrogen or antiandrogen effects to disrupt the development of puberty [20]. However, the correlation research of PFC exposure with the occurrence and development of PP is still in its infancy [19, 21, 22], and therefore extensive in-depth research and exploration is urgently required to clarify the exact response mechanism. The correlation analysis between PFCs and the clinical phenotype in girls with PP as well as the endogenous metabolites driven by PFCs will help to reveal the impact of PFCs on the occurrence and development of precocious puberty in girls and the preliminary mechanisms.

Based on this, the serum metabolic profiles of prepubertal, PP, PT, CPP, and adolescent girls were characterized by one-dimensional nuclear magnetic resonance hydrogen spectrum (1H-NMR) technology, and the metabolic differences and connections were analyzed by cross-metabolomics analyses, aiming to screen the specific biomarkers of CPP and PT. Furthermore, to reveal the effect of PFCs on the occurrence of CPP and PT, the metabolic modules driven by PFCs and clinical phenotype were identified by weighted gene co-expression network analysis (WGCNA).

The children were enrolled from the Department of Child Health, Women and Children’s Hospital, School of Medicine, Xiamen University. A total of 146 PP (including 30 CPP, 40 PT, and 76 unspecified PP) girls were enrolled at their first visit. The inclusion criteria of the patients are shown in detail in Fig. 1 according to the clinical guidelines and the related literatures [23,24,25]. In addition, 64 healthy girls were recruited as a control group for metabolic comparison, who were divided into 36 prepubertal and 28 adolescent girls based on their developmental status. Relevant clinical phenotypes were collected during the clinical examination. Morning fasting serum sample was collected from each girl through a clinical standard procedure and stood for 30 min, then centrifuged at 1000g for 10 min at 4 °C. The serum supernatant was transferred to a new centrifuge tube and stored at −80 °C until analysis.

Inclusion and screening process for study cohorts. *: Depending on the method of measurement. The detection method of basal LH and GnRH stimulation test in this study was immunochemiluminometric. LH: luteinizing hormone; FSH: follicular-stimulating hormone; BA: bone age (at the time of diagnosis); CA: chronologic age (at the time of diagnosis)

All serum samples were thawed at 4 °C, and 400 μL of serum was mixed with 200 μL of 60 mM phosphate buffer (pH 7.4, in 0.9% deuterated saline solution) and then vortexed for 10 s. After being centrifuged at 13,000g for 10 min at 4 °C, 550 μL of supernatant was transferred into a 5-mm NMR tube for 1H-NMR spectral acquisition.

The 1H-NMR spectra of serum samples were obtained on a 600-MHz Bruker Advance nuclear magnetic resonance (NMR) spectrometer (Bruker BioSpin, Germany) equipped with a triple resonance cryogenic probe operating at 600.13 MHz and 298.0 K. A typical water-suppressed Carr-Purcell-Meiboom-Gill (CPMG, [RD-90°-(τ-180°-τ)n-ACQ]) pulse sequence with a spectral width of 12,019.2 Hz, an acquisition time of 1.36 s, a relaxation delay of 4.0 s, a scan accumulation of 64 times, and a data point of 16 K was adopted to acquire 1H-NMR spectra.

Spectral processing was performed on MestReNova (version 14.1.1, Mestrelab Research S.L., Spain). All the free induction decays were zero-filled to 64 K data points and multiplied by an exponential function of 1.0 Hz line-broadening factor. The 1H-NMR spectra were manually phased, and baseline corrected, and then referenced to the doublet of endogenous lactate at δ1.33 after Fourier transformation. The spectral regions of δ4.70–δ5.17 and δ5.50–δ6.00 were removed to eliminate the interference of residual aquatic and urea signals. The remainder spectral regions (δ0.55–δ8.60) were integrally segmented into discrete regions of 0.002 ppm. To reduce the concentration difference between the samples, the obtained NMR spectral data were normalized to the total integrated area.

The serum PFCs of 40 PT and 30 CPP girls were measured by LC-MS/MS technique, and eleven kinds of PFCs, including perfluoro-n-octanoic acid (PFOA), potassium perfluoro-1-octanesulfonate (PFOS), perfluoro-n-butanoic acid (PFBA), perfluoro-n-undecanoic acid (PFUnDA), perfluoro-n-dodecanoic acid (PFDoDA), potassium perfluoro-1-butanesulfonate (PFBS), perfluoro-n-decanoic acid (PFDA), perfluoro-n-heptanoic acid (PFHpA), perfluoro-n-hexanoic acid (PFHA), potassium perfluoro-1-hexanesulfonate (TFHSA), and perfluoro-n-nonanoic acid (PFNA), were detected according to the references [26, 27]. The detailed procedures of serum sample preparation and PFC detection are shown in section S1 and Figure S1 in the Additional file 1.

The data are expressed as the means ± standard deviation (S.D.). All univariate statistical analysis was carried out using SPSS 20.0 software (SPSS Inc., Chicago, IL). According to distribution and variance homogeneity of indicator variables of data, Student’s t test or Mann-Whitney test were used for comparisons between two groups. If the variance was homogeneous, the p value was calculated directly; otherwise, Welch-Satterthwaite method was used. Differences were considered statistically significant when p < 0.05. Linear regression analysis was performed to correct the age factors.

The processed data of NMR were used for multivariate statistical analysis on SIMCA 14.1 software (Umetrics, Umea, Sweden) including principal component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA). The normalized dataset was scaled by unit variance, which makes each variable have the same variance. PCA is usually used for variable reduction and display of the relationship between samples such as whether there is clustering or outlier. The OPLS-DA models were applied to maximize the metabolic differences and extract the differential metabolites between the pair-wise groups. The NMR signals were assigned to individual metabolites with the referenced proton NMR peaks from Chenomx NMR Suite 8.1 (Chenomx Inc., EDBonton, AB, Canada) and confirmed by the public Human Metabolome Database (HMDB) (http://www.hmdb.ca/). In this study, the potential biomarkers were screened according to the following criteria: the value of variable importance for projection (VIP) > 1 of the metabolite and the p value after age correction (p-adj) < 0.05.

The comprehensive metabolic network was constructed by integrating all potential biomarkers identified from the present research using the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.ad.jp/kegg/), HMDB (http://www.hmdb.ca/) and MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/). Sensitivity and specificity of the potential biomarkers were analyzed by an exploratory receiver operating characteristic (ROC) curve based on random forest, and the area under ROC curve (AUC) and confidence interval (Cl) were determined correspondingly.

For the serum PFCs, after the raw data were transformed by natural logarithm, the Mann-Whitney test was utilized to analyze the statistical difference between PT and CPP girls. Spearman correlation analysis was performed to analyze the association of PFCs and clinical phenotype in PT and CPP girls, respectively.

In order to better understand the relationships of clinical phenotype-PFCs-metabolic characteristics, WGCNA were further used to identify the metabolite modules driven by clinical phenotype and PFCs. The protocol includes four aspects, namely network construction, module identification, the correlation between modules and features, and network visualization [28, 29]. Before analysis, the serum metabolic spectra (60 metabolites) were used as metabolite expression matrix, and clinical phenotype and PFCs of interest as trail matrix in the PT and CPP girls. Then the default WGCNA “step-by-step network construction” analysis was used to build modules. Firstly, the adjacencies between metabolites were calculated and constructed a topological overlap matrix (TOM). A hierarchical clustering tree was produced with the dissimilarity of the TOM and the modules were then selected by using the dynamic tree cut. Finally, the similar modules were merged by calculating the module eigenmetabolites (ME), clustering them and assigning a distance threshold (cut of 0.2). The detailed analysis process and network construction results could be found in section S2 and Figure S2 in the Additional file1.

According to the demographic and clinical characteristics (Table 1), the body mass index standard deviation score (BMISDS) of PP girls (0.48 ± 1.16) was significantly higher than that of both prepubertal (−0.32 ± 1.16) (p < 0.001) and adolescent girls (−0.42 ± 1.41) (p = 0.005), the BMISDS of CPP girls (0.70 ± 1.06) was significantly higher than that of both prepubertal (p < 0.001) and adolescent girls (p = 0.002), and the BMISDS of PT (0.33 ± 1.4) was also significantly higher than that of both prepubertal (p = 0.047) and adolescent (p = 0.0047). Obviously, the basal LH and FSH levels of PT girls (0.304 and 2.30 mIU/mL, respectively) were closer to prepubertal girls (0.256 and 2.46 mIU/mL, respectively), while those of CPP girls (2.73 and 5.01 mIU/mL, respectively) were closer to those of adolescent girls (2.93 and 5.85 mIU/mL, respectively), and estradiol and prolactin showed a similar trend. The serum levels of testosterone, dehydroepiandrosterone sulfate (DHEAS), free thyroxine (FT4), glucose, and urea/creatinine ratio were statistically significant between PT and CPP. The serum 25-hydroxyvitamin D (VD) level in CPP (61.86 nmol/L) was significantly lower than that of prepubertal (83.14 nmol/L) (p = 0.01) and PT (75.57 nmol/L) (p = 0.02).

A total of 60 metabolites were identified from the serum 1H-NMR spectra of PP, prepubertal, and adolescent girls (Additional file 1: Figure S3 and Table S3). No obvious sample separation was observed between PP and prepubertal and adolescent girls in the PCA score plots of the NMR data (Additional file 1: Figures S4A and S4B). The OPLS-DA highlighted and maximized the metabolic differences between PP and prepubertal and adolescent as shown in Fig. 2A1 and A2. The favorable model parameters, including R2 for indicating the explained variances of the original data and Q2 for indicating the predictive ability of the model, revealed the obvious metabolic differences between PP with prepubertal and adolescent (Additional file 1: Table S4), and the results were further externally cross-validated by the permutation tests (Additional file 1: Figure S5 and Table S4). According to the screening criteria (VIP >1 and p-adj < 0.05), a total of 16 metabolites were selected as the potential biomarkers of PP girls when compared with prepubertal girls and adolescent girls as demonstrated in Additional file 1: Table S5.

OPLS-DA score plots of serum samples. A1 The PP vs. prepubertal girls. A2 The PP vs. adolescent girls. B1 The CPP vs. prepubertal girls. B2 The CPP vs. adolescent girls. C1 The PT vs. prepubertal girls. C2 The PT vs. adolescent girls. The sample numbers of prepubertal, adolescent, PP, CPP, and PT girls were 36, 28, 146, 30, and 40, respectively. PPB: prepubertal; AD: adolescent

To clarify the respective pathogenesis of CPP and PT and differentiate the two classifications from the serum metabolome, the serum metabolic profile of CPP or PT girls was compared with prepubertal and adolescent, respectively. The OPLS-DA model revealed the obvious metabolic differences between CPP and prepubertal or adolescent girls (Fig. 2B1, B2 and Additional file 1: Table S4) and between PT and prepubertal or adolescent girls (Fig. 2C1, C2 and Additional file 1: Table S4) though no obvious separation was observed in the corresponding PCA score plots (Additional file 1: Figures S4C&S4D). The results were confirmed by the permutation tests (Additional file 1: Figure S5 and Table S4). According to the screening criteria, 23 metabolites were screened out as the potential biomarkers of CPP when compared with prepubertal girls and adolescent girls, and 25 metabolites were screened out as the potential biomarkers of PT when compared with prepubertal girls and adolescent girls as demonstrated in Table 2. Among them, formate, ethanol, and 3-hydroxybutyrate were both significantly upregulated in CPP and PT (Table 2).

The potential biomarkers of PP, CPP, and PT girls partially overlap (Fig. 3A), which is understandable because of their similar metabolic characteristics. A comparative analysis was performed to get their individual specific biomarkers. The results indicated that PP shared the potential biomarkers with CPP or PT, and eleven potential biomarkers, including glutamine, α-&β-glucose, dihydrothymine, methionine, hypoxanthine, isobutyrate, creatinine, valine, leucine, and phenylalanine, are specified to CPP, and eight potential biomarkers, including 1-methylhistidine, myo-inositol, N,N-dimethylglycine, very low-density lipoprotein (VLDL), glycerol, ornithine, asparagine, and glycine, are specified to PT (Fig. 3B). An exploratory ROC curve analysis based on random forest was used to evaluate the sensitivity and specificity of the specific biomarkers of the diseases (Fig. 4A). The AUC between CPP and prepubertal girls varied depending on the combinations of the eleven specific biomarkers, ranging from 0.66 to 0.73, and reached 0.721 (95% CI 0.553–0.874) when the eleven specific biomarkers were integrated (Fig. 4A1), while the AUC between CPP and prepubertal integrated adolescent girls was 0.646 (95% CI 0.526-0.757) at the combination of the eleven specific biomarkers (Fig. 4A2). The specific biomarkers for PT exhibit higher sensitivity. The AUC between PT and prepubertal girls was 0.97–0.98 and reached 0.972 (95% CI 0.904–1) at the combination of the eight specific biomarkers (Fig. 4A3). And it reached 0.822 (95% CI 0.698–0.92) between PT and prepubertal integrated adolescent girls when the eight specific biomarkers applied (Fig. 4A4). These results indicated the favorable sensitivity and specificity of the specific biomarkers of CPP and PT.

The potential biomarkers of PP, PT, and CPP. A The Venn plots of the potential biomarkers in the PP, PT, and CPP. B The bipartite graph of the potential biomarkers in the PP, PT, and CPP

An exploratory ROC curve analysis (A) and the validation PLS-DA models (B) based on the specific biomarkers. ROC curve analysis for the predictive power of specific biomarkers of CPP for distinguishing CPP from prepubertal girls (A1) and for distinguishing CPP from prepubertal integrated adolescent girls (A2). ROC curve analysis for the predictive power of specific biomarkers of PT for distinguishing PT from prepubertal girls (A3) and for distinguishing PT from prepubertal integrated adolescent girls (A4). The validation models were constructed to classify between CPP and prepubertal girls (B1), between CPP and prepubertal integrated adolescent girls (B2), between PT and prepubertal girls (B3), between PT and prepubertal integrated adolescent girls (B4). The sample numbers of prepubertal, adolescent, CPP and PT girls were 36, 28, 30, and 40, respectively. PPB: prepubertal; AD: adolescent

To further verify the superiority of disease-specific biomarkers in classification to the healthy girls, the PLS-DA models were reconstructed with the specific biomarkers as variables (Fig. 4B and Additional file 1: Table S4). The score plot showed a clear distinction between CPP and prepubertal girls (Fig. 4B1), and CPP could be clearly distinguished from prepubertal and adolescent girls with excellent predictive and explanatory power (Fig. 4B2 and Additional file 1: Table S4). A similar result was obtained in PT (Fig. 4B3, B4 and Additional file 1: Table S4).

To gain insight into the metabolic disorders of CPP and PT, the metabolic pathways were enriched based on the potential biomarkers through online database MetaboAnalyst 5.0, and the impact value was used to evaluate the crucial pathways involved in the occurrence of diseases. Based on the criterion of pathway impact > 0.1 and p < 0.05, six disturbed metabolic pathways were screened out from CPP, including aminoacyl-tRNA biosynthesis, valine, leucine, and isoleucine biosynthesis, phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, butanoate metabolism, and histidine metabolism (Additional file 1: Figure S6A), while seven metabolic pathways were significantly disordered in PT, including butanoate metabolism, synthesis and degradation of ketone bodies, glyoxylate and dicarboxylate metabolism, glycine, serine, and threonine metabolism, glycerolipid metabolism, histidine metabolism, and aminoacyl-tRNA biosynthesis (Additional file 1: Figure S6B). Furthermore, to better understand the process of diseases, the core metabolic network of CPP and PT were constructed to explain their individual pathogenesis based on the potential biomarkers via integrating the database of KEGG and HMDB. The links between hypothalamic-pituitary-gonadal-adrenal (HPGA) axis initiation and metabolism (including phenylalanine, tyrosine, and tryptophan biosynthesis, glycine, serine and threonine metabolism, glycolysis/gluconeogenesis, alanine, aspartate and glutamate metabolism, pyrimidine metabolism, aminoacyl-tRNA biosynthesis, and pyruvate metabolism and butanoate metabolism) were mainly shown in the core metabolic network of CPP (Fig. 5A). Metabolic pathway interconnections including glycerolipid metabolism, galactose metabolism, amino acid metabolism, pyruvate metabolism, butanoate metabolism, and pyrimidine metabolism were mainly shown in the core metabolic network of PT (Fig. 5B).

The core metabolic network of CPP (A) and PT (B). (HPG: hypothalamic-pituitary-gonadal; HPA: hypothalamic-pituitary-adrenal; GnRH: gonadotropin-releasing hormone; CRH: corticotropin releasing hormone; ACTH: adrenocorticotropic hormone; LH: luteinizing hormone; FSH: follicular-stimulating hormone; DHEAS: dehydroepiandrosterone sulfate; VLDL: very low-density lipoprotein.)

To understand the latent effects of PFCs on the occurrence of PP, serum PFCs in the PT and CPP girls were analyzed. The detailed detection results of PFCs are shown in Additional file 1: Table S6, which show the effectiveness of the detection method and the stability of the instrument. The levels of PFHPA and PFHA were statistically different between CPP and PT groups (p < 0.05), and the other nine PFCs were no statistical difference (Additional file 1: Table S7).

Spearman correlation analysis showed strong positive correlations among PFCs in the CPP girls. Estradiol was positively correlated with PFBA, TSH, and FSH. Prolactin was negatively correlated with PFBS and VD (Fig. 6a1 and Additional file 2: Table S8). In PT, the strong positive associations were also observed between PFCs, and prolactin was positively correlated with PFDUnDA, PFDA, and PFNA. VD was positively correlated with PFDA and PFNA, negatively correlated with PFBA (Fig. 6a2 and Additional file 2: Table S9).

Visualization plot of the association of clinical phenotype-PFCs-metabolic characteristics. A Spearman correlation analysis of PFCs and clinical phenotypes in the CPP (a1) and PT (a2) girls, red and blue represent positive and negative correlation, respectively (noted: only correlation values that meet p < 0.05 are presented on the heat map). B Module-trait associations of CPP (b1) and PT (b2), where each row corresponds to a module eigenmetabolite, column to a trait, respectively. Each cell contains the corresponding correlation and p value. C Visualization of module metabolites of CPP (c1) and PT (c2)

WGCNA is a systems biological method for understanding the correlated patterns between variables across different samples and has been widely used to find clusters or modules of metabolites. In this study, the metabolites were finally divided into eight and nine modules by WGCNA in CPP and PT, respectively, and the different modules were represented by different colors (the metabolites that do not belong to any module are classified as gray modules) (Additional file 1: Figure S2B). Furthermore, module-to-module correlation and cluster analysis are shown in Additional file 1: Figure S2C.

In module-trait heatmap, the driver module was determined based on the criterion |cor| > 0.30 and p < 0.05. In CPP, PFOA, PFNA, and PFDA mainly drive MEyellow, PFDoDA mainly drives MEbrown, DHEAS and VD mainly drive MEblack, FT4 mainly drives MEred, LH and FSH mainly drive MEbrown, and prolactin mainly drives MEblue (Fig. 6b1). In PT, PFOA mainly drives MEblack, TFHSA mainly drives MEpink, PFDoDA mainly drives MEyellow, and PFHpA mainly drives MEred, DHEAS mainly drives MEbrown, FSH mainly drives MEblack, BMISDS and body mass index (BMI) mainly drives MEgreen, FT4 mainly drives MEpink, TSH mainly drives MEturquoise, and estradiol mainly drives MEyellow (Fig. 6b2). The metabolites in each module are shown in Fig. 6C.

In this study, the levels of basal LH and FSH, estradiol, prolactin, testosterone, and DHEAS in PT and CPP indicated that the HPG axis is not activated in the PT girls but activated in the CPP girls, which was also confirmed by the GnRH stimulation test [30]. With activation of HPG axis, hypothalamus increasingly secretes gonadotropin GnRH, anterior pituitary gland increasingly secretes FSH and LH, and gonad increasingly secretes estradiol and testosterone. The higher level of estradiol (aromatized from testosterone) in girls before puberty is associated with an earlier thelarche, and it can affect and maintain cognitive function, regulate the sexual behavior and ovulation in the brain, and regulate higher order neural function [31, 32]. Behr et al. found that PFOA and PFOS can enhance estradiol-stimulated estrogen receptor β activity, and PFOS and PFBA can enhance dihydrotestosterone-stimulated androgen receptor activity [33]. In this study, PFBA was positively correlated with estradiol, PFBS was negatively correlated with prolactin in the CPP girls, and PFUnDA, PFDA, and PFNA were positively correlated with prolactin in the PT girls, which indicated that PFCs mainly caused developmental and reproductive toxicity by disrupting the body’s endocrine homeostasis, thus leading to an imbalance in the body’s steroid hormone secretion.

DHEAS is a stable marker for adrenal androgenic activity. At the biological level, DHEAS has effects on brain development, sexuality, mood and cognition, cardiovascular disease, stroke, and mortality [34]. Relevant longitudinal studies show that higher DHEAS level at the age of 8 predicts early menarche and doubles the risk of pubic hair development in girls [35]. In this study, the higher level of DHEAS in the CPP girls may be the result of a combination of adrenaline secretion and gonadal secretion, and it further leads to significant secondary sexual characteristics of the CPP girls. On the other hand, the increased level of DHEAS may potentially provide additional energy for the metabolic costs of early brain development, and also acts as a co-factor in promoting cortical maturation, thus leading to increased capacity for mentalizing and perspective-taking before the onset of reproductive maturation. This conclusion could be supported by the WGCNA analysis, where MEyellow in CPP and MEbrown in PT driven by DHEAS showed glucose aggregation. In addition, the significantly higher level of glucose in the CPP girls than in the PT girls also indicated a more active energy metabolism in CPP.

Interestingly, in this study, the BMI of CPP was higher than that of other groups. Studies have shown that the HPG axis initiates when girls’ body reaches a certain fat and/or protein mass [36]; therefore, the occurrence of CPP may be closely related to high BMI of girls. It is possible that obesity or overweight can promote the occurrence and development of CPP, which could be confirmed by the findings of many other scholars [37, 38]. In addition, the relationship between PFCs and VD in CPP and PT indicated PFCs partly disturbed osteogenesis [39], and also revealed new potential long-term PFC impacts on children.

CPP is primarily an early initiation of the HPGA axis, and hypothalamus generates a GnRH pulse that stimulates the pituitary gonadotropin secretion. Therefore, the corresponding core metabolic network mainly reflects the cause of activation of the HPGA axis and the metabolic disorder in the CPP girls. The link of clinical phenotype-PFCs-metabolic characteristics in the CPP girls indicated that PFCs may cause disturbance of metabolic network of CPP girls by disturbing endocrine homeostasis and/or directly affecting metabolite characteristics.

Tyrosine and phenylalanine are the precursor for catecholamines including tyramine, dopamine, epinephrine, and norepinephrine. The sympathetic nervous system neurotransmitters and catecholamines, especially norepinephrine, play an important role in the regulation of GnRH neurons [40]. The downregulated levels of phenylalanine and tyrosine in serum of the CPP girls indicated that the biosynthesis of phenylalanine, tyrosine, and tryptophan was disturbed, which keeps consistent with the previous metabolomics results [5]. The downregulation of phenylalanine and tyrosine levels may be due to the more consumption of norepinephrine for the activation of HPG axis, resulting in lower levels of its precursor substances. Such conclusion could be confirmed by the WGCNA analysis, where the modules MEblue and MEturquoise were primarily driven by prolactin and DHEAS (Fig. 6b1, c1).

Serine and glycine are connected through biosynthesis to provide necessary synthetic precursors of proteins, nucleic acids, and lipids. At the same time, serine homeostasis plays a vital role in maintaining brain energy metabolism [41]. The homeostasis of serine/glycine is essential for the proliferation of human primary muscle progenitor cells and efficient skeletal muscle regeneration [42]. Therefore, the decreased serine level in the CPP girls may be one of the reasons for the lifelong high short stature in adults, but further research is needed to confirm it. Creatine is essential in maintaining human growth, development, and health and can improve skin and bone health [43]. Creatinine was selected as a biomarker for predict CPP in Qi’s study [5], and it was regarded as one of the specific biomarkers for CPP in this study. Creatinine is linked to muscle mass, and the lower level of serum creatinine observed in the CPP girls could be a result of the excessive weight and relative low muscle/fat ratio caused by low physical activity. In the MEyellow and MEblue, creatine and creatinine are mainly related to amino acids and driven by multiple PFCs (PFOA, PFNA, PFDA), DHEAS, and prolactin (Fig. 6b1, c1), indicated that exposure to PFCs and/or endocrine disturbance will cause amino acid metabolism disturbance, thus affecting the bone growth of girls.

The increased levels of α-&β-glucose and lactose and decreased level of hypoxanthine in the CPP girls suggested the disturbed glycolysis/gluconeogenesis (energy metabolism). PFOA, PFDA, and PFNA were positively correlated with PFOS and mainly drove MEyellow containing α-&β-glucose. Studies have shown that the accumulation of PFCs, especially PFOS, contributes to the disorder of lipid and glucose metabolisms in children, but the mechanism is still in the initial stage of exploration [12, 44]. Compared with adolescent girls, the CPP girls had lower insulin sensitivity, glucose and lipid metabolism profile, and body composition, and the metabolic disturbance remained unchanged even after 1 year of GnRH treatment [45, 46]. The accumulation of PFCs may be a crucial reason for the decreased insulin sensitivity and high level of serum glucose in the CPP girls [47].

Glutamate and GABA are the principal excitatory and inhibitory neurotransmitters. Their interactions with GnRH neurons, including the regulation of GnRH gene and protein expression, hormone release, and modulation by estrogen, are critical to age-appropriate changes in reproductive function [48]. Glutamate is also crucial for bone growth and reconstruction [49]. In our study, the significantly decreased glutamate level in the CPP girls implies its dynamic change in the development of CPP. We speculated that the increased glutamate level stimulated the HPG axis and induced CPP in the early stage. However, the occurrence of CPP is often accompanied with rapid height growth, and subsequently the glutamate in bone and blood is excessively consumed at a certain stage of CPP.

In addition, the downregulated levels of isoleucine, alanine, valine, methionine, histidine, and α-ketoisovalerate in the CPP girls indicated the disturbed aminoacyl-tRNA biosynthesis. The disruption of aminoacyl-tRNA biosynthesis further supports our inference that body needs to consume large amounts of amino acids during the initial growth spurt and rapid bone maturation stage of the CPP girls. The significantly reduced level of VD in the CPP girls further indicated the affected bone growth of the CPP girls. It is possible that the lack of VD affected the absorption and deposition of calcium in the bone, thus affecting the health and growth of bones. Balance of aminoacyl-tRNA biosynthesis is closely related to bone health, and its disruption can lead to osteocyte protein synthesis dysfunction, marrow hypoplasia, and osteoporosis [50]. This may also be the reason for the short height of CPP girls at adult. However, whether appropriate supplementation of amino acids in diet can improve the adult height, or whether aminoacyl-tRNA biosynthesis can become a targeted therapeutic metabolic pathway for CPP girls is worthy of further research and discussion.

Purines and pyrimidines give prominent contributions in the development of the central nervous system, but the exact molecular mechanisms remain unclear [51]. PP girls have an increased risk of precocious sexual behavior and an increased prevalence of mental disorders such as depression and anxiety in adulthood [52, 53]. In this study, the downregulated serum levels of hypoxanthine, dihydrothymine, and uridine were observed in the CPP girls. Dihydrothymine belongs to MEblack jointly driven by PFOA, PFNA, DHEAS, FT4, and VD. This suggested that the disturbed pyrimidine metabolic pathway may be caused by the joint action of external environment and internal factors. CPP girls need to advance adaption to the cognitive, emotional, and changes of puberty prematurely, therefore, the sense of anxiety and depression increases significantly, resulting in a certain disorder of pyrimidine metabolism.

PT does not involve the early initiation of HPGA axis and keep normal growth of height and maturation of bone age. Its occurrence is mainly due to the breast enlargement caused by exogenous hormone intake. In Qi et al. study [5], upregulated levels of succinate and 1-methylhistidine in urine can effectively predict PT. In this study, it was found that most amino acids showed upregulated in PT girls, including the specific biomarker 1-methylhistidine and the potential biomarker succinate. Amino acids are not only the building blocks of proteins and an indispensable component of cells, but also play versatile roles in regulating cell metabolism, proliferation, differentiation, and growth by themselves or their derivatives. Their requirements vary at the various stages of children’s growth and development [54, 55]. Studies have shown that arginine and ornithine supplementation can promote the secretion of growth hormone and insulin growth factor-1 [56]. We speculated that the nutritional diets promote the high level of amino acid in the PT girls for providing energy for the growth and development. But the more relevant reason may be that exposure to PFCs disturbs the body’s amino acid metabolism, endocrine homeostasis, and vitamin level, which will affect the girl’s bone health and growth.

In this study, the metabolites in glycerolipid metabolic pathway and galactose metabolic pathway including N,N-dimethylglycine, ethanolamine, glycerol, dihydrothymine, and myo-inositol were upregulated in the PT girls. Studies have shown that thyroid dysfunction will affect the body’s glycerol metabolism and gluconeogenesis pathway [57, 58]. In this study, FT4 and TFHSA mainly drive MEpink containing glycerol, and PFHpA mainly drives MEred containing low-density lipoprotein, myo-inositol, phosphocholine, and glycerophosphorycholine, suggesting that the PFCs exposure and high level FT4 may be the main reason contributed to the disorder of serum glycerolipid metabolic and galactose metabolism in PT. However, this study cannot directly prove the causation between the elevated FT4 level and PFCs exposure in PT.

Interestingly, the disordered trend of pyruvate metabolism and butyrate metabolism kept almost consistent between CPP and PT (Fig. 5), indicating that pyruvate metabolism and butyrate metabolism disorder are possibly a common phenomenon of PP. In other words, simultaneously elevated levels in serum formate, ethanol, and 3-hydroxybutyrate may serve as the early diagnostic indicators for CPP and PT, i.e., precocious puberty in girls, but the stratification of PP still needs to be further determined based on the serum-specific biomarkers. Ethanol may act on bone remodeling including osteocyte apoptosis [59, 60]. Elevated level of 3-hydroxybutyrate will inhibit the differentiation and growth of (pre)-chondrocytes [61] and participate in the process of osteoporosis [62]. There is a strong correlation between PFOA and PFNA in CPP and PT. Ethanol and 3-hydroxybutyrate belong to MEred which were jointly driven by PFNA and FT4 in CPP. In PT, ethanol and 3-hydroxybutyrate belong to the MEblack which were jointly driven by PFOA and FT4. It is suggested that PFC exposure and/or fluctuations in FT4 level will lead to elevated ethanol and 3-hydroxybutyrate in PP and thus affects bone growth and development, but more evidence is needed to prove this point.

It should be pointed out that this study had several limitations. Firstly, the results were obtained from a smaller-sample-size cohort, and measurements in an expanded cohort are necessary to validate the disease-specific biomarkers and determine their universality. Secondly, clinical reports have revealed the gender-related differences of PP not only in symptoms but also in pathophysiology. Therefore, the present findings could not naturally extend to boys, and further study on the sex-specific metabolic characteristics of PP is needed to comprehensively understand the underlying pathogenic mechanisms of PP. Finally, the detailed mechanism of how PFCs lead to endocrine homeostasis imbalance still keeps unclear and needs further study.

In summary, based on the cross-metabolomics analyses, our study showed that formate, ethanol, and 3-hydroxybutyrate may serve as the early diagnostic indicators for the CPP and PT, i.e., precocious puberty in girls, and eleven CPP-specific biomarkers and eight PT-specific biomarkers in serum exhibited good sensitivity, which can facilitate the classification diagnosis of CPP and PT girls. The link of clinical phenotype-PFCs-metabolic characteristics in CPP and PT by WGCNA method revealed that PFC exposure is associated with endocrine homeostasis imbalance in the CPP and PT girls, and thus directly or indirectly drives metabolic changes and forms perturbations of the overall metabolic network. These findings will provide a potential diagnostic and stratification approach for the clinical diagnosis of precocious puberty in girls, and also raise social awareness that reducing exposure to PFC compounds may be an important strategy for preventing the occurrence and development of precocious puberty in girls.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

One-dimensional nuclear magnetic resonance hydrogen spectrum

Area under ROC curve

Body mass index

Body mass index standard deviation score

Confidence interval

Dehydroepiandrosterone sulfate

Follicular-stimulating hormone

Free thyroxine

Gonadotropin-releasing hormone

Human Metabolome Database

Hypothalamic-pituitary-gonadal

Hypothalamic-pituitary-gonadal-adrenal

Kyoto Encyclopedia of Genes and Genomes

Luteinizing hormone

Module eigenmetabolites

Nuclear magnetic resonance

Orthogonal partial least squares-discriminant analysis

p value after age correction

Principal component analysis

Perfluoro-n-butanoic acid

Potassium perfluoro-1-butanesulfonate

Perfluorinated compounds

Perfluoro-n-decanoic acid

Perfluoro-n-dodecanoic acid

Perfluoro-n-hexanoic acid

Perfluoro-n-heptanoic acid

Perfluoro-n-nonanoic acid

Perfluoro-n-octanoic acid

Potassium perfluoro-1-octanesulfonate

Perfluoro-n-undecanoic acid

Receiver operating characteristic

Potassium perfluoro-1-hexanesulfonate

Topological overlap matrix

Thyroid-stimulating hormone

25-hydroxyvitamin D

Variable importance for projection

Very low-density lipoprotein

Weighted gene co-expression network analysis

Cantas-Orsdemir S, Eugster EA. Update on central precocious puberty: from etiologies to outcomes. Expert Rev Endocrinol Metab. 2019;14(2):123–30.

Article CAS PubMed Google Scholar

Roberts SA, Kaiser UB. GENETICS IN ENDOCRINOLOGY: Genetic etiologies of central precocious puberty and the role of imprinted genes. Eur J Endocrinol. 2020;183(4):R107–17.

Article CAS PubMed PubMed Central Google Scholar

Mogensen SS, Aksglaede L, Mouritsen A, Sorensen K, Main KM, Gideon P, Juul A. diagnostic work-up of 449 consecutive girls who were referred to be evaluated for precocious puberty. J Clin Endocrinol Metab. 2011;96(5):1393–401.

Article CAS PubMed Google Scholar

Lee DS, Ryoo NY, Lee SH, Kim S, Kim JH. Basal luteinizing hormone and follicular stimulating hormone: is it sufficient for the diagnosis of precocious puberty in girls? Ann Pediatr Endocrinol Metab. 2013;18(4):196–201.

Article PubMed PubMed Central Google Scholar

Qi Y, Li P, Zhang Y, Cui L, Guo Z, Xie G, Su M, Li X, Zheng X, Qiu Y, et al. Urinary metabolite markers of precocious puberty. Mol Cell Proteomics. 2012;11(1):M111 011072.

Article PubMed Google Scholar

Yang L, Tang K, Qi Y, Ye H, Chen W, Zhang Y, Cao Z. Potential metabolic mechanism of girls’ central precocious puberty: a network analysis on urine metabonomics data. BMC Syst Biol. 2012;6(Suppl 3):S19.

Article PubMed PubMed Central Google Scholar

Ojo AF, Peng C, Ng JC. Assessing the human health risks of per- and polyfluoroalkyl substances: a need for greater focus on their interactions as mixtures. J Hazard Mater. 2021;407:124863.

Article CAS PubMed Google Scholar

Zhou W, Zhang L, Tong C, Fang F, Zhao S, Tian Y, Tao Y, Zhang J. Shanghai Birth Cohort S: plasma perfluoroalkyl and polyfluoroalkyl substances concentration and menstrual cycle characteristics in preconception women. Environ Health Perspect. 2017;125(6):067012.

Article PubMed PubMed Central Google Scholar

Ullah S, Ahmad S, Guo X, Ullah S, Ullah S, Nabi G, Wanghe K. A review of the endocrine disrupting effects of micro and nano plastic and their associated chemicals in mammals. Front Endocrinol (Lausanne). 2022;13:1084236.

Article PubMed Google Scholar

Rappazzo KM, Coffman E, Hines EP: Exposure to Perfluorinated Alkyl Substances and Health Outcomes in Children: A Systematic Review of the Epidemiologic Literature. Int J Environ Res Public Health. 2017;14(7):691.

Karwacka A, Zamkowska D, Radwan M, Jurewicz J. Exposure to modern, widespread environmental endocrine disrupting chemicals and their effect on the reproductive potential of women: an overview of current epidemiological evidence. Hum Fertil (Camb). 2019;22(1):2–25.

Article CAS PubMed Google Scholar

Salter DM, Wei W, Nahar PP, Marques E, Slitt AL. Perfluorooctanesulfonic Acid (PFOS) Thwarts the Beneficial Effects of Calorie Restriction and Metformin. Toxicol Sci. 2021;182(1):82–95.

Article CAS PubMed PubMed Central Google Scholar

Lee YJ, Jung HW, Kim HY, Choi YJ, Lee YA. Early-life exposure to per- and poly-fluorinated Alkyl substances and growth, adiposity, and puberty in children: a systematic review. Front Endocrinol (Lausanne). 2021;12:683297.

Article PubMed Google Scholar

Gore AC, Moore T, Groom MJ, Thompson LM. Prenatal exposure to an EDC mixture, NeuroMix: effects on brain, behavior, and stress responsiveness in rats. Toxics. 2022;10(3):122.

Article CAS PubMed PubMed Central Google Scholar

Qiu Z, Qu K, Luan F, Liu Y, Zhu Y, Yuan Y, Li H, Zhang H, Hai Y, Zhao C. Binding specificities of estrogen receptor with perfluorinated compounds: a cross species comparison. Environ Int. 2020;134:105284.

Article CAS PubMed Google Scholar

Godfrey A, Hooser B, Abdelmoneim A, Sepúlveda MS. Sex-specific endocrine-disrupting effects of three halogenated chemicals in Japanese medaka. J Appl Toxicol. 2019;39(8):1215–23.

Article CAS PubMed Google Scholar

Di Nisio A, Sabovic I, Valente U, Tescari S, Rocca MS, Guidolin D, Dall’Acqua S, Acquasaliente L, Pozzi N, Plebani M, et al. Endocrine disruption of androgenic activity by perfluoroalkyl substances: clinical and experimental evidence. J Clin Endocrinol Metab. 2019;104(4):1259–71.

Article PubMed Google Scholar

Ding N, Harlow SD, Randolph JF Jr, Loch-Caruso R, Park SK. Perfluoroalkyl and polyfluoroalkyl substances (PFAS) and their effects on the ovary. Hum Reprod Update. 2020;26(5):724–52.

Article CAS PubMed PubMed Central Google Scholar

Roepke TA, Sadlier NC. REPRODUCTIVE TOXICOLOGY: Impact of endocrine disruptors on neurons expressing GnRH or kisspeptin and pituitary gonadotropins. Reproduction. 2021;162(5):131–45.

Article Google Scholar

Kjeldsen LS, Bonefeld-Jørgensen EC. Perfluorinated compounds affect the function of sex hormone receptors. Environ Sci Pollut Res Int. 2013;20(11):8031–44.

Article CAS PubMed Google Scholar

Spaziani M, Tarantino C, Tahani N, Gianfrilli D, Sbardella E, Lenzi A, Radicioni AF. Hypothalamo-Pituitary axis and puberty. Mol Cell Endocrinol. 2021;520: 111094.

Article CAS PubMed Google Scholar

La Merrill MA, Vandenberg LN, Smith MT, Goodson W, Browne P, Patisaul HB, Guyton KZ, Kortenkamp A, Cogliano VJ, Woodruff TJ, et al. Consensus on the key characteristics of endocrine-disrupting chemicals as a basis for hazard identification. Nat Rev Endocrinol. 2020;16(1):45–57.

Article PubMed Google Scholar

Bangalore Krishna K, Fuqua JS, Rogol AD, Klein KO, Popovic J, Houk CP, Charmandari E, Lee PA, Freire AV, Ropelato MG, et al. Use of gonadotropin-releasing hormone analogs in children: update by an international consortium. Horm Res Paediatr. 2019;91(6):357–72.

Article CAS PubMed Google Scholar

Latronico AC, Brito VN, Carel JC. Causes, diagnosis, and treatment of central precocious puberty. Lancet Diabetes Endocrinol. 2016;4(3):265–74.

Article PubMed Google Scholar

Soriano-Guillen L, Argente J. Central precocious puberty, functional and tumor-related. Best Pract Res Clin Endocrinol Metab. 2019;33(3): 101262.

Article PubMed Google Scholar

Liu L, She J, Zhang X, Zhang J, Tian M, Huang Q, Shah Eqani SA, Shen H. Online background cleanup followed by high-performance liquid chromatography with tandem mass spectrometry for the analysis of perfluorinated compounds in human blood. J Sep Sci. 2015;38(2):247–53.

Article PubMed Google Scholar

Yusa V, Ye X, Calafat AM. Methods for the determination of biomarkers of exposure to emerging pollutants in human specimens. Trends Analyt Chem. 2012;38:129–42.

Article CAS PubMed PubMed Central Google Scholar

Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

Article PubMed PubMed Central Google Scholar

Liu W, Shi X, Dai T, Shen G, Feng J. Recoupled-STOCSY-based co-expression network analysis to extract phenotype-driven metabolite modules in NMR-based metabolomics dataset. Anal Chim Acta. 2022;1197:339528.

Article CAS PubMed Google Scholar

Vuralli D, Gonc EN, Ozon ZA, Alikasifoglu A. Adequacy of basal luteinizing hormone levels in the diagnosis of central precocious puberty. Turk Pediatri Ars. 2020;55(2):131–8.

PubMed PubMed Central Google Scholar

Luine VN. Estradiol and cognitive function: past, present and future. Horm Behav. 2014;66(4):602–18.

Article CAS PubMed PubMed Central Google Scholar

Pereira A, Corvalan C, Uauy R, Klein KO, Mericq V. Ultrasensitive estrogen levels at 7 years of age predict earlier thelarche: evidence from girls of the growth and obesity Chilean cohort. Eur J Endocrinol. 2015;173(6):835–42.

Article CAS PubMed Google Scholar

Behr AC, Lichtenstein D, Braeuning A, Lampen A, Buhrke T. Perfluoroalkylated substances (PFAS) affect neither estrogen and androgen receptor activity nor steroidogenesis in human cells in vitro. Toxicol Lett. 2018;291:51–60.

Article CAS PubMed Google Scholar

Campbell B. DHEAS and human development: an evolutionary perspective. Front Endocrinol (Lausanne). 2020;11:101.

Article PubMed Google Scholar

Thankamony A, Ong KK, Ahmed ML, Ness AR, Holly JM, Dunger DB. Higher levels of IGF-I and adrenal androgens at age 8 years are associated with earlier age at menarche in girls. J Clin Endocrinol Metab. 2012;97(5):E786-790.

Article CAS PubMed PubMed Central Google Scholar

Wahab F, Shahab M, Behr R. Hypothesis: Irisin is a metabolic trigger for the activation of the neurohormonal axis governing puberty onset. Med Hypotheses. 2016;95:1–4.

Article CAS PubMed Google Scholar

De Leonibus C, Marcovecchio ML, Chiavaroli V, de Giorgis T, Chiarelli F, Mohn A. Timing of puberty and physical growth in obese children: a longitudinal study in boys and girls. Pediatr Obes. 2014;9(4):292–9.

Article PubMed Google Scholar

Ullah R, Su Y, Shen Y, Li C, Xu X, Zhang J, Huang K, Rauf N, He Y, Cheng J, et al. Postnatal feeding with high-fat diet induces obesity and precocious puberty in C57BL/6J mouse pups: a novel model of obesity and puberty. Front Med. 2017;11(2):266–76.

Article PubMed Google Scholar

Liu S, Yang R, Yin N, Wang YL, Faiola F. Environmental and human relevant PFOS and PFOA doses alter human mesenchymal stem cell self-renewal, adipogenesis and osteogenesis. Ecotoxicol Environ Saf. 2019;169:564–72.

Article CAS PubMed Google Scholar

Bronzi CD, Orozco AS, Rodriguez D, Rastrilla AM, Sosa ZY, Casais M. Noradrenaline modulates the presence of gonadotropin-releasing hormone in ovary. The importance of its interrelation on the ovarian steroidogenesis and apoptosis on dioestrus II in rat. J Steroid Biochem Mol Biol. 2015;154:39–46.

Article CAS PubMed Google Scholar

Maugard M, Vigneron PA, Bolanos JP, Bonvento G. l-Serine links metabolism with neurotransmission. Prog Neurobiol. 2021;197:101896.

Article CAS PubMed Google Scholar

Gheller BJ, Blum JE, Lim EW, Handzlik MK, Hannah Fong EH, Ko AC, Khanna S, Gheller ME, Bender EL, Alexander MS, et al. Extracellular serine and glycine are required for mouse and human skeletal muscle stem and progenitor cell function. Mol Metab. 2021;43: 101106.

Article CAS PubMed Google Scholar

Wu G. Important roles of dietary taurine, creatine, carnosine, anserine and 4-hydroxyproline in human nutrition and health. Amino Acids. 2020;52(3):329–60.

Article CAS PubMed PubMed Central Google Scholar

Timmermann CA, Rossing LI, Grontved A, Ried-Larsen M, Dalgard C, Andersen LB, Grandjean P, Nielsen F, Svendsen KD, Scheike T, et al. Adiposity and glycemic control in children exposed to perfluorinated compounds. J Clin Endocrinol Metab. 2014;99(4):E608-614.

Article CAS PubMed Google Scholar

Hur JH, Park S, Jung MK, Kang SJ, Kwon A, Chae HW, Kim HS, Kim DH. Insulin resistance and bone age advancement in girls with central precocious puberty. Ann Pediatr Endocrinol Metab. 2017;22(3):176–82.

Article PubMed PubMed Central Google Scholar

Arcari AJ, Freire AV, Escobar ME, Ballerini MG, Ropelato MG, Bergada I, Gryngarten MG. One-year treatment with gonadotropin-releasing hormone analogues does not affect body mass index, insulin sensitivity or lipid profile in girls with central precocious puberty. J Pediatr Endocrinol Metab. 2019;32(2):181–6.

Article CAS PubMed Google Scholar

Margolis R, Sant KE. Associations between exposures to perfluoroalkyl substances and diabetes, hyperglycemia, or insulin resistance: a scoping review. J Xenobiot. 2021;11(3):115–29.

Article CAS PubMed PubMed Central Google Scholar

Maffucci JA, Gore AC. Chapter 2: hypothalamic neural systems controlling the female reproductive life cycle gonadotropin-releasing hormone, glutamate, and GABA. Int Rev Cell Mol Biol. 2009;274:69–127.

Article CAS PubMed PubMed Central Google Scholar

Skerry TM. The role of glutamate in the regulation of bone mass and architecture. J Musculoskelet Neuronal Interact. 2008;8(2):166–73.

CAS PubMed Google Scholar

Aleidi SM, Alnehmi EA, Alshaker M, Masood A, Benabdelkamel H, Al-Ansari MM, Abdel Rahman AM. A Distinctive Human Metabolomics Alteration Associated with Osteopenic and Osteoporotic Patients. Metabolites 2021;11(9):628.

Micheli V, Camici M, Tozzi MG, Ipata PL, Sestini S, Bertelli M, Pompucci G. Neurological disorders of purine and pyrimidine metabolism. Curr Top Med Chem. 2011;11(8):923–47.

Article CAS PubMed Google Scholar

Temelturk RD, Ilcioglu Ekici G, Beberoglu M, Siklar Z, Kilic BG. Managing precocious puberty: A necessity for psychiatric evaluation. Asian J Psychiatr. 2021;58:102617.

Article PubMed Google Scholar

Coban OG, Bedel A, Onder A, Adanir AS, Tuhan H, Parlak M. Psychiatric disorders, peer-victimization, and quality of life in girls with central precocious puberty. J Psychosom Res. 2021;143:110401.

Article PubMed Google Scholar

Gao H. Amino Acids in Reproductive Nutrition and Health. Adv Exp Med Biol. 2020;1265:111–31.

Article CAS PubMed PubMed Central Google Scholar

Pillai RR, Kurpad AV. Amino acid requirements in children and the elderly population. Br J Nutr. 2012;108(Suppl 2):S44-49.

Article CAS PubMed Google Scholar

Zajac A, Poprzecki S, Zebrowska A, Chalimoniuk M, Langfort J. Arginine and ornithine supplementation increases growth hormone and insulin-like growth factor-1 serum levels after heavy-resistance exercise in strength-trained athletes. J Strength Cond Res. 2010;24(4):1082–90.

Article PubMed Google Scholar

Yan Y, Niu Z, Sun C, Li P, Shen S, Liu S, Wu Y, Yun C, Jiao T, Jia S, et al. Hepatic thyroid hormone signalling modulates glucose homeostasis through the regulation of GLP-1 production via bile acid-mediated FXR antagonism. Nat Commun. 2022;13(1):6408.

Article CAS PubMed PubMed Central Google Scholar

Piras C, Pibiri M, Leoni VP, Balsamo A, Tronci L, Arisci N, Mariotti S, Atzori L. Analysis of metabolomics profile in hypothyroid patients before and after thyroid hormone replacement. J Endocrinol Invest. 2021;44(6):1309–19.

Article CAS PubMed Google Scholar

Maurel DB, Boisseau N, Benhamou CL, Jaffre C. Alcohol and bone: review of dose effects and mechanisms. Osteoporos Int. 2012;23(1):1–16.

Article CAS PubMed Google Scholar

Lin X, Zhu D, Wang K, Luo P, Rui G, Gao Y, Liu F, Yu H. Activation of aldehyde dehydrogenase 2 protects ethanol-induced osteonecrosis of the femoral head in rat model. Cell Prolif. 2022;55(6):e13252.

Article CAS PubMed PubMed Central Google Scholar

Wu Y, Teng Y, Zhang C, Pan Y, Zhang Q, Zhu X, Liu N, Su X, Lin J. The ketone body β-hydroxybutyrate alleviates CoCrMo alloy particles induced osteolysis by regulating NLRP3 inflammasome and osteoclast differentiation. J Nanobiotechnology. 2022;20(1):120.

Article CAS PubMed PubMed Central Google Scholar

Li X, Wang Y, Gao M, Bao B, Cao Y, Cheng F, Zhang L, Li Z, Shan J, Yao W. Metabolomics-driven of relationships among kidney, bone marrow and bone of rats with postmenopausal osteoporosis. Bone. 2022;156:116306.

Article PubMed Google Scholar

Download references

Not applicable.

This work was supported by the National Natural Science Foundation of China (Nos. 82072015 & 82103859), the Natural Science Foundation of Fujian Province of China (No. 2022J01062), the Guiding project of the Natural Science Foundation of Fujian Province of China (No. 2019D010), and the Young and Middle-aged Talent Cultivation Projects of Fujian Province of China (No. 2019-ZQNB-31).

Jinxia Wu and Jing Chen have contributed equally to this work.

Department of Electronic Science, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Siming District, 422 Siming South Road, Xiamen, 361005, Fujian, China

Jinxia Wu, Yajie Chang, Guiping Shen & Jianghua Feng

Department of Child Health, Women and Children’s Hospital, School of Medicine, Xiamen University, Xiamen, 361003, Fujian, China

Jing Chen, Rong Huang & Yanyan Lin

Department of Pediatrics, The First Affiliated Hospital of Bengbu Medical College, Anhui, Bengbu, 233000, China

Hongwei Zhu

Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Guangzhou, 510060, Guangdong, China

Lin Che

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

JXW, JC, and GPS designed the study, analyzed, and interpreted the data. JXW, LC, and YJC conducted the analytical experiments. JXW wrote the first draft of the manuscript. JHF and HWZ supervised the statistical analysis and experiments. JC, YJC, RH, and YYL organized samples collection. GPS, JHF, RH, and YYL contributed to the conceptualization of the study and the revision of the manuscript. All authors read and approved the final manuscript

Correspondence to Guiping Shen.

This study was performed under the guidance of the 1964 Declaration of Helsinki and its later amendments, and the procedure was approved by the Human Ethics Committees of Women and Children’s Hospital (KY-2016002), Xiamen University. All the families of the girls were informed of the aims of the study, and informed consent was written and obtained from the parents of all girls.

Not applicable.

The authors declare that they have no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Table S1. LC-MS/MS analytical mobility elution procedure. Table S2. Retention time, transitions and MS/MS conditions of the analytes. Table S3. The identified metabolites from 1H-NMR spectra of the girls’ serum samples. Table S4. Summary of the quality parameters of the multivariate statistical analysis. Table S5. The potential biomarkers in serum of PP. Table S6. Linear range, correlation coefficient, LOD and recovery of PFCs (n = 5). Table S7. Concentration and statistical analysis of serum PFCs in the CPP and PT groups. Figure S1. Sample preparation and detection procedure of serum PFCs. Figure S2. WGCNA related visualization diagram. Figure S3. Mean 1H-NMR spectra of serum from prepubertal, PP, PT, CPP and adolescent girls. Figure S4. PCA score plots of serum samples. Figure S5. Permutation test analysis to test the over-fitting of OPLS-DA model. Figure S6. The pathways enrichment analysis of the CPP (A) and PT (B) base on the potential biomarkers via MetaboAnalyst 5.0.

Table S8. Spearman correlation analysis between clinical phenotype and PFCs in CPP girls. Table S9. Spearman correlation analysis between clinical phenotype and PFCs in PT girls.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

Wu, J., Chen, J., Huang, R. et al. Metabolic characteristics and pathogenesis of precocious puberty in girls: the role of perfluorinated compounds. BMC Med 21, 323 (2023). https://doi.org/10.1186/s12916-023-03032-0

Download citation

Received: 19 February 2023

Accepted: 15 August 2023

Published: 25 August 2023

DOI: https://doi.org/10.1186/s12916-023-03032-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative