Next Article in Journal
Interaction between Omeprazole and Gliclazide in Relation to CYP2C19 Phenotype
Next Article in Special Issue
Influence of CYP2C9 Genetic Polymorphisms on the Pharmacokinetics of Losartan and Its Active Metabolite E-3174: A Systematic Review and Meta-Analysis
Previous Article in Journal
S100A6, Calumenin and Cytohesin 2 as Biomarkers for Cutaneous Involvement in Systemic Sclerosis Patients: A Case Control Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 5
10.3390/jpm11050369
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parathyroid Hormone Gene and Genes Involved in the Maintenance of Vitamin D Levels Association with Mandibular Retrognathism

by
Erika Calvano Küchler
1,
Caio Luiz Bitencourt Reis
2,
Guido Marañón-Vásquez
3,
Paulo Nelson-Filho
2,
Mírian Aiko Nakane Matsumoto
2,
Maria Bernadete Sasso Stuani
2,
Maria Angélica Hueb de Menezes Oliveira
4,
Peter Proff
1 and
Christian Kirschneck
1,*
1
Department of Orthodontics, University of Regensburg, Franz-Josef-Strauss-Allee 11, 93053 Regensburg, Germany
2
Department of Pediatric Dentistry, School of Dentistry of Ribeirão Preto, University of São Paulo, Av. do Café S/N, Ribeirão Preto 14040-904, Brazil
3
Department of Pediatric Dentistry and Orthodontics, School of Dentistry, Federal University of Rio de Janeiro, R. Prof. Rodolpho Paulo Rocco 325, Rio de Janeiro 21941-617, Brazil
4
School of Dentistry, University de Uberaba, R. Frei Paulino 30, Uberaba 38025-180, Brazil
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(5), 369; https://doi.org/10.3390/jpm11050369
Submission received: 29 March 2021 / Revised: 26 April 2021 / Accepted: 30 April 2021 / Published: 2 May 2021
(This article belongs to the Special Issue Clinical Utility of Pharmacogenetic Testing)
Browse Figures
Versions Notes

Abstract

:
In this study we evaluated whether single nucleotide polymorphisms (SNPs) in the genes encoding PTH, VDR, CYP24A1, and CYP27B1 were associated with mandibular retrognathism (MR). Samples from biologically-unrelated Brazilian patients receiving orthodontic treatment were included in this study. Pre-orthodontic lateral cephalograms were used to determine the phenotype. Patients with a retrognathic mandible were selected as cases and those with an orthognathic mandible were selected as controls. Genomic DNA was used for genotyping analysis of SNPs in PTH (rs694, rs6256, and rs307247), VDR (rs7975232), CYP24A1 (rs464653), and CYP27B1 (rs927650). Chi-squared or Fisher’s tests were used to compare genotype and allele distribution among groups. Haplotype analysis was performed for the SNPs in PTH. The established alpha was p < 0.05. Multifactor dimensionality reduction (MDR) was used to identify SNP–SNP interactions. A total of 48 (22 males and 26 females) MR and 43 (17 males and 26 females) controls were included. The linear mandibular and the angular measurements were statistically different between MR and controls (p < 0.05). In the genotype and allele distribution analysis, the SNPs rs694, rs307247, and rs464653 were associated with MR (p < 0.05). MDR analyses predicted the best interaction model for MR was rs694–rs927650, followed by rs307247–rs464653–rs927650. Some haplotypes in the PTH gene presented statistical significance. Our results suggest that SNPs in PTH, VDR, CYP24A1, and CYP27B1 genes are associated with the presence of mandibular retrognathism.

Graphical Abstract

1. Introduction

