216 Dental Journal (Majalah Kedokteran Gigi) 2021 December; 54(4): 216–220 Original article INTRODUCTION Class II skeletal malocclusion is often encountered in orthodontics and is one of the main reasons for patients to seek treatment. Previous studies have shown that the occurrence of class II malocclusion is quite high among populations around the world.1 A study in Denmark found a 25% occurrence of class II skeletal malocclusion in their population.2 In Asia, a study with people of Mongoloid ethnicity that showed 33.1% of this population had class II malocclusion3 and another study in India showed 30.1% class II malocclusion in the population.1 A study in Singapore with people of Chinese ethnicity showed a 23.1% occurrence of class II malocclusion4 and 70% occurrence was found in Iran.5 Rosenblum6 found out that in 103 class II skeletal malocclusion cases, 56.3% had maxillary prognathism and only 27% showed mandibular retrognathism. The differentiation of the osteoblast from its progenitor requires the activity of specific transcription factors which are expressed during development. Runt-related transcription factor-2 (Runx2) is a member Runt family of transcription factors that is essential for skeletal development.7-9 Runx2 is essential for osteoblast differentiation in endochondral and intramembranous ossification.10 The Runx gene consists of Runx1, Runx2 and Runx3. Runx2 plays an important role in osteogenesis, cartilage formation, cell migration, vascular bone invasion and tooth formation.8,11,12 Runx2 is essential Dental Journal (Majalah Kedokteran Gigi) p-ISSN: 1978-3728; e-ISSN: 2442-9740. Accredited No. 32a/E/KPT/2017. Open access under CC-BY-SA license. Available at https://e-journal.unair.ac.id/MKG/index DOI: 10.20473/j.djmkg.v54.i4.p216–220 Runx2 rs59983488 polymorphism in class II malocclusion in the Indonesian subpopulation Fadli Jazaldi1, Benny M. Soegiharto1, Astrid Dinda Hutabarat2, Noertami Soedarsono3 and Elza Ibrahim Auerkari3 1Department of Orthodontics, 2Orthodontic Postgraduate Program, 3Department of Oral Biology, Faculty of Dentistry, University of Indonesia, Jakarta, Indonesia ABSTRACT Background: Class II malocclusion is one of the main orthodontic issues for patients in seeking treatment. The prevalence of class II malocclusion varies in different populations. Variation in skeletal profile is mainly controlled internally by a regulatory gene. Runt-related transcription factor-2 (Runx2) plays a role in osteoblast differentiation and is highly expressed during development. Purpose: This study aimed to evaluate the relation of regulatory gene variation in the Runx2 promoter with class II malocclusion. Methods: DNA samples were acquired from 95 orthodontic patients in Jakarta, Indonesia, who were divided into two groups: class I skeletal malocclusion (control group) and class II malocclusion. A single nucleotide polymorphism was investigated using the polymerase chain reaction and restriction fragment length polymorphism techniques. The distribution of alleles was assessed using the Hardy-Weinberg test. The relationship between polymorphism and skeletal variation was assessed with the Chi-Square test and logistic regression. Results: The frequency distributions of genotypes and alleles were tested for Hardy-Weinberg equilibrium and found to be slightly deviated. There was an equal distribution of G and T alleles throughout class II and class I skeletal malocclusions and the Chi-Square test showed that this relationship was not significant (p=0.5). Conclusion: Runx2 rs59983488 polymorphism was found in the Indonesian subpopulation; however, an association between Runx2 rs59983488 polymorphism and class II skeletal malocclusion was not found. Keywords: class II; G330T; Indonesia; polymorphism; Runx2; rs59983488 Correspondence: Fadli Jazaldi, Department of Orthodontics, Faculty of Dentistry, University of Indonesia. Jl. Salemba Raya No.4, Jakarta Pusat, Jakarta 10430, Indonesia. Email: fadli.jazaldi@ui.ac.id mailto:fadli.jazaldi@ui.ac.id https://e-journal.unair.ac.id/MKG/index http://dx.doi.org/10.20473/j.djmkg.v54.i4.p216-220 217Jazaldi et al./Dent. J. (Majalah Kedokteran Gigi) 2021 December; 54(4): 216–220 in bone development and chondrogenesis as an upstream controller for several effector genes. Runx2 is required in the preosteoblast maturation process.10,11,13 Runx2 regulates matrix protein genes and Indian hedgehog transcription factor, an important regulator for chondrogenesis.11,14 Runx2 is upregulated by bone morphogenetic proteins (BMPs), fibroblast growth factors (FGFs) and retinoid acid. The Runx2 gene is located at 6p21.1 with 250 kb spans. Runx2 consists of two transcription factors, the distal P1 promoter and the proximal P2. Runx2 in osteoblasts was found mainly as transcripts from the P1 promoter. Both promoters are specifically expressed in osteoblast and chondrocytes. 8,11,15,16 Genetic variation of Runx2 has been investigated for its relation to skeletal variation and bone- related disease.2,17 Single nucleotide polymorphism (SNP) is the most common variation that occurs in the human genome and happens in more than 1% of the population.18 The association of Runx2 rs59983488 polymorphism with class II malocclusion was inconsistent in previous studies,16,19 therefore it is important to do further research to determine this relationship in this population. This study aimed to assess the relationship of Runx2 rs59983488 polymorphism to class II skeletal malocclusion in the Indonesian subpopulation. MATERIALS AND METHODS The research participants were orthodontic patients aged between 20 and 45 years from the University of Indonesia’s Dental Hospital. All study participants signed an informed consent document. Ethical clearance of this study was granted by the Ethical Committee, Faculty of Dentistry, University of Indonesia (no. 26/Ethical Approval/FKGUI/ VII/2015, Protocol number: 09370515). Age, sex, ethnicity and socioeconomic status were recorded. Malocclusion classification was based on sagittal relationship of Steiner’s Analysis measuring Sella-Nasion-A point (SNA), Sella- Nasion-B point (SNB) and A point-Nasion-B Point (ANB). Study participants included 48 class II and 52 class I malocclusion patients. Class I subjects were recruited with SNA 82°±2°, SNB 78°±2°, ANB 2°±2°. Class II subjects were recruited as: class II with mandibular retrusion (SNA 82°±2°, SNB <76°, ANB >4°, Overjet >4 mm), class II with maxillary protrusion (SNA >84°, SNB 78°±2°, ANB >4°, Overjet >4 mm), and class II with maxillary protrusion and mandibular retrusion (SNA >84°, SNB <76°, ANB >4°, Overjet >4 mm). All lateral cephalometric images were traced the by same operator (F.J), who is an experienced orthodontist. To reduce the effect of operator measurement error, the tracing was analysed three separate times at two weeks intervals with single blinding. Genomic DNA was collected from the peripheral blood of all 95 subjects. DNA extraction was performed by using 3 mL peripheral blood and 9 mL red blood cell lysis solution (1.45 M NH3Cl, 5 mM anhydrous EDTA, and 0.1 KHCO3, [InstaGene Matrix, Bio-Rad, California, United State of America]) incubated at room temperature for 10 min. The mixture was then centrifuged at 1500 rpm for 10 min and the supernatant was removed. 2 mL cell lysis solution was added and incubated at 37℃ for 60 min. 1.3 ml protein precipitation (PP) solution [Promega, Wisconsin, United State of America] was added and then vortexed and centrifuged at 40℃ at 3000 rpm for 5 min. The supernatant was transferred into a new tube containing 2.3 mL isopropanol, and then inverted to dry. This process was repeated until the DNA turned white. After that, the DNA was rehydrated with 200–300 µL Tris-HCl EDTA, incubated in 37℃ for 2 h and stored at -20℃. Runx2 rs59983488 G330T polymorphism was detected by using the polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) method. PCR amplification was done by using the forward primer 5’-AAA GCA AAG GAG GTT GAC CGG-3’ and reverse primer 5’-CCC TGC CCT TCT TTC TCT CTC-3’ with Perkin Elmer Gene Amp PCR system 9700 [Kumamoto, Japan]. Each reagent contained 10 µL buffer solution [Promega, Wisconsin, United State of America], 2 µL forward and reverse primers, 7.5 µL ddH2O, and 0.5 µL genomic DNA. The conditions for this PCR were 94℃ for 6 min followed by 35 cycles of the amplification denaturation phase at 94℃ for 60 s, the annealing phase at 62℃ for 30 s, elongation at 72℃ for 30 s and final elongation at 72℃ for 5 min. The PCR product was digested with BSaJ1 enzyme restriction [New England Biolabs (NEB), Massachusetts, United States of America]. The mixture consisted of 1 µL restriction enzyme BSaJ1 (10U/µL), 2 µL enzyme buffer and 18 µL ddH2O was added to 10 µL PCR product. The mixture was incubated at 60℃ for 4 h and inactivated by incubating at 80℃ for 20 min. The restriction products were subjected to electrophoresis on 2% agarose gel [Promega, Wisconsin, United State of America] in 1x TAE buffer solution (0.04 M Tris–acetate, 0.002 M EDTA, pH 8.0) at 80 V, 400 mA, for 60 min, resulting in three possible genotypes: GG, GT and TT. The restriction products were visualised with GelDoc 2000 [Bio-Rad, California, United State of America] and showed qualitative data for the G allele of 205 bp and 20 bp and the T allele of 225 bp. This study used the Hardy-Weinberg equilibrium to evaluate genetic variation that occurred in the population. Intra-observer reliability was evaluated using the intraclass correlation coefficient (ICC). Chi-Square and logistic regression tests were used to assess the relation of Runx2 rs59983488 polymorphism with class II skeletal malocclusion. Statistical analysis was performed using the Program for Social Science (IBM SPSS) version 22 [Chicago, United State of America]. RESULTS A total of 95 orthodontic patients, consisting of 48 class II skeletal malocclusion cases and 52 class I skeletal Dental Journal (Majalah Kedokteran Gigi) p-ISSN: 1978-3728; e-ISSN: 2442-9740. Accredited No. 32a/E/KPT/2017. Open access under CC-BY-SA license. Available at https://e-journal.unair.ac.id/MKG/index DOI: 10.20473/j.djmkg.v54.i4.p216–220 https://e-journal.unair.ac.id/MKG/index http://dx.doi.org/10.20473/j.djmkg.v54.i4.p216-220 218 Jazaldi et al./Dent. J. (Majalah Kedokteran Gigi) 2021 December; 54(4): 216–220 malocclusion cases, were recruited at the Dental Hospital, Faculty of Dentistry, University of Indonesia. The reliability of the cephalometric measurement was good (ICC > 0.820). Genotyping using the PCR-RFLP methods showed G allele cuts in 205 bp and 20 bp and T allele cuts in 225 bp (Figure 1). Therefore, the GG genotype (wildtype) showed two bands, the GT genotype (heterozygous) showed three bands and the TT genotype (mutant) showed a single band. The distributions of genotypes and alleles are shown in Table 1. This population was tested against Hardy-Weinberg equilibrium and showed a slight deviation (p=0.04959). The distributions of the genotype and allele MOI (class I skeletal malocclusion) and MOII (class II skeletal malocclusion) were assessed with the Chi-Square test to identify their relationship within each group and were not found to be significant. The distributions of genotype and allele G and T were equal in class II skeletal malocclusion and the control (class I skeletal malocclusion). A logistic regression test was used and found significant (p<0.05) with a slight reduction of the odds ratio 0.908 (95% Confidence Interval :0.51–1.617). DISCUSSION Runx2 rs59983488 shows bases alteration from G to T at -330 in the P1 region. Genotyping using PCR-RFLP methods showed GG genotype cuts in 205 bp and 20 bp, GT genotype cuts in 225 bp, 205 bp and 20 bp, and TT genotype cuts in 225 bp (Figure 1). The 20 bp fragment is hard to see in Figure 1 because traditional agarose gels are most effective in separating DNA fragments between 100 bp and 25 kb. Agarose gel is popular for separating moderate- to large-sized nucleic acids which have a wide range of separations. There are few alternate methods for separating small DNA fragments. Polyacrylamide gel can be used as an alternative since this gel has a higher concentration resulting in better resolution when run vertically.20 Another method to separate smaller DNA fragments would be to use agarose gel doped with graphene oxide.21 The genotype and allotype frequency distributions in this study are slightly deviated from the Hardy-Weinberg equilibrium (HWE) (p<0.05). This finding is consistent with previous Runx2 studies in Indonesian and Chinese populations.22,23 Mokhtar et al.24 found different results Table 1. Distributions of genotypes and alleles of Runx2 T330G Genotype MOI MOII P n % n % TT 7 13.5 6 13.9 0.831 GT 32 61.5 24 55.9 GG 13 25 13 30.2 T allele 46 44.2 36 41.8 0.5 G allele 58 55.8 50 58.2 Figure 1. Genotyping Runx2 rs59983488; A and D. GT genotype (225 bp, 205 bp, 20 bp); B and F. TT genotype (225 bp); C. GG genotype (205 bp, 20 bp); Lad. Ladder 50 bp; G. positive control; H. negative control. MOI: Class I skeletal malocclusion; MOII: Class II skeletal malocclusion. Chi-Square test, p<0.05 was considered significant. Logistic regression test p<0.05, OR 0.908, Confidence Interval 95% 0.51–1.617 Dental Journal (Majalah Kedokteran Gigi) p-ISSN: 1978-3728; e-ISSN: 2442-9740. Accredited No. 32a/E/KPT/2017. Open access under CC-BY-SA license. Available at https://e-journal.unair.ac.id/MKG/index DOI: 10.20473/j.djmkg.v54.i4.p216–220 https://e-journal.unair.ac.id/MKG/index http://dx.doi.org/10.20473/j.djmkg.v54.i4.p216-220 219Jazaldi et al./Dent. J. (Majalah Kedokteran Gigi) 2021 December; 54(4): 216–220 in the Malaysian population that showed consistency of genotypes and alleles with the HWE. Another study of Runx2 polymorphism in Indonesia has been done in relation to osteoporosis, which had similar findings to ours, but showing a higher deviation in HWE value.23 This study involved a different population from that study, although they shared the same location, and therefore the populations had similar characteristics. The reason for this deviation was not clear. We suspect it was due to population stratification and a limited number of study participants.25-27 The population in this study was relatively small – 95 subjects in total, with 13 having the TT genotype (13.7%) and only 6 (0.06%) in the class II skeletal malocclusion group. In this study, an association between Runx2 rs59983488 polymorphism and class II skeletal malocclusion was not found. The T and G alleles were spread evenly in both groups of class I and class II malocclusion. Küchler et al.19 found a different result in the Caucasian population. He studied two sites of Runx2 promoter and found an association between rs59983488 polymorphism and maxillary protrusion, although he did not find any association between rs1200425 and either maxillary protrusion or class II malocclusion. The study also found a strong interaction between Runx2 rs59983488 SNP and SMAD6 rs3934908 SNP in the class II malocclusion model.19 Another study in Malaysia investigated the relation of Runx2 rs6930053 polymorphism with class II malocclusion and found a significant difference in both allele and genotype.24 Bigger sample size and additional maxillary measurement are recommended to evaluate any association of Runx2 SNP with class II skeletal malocclusion and maxillary growth, and the possible mechanism for this. Runx2 polymorphism at various sites has been widely investigated in its relation to different phenotypes of skeletal growth (Table 2). Chang et al.28 reported that there was an association between Runx2 polymorphism and ossification of the posterior longitudinal ligament. Auerkari et al.23 and Bustamante et al.16 reported no association between Runx2 rs59983488 polymorphism and osteoporosis. There were also studies that reported on the relationship between Runx2 polymorphism and non-syndromic cleft lip and palate.22,29 There are still inconsistencies regarding the relationship of polymorphism in Runx2 with skeletal phenotypes, and further research is needed to assess the role of Runx2 in skeletal growth and development. The level of Runx2 expression and its activity fluctuate during the developmental process. The activity of the Runx2 protein is regulated by various post-translational modifications, such as phosphorylation, methylation and acetylation or ubiquitination. BMP-2 activity was reported to stimulate Runx2 acetylation.11 Post-translational regulation by phosphorylation and acetylation is essential for FGF-induced cranial development, and this has been investigated for premature fusion of cranial sutures and the midfacial hypoplastic changes that occur because of lack of growth at the fused sutures. Enhanced stability of Runx2 by post-translational regulation could offer a methodology for bone regeneration therapy.30 One study shows that epigenetic histone modifications on receptors for vitamin D (VDR) and bone cell differentiation affect the binding of Runx2 to the genome, modify and restrict patterns of gene expression and alter cellular response to the vitamin D hormone.31 Yan et al.32 found that METTL3 could directly induce m6A methylation of Runx2 mRNA to enhance its cellular stability and indirectly upregulate the cellular level of Runx2 by m6A methylation of pre-miR-320. Ling et al.33 reported that nicotinamide phosphoribosyltransferase (NAMPT) promotes osteogenesis through the epigenetic regulation of Runx2 expression, thus upregulating Runx2 expression. Future studies involving Runx2 post- translational regulation modification should be carried out to investigate its relation with maxillary growth and class II skeletal malocclusion. Table 2. Runx2 polymorphism studies RS P phenotype Reference rs967588 rs16873379 rs3749863 rs6908650 rs1321075 rs1406846 rs2677108 C>T T>C A>C G>C A>C T>A C>T Not significant Significant Not significant Not significant Not significant Significant Significant Ossification posterior longitudinal ligament Chang et al., 201728 rs59983488 G>T Not significant Osteoporosis Auerkari et al., 201423 deletion 17A>11A Significant Bone fracture Morrison et al., 201334 rs6930053 Significant Class II malocclusion Mokhtar et al., 201824 rs545239 rs1200425 rs16873396 Significant Non-syndromic cleft lip and palate Wu al., 201222 rs1934328 Significant Non-syndromic cleft lip and palate Jung et al.,201429 rs59983488 rs1200425 rs59983488 rs59983488 G>T G>A G>T G>T Significant Not significant Not Significant Not Significant Maxillary protrusion Osteoporosis Class II malocclusion Küchler et al., 202019 Bustamante et al., 200716 This study Dental Journal (Majalah Kedokteran Gigi) p-ISSN: 1978-3728; e-ISSN: 2442-9740. Accredited No. 32a/E/KPT/2017. Open access under CC-BY-SA license. Available at https://e-journal.unair.ac.id/MKG/index DOI: 10.20473/j.djmkg.v54.i4.p216–220 https://e-journal.unair.ac.id/MKG/index http://dx.doi.org/10.20473/j.djmkg.v54.i4.p216-220 220 Jazaldi et al./Dent. J. (Majalah Kedokteran Gigi) 2021 December; 54(4): 216–220 In conclusion, Runx2 rs59983488 polymorphism was detected in the Indonesian subpopulation. However, a relationship between the Runx2 polymorphism promoter and class II skeletal malocclusion was not found. A larger sample of this population and an additional maxillary component of measurement are recommended to further assess the link between Runx2 rs59983488 polymorphism and class II skeletal malocclusion and maxillary growth. ACKNOWLEDGEMENT This study was supported by Post Graduate Grant 2015, DPRM, University of Indonesia. REFERENCES 1. Joshi N, Hamdan AM, Fakhouri WD. Skeletal malocclusion: a developmental disorder with a life-long morbidity. J Clin Med Res. 2014; 6(6): 399–408. 2. Proffit WR, Fields HW, Larson B, Sarver DM. Contemporary orthodontics. 6th ed. St. Louis: Mosby; 2018. p. 685–715. 3. Wahab RMA, Idris H, Yacob H, Ariffin SHZ. Cephalometric and malocclusion analysis of Kadazan dusun ethnic orthodontic patients. Sains Malaysiana. 2013; 42: 25–32. 4. Chan GXL, Tan ELY, Chew MT, Wong HC, Foong KWC, Yow M. Prevalence of Class I, II and III skeletal relationships and its association with dental anomalies in an ethnic Chinese orthodontic population. Proc Singapore Healthc. 2021; : 201010582110007. 5. Oshagh M, Ghaderi F, Pakshir HR, Baghmollai AM. Prevalence of malocclusions in school-age children attending the orthodontics department of Shiraz University of Medical Sciences. East Mediterr Heal J. 2010; 16(12): 1245–50. 6. Rosenblum RE. Class II malocclusion: mandibular retrusion or maxillary protrusion? Angle Orthod. 1995; 65(1): 49–62. 7. Takarada T, Hinoi E, Nakazato R, Ochi H, Xu C, Tsuchikane A, Takeda S, Karsenty G, Abe T, Kiyonari H, Yoneda Y. An analysis of skeletal development in osteoblast-specific and chondrocyte-specific runt-related transcription factor-2 (Runx2) knockout mice. J Bone Miner Res. 2013; 28(10): 2064–9. 8. Komori T. Roles of Runx2 in skeletal development. In: Groner Y, Ito Y, Liu P, Neil J, Speck N, van Wijnen A, editors. RUNX Proteins in Development and Cancer Advances in Experimental Medicine and Biology. Singapore: Springer; 2017. p. 83–93. 9. Lam K, Zhang D-E. RUNX/CBF transcription factors☆. In: Reference Module in Biomedical Sciences. Elsevier; 2015. p. 1–5. 10. Long F. Building strong bones: molecular regulation of the osteoblast lineage. Nat Rev Mol Cell Biol. 2012; 13(1): 27–38. 11. Wysokinski D, Pawlowska E, Blasiak J. RUNX2: A master bone growth regulator that may be involved in the DNA damage response. DNA Cell Biol. 2015; 34(5): 305–15. 12. Bruderer M, Richards R, Alini M, Stoddart M. Role and regulation of RUNX2 in osteogenesis. Eur Cells Mater. 2014; 28: 269–86. 13. Kawane T, Qin X, Jiang Q, Miyazaki T, Komori H, Yoshida CA, Matsuura-Kawata VK dos S, Sakane C, Matsuo Y, Nagai K, Maeno T, Date Y, Nishimura R, Komori T. Runx2 is required for the proliferation of osteoblast progenitors and induces proliferation by regulating Fgfr2 and Fgfr3. Sci Rep. 2018; 8(1): 13551. 14. Ono K, Hata K, Nakamura E, Ishihara S, Kobayashi S, Nakanishi M, Yoshida M, Takahata Y, Murakami T, Takenoshita S, Komori T, Nishimura R, Yoneda T. Dmrt2 promotes transition of endochondral bone formation by linking Sox9 and Runx2. Commun Biol. 2021; 4(1): 326. 15. Komori T. The functions of Runx family transcription factors and Cbfb in skeletal development. Oral Sci Int. 2015; 12(1): 1–4. 16. Bustamante M, Nogués X, Águeda L, Jurado S, Wesselius A, Cáceres E, Carreras R, Ciria M, Mellibovsky L, Balcells S, Díez- Pérez A, Grinberg D. Promoter 2 -1025 T/C polymorphism in the RUNX2 gene is associated with femoral neck BMD in Spanish postmenopausal women. Calcif Tissue Int. 2007; 81(4): 327–32. 17. Carlson DS, Buschang PH. Craniofacial growth and development: Developing a perspective. In: Graber LW, Vanarsdall RL, Vig KWL, Huang GJ, editors. Orthodontics: Current principles and techniques. 6th ed. St. Louis: Mosby Elsevier; 2017. p. 1–30. 18. Alberts B, Hopkin K, Johnson A, Morgan D, Raff M, Roberts K, Walter P. Essential cell biology. 5th ed. UK: W. W. Norton; 2018. p. 864. 19. Küchler EC, Reis CLB, Carelli J, Scariot R, Nelson‐Filho P, Coletta RD, Paza AO, Matsumoto MAN, Proff P, Kirschneck C. Potential interactions among single nucleotide polymorphisms in bone‐ and cartilage‐related genes in skeletal malocclusions. Orthod Craniofac Res. 2021; 24(2): 277–87. 20. L e e PY, Cost u mbrado J, Hsu C-Y, K i m Y H. Aga rose gel electrophoresis for the separation of DNA fragments. J Vis Exp. 2012; (62): 3923. 21. Li J, Yang Y, Mao Z, Huang W, Qiu T, Wu Q. Enhanced resolution of DNA separation using agarose gel electrophoresis doped with graphene oxide. Nanoscale Res Lett. 2016; 11(1): 404. 22. Wu T, Fallin MD, Shi M, Ruczinski I, Liang KY, Hetmanski JB, Wang H, Ingersoll RG, Huang S, Ye X, Wu-Chou Y-H, Chen PK, Jabs EW, Shi B, Redett R, Scott AF, Murray JC, Marazita ML, Munger RG, Beaty TH. Evidence of gene-environment interaction for the RUNX2 gene and environmental tobacco smoke in controlling the risk of cleft lip with/without cleft palate. Birth Defects Res Part A Clin Mol Teratol. 2012; 94(2): 76–83. 23. Auerkari EI, Suryandari DA, Umami SS, Kusdhany LS, Siregar TWA, Rahardjo TBW, Talbot C, Hogervorst E. Gene promoter polymorphism of RUNX2 and risk of osteoporosis in postmenopausal Indonesian women. SAGE Open Med. 2014; 2: 205031211453157. 24. Mokhtar KI, Abu Bakar N, Kharuddin AF. Determination of RUNX2 single nucleotide polymorphism rs6930053 in Class I, II And III malocclusions. IIUM Med J Malaysia. 2017; 16(2): 1060. 25. Salanti G, Amountza G, Ntzani EE, Ioannidis JPA. Hardy–Weinberg equilibrium in genetic association studies: an empirical evaluation of reporting, deviations, and power. Eur J Hum Genet. 2005; 13(7): 840–8. 26. Chen B, Cole JW, Grond-Ginsbach C. Departure from Hardy Weinberg equilibrium and genotyping error. Front Genet. 2017; 8: 167. 27. Zintzaras E. Impact of Hardy–Weinberg equilibrium deviation on allele-based risk effect of genetic association studies and meta- analysis. Eur J Epidemiol. 2010; 25(8): 553–60. 28. Chang F, Li L, Gao G, Ding S, Yang J, Zhang T, Zuo G. Role of Runx2 polymorphisms in risk and prognosis of ossification of posterior longitudinal ligament. J Clin Lab Anal. 2017; 31(4): e22068. 29. Jung SH, Lee A-Y, Park JW, Baek S-H, Kim YH. Investigation of parental transmission of RUNX2 single nucleotide polymorphism and its association with nonsyndromic cleft lip with or without palate. Cleft Palate-Craniofacial J. 2014; 51(2): 234–9. 30. Kim W-J, Shin H-L, Kim B-S, Kim H-J, Ryoo H-M. RUNX2- modifying enzymes: therapeutic targets for bone diseases. Exp Mol Med. 2020; 52(8): 1178–84. 31. Pike JW, Meyer MB, St. John HC, Benkusky NA. Epigenetic histone modifications and master regulators as determinants of context dependent nuclear receptor activity in bone cells. Bone. 2015; 81: 757–64. 32. Yan G, Yuan Y, He M, Gong R, Lei H, Zhou H, Wang W, Du W, Ma T, Liu S, Xu Z, Gao M, Yu M, Bian Y, Pang P, Li X, Yu S, Yang F, Cai B, Yang L. m6A methylation of precursor-miR-320/RUNX2 controls osteogenic potential of bone marrow-derived mesenchymal stem cells. Mol Ther - Nucleic Acids. 2020; 19: 421–36. 33. Ling M, Huang P, Islam S, Heruth DP, Li X, Zhang LQ, Li D-Y, Hu Z, Ye SQ. Epigenetic regulation of Runx2 transcription and osteoblast differentiation by nicotinamide phosphoribosyltransferase. Cell Biosci. 2017; 7(1): 27. 34. Morrison NA, Stephens AS, Osato M, Pasco JA, Fozzard N, Stein GS, Polly P, Griffiths LR, Nicholson GC. Polyalanine repeat polymorphism in RUNX2 is associated with site-specific fracture in post-menopausal females. Grant SFA, editor. PLoS One. 2013; 8(9): e72740. Dental Journal (Majalah Kedokteran Gigi) p-ISSN: 1978-3728; e-ISSN: 2442-9740. Accredited No. 32a/E/KPT/2017. Open access under CC-BY-SA license. Available at https://e-journal.unair.ac.id/MKG/index DOI: 10.20473/j.djmkg.v54.i4.p216–220 https://e-journal.unair.ac.id/MKG/index http://dx.doi.org/10.20473/j.djmkg.v54.i4.p216-220