Skip to main content
  • Research article
  • Open access
  • Published:

Associations among IGF-1, IGF2, IGF-1R, IGF-2R, IGFBP-3, insulin genetic polymorphisms and central precocious puberty in girls

Abstract

Background

Insulin and insulin-like growth factor (IGF)-1 coupled with growth hormone helps control timing of sexual maturation. Mutations and variants in multiple genes are associated with development or reduced risk of central precocious puberty (CPP).

Methods

We assessed single nucleotide polymorphisms (SNPs) in the IGF-1, IGF-2, IGF-3, IGF-1 receptor (IGF1R), IGF-2 receptor (IGF2R), and IGF -binding protein 3 (IGFBP-3) genes, and their association with demographics and metabolic proteins in girls with CPP. Z-scores of height, weight, and body mass index (BMI) were calculated with the WHO reference growth standards for children.

Results

IGF-1 serum levels of CPP group exhibited a higher correlation with bone age, z-scores of height and weight, and luteinizing hormone (LH) than those of control group, regardless of BMI adjustment. In the CPP group, height was associated with IGF-2(3580), an adenine to guanine (A/G) SNP at position + 3580. BMI in the CPP group was associated with IGF-2(3580), IGF1R, and the combinations of [IGF-2(3580) + IGF2R], and [IGF-2(3580) + IGFBP-3]. Body weight in the CPP group was associated with the combination of [IGF-2(3580) + IGFBP-3] (p = 0.024). Weight and BMI were significantly associated with the combination of [IGF-2(3580) + IGF2R + IGFBP-3] in the CPP group. These associations were not significantly associated with z-scores of weight, height, or BMI. The distribution of these genotypes, haplotypes, and allele frequencies were similar between control and CPP groups.

Conclusions

These known SNPs of these IGF-1 axis genes appear to play minor roles in the risk for development of CPP.

Peer Review reports

Background

Precocious puberty (PP) is defined as early development of puberty in females and males (before 8 years and 9 years of age, respectively). The two types of PP, central precocious puberty (CPP) and pseudo or peripheral precocious puberty (PPP) differ in their etiology. Secondary sex characteristics of PPP develop early due to excessive hormonal secretion from reproductive tumors or adrenal hyperplasia, and PPP is considered gonadotropin-independent. In contrast, CPP is gonadotropin-dependent and involves the premature maturation of the hypothalamic-pituitary-gonadal axis which induces the early secretion of testosterone from boys’ testes and estrogens from girls’ ovaries. Although the cause of CPP is not elucidated in most cases, CPP can be associated with obesity, brain structural abnormalities, head injuries, and products that contain compounds which mimic hormones including cosmetic products, some insecticides, and some foods [1]. CPP causes early onset of menarche and initial breast development in girls, psychosocial challenges, and early epiphyseal fusion in bones which reduces further growth spurts and thus may lower final height. Treatment has the goal of preserving growth and height potential and can involve gonadotropin-releasing hormone (GnRH) analogs [2]. Treatment for CPP did not significantly affect the risk of cancer death, obesity, or metabolic disorders in 30- to 50-year-old women who had had CPP [3].

Environmental and metabolic factors may influence the initiation of puberty in 20 to 50% of cases [4]. Poor nutrition may delay puberty while obesity promotes earlier initiation of puberty in girls [5,6,7,8,9].

The hormones, insulin, insulin-like growth factor-1 (IGF-1), and growth hormone (GH) are linked to precocious puberty [10,11,12]. GH and the IGF signaling pathways play major roles in regulating endocrine secretions involved in growth and sex maturation. The complex IGF network involves several growth factors (IGF-1, IGF-2), high affinity insulin-like growth factor binding proteins (IGFBP3, IGFALS) and cell surface receptors (IGF-1R). Girls with CPP have higher insulin and IGF-1 blood levels than girls without PP [11,12,13]. In animal studies, insulin or IGF-1 stimulation augments GnRH secretion, confirming that insulin receptors (IRs) and IGF-IR are present on GnRH neurons. Furthermore, puberty is delayed in female mice lacking IGF-1R, but not in males [14], suggesting a role of IGF-1 signaling in timing of puberty in females.

Mutations in the MKRN3 gene and the KISS gene are highly associated with CPP [15, 16]. However, the influence of the IGF axis genes on puberty in humans, especially CPP, is unclear. As tangential evidence, girls diagnosed with PP have a higher risk for developing breast cancer than girls without PP [17]. Abnormalities in the IGF signaling pathways affect progression of breast cancer [18].

The hypothesis of this study is that polymorphisms in one or more of the IGF genes may influence sex hormonal changes and affect the development of precocious puberty. Since members of the IGF family are involved in the onset of puberty, we aimed to identify polymorphisms in insulin (INS), IGF-1, IGF-2, IGF-1 receptor (IGF1R), IGR2R, and IGF binding protein 3 (IGFBP-3) that may alter the risk for development of central precocious puberty.

Methods

Study population

This two-cohort study assessed ≤8 year old girls who showed development of secondary sex characteristics and sought care in the genetic/metabolic outpatient department at the Hospital of Chung Shan Medical University. The study included a total of 489 girls with 264 girls in the CPP group, and 225 girls in the early puberty (EP) control group. During the study’s implementation period, all girls received treatment and consultation. All girls were examined for bone age (BA), weight, and height; their blood samples were collected for estradiol (E2), luteinizing hormone (LH), and luteinizing hormone-releasing hormone (LHRH) analysis; and their BA/CA (bone age/chronological age) ratio and BMI were calculated during their first visit. Follow-up visits occurred every 3 months.