Skeletal malocclusions are a set of human craniofacial morphologic characteristics that result in an improper skeletal relationship of the jaws and specific facial patterns. There is considerable evidence from family and twin studies suggesting that genes play a significant role in the etiology of skeletal malocclusions [1]. Familial occurrence of skeletal Class II malocclusion has been reported in some studies including twin and triplet studies design and in family pedigrees studies design [1]. A clinical and cephalometric study evaluating 114 patients (48 twin pairs and six sets of triplets) performed an intra- and inter-pair comparison to determine concordance/discordance rates for monozygotic and dizygotic twins. The monozygotic twin pairs presented 100% concordance, while almost 90% of the dizygotic twin pairs were discordant for the Class II division 2 malocclusion [2]. A study evaluating skeletal and dental patterns demonstrated that hereditary factors were responsible for 40% of the variations resulting in malocclusion and the genetic component was higher for skeletal than for dental patterns [3]. In a classical study, cephalometric radiographs of several families were analyzed to evaluate facial morphologic differences among family members. The authors observed a high correlation between parents and their offspring and concluded that a strong familial tendency for skeletal malocclusions exists, supporting the impact of genetic background on face morphology [4]. More recently, many studies in different populations have been performed and demonstrated that genes involved in a variety of functions are associated with skeletal malocclusions and morphological patterns of the face [5,6]. Mandibular retrognathism (MR) is a common skeletal malocclusion in humans. It refers to a retruded position of the mandible as a consequence of an anomaly of the skeletal jaw–cranial base relationship. This unfavorable positional relationship of developing jaws is diagnosed through image exams [7].
Single nucleotide polymorphisms (SNPs) are a type of genetic variant involving variation of a single base pair in the genome. This is the most common type of genetic variation in humans, and it has been stated that these variants could explain differences in individual predisposition to present complex traits [8]. Previous studies have shown that MR was associated with SNPs in various genes such as Myosin 1H (MYO1H) [9], Matrilin 1 (MATN1) [7], ADAM Metallopeptidase with Thrombospondin Type 1 Motif 9 (ADAMTS9) [10], and Bone Morphogenic Protein 2 (BMP2) [11] genes. In addition, SNPs in the gene encoding growth hormone receptor (GHR) have been associated with variations in the mandibular pattern including prognathism [5,12], supporting that genes coding for hormones, hormone receptors, hormone precursors, and molecules involved in hormonal synthesis could also be involved in the etiology of other phenotypes.
Two important hormones that play a crucial role in bone development are parathyroid hormone (PTH) and vitamin D (a secosteroid hormone). These are major regulators of mineral metabolism involved in calcium and phosphate homeostasis as well as in bone growth and development [13]. The biological actions of vitamin D are exerted by binding to the nuclear vitamin D receptor (VDR) [14]. VDR is expressed in the parathyroid glands acting as sensors for the detection and maintenance of adequate vitamin D levels, regulating PTH synthesis and release [15], which among other tissues affects the periodontal ligament [16]. Additionally, other molecules such as the vitamin D 24-hydroxylase (CYP24A1) and 1-hydroxylase (CYP27B1) enzymes participate in related processes, being considered as pivotal determinants of the local concentration of active vitamin D [17].
Nutrition and biomechanical factors can affect the facial pattern, however, the most important factor related with craniofacial growth and development seem to be hormones, genetic, and molecular mechanisms as well as the interplay among them [18]. There is some evidence from animal model studies that PTH is involved in mandible and mandibular condyle development and mineralization [19,20]. Additionally, a study in rats using cephalometric analysis to investigate the effects of a low calcium and vitamin D-deficient diet on craniofacial morphology and growth observed that vitamin D-free and low calcium diet affects mandible development [21]. However, the role of PTH, vitamin D, and genes associated with these two hormones in the human MR is still unclear. Therefore, in the present study, we evaluated if SNPs in genes encoding PTH, VDR, CYP24A1, and CYP27B1 as well as the interplay among them were associated with MR. A tightly controlled connection between vitamin D, serum calcium, and genes involved in the maintenance of vitamin D levels orchestrates mineral homeostasis and development.

2. Materials and Methods

2.1. Sample

