- Research article
- Open Access
Association between MKRN3 and LIN28B polymorphisms and precocious puberty
BMC Geneticsvolume 19, Article number: 47 (2018)
The present study aimed to investigate the association between MKRN3 and LIN28B gene polymorphisms and precocious puberty in Korean boys and girls.
Children 7 to 9 years of age in 2011 to 2012 who were part of the Ewha Birth & Growth Cohort Study were recruited for this study. A total of 103 girls and 70 boys were included in the analyses. Seven girls and 26 boys were identified to have precocious puberty. Among four single nucleotide polymorphisms (SNPs) of MKRN3 and two SNPs of LIN28B examined, three SNPs (rs2239669, rs6576457, and rs12441827) showed significant associations with precocious puberty in additive models in boys but no significance was found in any SNPs in girls. From the logistic regression analysis, boys with TT alleles in rs12441827 had about a four-times greater risk for precocious puberty when compared to C allele carriers (OR = 3.95, 95% CI = 1.27–12.32 in model 1). eQTL analysis revealed that SNPs of statistical significance from our study did not show the variation in expression profiles nor found in the database.
This study supports the impact of MKRN3 SNP rs12441827 on precocious puberty in Korean boys. The results add a further aspect to genetic association in precocious puberty along with complex interactions of environmental, nutritional and socioeconomic factors.
Puberty is a complex developmental process required for reproductive competency. The onset of puberty is dependent on the activation of hypothalamus-pituitary-gonadal axis. Gonadotropin-releasing hormone (GnRH) plays a crucial role in reproduction and sexual development by leading to secretion of luteinizing hormone and follicle stimulating hormone, which in turn stimulates gonadal secretion of testosterone and estradiol . Onset of puberty is regulated by complex interactions of environmental, genetic, nutritional, and socioeconomic factors [2, 3].
Precocious puberty is clinically defined by the development of secondary sexual characteristics before the age of 8 years in girls and 9 years in boys . The prevalence of precocious puberty is substantially increasing in Western countries and Asia in both girls and boys . The overall incidence of precocious puberty in Korea was 15.3 per 100,000 girls, and 0.6 per 100,000 boys between the years 2004 to 2010 . Precocious puberty is associated with many problems including early bone maturation, reduced adult height, and increased risk for the development of psychological problems . Pubertal timing varies within the population and differs between boys and girls . There is evidence of the influence of genetic factors on pubertal timing [9,10,11,12]. Recently, genetic variation in and near the makorin RING finger protein 3 (MKRN3) was reported to be associated with cases of familial central precocious puberty as well as in non-familial central precocious puberty [13, 14]. Genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) in or near the Lin-28 homolog B (LIN28B) gene as linked to age at menarche. Hence, it is possible that menarche related SNPs may be related to precocious puberty [15, 16].
Although the functions of these genes are not completely known, MKRN3 and LIN28B are regarded as important genes for the onset of puberty. However, few studies have been focused on these genes concerning precocious puberty in Korea and especially in boys. Therefore, we aimed to investigate the correlation between MKRN3 and LIN28B gene polymorphisms and precocious puberty in Korean girls and boys.
As part of the Ewha Birth & Growth Cohort Study, we conducted the present investigation to assess the effects of gene polymorphisms on precocious puberty. The Ewha Birth & Growth Cohort Study involves a prospective birth cohort established in 2001–2006. A detailed explanation of this cohort has been provided in previous studies [17, 18]. Briefly, we recruited pregnant women at gestational age 24–28 weeks who sought prenatal care at the Ewha Womans University Mokdong Hospital, Seoul, Korea. A total of 940 women participated and their children have been followed-up since 2005. The physical exam was performed every year starting 2011 in girls 7 and 8 years old and in boy 7 to 9 years old. A total of 344 parents were telephoned and asked to participate in the study. Detailed follow-up data were obtained from July to August of 2011 and peripheral blood samples were collected in EDTA-coated tubes for genotype analysis during the follow-up visits . The information collected included age, gender, height, body weight, and Tanner stage. Subjects were excluded if data on Tanner stage, height, or weight was missing. Breast and pubic hair development in girls and genital enlargement and pubic hair development in boys were evaluated by pediatrician in accordance with the Tanner stage system. Testes size was classified per Tanner staging, which stage 1 refers to testes smaller than 1.5 mL or long axis < 2.5 cm, stage 2 refers to size 1.6–6 mL or long axis 2.5–3.2 cm, and stage 3 refers to size 6–12 mL or long axis 3.3–4.0 cm .
This cohort study was performed in Ewha Womans University Mokdong Hospital. The study protocol was approved by ethics committee of the Ewha Womans University Mokdong Hospital Institutional Review Board, and written consent was obtained from the parents of each subject.
Genomic DNA was obtained from peripheral blood collected at the follow-up visit. Genomic DNA was isolated from whole blood sample using the QIAamp DNA Blood Mini Kit (QIAGEN GmbH, Hilden, Germany) following the manufacturer’s protocol. Four SNPs of the MKRN3 gene (rs2239669, rs12441827, rs657457, and rs12148769) and two SNPs of the LIN28B gene (rs314280 and rs314276) were selected based on the data from previous studies [16, 20,21,22] and considering the minor allele frequency (MAF) of 20% or above in combined Chinese and Japanese (CHB + JPT) panels from Hapmap v3 (release #2). Primer sequences of these SNPs are provided in Table 1. The genotypes of five SNPs were screened using single base primer extension assay using ABI PRISM SNaPshot Multiplex kits (ABI, Foster City, CA, USA) according to manufacturer’s protocol. LIN28B rs314276 was analyzed by sequencing due to the failure of SNaPshot assay. Genotyping results have been deposited in the TreeBASE database (Study Accession URL: http://purl.org/phylo/treebase/phylows/study/TB2:S22535).
Univariate analysis was conducted in order to identify the association of precocious puberty with factors including age, weight, height, and body mass index (BMI). Continuous variables were compared by the Independent-samples t-test. The chi-square test was done to compare categorical variables. Tests of normality were performed and determined by Shapiro-Wilk’s p value. Variables without normal distribution were analyzed by nonparametric method such as Mann-Whitney test. Each genotype of examined polymorphisms was assessed by unconditional logistic regression analyses under the additive (major allele homozygotes vs. heterozygotes or minor allele homozygotes vs. heterozygotes) models, dominant (major allele homozygotes vs. heterozygotes + minor allele homozygotes), and recessive (major allele homozygotes + heterozygotes vs. minor allele homozygotes) models, respectively. Pseudo-R-square statistics of Nagelkerke R Square were analyzed and Hosmer–Lemeshow test was used to assess the goodness of fit of multivariate model.
All statistical tests were two-sided and p-values < 0.05 were considered to be statistically significant. All analyses were performed with the Package for Social Sciences Version 23.0 for Windows (SPSS 23.0 K, SPSS INC, Chicago, IL, USA).
A total of 103 girls and 70 boys were included in the analyses. There were 96 girls (93.2%) in Tanner stage 1 at the age of 7 and 8 years, 7 girls in stage 2 at the age of 7 years (6.8%), and none in stage 3. Among boys, 44 (62.9%) were found to be in stage 1 at the age of 7 to 9 years, 21 (30.0%) in stage 2 at the age of 7 and 8 years, and 5 (7.1%) in stage 3 at the age of 7 and 8 years. Since Tanner stage 2 or above before the age of 8 years in girls and 9 years in boys represents precocious puberty, 7 girls and 26 boys were identified to have a precocious puberty. The mean age of the study participants was 7.3 years in girls and 7.9 years in boys. From univariate analysis, weight was a significant factor for the precocious puberty in girls but not in boys (Table 2).
Haploview analyses indicated high linkage disequilibrium (LD) between rs2239669 and rs6576457 (r2 = 0.97 in girls and r2 = 1 in boys) and rs314280 and rs314276 (r2 = 0.97 in girls and r2 = 0.96 in boys) in this population.
Additive model analysis revealed significant associations between three SNPs in MKRN3 (rs2239669, rs6576457, and rs12441827) and precocious puberty in boys (OR = 7.73, 9.27, and 0.17, 95% CI = 1.54–38.67, 2.89–45.48, and 0.08–0.37, respectively). Also, rs12441827 showed statistical significance in dominant models (OR = 0.18, 95% CI = 0.06–0.51) and rs2239669 and rs6576457 were in significant association with precocious puberty in recessive models (OR = 10.75 and 12.9, 95% CI = 1.78–98.15 and 1.46–114.40, respectively) in boys. In contrast, none of the SNPs showed statistical significance with precocious puberty in girls. Also, SNPs in the LIN28B gene were not associated with precocious puberty in both boys and girls (Table 3).
Multivariate stepwise regression model was used to clarify independent association between three variants ((rs2239669, rs6576457, and rs12441827) and precocious puberty in boys. The results confirmed that rs12441827 SNP was strong independent factor associated with precocious puberty. This indicated boys with TT alleles in rs12441827 have about four-times greater risk for precocious puberty when compared to boys carrying C allele (OR = 3.95, 95% CI = 1.27–12.32 in model 1) (Table 4).
In order to understand the effects of chosen SNPs on MKRN3 gene expression, we examined the expression profiling using GTEx datasets . eQTL analysis was performed on Genotype-Tissue Expression (GTEx) dataset which contains RNA sequencing and whole genome sequencing data from non-diseased human tissues across 960 donors . Unfortunately, three SNPs (rs12441827, rs6576457, rs2239669) were not found in GTEx datasets. The SNP rs12441827 was recorded as a significant expression quantitative trait loci (eQTL) with MKRN3 transcript (p = 5.1 × 10− 7), although we could not find a significant difference of the chosen genotype frequency with precocious puberty in our cohort. The expression profiles for LIN28B revealed significant regional expression differences showing lower expression with variant type alleles in rs314276 and wild type alleles in rs314280 in pituitary (p = 9.4 × 10− 2 and p = 7.6 × 10− 8, respectively). However, both SNPs of LIN28B were not shown to be associated with precocious puberty in our study.
This study aimed to investigate the association between MKRN3 and LIN28B gene polymorphisms and precocious puberty in Korean children by analyzing six SNPs. From the multivariate regression analysis, a significant association between rs12441827 mutation within the MKRN3 gene and a decreased risk for precocious puberty in boys was observed.
MKRN3 is an intronless gene located on chromosome 15q11.2, in the Prader–Willi syndrome critical region . The role of MKRN3 in pubertal initiation was first described in a genetic study of several families with precocious puberty in 2013. In this study, the authors researched 40 members of 15 families with precocious puberty from several ethnic groups (12 Brazilian families, 2 American families, and one Belgian familys) and demonstrated the effect of mutations of MKRN3 for precocious puberty. Settas et al. described a novel heterozygous missense mutation (p.Cys340Gly) in MKRN3 in two Greek siblings, a girl with central precocious puberty and a boy with early puberty . These studies further expanded the MKRN3 mutational spectrum in familial precocious puberty on top of the findings that MKRN3 gene is expressed from paternally inherited copy .
We studied MKRN3 polymorphisms in unrelated children (103 girls and 70 boys) and identified an association with precocious puberty in boys but not in girls. This could be due to the lower incidence rate of precocious puberty in girls in our cohort compared to previous studies reporting higher frequency in girls than in boys [6, 28, 29]. Another reason that we could not find an association with SNPs in MKRN3 and precocious puberty in girls could be related to the lack of information on family history. MKRN3 is a paternally expressed genes but we were not able to consider whether the SNP was inherited from the father or not.
Although the function of MKRN3 gene is not completely understood, it is known to lead the inhibition of factors promoting the pulsatile GnRH secretion . MKRN3 expression of male and female mice are comparable in hypothalamic arcuate nucleus which indicates prepubertal inhibitory tonus of MKRN3 in both sexes, but allows the possibility of gender-specific action [1, 31]. Girls begin puberty at a younger age than boys and are linked to higher rate of central precocious puberty while boys have a higher incidence of delayed puberty . It’s likely to be indicated that the inhibition of GnRH during childhood is weaker in girls than in boys allowing girls to be more sensitive to disturbances in pubertal onset.
The association of MKRN3 polymorphism and precocious puberty was studied previously and the author reported low frequency of MKRN3 mutations in central precocious puberty in Korean girls . Bessa et al. demonstrated a high frequency of MKRN3 mutations in boys when compared to previously reported female data (40% versus 6.4%, respectively, p < 0.001), which is consistent with our findings . The observed MAF of rs2239669 was 0.16 in a study done in Italian girls with idiopathic CPP, which is lower than the reported MAF of European population . Another study done in France reported higher frequency of mutations in MKRN3 coding region . These results show MKRN3 mutations are common in familial idiopathic CPP and lower in sporadic cases. Although our study included small number of children, the result is meaningful as this is the first study to investigate the association of MKRN3 and precocious puberty in boys in Korea. Also, we described novel finding that rs12441827 is associated with precocious puberty, which is the first report regarding this SNP.
We analyzed two SNPs (rs314276 and rs314280) in LIN28B and found no significant association with precocious puberty in both girls and boys in all models. Some studies reported rs314280 and rs314276 were significantly associated with precocious puberty [15, 16] while others failed to show significance in those SNPs [36, 37]. This inconsistency could be explained by the small sample size or the ethnic differences. Certainly, as other SNPs within LIN28B need to be investigated including non-coding regions. Given the finding of a study in mice that genetic modification of the Lin28-Let-7 pathway by over-expression of Lin28B increased body size and delayed the onset of puberty , involvement of LIN28B in the regulation of pubertal timing in human still warrants further studies in different patient groups and populations.
One of the limitation of the present study is the low statistical power due to the small number of samples, as calculated adequate sample size was 379 in each group [6, 13]. This warrants further investigation with larger sample size. Another limitation is that GnRH stimulation test was not performed. The diagnosis of precocious puberty was done by pediatrician through physical exams including breast development, size of the testes and penis, examination of pubic hair. Furthermore, for a comprehensive relationship between SNPs and precocious puberty, SNPs in coding region also needed to be analyzed in a near future.
In conclusion, our study has found evidence that rs12441827 of MKRN3 is associated with precocious puberty in Korean boys. The exact mechanisms involved in pubertal development are still not well understood. Therefore, functional researches and studies with large-scale samples with MKRN3 and IN28B mutations are warranted to confirm our findings.
Body mass index
Expression quantitative trait loci
Lin-28 homolog B
Minor allele frequency
Makorin RING finger protein 3
Single nucleotide polymorphism
Abreu AP, Macedo DB, Brito VN, Kaiser UB, Latronico AC. A new pathway in the control of the initiation of puberty: the MKRN3 gene. J Mol Endocrinol. 2015;54(3):R131–9.
Bulcao Macedo D, Nahime Brito V, Latronico AC. New causes of central precocious puberty: the role of genetic factors. Neuroendocrinology. 2014;100(1):1–8.
Shin YL. An update on the genetic causes of central precocious puberty. Ann Pediatr Endocrinol Metab. 2016;21(2):66–9.
Fuqua JS. Treatment and outcomes of precocious puberty: an update. J Clin Endocrinol Metab. 2013;98(6):2198–207.
Sorensen K, Mouritsen A, Aksglaede L, Hagen CP, Mogensen SS, Juul A. Recent secular trends in pubertal timing: implications for evaluation and diagnosis of precocious puberty. Horm Res Paediatr. 2012;77(3):137–45.
Kim SH, Huh K, Won S, Lee KW, Park MJ. A significant increase in the incidence of central precocious puberty among Korean girls from 2004 to 2010. PLoS One. 2015;10(11):e0141844.
Wojniusz S, Callens N, Sutterlin S, Andersson S, De Schepper J, Gies I, Vanbesien J, De Waele K, Van Aken S, Craen M, et al. Cognitive, emotional, and psychosocial functioning of girls treated with pharmacological puberty blockage for idiopathic central precocious puberty. Front Psychol. 2016;7:1053.
Cousminer DL, Widen E, Palmert MR. The genetics of pubertal timing in the general population: recent advances and evidence for sex-specificity. Curr Opin Endocrinol Diabetes Obes. 2016;23(1):57–65.
Palmert MR, Hirschhorn JN. Genetic approaches to stature, pubertal timing, and other complex traits. Mol Genet Metab. 2003;80(1–2):1–10.
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(3):188–92.
Parent AS, Matagne V, Westphal M, Heger S, Ojeda S, Jung H. Gene expression profiling of hypothalamic hamartomas: a search for genes associated with central precocious puberty. Horm Res. 2008;69(2):114–23.
Luan X, Yu H, Wei X, Zhou Y, Wang W, Li P, Gan X, Wei D, Xiao J. GPR54 polymorphisms in Chinese girls with central precocious puberty. Neuroendocrinology. 2007;86(2):77–83.
Simon D, Ba I, Mekhail N, Ecosse E, Paulsen A, Zenaty D, Houang M, Jesuran Perelroizen M, de Filippo GP, Salerno M, et al. Mutations in the maternally imprinted gene MKRN3 are common in familial central precocious puberty. Eur J Endocrinol. 2016;174(1):1–8.
Ortiz-Cabrera NV, Riveiro-Alvarez R, Lopez-Martinez MA, Perez-Segura P, Aragon-Gomez I, Trujillo-Tiebas MJ, Soriano-Guillen L. Clinical exome sequencing reveals MKRN3 pathogenic variants in familial and nonfamilial idiopathic central precocious puberty. Horm Res Paediatr. 2017;87(2):88–94.
Hu Z, Chen R, Cai C. Association of genetic polymorphisms around the LIN28B gene and idiopathic central precocious puberty risks among Chinese girls. Pediatr Res. 2016;80(4):521–5.
Ong KK, Elks CE, Li S, Zhao JH, Luan JA, Andersen LB, Bingham SA, Brage S, Smith GD, Ekelund U, et al. Genetic variation in LIN28B is associated with the timing of puberty. Nat Genet. 2009;41(6):729–33.
Min JW, Kong KA, Park BH, Hong JH, Park EA, Cho SJ, Ha EH, Park H. Effect of postnatal catch-up growth on blood pressure in children at 3 years of age. J Hum Hypertens. 2007;21(11):868–74.
Lee HA, Kim YJ, Lee H, Gwak HS, Park EA, Cho SJ, Oh SY, Ha EH, Kim HS, Park H. Association of vitamin D concentrations with adiposity indices among preadolescent children in Korea. J Pediatr Endocrinol Metab. 2013;26(9–10):849–54.
Lee HA, Kim YJ, Lee H, Gwak HS, Hong YS, Kim HS, Park EA, Cho SJ, Ha EH, Park H. The preventive effect of breast-feeding for longer than 6 months on early pubertal development among children aged 7-9 years in Korea. Public Health Nutr. 2015;18(18):3300–7.
Perry JRB, Day F, Elks CE, Sulem P, Thompson DJ, Ferreira T, He C, Chasman DI, Esko T, Thorleifsson G, et al. Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature. 2014;514(7520):92–7.
Schreiner F, Gohlke B, Hamm M, Korsch E, Woelfle J. MKRN3 mutations in familial central precocious puberty. Horm Res Paediatr. 2014;82(2):122–6.
He C, Kraft P, Chen C, Buring JE, Pare G, Hankinson SE, Chanock SJ, Ridker PM, Hunter DJ, Chasman DI. Genome-wide association studies identify loci associated with age at menarche and age at natural menopause. Nat Genet. 2009;41(6):724–8.
Carithers LJ, Ardlie K, Barcus M, Branton PA, Britton A, Buia SA, Compton CC, DeLuca DS, Peter-Demchok J, Gelfand ET, et al. A novel approach to high-quality postmortem tissue procurement: the GTEx project. Biopreserv Biobanking. 2015;13(5):311–9.
GTEx Consortium. Human genomics. The genotype-tissue expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348(6235):648–60.
Jong MT, Gray TA, Ji Y, Glenn CC, Saitoh S, Driscoll DJ, Nicholls RD. A novel imprinted gene, encoding a RING zinc-finger protein, and overlapping antisense transcript in the Prader-Willi syndrome critical region. Hum Mol Genet. 1999;8(5):783–93.
Settas N, Dacou-Voutetakis C, Karantza M, Kanaka-Gantenbein C, Chrousos GP, Voutetakis A. Central precocious puberty in a girl and early puberty in her brother caused by a novel mutation in the MKRN3 gene. J Clin Endocrinol Metab. 2014;99(4):E647–51.
Macedo DB, Abreu AP, Reis AC, Montenegro LR, Dauber A, Beneduzzi D, Cukier P, Silveira LF, Teles MG, Carroll RS, 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(6):E1097–103.
Teilmann G, Pedersen CB, Jensen TK, Skakkebaek NE, Juul A. Prevalence and incidence of precocious pubertal development in Denmark: an epidemiologic study based on national registries. Pediatrics. 2005;116(6):1323–8.
Brito VN, Latronico AC, Arnhold IJ, Mendonca BB. Update on the etiology, diagnosis and therapeutic management of sexual precocity. Arq Bras Endocrinol Metabol. 2008;52(1):18–31.
Grandone A, Cantelmi G, Cirillo G, Marzuillo P, Luongo C, Miraglia del Giudice E, Perrone L. A case of familial central precocious puberty caused by a novel mutation in the makorin RING finger protein 3 gene. BMC Endocr Disord. 2015;15:60.
Abreu AP, Dauber A, Macedo DB, Noel SD, Brito VN, Gill JC, Cukier P, Thompson IR, Navarro VM, Gagliardi PC, et al. Central precocious puberty caused by mutations in the imprinted gene MKRN3. N Engl J Med. 2013;368(26):2467–75.
Bessa DS, Macedo DB, Brito VN, Franca MM, Montenegro LR, Cunha-Silva M, Silveira LG, Hummel T, Bergada I, Braslavsky D, et al. High frequency of MKRN3 mutations in male central precocious puberty previously classified as idiopathic. Neuroendocrinology. 2017;105(1):17–25.
Lee HS, Jin HS, Shim YS, Jeong HR, Kwon E, Choi V, Kim MC, Chung IS, Jeong SY, Hwang JS. Low frequency of MKRN3 mutations in central precocious puberty among Korean girls. Horm Metab Res. 2016;48(2):118–22.
Grandone A, Capristo C, Cirillo G, Sasso M, Umano GR, Mariani M, Miraglia Del Giudice E, Perrone L. Molecular screening of MKRN3, DLK1, and KCNK9 genes in girls with idiopathic central precocious puberty. Horm Res Paediatr. 2017;88(3–4):194–200.
Bianco SD. A potential mechanism for the sexual dimorphism in the onset of puberty and incidence of idiopathic central precocious puberty in children: sex-specific kisspeptin as an integrator of puberty signals. Front Endocrinol. 2012;3:149.
Park SW, Lee ST, Sohn YB, Cho SY, Kim SH, Kim SJ, Kim CH, Ko AR, Paik KH, Kim JW, et al. LIN28B polymorphisms are associated with central precocious puberty and early puberty in girls. Korean J Pediatr. 2012;55(10):388–92.
Tommiska J, Sorensen K, Aksglaede L, Koivu R, Puhakka L, Juul A, Raivio T. LIN28B, LIN28A, KISS1, and KISS1R in idiopathic central precocious puberty. BMC Res Notes. 2011;4:363.
Zhu H, Shah S, Shyh-Chang N, Shinoda G, Einhorn WS, Viswanathan SR, Takeuchi A, Grasemann C, Rinn JL, Lopez MF, et al. Lin28a transgenic mice manifest size and puberty phenotypes identified in human genetic association studies. Nat Genet. 2010;42(7):626–30.
This work was supported by the National Research Foundation of Korea Grant funded by the Korean Government (No. 2010–0026225). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Genotyping results have been deposited in the TreeBASE database (Study Accession URL: http://purl.org/phylo/treebase/phylows/study/TB2:S22535).
Ethics approval and consent to participate
The study protocol was approved by ethics committee of the Ewha Womans University Mokdong Hospital Institutional Review Board, and written consent was obtained from the parents of each subject.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.