Girls were diagnosed with CPP if the girls’ bone age examination results were greater than their age growth, estradiol was ≥10 pg/ml, and highest LHRH and LH values were ≥ 10 ml U/ml [19]. Girls with CPP received consultation, health education, and treatment from a genetic/metabolic counselor. Eligibility criteria for the CPP group included continuing outpatient treatment at this hospital after outpatient examination and diagnosis with CPP.

Exclusion criteria for the CPP group consisted of refusal to sign a consent form or having a disease that may have caused CPP, such as a chromosome anomaly; a noncancerous tumor in the brain or pituitary gland, brain injury; an infection in the brain (e.g. meningitis); radiation or chemotherapy for cancer treatment.

If the girls’ bone age, estradiol, LH, and LHRH test results did not fully meet the diagnostic standard in the first visit and they only had emergent secondary sex characteristics, they were considered to have EP. The girls with EP received consultation and health education from a genetic/metabolic counselor and continued to receive outpatient follow-up. Exclusion criteria for the EP group were a GnRH homolog treatment, failure to continue to receive follow-ups, refusal of family members to sign the participation consent forms, or the aforementioned diseases or treatments in the exclusion criteria for CPP.

Ethical considerations

The study protocol was approved by the Human Investigations Committee of the Hospital of Chung Shan Medical University before the study started. Because all participating girls were 8 years old or younger, the parents or guardians of the participating girls uniformly provided signed informed consent.

Bone age, body mass index (BMI)

Left-hand X-rays were performed on all subjects, and bone age (BA) was determined using the method of Greulich and Pyle [20]. Three replicate measurements of standing height were made using a wall-mounted stadiometer. BMI was calculated by dividing body weight (kg) by the square of height in meters (m2). Tanner stage standards were used to assess breast and pubic hair development [21, 22].

Blood and data collection

After girls fasted eight-hours, blood specimens were collected for measurement of basal E2 (pmol/L), follicle stimulating hormone (FSH, U/L), LH (U/L), growth hormone (GH), IGF-1, and insulin-like growth factor binding protein 3 (IGFBP-3).

GnRH test

After the subjects had fasted overnight, venous access was secured with a three-way stopcock and heparinized saline for the GnRH test. After baseline blood samples (2 mL) were drawn for LH & FSH or estradiol measurements, GnRH (range for children: 2.5 mcg/kg to 100 mcg/kg) was administered intravenously as a bolus. After 20 min and 60 min, blood samples (2 mL) for LH & FSH assessments were drawn and serum was harvested [23, 24].

Immunoassays

Serum E2 levels were measured using a commercial radioimmunoassay kit (Diagnostic Systems Laboratories). Sensitivity was 2.2 pg/ml, with intra- and inter-assay coefficients of variation (CV) at 7.5 and 9.3%, respectively. Serum FSH and LH levels were measured by enzyme immunoassay (FSH: Abbott Laboratories, Rome, Italy; LH: Dade Behring, Milan, Italy). Sensitivity for both assays was 0.2 mIU/ml. Intra- and inter-assay CVs were 4.7 and 8.9%, respectively, for FSH, and 3.1 and 4.0%, respectively, for LH.

Quantification of serum GH, IGF-1 and IGFBP-3 levels

After acid extraction, serum GH, IGF-I and IGFBP-3 levels were measured using the commercial radioimmunoassay (RIA) kits (Diagnostic Systems Laboratories, Webster, TX, USA). The sensitivities were 0.01 ng/mL, 0.9 mg/L and 0.01 ng/mL for GH, IGF-1 and IGFBP-3, respectively. The intra- and interassay CVs were 5.3 and 5.7% for GH; 7.2 and 9.8% for IGF-1; and 5.8 and 8.2% for IGFBP-3.

Genomic DNA extraction

Venous blood from each subject was drawn into Vacutainer tubes containing EDTA and stored at 4 °C. Genomic DNA was extracted by QIAamp DNA Blood Mini Kits (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. DNA was dissolved in Tris-EDTA buffer (10 mM Tris (pH 7.8) and 1 mM EDTA) and then quantitated by a measurement of optical density at 260 nm. The final preparation was stored at − 20 °C and used as a template for PCR.

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP)

The IGF-1R, IGF-2, IGF-2R and INS gene polymorphisms were determined by PCR-RFLP assay. The sequences of primers used to amplify the related genotype and restriction enzyme for digestion as well as PCR products after digestion were listed in Table 1. PCR was performed in a 10 μL volume containing 100 ng DNA template, 1.0 μL 10´ PCR buffer (Invitrogen, Carslbad, CA, USA), 0.25 U Taq DNA polymerase (Invitrogen), 0.2 mM dNTPs (Promega, Madison, WI, USA), and 200 nM primer (MDBio Inc., Taipei, Taiwan). The PCR products of gene polymorphisms were subjected to enzymatic digestion by incubation with related restriction enzyme for 4 h at 37 °C and were then submitted to electrophoresis in 3% agarose gels.

Table 1 Sequences of primers used to amplify related genotypes, restriction enzymes for digestion, and sizes of PCR products after digestion

Real-time PCR

The IGF-1 + 1770, + 6093, and IGFBP-3 -202 genes polymorphisms were determined by real-time PCR assay. Real-time PCR based on VIC (green) and FAM (blue) fluorescent dyes were applied for accurate quantification of the target sequence. The sequences of IGF-1 + 1770, + 6093, and IGFBP-3 -202-specific primers and PCR conditions were listed in Table 1. Strength of fluorescence for each sample was detected in each reaction cycle and plotted the fluorescent values against cycle number. The quantity and polymorphisms of each gene were observed by the intensity and color of fluorescence.

Statistical analysis

Subjects’ demographics and characteristics data were represented as mean ± standard deviations (SD). The height, weight, and BMI values were used to calculate the z-scores with the WHO Child Growth Standards for subjects 5 years and younger [25] and the WHO 2007 reference for ages 5 years to 19 years [26]. The comparison of the subjects’ demographic and characteristic data between groups was performed using Mann-Whitney U test because the data were not normally distributed. Moreover, a Kruskall-Wallis test was conducted for to assess associations between two or more than two types of genotypes or combinations of genotypes and the subjects’ demographic and characteristics data. A Spearman correlation analysis was applied for identifying the correlation between serum IGF-1 and IGFBP-3 levels with subjects’ demographic and characteristic data; the coefficient of correlation r was calculated. A Pearson’s correlation was also applied for the correlation analysis with adjusted subjects’ BMI. All statistical analyses were carried out with IBM SPSS statistical software version 22 for Windows (IBM Corp., Armonk, NY, USA).

Results

This study enrolled a total of 489 girls who were classified into the CPP (n = 264) and EP control (n = 225) groups. Demographics and characteristics of CPP and control groups are summarized in Table 2. The mean chronological age of all subjects was 8.61 yrs. (SD = 1.36) and was not significantly different between groups. There were significant differences between the CPP and control groups in most demographics and characteristics, except for GH levels (Table 2).

Table 2 Demographic and clinical characteristics of CPP and control groups

We assessed the correlation of serum IGF-1 (Table 3) and IGFBP-3 (Table 4) levels with subjects’ characteristics. The serum IGF-1 level in the total population was positively correlated with most characteristics, especially bone age and LH. The IGF-1 levels in the CPP group showed a higher correlation with bone age, z-scores of height, z-scores of weight, z-scores of BMI, and LH compared to the control group. After adjusting for subjects’ BMI which differed significantly between the control and CPP groups, the IGF-1 levels in the CPP group also showed a higher correlation with bone age, CPP onset age, z-scores of height and weight, and LH levels compared to the control group (Table 3). The correlations between IGF-1 levels and the different demographic and pathological features after adjusting for bone age are summarized in Table 5.

Table 3 Correlations of IGF-1 levels with demographic and pathological features with and without BMI adjustment by group
Table 4 Correlations of IGFBP-3 levels with demographic and pathological features with and without BMI adjustment by group
Table 5 Correlations between IGF-1 levels and demographic and pathological features after adjusting for bone age by group

The serum IGFBP-3 level of the total population was significantly and positively correlated to bone age, BA/CA ratio, and FSH (Table 4). The serum IGFBP-3 level was positively correlated to FSH in the control group but not in the CPP group (Table 4), both prior to and after adjusting for the subjects’ BMI (Table 4).

There were no significant differences in the distribution of genotypes and allele frequency of IGF1R, IGF-1(6093), IGF-1(1770), IGF-2(3123), IGF2R, IGF-2(3580), IGFBP-3, and insulin between the control and CPP groups (Additional file 1: Table S1). There were also no significant differences in genotype distribution of combinations of two SNPs in two distinct genes between the CPP and control groups (Additional file 1: Tables S2–1; S2–2).

We evaluated the association of demographic and pathological features with SNP genotype in both the control and CPP groups. In the control group, bone age was significantly associated with IGF1R (p = 0.033) and IGF-2(3580) (p = 0.046); E2 levels were significantly associated with insulin (p = 0.019); and the LH levels were significantly associated with IGF-1(1770) (p = 0.020). IGF-1 levels were significantly associated with IGFBP-3 (p = 0,014) (Fig. 1; Additional file 1: Table S3–1). IGF-1 levels were significantly associated with IGFBP-3 in the CPP group (p = 0.038) (Fig. 2; Additional file 1: Table S3–2).

Fig. 1
figure 1

Data were presented as mean +-SD given SNP genotype (a. Bone age v.s IGF1R; b. LH v.s IGF-1(1770); c. Bone age v.s IGF-2(3580); d. IGF-1 v.s IGFBP-3(AA, AC, CC); e. IGF-1 v.s IGFBP-3(AC+CC, AA); f. E2 v.s Insulin)

Fig. 2
figure 2

Demographic and pathological features significantly associated with SNP genotype in CPP group. Data were presented as mean ± SD given SNP genotype and p-value was presented form difference among genotypes

We compared the association of demographic and pathological features with a combination of two genotypes in the control (Fig. 3; Additional file 1: Table S4–1) and CPP (Fig. 4; Additional file 1: Table S4–2) groups. For the control group, the serum IGF-1 level was associated with the combinations of [IGF-1(1770) + IGFBP-3] (p = 0.038), [IGF-1(6093) + IGFBP-3] (p = 0.013), and [IGF-2(3580) + IGFBP-3] (p = 0.036). The z-scores of BMI were associated with the combination, [IGF-2(3123) + IGF-2R] (p = 0.012) (Fig. 3). For CPP group, the age of onset was shown associated with combination [IGF-2(6093) + IGFBP-3] (p = 0.039) (Fig. 4). However, none of the other demographic or pathological features were significantly associated with the combination of two SNP genotypes (Additional file 1: Table S4–2).

Fig. 3
figure 3

Demographic and pathological features significantly associated with a combination of two SNP genotypes in the control group. Data were presented as mean ± SD for a given two combination of SNP genotypes. (a. IGF-1 v.s IGF-2(3580)(GG), IGF-2R(GG); b. Z score of weight v.s IGF-2 (3123)(AA), IGF2R(GG); c. IGF-1 v.s IGF-2 (1770)(TT), IGFBP-3(AA); d. IGF-1 v.s IGF-1(6093)(GG), IGFBP-3(AA); e. IGF-1 v.s IGF-2(3580)(GG), IGFBP-3(AA); f. IGF-1 v.s IGF-1(3123(AA), IGFBP-3(AA))

Fig. 4
figure 4

Demographic and pathological features significantly associated with a combination of two SNP genotypes in the CPP group. Data were presented as mean ± SD for a given combination of two SNP genotypes. p-value was presented for differences among genotypes

Furthermore, we compared the association between demographic and pathological features and a combination of the IGFBP-3 genotypes with two additional genes, in the control (Additional file 1: Table S5–1) and CPP (Additional file 1: Table S5–2) groups. In the control group, there was a significant association between IGF-1 levels and the combination, [IGF-1(6093) + IGF1R + IGFBP-3] (p = 0.026), and [IGF-2(3580) + IGF2R + IGFBP-3] (p = 0.010). The z-scores of weight was associated with the gene combination, [IGF-2(3123) + IGF2R + IGFBP-3] (p = 0.035) (Additional file 1: Table S5–1; Fig. 5). The CPP group showed a significant association between age of onset and the combination [IGF-2(3580) + IGF2R + IGFBP-3] (p = 0.020; Fig. 6). However, none of the other demographic or pathological features were significantly associated with the combination of three SNP genotypes (Additional file 1: Table S5–2).

Fig. 5
figure 5

Demographic and pathological features significantly associated with a combination of three SNP genotypes in the control group. Data were presented as mean ± SD for the given combination of three SNP genotypes. (a. IGF-1 v.s IGF-1(6093)(GG), IGF-1R(GG), IGFBP-3(AA); b. IGF-1 v.s IGF-2(3580)(GG), IGF-2R(GG), IGFBP-3(AA); c. Z score of weight v.s IGF-2(3123)(AA), IGF-2R(GG), IGFBP-3(AA))

Fig. 6
figure 6

Demographic and pathological features significantly associated with a combination of three SNP genotypes in the CPP group. Data were presented as mean ± SD for a given combination of three SNP genotypes. p-value was presented for difference among genotypes

Discussion

Our findings showed that the IGF-1 serum levels of the CPP group exhibited a higher correlation with bone age, z-scores of height, z-scores of weight, and LH than those of the control group, regardless of adjustment for BMI. The height of the CPP group, but not the z-score of the height, was associated with IGF-2(3580). The BMI of the CPP group, but not the z-score of the BMI, was associated with IGF-2(3580), IGF1R, and the combinations of [IGF-2(3580) + IGF2R], and [IGF-2(3580) + IGFBP-3]. The weight of the CPP group, but not the z-score of the weight, was associated with the combination of [IGF-2(3580) + IGFBP-3] (p = 0.024). The CPP group showed a significant association between weight and BMI with the combination of [IGF-2(3580) + IGF2R + IGFBP-3]. The distribution of genotypes and allele frequency in IGF1R, IGF-1(6093), IGF-1(1770), IGF-2(3123), IGF2R, IGF-2(3580), IGFBP-3, and insulin between control and CPP groups were similar. Likewise, the distributions of the gene SNP combinations were similar.

A number of studies have previously investigated the relationship between various IGF family SNPs and variables such as age, BMI, and weight. The IGF family plays an important role in stimulating skeletal growth, cell differentiation and metabolism, and has been shown to influence body composition [27]. IGF-2 has been reported to play a role in fetal development, while IGF-1 is expressed after birth [28]. A polymorphism in the IGF-1 promoter was reported to be associated with IGF-1 serum levels, birth weight and body height in girls with CPP, as well as in adults. This polymorphism was also shown to be associated with higher body weight, BMI, fat mass, and waist circumference in young subjects [27]. A SNP in IGF-1R was shown to influence free IGF-1 plasma concentrations. This A → G variant was predicted to generate a silent mutation, E1013E, and was associated with higher IGF-1 concentrations in Italian adults [29]. Homozygosity for the A variant was associated with the lowest mean IGF-1 concentration, whereas heterozygosity was associated with a slightly higher mean IGF-1 concentration. Homozygosity for the G variant was associated with the highest mean IGF-1 concentration. These data suggested that IGF-1 concentrations were influenced by the IGF-1R genotype at codon 1013 [29].

The AA, AG, and GG have been reported there have a dramatic difference in birth weight standard deviation scores (SDSs) in neonatal + 3123/ApaI genotypes of IGF-2 [30]. AA homozygotes had a mean birth weight SDS of 0.18 lower than that of GG homozygotes (p = 0.01), while heterozygotes showed an intermediate mean value. In contrast, there was no significant difference in birth weight SDSs between the AA, AG and GG maternal + 3123/ApaI genotypes. There was also no significant difference in birth weight SDSs between the AA, AG, and GG + 3580/MspI genotypes in both neonatal and maternal samples. However, the association between IGF-2 polymorphisms and size at birth remains controversial [30]. Analysis of the + 3123/ApaI genotype in 693 Hertfordshire adults showed that birth weight was highest for the GG genotype but the differences were not statistically significant [31]. Other data showed that the IGF-2 genotype was not significantly associated with BMI and/or birth weight in 294 healthy volunteers, but there was a statistically significant correlation between birth weight and BMI in subjects with the GG genotype whose birth weight was higher than 3.5 kg [32]. One reason for these different findings regarding the association between + 3123/ApaI polymorphism and size at birth may be differences in race and size of the study population, or differences in handling of the somatoscopic characteristics [32].

A recent study found a significant difference between the experimental and control groups in the distribution frequency of the IGF-2 + 3580 polymorphism. Additionally, multiple regression model analyses showed that the presence of the IGF-2R AA or AG genotypes may exert a protective effect against hepatitis C [odds ratio (OR) = 0.35, 95% confidence interval (CI) = 0.15–0.82]. The combination of IGF-2 + 3580 AA genotype and IGF-2R GG genotype may be associated with a significantly lower risk of HCC (OR = 0.20, 95% CI = 0.05–0.87). No polymorphisms of any IGF genes were associated with liver-related clinicopathological markers in serum [33]. Comparison of allele frequencies between the premature pubarche (PP), hyperandrogenism (HA), and healthy control subjects showed a significantly higher frequency of the G allele in the PP group compared to the other groups, (P = .0781). However, allele frequencies were comparable in the HA and the healthy control subjects [34].

IGF-1 levels are significantly higher in girls of any age undergoing puberty than those in prepubertal girls. Our results suggest that the known variants in the genes of the IGF-1 axis play minor roles in the timing of elevation of IGF-1 and IGFBP-3 protein levels during puberty, similar to the role that SNPs in IGF-1 axis genes play in breast cancer [35]. In a similar manner, although the CPP group had significantly higher leptin levels than the control group, the differences in timing and expression level of leptin could not be explained by single nucleotide polymorphisms in either leptin or the leptin receptor [36].

The similar distribution of SNPs of the genes in the IGF-1 axis between CPP and control groups were in contrast to the skewed distribution of mutations and SNPs of four genes that were more prevalent in patients with CPP: the autosomal dominant GPR54 R386P mutation [37, 38], several polymorphisms (55,648,184; 55,648,186) in the KISS1 gene [39], the intron 4 (TTTA)13 repeat in the cytochrome P450 19A1 gene CYP19A1 gene [40], and a haplotype in the 5′ promoter region of the LH β gene [41]. In addition, eight unrelated girls with CPP at age 6 had one of five novel heterozygous loss-of-function mutations in the makeorin ring finger 3 (MKRN3), which normally suppresses or delays GnRH secretion [42]. Conversely, the mutations and SNPs of four genes were less prevalent in patients with CPP than in controls: the polymorphism (55648176) in the KISS1 gene [39], the AC haplotype of Lin28B in two positions (SNPs rs4946651, RS369065) [43], and the cytochrome P450 CYP1B1 Eco571 variant (V432 L) [44]. In addition, leptin levels were significantly higher in the CPP group than those in the control group but SNPs in either the leptin receptor or leptin genes were not able to explain the differences [36]. Taken together, these studies indicate that multiple genes can influence the onset of puberty, and specific genotypes and haplotypes can increase the risk for development of CPP [39,40,41,42,43].

In our present study, we found no significant association between the SNPs evaluated and z-scores of height, weight, or BMI in either the EP or CPP groups. However, our data showed that the bone ages of subjects in the IGF-1R + 1013 (AG) and IGF-2 + 3580 (AG + AA) groups were more advanced in the EP group. This could possibly be because although the girls did not appear to have entered puberty, their bone age had already acquired the characteristics of puberty. It is possible that differences between our data and previous studies could be due to ethnic or other demographic differences in our study population. Although our data did not directly prove that IGF-1R and IGF-2 + 3580 were related to precocious puberty in girls, our results showed that the IGF-1R G variant and the IGF-2 + 3580 A variant were associated with CPP. In addition, we also believe that the interaction between IGF-I and IGF-II polymorphisms could play an important role, and warrants further investigation.

Some important limitations of this study were 1) healthy subjects were not included, 2) IGF-2 levels were not measured, and 3) the role of other genetic pathways which could play a role in CPP were not investigated.

Conclusion

In conclusion, specific genotypes from several genes (GPR54, KISS, CYP19A1, and Lin28B) can accelerate or slow the onset of puberty and have been associated with higher or lower prevalence in girls with CPP. The IGF-1 protein levels coupled with human GH levels affect timing of menarches [45,46,47]. This study showed that single SNPs of the genes in the IGF-1 axis (IGF-1(6093), IGF-1(1770), IGF1R, IGF-2(3123), IGF-2(3580), IGF2R, IGFBP-3(− 202)) did not appear to exert a significant role in the risk for development of CPP. However, several combinations were significantly associated with higher IGF-1 blood levels. Whether epigenetic modulation of the genes in the IGF axis plays a more prominent role in the risk for CPP than SNPs will require further research. Alternatively, other genes [37,38,39,40] or environmental factors [48, 49]) appear to play a more prominent role in triggering the development of CPP.

Abbreviations

BA:

Bone age

BMI:

Body mass index

CPP:

Central precocious puberty

CV:

Coefficients of variation

E2:

Estradiol

EP:

Early puberty

FSH:

Follicle stimulating hormone

GH:

Growth hormone

GnRH:

Gonadotropin-releasing hormone

IGF-1:

Insulin-like growth factor-1

IGF-1R:

Insulin-like growth factor 1 receptor

IGFALS:

Insulin-like growth factor binding protein, acid labile subunit

IGFPB3:

Insulin-like growth factor-binding protein 3

IRs:

Insulin receptors

LH:

Luteinizing hormone

LHRH:

Luteinizing hormone-releasing hormone

MKRN3:

Makeorin ring finger 3

PCR:

Polymerase chain reaction

PP:

Precocious puberty

PPP:

Pseudo or peripheral precocious puberty

RFLP:

Restriction fragment length polymorphism

RIA:

Radioimmunoassay

SD:

Standard deviation

References

  1. Berberoglu M. Precocious puberty and normal variant puberty: definition, etiology, diagnosis and current management. J Clin Res Pediatr Endocrinol. 2009;1:164–74.

    Article  Google Scholar 

  2. Chen M, Eugster EA. Central precocious puberty: update on diagnosis and treatment. Paediatr Drugs. 2015;17:273–81.

    Article  Google Scholar 

  3. Lazar L, Lebenthal Y, Yackobovitch-Gavan M, Shalitin S, de Vries L, et al. Treated and untreated women with idiopathic precocious puberty: BMI evolution, metabolic outcome, and general health between third and fifth decades. J Clin Endocrinol Metab. 2015;100:1445–51.

    Article  CAS  Google Scholar 

  4. Palmert MR, Boepple PA. Variation in the timing of puberty: clinical spectrum and genetic investigation. J Clin Endocrinol Metab. 2001;86:2364–8.

    Article  CAS  Google Scholar 

  5. Euling SY, Herman-Giddens ME, Lee PA, Selevan SG, Juul A, et al. Examination of US puberty-timing data from 1940 to 1994 for secular trends: panel findings. Pediatrics. 2008;121(Suppl 3):S172–91.

    Article  Google Scholar 

  6. Biro FM, Khoury P, Morrison JA. Influence of obesity on timing of puberty. Int J Androl. 2006;29:272–7. discussion 286-90

    Article  Google Scholar 

  7. Slyper AH. The pubertal timing controversy in the USA, and a review of possible causative factors for the advance in timing of onset of puberty. Clin Endocrinol. 2006;65:1–8.

    Article  Google Scholar 

  8. Kaplowitz PB. Link between body fat and the timing of puberty. Pediatrics. 2008;121(Suppl 3):S208–17.

    Article  Google Scholar 

  9. Kaplowitz PB, Slora EJ, Wasserman RC, Pedlow SE, Herman-Giddens ME. Earlier onset of puberty in girls: relation to increased body mass index and race. Pediatrics. 2001;108:347–53.

    Article  CAS  Google Scholar 

  10. Veldhuis JD, Roemmich JN, Richmond EJ, Bowers CY. Somatotropic and gonadotropic axes linkages in infancy, childhood, and the puberty-adult transition. Endocr Rev. 2006;27:101–40.

    Article  CAS  Google Scholar 

  11. Anderson CA, Zhu G, Falchi M, van den Berg SM, Treloar SA, et al. A genome-wide linkage scan for age at menarche in three populations of European descent. J Clin Endocrinol Metab. 2008;93:3965–70.

    Article  CAS  Google Scholar 

  12. Veldhuis JD, Roemmich JN, Richmond EJ, Rogol AD, Lovejoy JC, et al. Endocrine control of body composition in infancy, childhood, and puberty. Endocr Rev. 2005;26:114–46.

    Article  CAS  Google Scholar 

  13. Renehan AG, Frystyk J, Flyvbjerg A. Obesity and cancer risk: the role of the insulin-IGF axis. Trends Endocrinol Metab. 2006;17:328–36.

    Article  CAS  Google Scholar 

  14. Divall SA, Williams TR, Carver SE, Koch L, Bruning JC, et al. Divergent roles of growth factors in the GnRH regulation of puberty in mice. J Clin Invest. 2010;120:2900–9.

    Article  CAS  Google Scholar 

  15. Silveira LG, Noel SD, Silveira-Neto AP, Abreu AP, Brito VN, et al. Mutations of the KISS1 gene in disorders of puberty. J Clin Endocrinol Metab. 2010;95:2276–80.

    Article  CAS  Google Scholar 

  16. Abreu AP, Dauber A, Macedo DB, Noel SD, Brito VN, et al. Central precocious puberty caused by mutations in the imprinted gene MKRN3. N Engl J Med. 2013;368:2467–75.

    Article  CAS  Google Scholar 

  17. Hamilton AS, Mack TM. Puberty and genetic susceptibility to breast cancer in a case-control study in twins. N Engl J Med. 2003;348:2313–22.

    Article  Google Scholar 

  18. Mu L, Tuck D, Katsaros D, Lu L, Schulz V, et al. Favorable outcome associated with an IGF-1 ligand signature in breast cancer. Breast Cancer Res Treat. 2012;133:321–31.

    Article  CAS  Google Scholar 

  19. Garibaldi L, Chemaitilly W. Chapter 556, Disorders of pubertal development. In: Kliegman RM, BMD S, editors. Nelson textbook of pediatrics. 2. 19th ed. Philadelphia: Elsevier Saunders; 2011. p. 1886–1894.e1883.

    Google Scholar 

  20. Greulich W, Pyle S. Radiographic Altas of skeletal development of hand and wrist. Stanford: Stanford University Press; 1999.

    Google Scholar 

  21. Tanner JM, Whitehouse RH, Takaishi M. Standards from birth to maturity for height, weight, height velocity, and weight velocity: British children, 1965. I. Arch Dis Child. 1996;41:454–71.

    Article  Google Scholar 

  22. Tanner JM, Whitehouse RH, Takaishi M. Standards from birth to maturity for height, weight, height velocity, and weight velocity: British children, 1965. II. Arch Dis Child. 1966;41:613–35.

    Article  CAS  Google Scholar 

  23. Wu FC, Butler GE, Kelnar CJ, Sellar RE. Patterns of pulsatile luteinizing hormone secretion before and during the onset of puberty in boys: a study using an immunoradiometric assay. J Clin Endocrinol Metab. 1990;70:629–37.

    Article  CAS  Google Scholar 

  24. Yen SS, VandenBerg G, Rebar R, Ehara Y. Variation of pituitary responsiveness to synthetic LRF during different phases of the menstrual cycle. J Clin Endocrinol Metab. 1972;35:931–4.

    Article  CAS  Google Scholar 

  25. World health Organization. WHO child growth standards length/height-for-age, weight-for-age, weight-for-length, weight-for-height and body mass index-for-age methods and development. Geneva, Switzerland: World Health Organization,; 2007.

  26. de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, et al. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007;85:660–7.

    Article  Google Scholar 

  27. Voorhoeve PG, van Rossum EFC, te Velde SJ, Koper JW, Kemper HC, et al. Association between an IGF-I gene polymorphism and body fatness: differences between generations. Eur J Endocrinol. 2006;154:379–88.

    Article  CAS  Google Scholar 

  28. Randhawa R, Cohen P. The role of the insulin-like growth factor system in prenatal growth. Mol Genet Metab. 2005;86:84–90.

    Article  CAS  Google Scholar 

  29. Bonafè M, Barbieri M, Marchegiani F, Olivieri F, Ragno E, et al. Polymorphic variants of insulin-like growth factor I (IGF-I) receptor and phosphoinositide 3-kinase genes affect IGF-I plasma levels and human longevity: cues for an evolutionarily conserved mechanism of life span control. J Clin Endocrinol Metab. 2003;88:3299–304.

    Article  Google Scholar 

  30. Kaku K, Osada H, Seki K, Sekiya S. Insulin-like growth factor 2 (IGF2) and IGF2 receptor gene variants are associated with fetal growth. Acta Paediatr. 2007;96:363–7.

    Article  CAS  Google Scholar 

  31. Sayer AA, Syddall H, O’Dell SD, Chen XH, Briggs PJ, et al. Polymorphism of the IGF2 gene, birth weight and grip strength in adult men. Age Ageing. 2002;31:468–70.

    Article  Google Scholar 

  32. Gomes MV, Soares MR, Pasqualim-Neto A, Marcondes CR, Lobo RB, et al. Association between birth weight, body mass index and IGF2/ApaI polymorphism. Growth Hormon IGF Res. 2005;15:360–2.

    Article  CAS  Google Scholar 

  33. Weng CJ, Hsieh YH, Tsai CM, Chu YH, Ueng KC, et al. Relationship of insulin-like growth factors system gene polymorphisms with the susceptibility and pathological development of hepatocellular carcinoma. Ann Surg Oncol. 2010;17:1808–15.

    Article  Google Scholar 

  34. Roldan MB, White C, Witchel SF. Association of the GAA1013→GAG polymorphism of the insulin-like growth factor-1 receptor (IGF1R) gene with premature pubarche. Fertil Steril. 2007;88:410–7.

    Article  CAS  Google Scholar 

  35. Canzian F, McKay JD, Cleveland RJ, Dossus L, Biessy C, et al. Polymorphisms of genes coding for insulin-like growth factor 1 and its major binding proteins, circulating levels of IGF-I and IGFBP-3 and breast cancer risk: results from the EPIC study. Br J Cancer. 2006;94:299–307.

    Article  CAS  Google Scholar 

  36. Su PH, Yang SF, Yu JS, Chen SJ, Chen JY. Study of leptin levels and gene polymorphisms in patients with central precocious puberty. Pediatr Res. 2012;71:361–7.

    Article  CAS  Google Scholar 

  37. Teles MG, Bianco SD, Brito VN, Trarbach EB, Kuohung W, et al. A GPR54-activating mutation in a patient with central precocious puberty. N Engl J Med. 2008;358:709–15.

    Article  CAS  Google Scholar 

  38. Luan X, Yu H, Wei X, Zhou Y, Wang W, et al. GPR54 polymorphisms in Chinese girls with central precocious puberty. Neuroendocrinology. 2007;86:77–83.

    Article  CAS  Google Scholar 

  39. Rhie YJ, Lee KH, Ko JM, Lee WJ, Kim JH, et al. KISS1 gene polymorphisms in Korean girls with central precocious puberty. J Korean Med Sci. 2014;29:1120–5.

    Article  CAS  Google Scholar 

  40. Lee HS, Kim KH, Hwang JS. Association of aromatase (TTTA)n repeat polymorphisms with central precocious puberty in girls. Clin Endocrinol. 2014;81:395–400.

    Article  CAS  Google Scholar 

  41. Zhao Y, Chen T, Zhou Y, Li K, Xiao J. An association study between the genetic polymorphisms within GnRHI, LHbeta and FSHbeta genes and central precocious puberty in Chinese girls. Neurosci Lett. 2010;486:188–92.

    Article  CAS  Google Scholar 

  42. Macedo DB, Abreu AP, Reis AC, Montenegro LR, Dauber A, et al. Central precocious puberty that appears to be sporadic caused by paternally inherited mutations in the imprinted gene makorin ring finger 3. J Clin Endocrinol Metab. 2014;99:E1097–103.

    Article  CAS  Google Scholar 

  43. Park SW, Lee ST, Sohn YB, Cho SY, Kim SH, et al. LIN28B polymorphisms are associated with central precocious puberty and early puberty in girls. Korean J Pediatr. 2012;55:388–92.

    Article  CAS  Google Scholar 

  44. Matsuzaki CN, Junior JM, Damiani D, de Azevedo Neto RS, Carvalho KC, et al. Are CYP1A1, CYP17 and CYP1B1 mutation genes involved on girls with precocious puberty? A pilot study. Eur J Obstet Gynecol Reprod Biol. 2014;181:140–4.

    Article  CAS  Google Scholar 

  45. Gamba M, Pralong FP. Control of GnRH neuronal activity by metabolic factors: the role of leptin and insulin. Mol Cell Endocrinol. 2006;254-255:133–9.

    Article  CAS  Google Scholar 

  46. Sorensen K, Mouritsen A, Mogensen SS, Aksglaede L, Juul A. Insulin sensitivity and lipid profiles in girls with central precocious puberty before and during gonadal suppression. J Clin Endocrinol Metab. 2010;95:3736–44.

    Article  Google Scholar 

  47. Sorensen K, Aksglaede L, Petersen JH, Andersson AM, Juul A. Serum IGF1 and insulin levels in girls with normal and precocious puberty. Eur J Endocrinol. 2012;166:903–10.

    Article  Google Scholar 

  48. Wolff MS, Teitelbaum SL, Pinney SM, Windham G, Liao L, et al. Investigation of relationships between urinary biomarkers of phytoestrogens, phthalates, and phenols and pubertal stages in girls. Environ Health Perspect. 2010;118:1039–46.

    Article  CAS  Google Scholar 

  49. Wolff MS, Teitelbaum SL, McGovern K, Pinney SM, Windham GC, et al. Environmental phenols and pubertal development in girls. Environ Int. 2015;84:174–80.

    Article  CAS  Google Scholar 

Download references

Funding

This study was supported by Chung Shan Medical University Hospital (Grant No.: CSH-2015-C-017 and CSH-2016-C-029).

Availability of data and materials

All deidentified patient data were collected prospectively from the two cohorts and are available in the main text of the study report only. Because of the extensive amount of data from 489 participants, the authors are not sharing the complete deidentified raw data from the supplementary materials in online sources.

Author information

Authors and Affiliations

Authors

Contributions

CH participated in the design of the study, performed in the acquisition of data and analysis and interpretation of data, helped to draft the manuscript. YS participated in the design of the study, performed in the acquisition of data and analysis and interpretation of data, helped to draft the manuscript. WS participated in the design of the study, performed in the acquisition of data and analysis and interpretation of data. SP participated in the design of the study, performed in the acquisition of data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pen-Hua Su.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Human Investigations Committee of the Hospital of Chung Shan Medical University before the study started. Because all participating girls were 8 years old or younger, the parents or guardians of the participating girls uniformly provided signed informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Table S1. Distribution of SNPs in IGF1R, IGF-1(6093), IGF-1(1770), IGF-2(3123), IGF2R, IGF-2(3580), IGFBP-3, and insulin in CPP and control groups. Table S2–1 Summary of genotype distribution of two SNP combinations in two distinct genes by group. Table S2–2. (continued) Summary of genotype distribution in two SNP combinations in two distinct genes by group. Table S3–1. Associations between demographic and pathological features and SNP genotypes in control group. Table S3–2. Associations between demographic and pathological features and SNP genotypes in CPP group. Table S4–1. Comparison of associations between demographic and pathological features with two SNP genotype combinations in the control group. Table S4–2. Associations between demographic and pathological features and two SNP genotype combinations in the CPP group. Table S5–1. Associations between demographic and pathological features and combinations of IGFBP-3 and two additional genes in the control group. Table S5–2. Associations between demographic and pathological features and combinations of IGFBP-3 and two additional genes in the CPP group. Table S6. Summary of the power for given control and CPP groups. (DOCX 144 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chang, HP., Yang, SF., Wang, SL. et al. Associations among IGF-1, IGF2, IGF-1R, IGF-2R, IGFBP-3, insulin genetic polymorphisms and central precocious puberty in girls. BMC Endocr Disord 18, 66 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12902-018-0271-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12902-018-0271-1

Keywords