This nested case-control study was previously approved by the Human Ethics Committee of the University of São Paulo—Ribeirão Preto Dental School (# 01451418.3.0000.5419). Informed consent/assent was obtained from all participants and/or their legal guardians where age-appropriate. This project was performed according to the Helsinki Declaration and its amendments. The Strengthening the Reporting of Genetic Association study (STREGA) statement checklist [22] was followed to develop and report the results of this study.
Patients undergoing orthodontic treatment at the University of São Paulo were recruited and consecutively included in this study from 2015 to 2017. Patients with syndromes, congenital alterations, hormonal, and/or metabolic disorders or those with previous orthodontic and/or orthopedic treatments were not included. None of the patients reported history of vitamin D deficiency. All patients included were Brazilian biologically-unrelated and self-reported as Caucasian. Additionally, none of the patients were using vitamin D supplementation and all patients were from the Ribeirão Preto area, which presents a high total sunshine hour and UVB incidence over the entire year.

2.2. Phenotypes Definition

Pre-orthodontic lateral cephalograms with the mandible in centric relationship were used and digital cephalometric tracings performed by a calibrated orthodontist using the software Dolphin Imaging Version 8.0 (Dolphin Imaging, Chatsworth, CA, USA), as demonstrated in Figure 1. The following landmarks were used to determine the phenotype: point A, point B, sella (S), and nasion (N) and, therefore, the angular measurements SNB and ANB were calculated. Patients having a retrognathic mandible (SNB < 78°) were selected as cases, and those with an orthognathic mandible (SNB = 78°–82°) were selected as controls. Patients with mandibular prognathism (SNB > 82°) were excluded.
Additionally, the linear measurements associated with mandibular size (mandibular length, Co-Gn; length of mandibular base, Go-Pg; and mandibular ramus height, Co-Go) were measured in millimeters (mm) and compared between the patients with retrognathic mandible and the patients with orthognathic mandible.

2.3. Allelic Discrimination

Genomic DNA extracted from saliva was used for genotyping analysis. Briefly, for saliva collection, saline mouth solution to rinse in the mouth for 60 s was used. Therefore, the genomic DNA was extracted from buccal epithelial cells from saliva samples as previously described [23]. Quantification of the concentration and purity of the DNA was determined by spectrophotometry (Nanodrop 1000; Thermo Scientific, Wilmington, DE, USA).
Six SNPs were evaluated in the present study and are reported in Table 1. The genotyping was blindly performed using the Taqman™ method for real-time PCR in the StepOnePlusTM sequence detection system, Applied Biosystems™ (Foster City, CA, USA) or in the Mastercycler® ep realplex-S thermocycler, Eppendorf AG (Hamburg, Germany). Additionally, 10% of the sample was genotyped twice and an agreement of 100% was observed. The reaction was previously described in [24].

2.4. Statistical Analysis

The Chi-squared test was used to estimate the Hardy–Weinberg equilibrium (https://wpcalc.com/en/equilibrium-hardy-weinberg/) (accessed on 2 April 2021).
Chi-squared or Fisher’s exact tests were used to compare gender, genotypic, and allelic distribution among groups. The Mann–Whitney U test was applied for the comparison of continuous data among the groups after evaluating the normality by the Shapiro–Wilk test, and Levene’s test was used to assess the homogeneity of variance. These analyses were performed by GraphPad Prism version 7.0 for Windows (GraphPad Software, San Diego, CA, USA).
Haplotype analysis was performed for the SNPs in PTH by PLINK version 1.06 (https://zzz.bwh.harvard.edu/plink/ld.shtml) (accessed on 2 April 2021). The established alpha for these analyses was p < 0.05.
Multifactor dimensionality reduction (MDR) [31] was done to identify SNP-SNP interactions using gender and age as co-variables as previously reported [11]. Two specific software programs (Multifactor Dimensionality Reduction 3.0.2 [31], and MDR Permutation Testing Module 1.0 beta 2 [32], available in sourceforge.net/projects/mdr/files (accessed on 2 April 2021)), were used to perform a 10-fold cross-validation consistency (CVC), testing balancing accuracy (TBA), and the 1000 permutation test to determine the statistical significance of the models. Models with the cross-validation consistency of 9/10 or 10/10 and the TBA > 0.55 and p ≤ 0.05 were considered as the best models. Entropy values were calculated according to Jakulin and Bratko [33], and MDR created dendrograms and interaction graphs using these values. MDR analysis can be performed freely in a JAVA language.

3. Results

A total of 91 patients were included. Figure 2 presents the flow diagram of the study participants and genotype success rate for each SNP. The sample characteristics and mandibular parameters according to the groups (MR and control) are shown in Table 2. There was no significant difference between the groups according to age and gender distribution (p > 0.05). The participants in the MR group had significantly lower linear measurements in the mandible than the control group (p < 0.05). Furthermore, the SNB angle was significantly lower and the ANB angle was significantly higher in the retrognathic subjects (p < 0.001).
All SNPs assessed were within the Hardy–Weinberg equilibrium. Genotype and allele distributions are demonstrated in Table 3. The SNPs rs694 and rs307247 in PTH and the SNP rs464653 in CYP27B1 were significantly associated with mandibular retrognathism (p < 0.05).
The best models for MDR analyses were rs694 (PTH), rs927650 (CYP24A1), rs307247 (PTH), rs464653 (CYP27B1), rs927650 (CYP24A1) and rs694 (PTH), rs6256 (PTH), rs307247 (PTH), rs7975232 (VDR), rs464653 (CYP27B1), and rs927650 (CYP24A1) with a cross-validation consistency of 10 out of 10. Table 4 shows the MDR-predicted interaction models.
Figure 3 shows the interactions between SNPs (dendrogram and interaction map). The strongest synergism interaction effect was between rs307247 (PTH) and rs927650 (CYP24A1), followed by the interaction between rs6256 (PTH) and rs464653 (CYP27B1).
The haplotype analysis for the SNPs in the PTH gene is presented in Table 5. Some haplotypes were associated with mandibular retrognathism.

4. Discussion

Disorders of the face and dental jaws such as skeletal malocclusions are very common developmental disorders in all ethnic populations [34]. Skeletal malocclusion affects dental and facial tissues [34]. Mandibular retrognathism is a condition that not only affects facial aesthetics, but is also associated with problems such as temporomandibular disorders [35] and alterations in the respiratory pattern and normal sleep [34,36]. The prevalence of skeletal malocclusion ranges in different ethnic groups [34]. In the orthodontic population studied here with evaluated Brazilian patients, skeletal class II malocclusion affects about 30% of the sample [11,37] and mandibular retrognathism affects 35% of the sample [8]. Retrognathism is a much more frequent condition than mandibular prognathism; however, its etiology has been the subject of only a few studies. Most studies on the genetic background of skeletal malocclusion have focused on prognathism and Class III malocclusion, which is less frequent in the general population. Studies on genes involved in skeletal Class II malocclusion, retrognathism, and micrognathia are rarer [2] and only a few genes and SNPs have been evaluated thus far [7,9,10,11,38,39]. Therefore, the present study aimed to explore the role of some SNPs as well as their interaction in the mandibular retrognathism phenotype in humans.
Sequence variation in human genes is largely confined to SNPs and is valuable in tests of association with common traits such as retrognathism. Uncovering the SNPs and genes responsible for the regulation of facial morphology is not a trivial task. Human facial development is a complex multistep process, implicating several signaling cascades of factors. The mechanisms involved in this process include the expression of innumerous genes and protein translation. These events are precisely timed and are under hormonal control [40]. Therefore, in the present study, we evaluated whether common SNPs in genes involved in hormonal synthesis and metabolisms were associated with mandibular retrognathism.
PTH is an 84-amino acid peptide hormone synthesized in the cells of the parathyroid glands. This hormone is a major mediator of bone remodeling and plays a crucial role in calcium homeostasis, showing several effects on the bone remodeling process, resulting in anabolic activity (bone formation) and catabolic activity (bone resorption) [41,42]. PTH promotes calcium release at the bone level, in which a hypocalcemic signal will lead to a higher release and synthesis of PTH, restoring the serum calcium to normal [43]. In our study, two SNPs in the gene encoding PTH were associated with retrognathism in the univariate analysis: the intronic SNP rs694 and the SNP rs307247, which is located in a 3’ untranslated region (3’UTR). One important aspect in genetic association studies is that the majority of traits were associated with non-coding regions (intronic and intergenic) called regulatory SNPs [44]. Intronic variants can impact alternative splicing by interfering with splice site recognition [45]. Additionally, 3’UTRs can modify gene expression by controlling the mRNA nuclear export, cytoplasmic localization and stability, or by affecting translational efficiency. These gene fragments are targeted by microRNA as well as regulatory molecules [46]. However, one aspect to be highlighted for future studies designs is that a redundancy was observed between rs694 and rs307247 in the interaction map.
In our study, we also evaluated the SNP rs6256 located in exon 3 of the PTH gene. This nonsynonymous variant might contribute to the altered gene expression. A previous study demonstrated that serum PTH levels were higher in individuals carrying the rs6256 AA genotype [26]. We observed that this SNP was associated with MR in the MDR analysis and in the haplotype analysis.
It is well known that vitamin D requires the involvement of several key proteins in a closely regulated process. The vitamin D produced in the skin, in response to UVB light or from exogenous supplements, is sequentially hydroxylated into the active metabolite 1,25-Dihydroxyvitamin D (1,25(OH)2D, the biologically active form of vitamin D) in the kidney and other tissues via the enzyme CYP27B1 [47]. Additionally, 1,25(OH)2D concentrations are regulated by CYP24A1 [48]. In the genotype distribution as well as in the MDR analysis, the intronic SNP rs464653 in CYP27B1 was associated with mandibular retrognathism. A strong interaction and synergism were observed between the SNP rs307247 in PTH and rs927650 in CYP24A1. Parathyroid glands contain CYP27B1 and CYP24A1 and the expression levels of both these enzymes are transcriptionally regulated by PTH [49]. Therefore, our results suggest that SNPs in genes involved in this highly regulated interaction are involved in mandibular morphology in humans.
A study evaluating animals fed with a low calcium and vitamin D-deficient diet observed that this diet caused alterations in craniofacial morphology including reduced mandibular dimensions [21]. The biological effects of vitamin D are mediated by binding to its receptor, a member of the nuclear receptor superfamily, encoded by the gene VDR [14]. Although the studied SNP rs7975232 in VDR seems to only have a small role in the mandibular retrognathism phenotype, this result should be interpreted with caution. In our analysis, only one SNP was evaluated in VDR, which is a gene with some well-known SNPs. Therefore, the coverage of the gene was a limitation of our study. It is possible that other SNPs in VDR as well as their interaction are also involved in the retrognathism phenotype. Additionally, the generalizability of our results to populations of other ancestries is unknown. Further studies are necessary to evaluate the role of these SNPs in skeletal malocclusions in different populations.

5. Conclusions

Briefly, MR is a polygenic trait. SNPs in coding and regulatory regions can lead to different gene activities and possible contributions to interindividual variability including variability in facial morphology. Our results support that genes involved in the maintenance of vitamin D levels are involved in the etiology of human MR.

Author Contributions

Conceptualization, E.C.K., P.N.-F., M.A.N.M., P.P. and C.K.; Resources, E.C.K., P.P., P.N-.F. and C.K.; Supervision, M.A.N.M. and M.B.S.S.; Project, E.C.K. and C.K.; Sample collection, G.M.-V., M.A.N.M. and M.B.S.S.; Phenotype analysis, G.M.-V., M.A.N.M. and M.B.S.S.; Statistical analysis, M.A.H.d.M.O. and C.L.B.R.; Data interpretation, E.C.K., C.L.B.R., G.M.-V. and C.K.; Funding acquisition E.C.K. and C.K.; Writing—original draft preparation, E.C.K., C.L.B.R., G.M.-V. and C.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brasil (CAPES)-Finance Code 001 and Alexander-von-Humboldt-Foundation (Küchler/Kirschneck accepted on July 4th, 2019). This work was also supported by the São Paulo Research Foundation (FAPESP) (2015/06866-5) and individual scholarships (CAPES and CNPq).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Human Ethics Committee of the University of São Paulo—Ribeirão Preto Dental School (# 01451418.3.0000.5419).

Informed Consent Statement

Informed consent/assent was obtained from all participants and/or their legal guardians where age-appropriate.

Data Availability Statement

The data generated in the present study are included within the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mossey, P.A. The heritability of malocclusion: Part 2. The influence of genetics in malocclusion. Br. J. Orthod. 1999, 26, 195–203. [Google Scholar] [CrossRef]
  2. Markovic, M.D. At the crossroads of oral facial genetics. Eur. J. Orthod. 1992, 14, 469–481. [Google Scholar] [CrossRef] [PubMed]
  3. Lundström, A. Nature versus nurture in dento-facial variation. Eur. J. Orthod. 1994, 6, 77–91. [Google Scholar] [CrossRef] [PubMed]
  4. Nakasima, A.; Ichinose, M.; Nakata, S.; Takahama, Y. Hereditary factors in the craniofacial morphology of Angle’s Class II and Class III malocclusions. Am. J. Orthod. 1982, 82, 1–26. [Google Scholar] [CrossRef]
  5. Uribe, L.M.M.; Miller, S.F. Genetics of the dentofacial variation in human malocclusion. Orthod. Craniofac. Res. 2015, 18, 91–99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Dehesa-Santos, A.; Iber-Diaz, P.; Iglesias-Linares, A. Genetic factors contributing to skeletal class III malocclusion: A systematic review and meta-analysis. Clin. Oral Investig. 2021, 25, 1587–1612. [Google Scholar] [CrossRef] [PubMed]
  7. Balkhande, P.B.; Lakkakula, B.V.K.S.; Chitharanjan, A.B. Relationship between matrilin-1 gene polymorphisms and mandibular retrognathism. Am. J. Orthod. Dentofac. Orthop. 2018, 153, 255–261. [Google Scholar] [CrossRef]
  8. Teama, S. DNA polymorphisms: DNA-based molecular markers and their application in medicine. In Genetic Diversity and Disease Susceptibility, 1st ed.; Liu, Y., Ed.; IntechOpen: London, UK, 2018; Volume 25. [Google Scholar]
  9. Arun, R.M.; Lakkakula, B.V.; Chitharanjan, A.B. Role of myosin 1H gene p morphisms in mandibular retrognathism. Am. J. Orthod. Dentofac. Orthop. 2016, 149, 699–704. [Google Scholar] [CrossRef]
  10. Wang, C.; Ni, Z.; Cai, Y.; Zhou, Y.; Chen, W. Association of Polymorphism rs67920064 in ADAMTS9 Gene with Mandibular Retrognathism in a Chinese Population. Med. Sci. Monit. 2020, 30, e925965. [Google Scholar]
  11. Küchler, E.C.; Reis, C.L.B.; Carelli, J.; Scariot, R.; Nelson-Filho, P.; Coletta, R.D.; Paza, A.O.; Matsumoto, M.A.N.; Proff, P.; Kirschneck, C. Potential interactions among single nucleotide polymorphisms in bone- and cartilage-related genes in skeletal malocclusions. Orthod. Craniofac. Res. 2021, 24, 277–287. [Google Scholar] [CrossRef]
  12. Nascimento, M.A.; Oliveira, D.S.B.; Reis, C.L.B.; Wambier, L.M.; Horta, K.C.; Romano, F.L.; Bezerra da Silva, L.A.; Raquel da Silva, A.B.; Nelson-Filho, P.; Kuchler, E.C. Association between P561T polymorphism in growth hormone receptor gene and mandibular prognathism: Systematic review and metaanalysis. Rev. Científica Do CRO-RJ 2019, 4, 2–11. [Google Scholar]
  13. Khundmiri, S.J.; Murray, R.; Lederer, E. PTH and Vitamin D. Compr. Physiol. 2016, 15, 561–601. [Google Scholar]
  14. Van der Velden, U.; Kuzmanova, D.; Chapple, I.L.C. Micronutritional approaches to periodontal therapy. J. Clin. Periodontol. 2011, 38, 142–158. [Google Scholar] [CrossRef]
  15. Beckerman, P.; Silver, J. Vitamin D and the parathyroid. Am. J. Med. Sci. 1999, 317, 363–369. [Google Scholar]
  16. Wolf, M.; Lossdörfer, S.; Marciniak, J.; Römer, P.; Kirschneck, C.; Craveiro, R.; Deschner, J.; Jäger, A. CD8+ T cells mediate the regenerative PTH effect in hPDL cells via Wnt10b signaling. Innate Immun. 2016, 22, 674–681. [Google Scholar] [CrossRef] [Green Version]
  17. Zhalehjoo, N.; Shakiba, Y.; Panjehpour, M. Gene expression profiles of CYP24A1 and CYP27B1 in malignant and normal breast tissues. Mol. Med. Rep. 2017, 15, 467–473. [Google Scholar] [CrossRef] [Green Version]
  18. Williams, S.E.; Slice, D.E. Regional shape change in adult facial bone curvature with age. Am. J. Phys. Anthropol. 2010, 143, 437–447. [Google Scholar] [CrossRef] [PubMed]
  19. Tsutsui, T.W.; Riminucci, M.; Holmbeck, K.; Bianco, P.; Robey, P.G. Development of craniofacial structures in transgenic mice with constitutively active PTH/PTHrP receptor. Bone 2008, 42, 321–331. [Google Scholar] [CrossRef] [Green Version]
  20. Dutra, E.H.; O’Brien, M.H.; Gutierrez, T.; Lima, A.; Nanda, R.; Yadav, S. PTH [1-34]-induced alterations predispose the mandibular condylar cartilage to mineralization. Orthod. Craniofacial Res. 2017, 20, 162–166. [Google Scholar] [CrossRef]
  21. Engström, C.; Linde, A.; Thilander, B. Craniofacial morphology and growth in the rat. Cephalometric analysis of the effects of a low calcium and vitamin D-deficient diet. J. Anat. 1982, 134 Pt 2, 299–314. [Google Scholar]
  22. Little, J.; Higgins, J.P.; Ioannidis, J.P.; Moher, D.; Gagnon, F.; Von Elm, E.; Khoury, M.J.; Cohen, B.; Davey-Smith, G.; Grimshaw, J.; et al. STrengthening the REporting of Genetic Association Studies (STREGA)—an extension of the STROBE statement. Genet. Epidemiol. Off. Publ. Int. Genet. Epidemiol. Soc. 2009, 33, 581–598. [Google Scholar]
  23. Küchler, E.C.; Tannure, P.N.; Falagan-Lotsch, P.; Lopes, T.S.; Granjeiro, J.M.; Amorim, L.M.F. Buccal cells DNA extraction to obtain high quality human genomic DNA suitable for polymorphism genotyping by PCR-RFLP and Real-Time PCR. J. Appl. Oral Sci. 2012, 20, 467–471. [Google Scholar] [CrossRef] [Green Version]
  24. Cunha, A.S.; Dos Santos, L.V.; Marañón-Vásquez, G.A.; Kirschneck, C.; Gerber, J.T.; Stuani, M.B.; Matsumoto, M.A.N.; Vieira, A.R.; Scariot, R.; Küchler, E.C. Genetic variants in tooth agenesis–related genes might be also involved in tooth size variations. Clin. Oral Investig. 2021, 25, 1307–1318. [Google Scholar] [CrossRef]
  25. Ertl, D.A.; Stary, S.; Streubel, B.; Raimann, A.; Haeusler, G. A novel homozygous mutation in the parathyroid hormone gene (PTH) in a girl with isolated hypoparathyroidism. Bone 2012, 51, 629–632. [Google Scholar] [CrossRef] [PubMed]
  26. Kanzawa, M.; Sugimoto, T.; Kobayashi, T.; Kobayashi, A.; Chihara, K. Parathyroid hormone gene polymorphisms in primary hyperparathyroidism. Clin. Endocrinol. 1999, 50, 583–588. [Google Scholar] [CrossRef] [PubMed]
  27. Tenne, M.; McGuigan, F.; Jansson, L.; Gerdhem, P.; Obrant, K.J.; Luthman, H.; Åkesson, K. Genetic variation in the PTH pathway and bone phenotypes in elderly women: Evaluation of PTH, PTHLH, PTHR1 and PTHR2 genes. Bone 2008, 42, 719–727. [Google Scholar] [CrossRef]
  28. Al-Ghafari, A.B.; Balamash, K.S.; Al Doghaither, H.A. Relationship Between Serum Vitamin D and Calcium Levels and Vitamin D Receptor Gene Polymorphisms in Colorectal Cancer. BioMed Res. Int. 2019, 2019, 8571541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Zeljic, K.; Supic, G.; Stamenkovic, R.M.; Jovic, N.; Kozomara, R.; Magic, Z. Vitamin D receptor, CYP27B1 and CYP24A1 genes polymorphisms association with oral cancer risk and survival. J. Oral Pathol. Med. 2012, 41, 779–787. [Google Scholar] [CrossRef]
  30. Hibler, E.A.; Klimentidis, Y.C.; Jurutka, P.W.; Kohler, L.N.; Lance, P.; Roe, D.J.; Jacobs, E.T. CYP24A1 and CYP27B1 polymorphisms, concentrations of vitamin D metabolites, and odds of colorectal adenoma recurrence. Nutr. Cancer 2015, 67, 1131–1141. [Google Scholar] [CrossRef] [Green Version]
  31. Ritchie, M.D.; Hahn, L.W.; Roodi, N.; Bailey, L.R.; Dupont, W.D.; Parl, F.F.; Moore, J.H. Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am. J. Human Genet. 2001, 69, 138–147. [Google Scholar] [CrossRef] [Green Version]
  32. Pattin, K.A.; White, B.C.; Barney, N.; Gui, J.; Nelson, H.H.; Kelsey, K.T.; Moore, J.H. A computationally efficient hypothesis testing method for epistasis analysis using multifactor dimensionality reduction. Genet. Epidemiol. Off. Publ. Int. Genet. Epidemiol. Soc. 2009, 33, 87–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jakulin, A.; Bratko, I. Quantifying and visualizing attribute interactions. arXiv 2004, arXiv:cs/0308002v3. [Google Scholar]
  34. Joshi, N.; Hamdan, A.M.; Fakhouri, W.D. Skeletal malocclusion: A developmental disorder with a life-long morbidity. J. Clin. Med. Res. 2014, 6, 399–408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Miller, J.R.; Burgess, J.A.; Critchlow, C.W. Association between mandibular retrognathia and TMJ disorders in adult females. J. Public Health Dent. 2004, 64, 157–163. [Google Scholar] [CrossRef] [PubMed]
  36. Solow, B.; Siersbæk-Nielsen, S.; Greve, E. Airway adequacy, head posture, and craniofacial morphology. Am. J. Orthod. 1984, 86, 214–223. [Google Scholar] [CrossRef]
  37. Costa, A.M.G.; Trevizan, M.; Matsumoto, M.A.N.; da Silva, R.A.B.; da Silva, L.A.B.; Horta, K.C.; Romano, F.L.; Nelson-Filho, P.; Küchler, E.C. Association between tooth agenesis and skeletal malocclusions. J. Oral Maxillofac. Res. 2017, 8, e3. [Google Scholar] [CrossRef] [Green Version]
  38. Gutierrez, S.J.; Gomez, M.; Rey, J.A.; Ochoa, M.; Gutierrez, S.M.; Prieto, J.C. Polymorphisms of the noggin gene and mandibular micrognathia: A first approximation. Acta Odontol Lat. 2010, 23, 13–19. [Google Scholar]
  39. da Fontoura, C.G.; Miller, S.F.; Wehby, G.L.; Amendt, B.A.; Holton, N.E.; Southard, T.E.; Moreno Uribe, L.M. Candidate gene analyses of skeletal variation in malocclusion. J. Dent. Res. 2015, 94, 913–920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Barash, M.; Bayer, P.E.; van Daal, A. Identification of the Single Nucleotide Polymorphisms Affecting Normal Phenotypic Variability in Human Craniofacial Morphology Using Candidate Gene Approach. J. Genet. Genome Res. 2018, 5, 041. [Google Scholar]
  41. Barros, S.P.; Silva, M.A.D.; Somerman, M.J.; Nociti, F.H., Jr. Parathyroid hormone protects against periodontitis-associated bone loss. J. Dent. Res. 2013, 82, 791–795. [Google Scholar] [CrossRef]
  42. Tokunaga, K.; Seto, H.; Ohba, H.; Mihara, C.; Hama, H.; Horibe, M.; Yoneda, S.; Nagata, T. Topical and intermittent application of parathyroid hormone recovers alveolar bone loss in rat experimental periodontitis. J. Periodontal Res. 2011, 46, 655–662. [Google Scholar] [CrossRef]
  43. Silva, B.C.; Costa, A.G.; Cusano, N.E.; Kousteni, S.; Bilexikian, J.P. Catabolic and anabolic actions of parathyroid hormone on the skeleton. J. Endocrinol. Investig. 2011, 34, 801–810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hindorff, L.A.; Sethupathy, P.; Junkins, H.A.; Ramos, E.M.; Mehta, J.P.; Collins, F.S.; Manolio, T.A. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc. Natl. Acad. Sci. USA 2009, 106, 9362–9367. [Google Scholar] [CrossRef] [Green Version]
  45. Lin, H.; Hargreaves, K.A.; Li, R.; Reiter, J.L.; Wang, Y.; Mort, M.; Cooper, D.N.; Zhou, Y.; Zhang, C.; Eadon, M.T.; et al. RegSNPs-intron: A computational framework for predicting pathogenic impact of intronic single nucleotide variants. Genome Biol. 2019, 20, 1–16. [Google Scholar] [CrossRef]
  46. Hughes, T.A. Regulation of gene expression by alternative untranslated regions. Trends Genet. 2006, 22, 119–122. [Google Scholar] [CrossRef]
  47. Lehmann, B.; Meurer, M. Vitamin D metabolism. Dermatol. Ther. 2010, 23, 2–12. [Google Scholar] [CrossRef] [PubMed]
  48. Tashiro, K.; Abe, T.; Oue, N.; Yasui, W.; Ryoji, M. Characterization of vitamin D-mediated induction of the CYP 24 transcription. Mol. Cell. Endocrinol. 2004, 226, 27–32. [Google Scholar] [CrossRef] [PubMed]
  49. Segersten, U.; Björklund, P.; Hellman, P.; Åkerström, G.; Westin, G. Potentiating effects of nonactive/active vitamin D analogues and ketoconazole in parathyroid cells. Clin. Endocrinol. 2007, 66, 399–404. [Google Scholar] [CrossRef] [PubMed]
Jpm 11 00369 g001 550
Figure 1. Cephalometric tracings. (A) Exemplification of retrognathism (SNB > 78°). (B) Exemplification of orthognathic mandible (SNB ranging from 78° to 82°).
Figure 1. Cephalometric tracings. (A) Exemplification of retrognathism (SNB > 78°). (B) Exemplification of orthognathic mandible (SNB ranging from 78° to 82°).
Jpm 11 00369 g001
Jpm 11 00369 g002 550
Figure 2. Flow chart of patient selection.
Figure 2. Flow chart of patient selection.
Jpm 11 00369 g002
Jpm 11 00369 g003 550
Figure 3. Interaction analysis for the studied SNPs by MDR. (A) Interaction dendrogram. Dendrogram shows keys connecting the SNPs. Strongly interacting SNPs are connected by keys farther to the right, whereas the weakly interacting SNPs are connected by keys farther to the left. (B) Interaction entropy graph. Each SNP are represented by a node and percentages inside the nodes are their entropy values. The percentages between the nodes indicate the entropy values resulting from the combination between SNPs.
Figure 3. Interaction analysis for the studied SNPs by MDR. (A) Interaction dendrogram. Dendrogram shows keys connecting the SNPs. Strongly interacting SNPs are connected by keys farther to the right, whereas the weakly interacting SNPs are connected by keys farther to the left. (B) Interaction entropy graph. Each SNP are represented by a node and percentages inside the nodes are their entropy values. The percentages between the nodes indicate the entropy values resulting from the combination between SNPs.
Jpm 11 00369 g003
Table
Table 1. Characteristics of the studied SNPs.
Table 1. Characteristics of the studied SNPs.
GeneSNPBase
Change #		Functional
Consequence #		Clinical
Significance †,‡		Biological Effects (Reference)
PTHrs694C > TIntron VariantBenignLow levels of PTH serum [25]
rs6256 *G > TStop Gained *BenignLow levels of PTH serum [26]
rs307247G > A3′ untranslated regionBenignLow levels of PTH serum [27]
VDRrs7975232 #C > AIntron VariantBenignHigh levels of Vitamin-D serum [28]
CYP27B1rs464653A > GIntron VariantBenignOral Neoplasm [29]
CYP24A1rs927650C > TIntron VariantUncertain SignificanceHigh levels of Vitamin-D serum [30]
Note: * (Arg) > (Ter). # known as Apal. Information was obtained in ncbi.nlm.nih.gov/snp/ (#) (accessed on 22 December 2020), ncbi.nlm.nih.gov/CBBresearch/Lu/Demo/LitVar/# (†) (accessed on 22 December 2020) and varsome.com (‡) (accessed on 22 December 2020).
Table
Table 2. Comparison of cephalometric variables between mandibular retrognathism and control groups.
Table 2. Comparison of cephalometric variables between mandibular retrognathism and control groups.
VariablesMRControlp-Value
Gender n (%)
Male22 (45.8)17 (39.5)0.672
Female26 (54.2)26 (60.1)
Age
Median (IQR)12.0 (4.0)12.0 (4.5)0.809
SNB (°)
Median (IQR)76.1 (3.5)80.0 (2.0)<0.001 *
ANB (°)
Mean (SD)4.2 (2.3)2.4 (2.3)<0.001 *
Co-Gn (mm)
Mean (SD)110 (8.0)115 (9.4)0.008 *
Go-Pg (mm)
Median (IQR)64.4 (6.9)68.2 (5.5)0.014 *
Co-Go (mm)
Median (IQR)54.0 (8.8)56.8 (7.4)0.044 *
Note: IQR means interquartile range; SD means standard deviation. MR means mandibular retrognathism. * means statistically significant difference (p < 0.05).
Table
Table 3. Genotype and allele distribution among the groups.
Table 3. Genotype and allele distribution among the groups.
GeneSNP Frequency-n (%)p-ValueOR (CI 95%)
Genotype/AlleleControlMR
PTHrs694GenotypeTT4 (10.8)19 (44.2)Reference-
CT22 (59.5)19 (44.2)0.004 *0.18 (0.06–0.59)
CC11 (29.7)5 (11.6)0.001 *0.09 (0.02–0.44)
AlleleT30 (40.5)57 (59.4)Reference-
C44 (59.5)39 (40.6)0.014 *0.46 (0.25–0.87)
rs6256GenotypeGG30 (81.1)32 (76.2)Reference
GT7 (18.9)9 (21.4)0.7401.20 (0.39–3.39)
TT0 (0.0)1 (2.4)>0.999-
AlleleG67 (90.5)73 (86.9)Reference-
T7 (9.5)11 (13.1)0.4721.44 (0.51–3.71)
rs307247GenotypeGG12 (33.3)26 (60.5)Reference
AG14 (38.9)12 (27.9)0.0740.39 (0.14–1.16)
AA10 (27.8)5 (11.6)0.019 *0.23 (0.06–0.78)
AlleleG38 (52.8)64 (75.0)Reference
A34 (47.2)22 (25.0)0.003 *0.37 (0.19–0.73)
VDRrs7975232GenotypeAA14 (35.0)18 (41.9)Reference
AC17 (42.5)21 (48.8)0.9340.96 (0.37–2.44)
CC9 (22.5)4 (9.3)0.1210.34 (0.10–1.31)
AlleleA45 (56.2)57 (66.3)Reference-
C35 (43.8)29 (33.7)0.1840.65 (0.34–1.20)
CYP27B1rs464653GenotypeAA13 (33.3)25 (55.6)Reference-
AG21 (53.9)15 (33.3)0.037 *0.37 (0.13–0.93)
GG5 (12.8)5 (11.1)0.3580.52 (0.12–2.09)
AlleleA47 (60.3)45 (64.3)Reference-
G31 (39.7)25 (35.7)0.6130.84 (0.44–1.68)
CYP24A1rs927650GenotypeCC17 (44.7)15 (34.1)Reference-
CT16 (42.1)26 (59.1)0.1971.84 (0.75–4.69)
TT5 (13.2)3 (6.8)0.6330.68 (0.16–3.45)
AlleleC50 (65.8)56 (63.6)Reference-
T26 (34.2)32 (36.4)0.7731.09 (0.57–2.12)
Note: * means statistically significant difference (p < 0.05).
Table
Table 4. Summary of MDR analysis results.
Table 4. Summary of MDR analysis results.
Locus NumberBest CombinationCVC #TBA †p-Value ‡
2rs694 (PTH), rs927650 (CYP24A1)10/100.67420.040 *
3rs307247 (PTH), rs464653 (CYP27B1), rs927650 (CYP24A1)10/100.7651<0.001 *
4rs307247 (PTH), rs7975232 (VDR), rs464653 (CYP27B1), rs927650 (CYP24A1)8/100.70160.009 *
5rs694 (PTH), rs307247 (PTH), rs7975232 (VDR), rs464653 (CYP27B1), rs927650 (CYP24A1)8/100.68320.026 *
6rs694 (PTH), rs6256 (PTH), rs307247 (PTH), rs7975232 (VDR), rs464653 (CYP27B1), rs927650 (CYP24A1)10/100.70850.008 *
Note: * means statistically significant difference (p < 0.05). # The 10-fold CVC indicates a prediction error. † TBA values show the proportion of individuals correctly classified as the case or control. ‡ p-values were based on the 1000 permutations test by a specifically developed software for MDR analysis. More details can be consulted in Pattin et al. (2009) [32].
Table
Table 5. Haplotype association analysis of SNPs in the PTH gene.
Table 5. Haplotype association analysis of SNPs in the PTH gene.
SNPSHaplotypeMRControlp-Value
rs694, rs6256, rs307247C-G-A0.250.470.004 *
C-T-G0.010.030.510
T-T-G0.110.060.231
C-G-G0.060.090.530
T-G-G0.540.340.010 *
rs694, rs6256T-T0.130.090.472
C-G0.330.590.001 *
T-G0.530.310.004 *
rs694, rs307247C-A0.250.470.004 *
C-G0.080.120.365
T-G0.660.400.001 *
rs307247, rs6256A-G0.260.470.006 *
G-T0.130.090.510
G-G0.600.430.027 *
Note: * means statistically significant difference.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Küchler, E.C.; Reis, C.L.B.; Marañón-Vásquez, G.; Nelson-Filho, P.; Matsumoto, M.A.N.; Stuani, M.B.S.; Oliveira, M.A.H.d.M.; Proff, P.; Kirschneck, C. Parathyroid Hormone Gene and Genes Involved in the Maintenance of Vitamin D Levels Association with Mandibular Retrognathism. J. Pers. Med. 2021, 11, 369. https://doi.org/10.3390/jpm11050369

AMA Style

Küchler EC, Reis CLB, Marañón-Vásquez G, Nelson-Filho P, Matsumoto MAN, Stuani MBS, Oliveira MAHdM, Proff P, Kirschneck C. Parathyroid Hormone Gene and Genes Involved in the Maintenance of Vitamin D Levels Association with Mandibular Retrognathism. Journal of Personalized Medicine. 2021; 11(5):369. https://doi.org/10.3390/jpm11050369

Chicago/Turabian Style

Küchler, Erika Calvano, Caio Luiz Bitencourt Reis, Guido Marañón-Vásquez, Paulo Nelson-Filho, Mírian Aiko Nakane Matsumoto, Maria Bernadete Sasso Stuani, Maria Angélica Hueb de Menezes Oliveira, Peter Proff, and Christian Kirschneck. 2021. "Parathyroid Hormone Gene and Genes Involved in the Maintenance of Vitamin D Levels Association with Mandibular Retrognathism" Journal of Personalized Medicine 11, no. 5: 369. https://doi.org/10.3390/jpm11050369

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop