Next Article in Journal
Pomegranate: A Source of Multifunctional Bioactive Compounds Potentially Beneficial in Alzheimer’s Disease
Next Article in Special Issue
Impact of Pharmacogenomics in Clinical Practice
Previous Article in Journal
Amlexanox: Readthrough Induction and Nonsense-Mediated mRNA Decay Inhibition in a Charcot–Marie–Tooth Model of hiPSCs-Derived Neuronal Cells Harboring a Nonsense Mutation in GDAP1 Gene
Previous Article in Special Issue
CYP2D6 Genotype and Pharmacovigilance Impact on Autism Spectrum Disorder: A Naturalistic Study with Extreme Phenotype Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis

1
Graduate School of Pharmaceutical Sciences, College of Pharmacy, Ewha Womans University, 52 Ewhayeodae-gil, Seodaemun-gu, Seoul 03760, Republic of Korea
2
Department of Oral and Maxillofacial Surgery, School of Medicine, Ewha Womans University Medical Center, Mokdong Hospital, 1071 Anyangcheon-ro, Yangcheon-gu, Seoul 07985, Republic of Korea
3
Institute of Pharmaceutical Science and Technology, College of Pharmacy, Hanyang University, 55 Hanyangdaehak-ro, Sangnok-gu, Ansan 15588, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Pharmaceuticals 2023, 16(7), 1035; https://doi.org/10.3390/ph16071035
Submission received: 30 June 2023 / Revised: 18 July 2023 / Accepted: 20 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Pharmacogenomics - A Genetic Approach to Drug Therapy and Development)

Abstract

:
Bisphosphonate-related osteonecrosis of the jaw (BRONJ) is a rare but severe adverse effect that can occur as a result of bisphosphonate treatment. This study aimed to examine the relationship between PPARγ and PPARGC1A polymorphisms and the BRONJ development in female osteoporosis patients undergoing bisphosphonate treatment. We prospectively conducted this nested case–control study at the Ewha Womans University Mokdong Hospital between 2014 and 2018. We assessed five single-nucleotide polymorphisms (SNPs) of PPARγ and six SNPs of PPARGC1A and performed a multivariable logistic regression analysis to determine the independent risk factors for developing BRONJ. There were a total of 123 patients included in this study and 56 patients (45.5%) developed BRONJ. In the univariate analysis, PPARGC1A rs2946385 and rs10020457 polymorphisms were significantly associated with BRONJ (p = 0.034, p = 0.020, respectively), although the results were not statistically significant in the multivariable analysis. Patients with the combined genotypes of GG in both PPARγ rs1151999 and PPARGC1A rs2946385 showed a 3.03-fold higher risk of BRONJ compared to individuals with other genotype combinations after adjusting for confounders (95% confidence interval (CI): 1.01–9.11). Old age (≥70 years) and duration of bisphosphonate use (≥60 months) increased the risk of BRONJ. The area under the receiver operating characteristic curve for the predicted probability was 0.78 (95% CI: 0.69–0.87, p < 0.001), demonstrating a satisfactory level of discriminatory power. Our study elucidated that PPARγ and PPARGC1A polymorphisms were interactively associated with BRONJ development. These results have potential implications for tailoring personalized treatments for females undergoing bisphosphonate therapy for osteoporosis.

1. Introduction

Bisphosphonates are widely used to treat bone diseases, such as osteoporosis and cancer-induced bone metastasis. They impair osteoclast activity, thereby increasing bone mass [1]. While bisphosphonates are generally well tolerated, they can cause renal toxicity, gastric ulcers, and osteonecrosis of the jaw (ONJ) [2].
Bisphosphonate-related ONJ (BRONJ) is a rare but severe adverse effect related to bisphosphonate therapy [3]. Clinical BRONJ manifestations include the presence of exposed bone in the maxillofacial area that persists for more than eight weeks in patients who have been administered bisphosphonates, particularly without any prior history of radiation therapy [4]. Since the first case reported by Marx in 2003, many studies have been released to distinguish the pathophysiology and risk factors of BRONJ [5]. Nevertheless, the precise underlying cause remains unidentified. Some intriguing pathogenic mechanisms have been proposed, such as the disturbance of the normal bone turnover cycle, the hindrance of the cellular wound healing process, and impaired angiogenesis [6,7,8]. In terms of risk factors, several have been reported in association with BRONJ development, including underlying malignant disease, the intravenous use of bisphosphonates, prolonged or high-dose bisphosphonate therapy, concomitant use of certain medications, dental infections, and jawbone surgical procedures [9]. In addition to environmental factors, genetic predisposition has been identified as a risk factor [10,11]. Several studies have reported BRONJ-associated genes such as CYP2C8, RBMS3, SIRT1, COL1A1, RANK, CYP19A1, VEGFA, and PPARγ through genome-centered or single-nucleotide polymorphism (SNP) analyses [12,13,14,15,16].
Peroxisome proliferator-activated receptor gamma (PPARγ), a nuclear receptor of transcription factors, regulates various physiological processes, such as adipogenesis, glucose and lipid metabolism, inflammation, angiogenesis, and bone remodeling [17,18,19,20]. In bone remodeling, PPARγ inhibits the differentiation of mesenchymal stem cells into osteoblasts and promotes the differentiation of hematopoietic stem cells into osteoclasts. PPARγ modulates transcription by recruiting coactivators [18]. One specific coactivator of PPARγ in bone homeostasis is PPARγ coactivator 1 alpha (PPARGC1A). Recent studies have reported the associations between PPARGC1A and the regulation of bone homeostasis [21,22]. Regarding these findings, genetic variants of PPARγ and PPARGC1A may affect bisphosphonate’s mechanism of action in bone remodeling, contributing to the development of BRONJ.
The association between PPARγ polymorphisms and the onset of BRONJ has been demonstrated in previous research [16,23,24]. However, a meta-analysis found no significant relationship between this gene and the BRONJ with low-quality evidence [10]. The majority of studies concentrated on patients with multiple myeloma or solid cancers, leaving a research gap regarding osteoporosis patients. Moreover, the impact of PPARGC1A polymorphisms on BRONJ has not been investigated. Our study objective was to examine the relationship between PPARγ and PPARGC1A polymorphisms and the development of BRONJ in patients with osteoporosis patients undergoing bisphosphonate treatment.

2. Results

Of the 147 patients who were initially screened, 24 were excluded as 2 were male, 20 used bisphosphonates for cancer, and 2 had insufficient clinical data. Consequently, a total of 123 patients were eligible for this study, among whom 56 patients (45.5%) developed BRONJ. The mean age and duration of bisphosphonate use were 72.9 ± 9.5 years and 54.7 ± 79.1 months, respectively. The relationship between baseline characteristics and BRONJ is demonstrated in Table 1. Age was statistically significantly associated with BRONJ, as individuals aged over 70 years were at a higher risk compared to others (p = 0.003). The BRONJ group had a higher hypertension prevalence compared to the non-BRONJ group (62.5% vs. 41.8%, p = 0.022). No significant differences were observed in the types and routes of bisphosphonate administration between the cases and controls. However, there was a trend towards an increased risk of BRONJ with a bisphosphonate use duration of 60 months or longer (p = 0.001).
Table 2 displays the genotypes and minor allele frequency (MAF) of the analyzed SNPs in PPARγ and PPARGC1A. Among the PPARγ SNPs studied, there were no significant SNPs showing an association with BRONJ. However, when examining PPARGC1A, rs2946385 and rs10020457 were significantly associated with the occurrence of BRONJ. Additionally, when considering the combination of PPARγ/PPARGC1A SNPs, rs1151999/rs294638 and rs1151999/rs10020457 displayed a significant association with BRONJ.
As shown in Table 3, in the univariate analysis, patients with the GG genotype of PPARGC1A rs2946385 had a higher frequency of BRONJ development than individuals with the T allele (Odds ratio (OR) 2.18; 95% confidence interval (CI): 1.06–4.50). For PPARGC1A rs10020457, individuals with the A allele had a 3.0-fold increased risk of developing BRONJ compared to those with GG genotypes (95% CI: 1.15–7.54). We conducted the multivariable analyses using three models to evaluate the independent and interactive associations between the SNPs and BRONJ. Model I included demographic and genetic factors that exhibited a p < 0.1 in the univariate analysis. Only age ≥70 years and a duration of bisphosphonate use ≥60 months were statistically significantly associated with BRONJ after adjusting for the covariates.
To assess the interactive effect of combined genotypes, we constructed additional models. Model II involved the combination of PPARγ rs1151999 GG and PPARGC1A rs2946385 GG genotypes, while Model III involved PPARγ rs1151999 GT and PPARGC1A rs10020457 GA genotypes. In Model II, patients with the combined genotypes of GG in both PPARγ rs1151999 and PPARGC1A rs2946385 showed a higher risk of developing BRONJ compared to patients with other genotype combinations (adjusted odds ratio (AOR): 3.03, 95% CI: 1.01–9.11). However, the effect of combined genotypes in Model III was not significant in the multivariable analysis. Regardless of the models, age (≥70 years) and duration of bisphosphonate use (≥60 months) remained significant factors. The result of the Hosmer–Lemeshow test of Model II indicated a reasonably good fit to the observed data (χ2 = 3.998, 7 degrees of freedom, and p = 0.780). The area under the receiver operating curve (AUROC) for the predicted probability of the multivariable logistic regression for Model II was 0.78 (95% CI: 0.69–0.87, p < 0.001; Figure 1), demonstrating a satisfactory level of discriminatory power.

3. Discussion

Previous studies have demonstrated the relationship between PPARγ and BRONJ; however, conflicting results were found. Furthermore, these studies were confined to cancer patients, leaving a research gap concerning patients with osteoporosis. In this study, our objective was to investigate the association between PPARγ and PPARGC1A polymorphisms and the development of BRONJ in osteoporosis patients undergoing bisphosphonate treatment. This study indicated that the combined effect of PPARγ rs1151999 GG and PPARGC1A rs2946385 GG genotypes was related to BRONJ, along with age and treatment duration after adjusting for confounders. The AUROC for predicted probability revealed good model performance.
Although the individual SNPs of PPARγ rs1151999 and PPARGC1A rs2946385 were not statistically significant in the multivariable analysis, a significant interactive effect was observed when these SNPs were analyzed in combination. PPARGC1A, a coactivator of PPARγ, interacts with factors involved in RNA splicing and transcript elongation, suggesting a potential role in mRNA maturation [25]. A study conducted by Ruchat et al. reported independent and interactive associations between PPARγ and PPARGC1A with insulin and glucose homeostasis [26]. Subjects who had the reference allele in PPARGC1A and the variant allele of PPARγ exhibited higher fasting insulin levels, insulin resistance, and insulin area under the curve compared to those with the other combinations. An interactive association could occur when both the transcriptional factor and its coactivator possess the polymorphism, resulting in maximal and synergistic effects on the development of BRONJ.
PPARγ, located on chromosome 3p25.3, plays a vital role in maintaining skeletal homeostasis. Prolonged use of thiazolidinediones, which are synthetic agents that activate PPARγ, has been correlated with an increased incidence of fracture in individuals with diabetes mellitus [27,28]. The study findings revealed that patients with osteoporosis exhibited an increased number of adipose cells in the bone marrow. This was observed as a result of the stimulation of PPARγ, which led to the differentiation of the stem cells into adipocytes. The association between increased adipose cells and osteoporosis can be attributed to the inhibitory effect of PPARγ on bone formation [29,30]. Hence, polymorphisms in the PPARγ gene are often linked to a higher susceptibility to skeletal disorders. A cohort study conducted in Denmark found that genetic variation increased the risk of vertebral fracture [31]. A pilot study conducted by Dragojevič et al. also indicated that the PPARγ polymorphism could play a role in the incidence of non-traumatic hip fractures among older individuals [32].
Regarding BRONJ, several studies have been conducted on the PPARγ rs1152003 polymorphism. However, the recently published meta-study did not find a significant association between the polymorphism of rs1152003 and the incidence of BRONJ [10]. Similarly, our study observed no significant difference in the association of PPARγ rs1152003 with BRONJ. With respect to PPARγ rs1151999, it did not show a significant association with BRONJ, whereas its combination with either PPARGC1A rs2946385 or rs10020457 demonstrated a significant correlation with BRONJ in the univariate analysis. Although rs1151999 is an intron variant, intron regions may impact mRNA splicing, thereby modifying protein expression or activity [33,34]. Kiel et al. showed that a haplotype containing the G allele of rs1151999 was associated with decreased femoral neck bone mineral density [35]. This finding is consistent with our study, which showed that the GG genotype was associated with an increased risk of BRONJ, likely due to decreased bone mineral density. In addition to bone-related effects, the polymorphism of this gene has been shown to impact other systems significantly. Cirelli et al. demonstrated that the wild-type allele of rs1151999 significantly reduced the risk of hyperglycemia, hyperlipidemia, and obesity even after adjustments for age, sex, and smoking status in a Brazilian population [36]. Moreover, the rs1151999 SNP was associated with protection against Alzheimer’s disease [37]. While this potential mechanism is intriguing, there is no evidence that rs1151999 directly affects PPARγ mRNA expression levels. Further research is needed for a more comprehensive understanding of this SNP.
PPARGC1A is a transcriptional coactivator that interacts with PPARγ, allowing it to interact with multiple transcription factors. While the role of PPARGC1A in energy metabolism is well established, its impact on bone in humans has been barely studied. However, recent in vitro and animal studies have shown the involvement of PPARGC1A in bone-related processes. In in vitro studies, it positively influenced SIRT3 activity, thereby promoting osteogenic differentiation [38] and increasing the expression of essential genes related to osteocyte function [39]. In aged mice, a PPARGC1A deficiency has negative effects on bone mass and strength [40]. PPARGC1A-deficient mice exhibited impaired bone structure and lower trabecular thickness, leading to a significant decline in bending strength. Consistently, Yu et al. demonstrated that PPARGC1A regulated the balance between bone and fat in a mouse model of osteoporosis. Loss of function of this gene in skeletal stem cells inhibited bone formation and facilitated the accumulation of adipose tissue in the bone marrow, while the induction alleviated bone loss [21].
PPARGC1A also plays a crucial role in the angiogenesis and migration of associated endothelial cells [41]. The growth of blood vessels is crucial for initiating and sustaining the wound healing process. Bi et al. demonstrated that bisphosphonate was the primary cause of BRONJ-like disease in mice, mediated in part by its ability to suppress osseous angiogenesis [42]. Consequently, the anti-angiogenic effect has been considered as a possible mechanism underlying BRONJ in humans. Arduino et al. investigated the polymorphism of the vascular endothelial growth factor (VEGF) gene in relation to female patients who developed BRONJ in Italy, and the results suggested a possible haplotype effect of VEGF polymorphism expression [15]. In a Korean cohort, VEGF polymorphisms showed a significant association with an increased risk of BRONJ [43]. These findings imply that the PPARGC1A polymorphisms can contribute to the pathogenesis of BRONJ through multiple mechanisms.
PPARGC1A rs2946385, a synonymous variant, was a significant risk factor in our study. While the specific function of this SNP in maintaining bone homeostasis has not been thoroughly characterized, it has been investigated in relation to type 2 diabetes mellitus (T2DM). A study conducted on Indian patients indicated that the variant allele was associated with increased susceptibility to T2DM (OR 1.72, 95% CI: 1.05–2.82) [44]. This finding could be interpreted that the polymorphism resulted in decreased mRNA expression of PPARGC1A, which eventually influenced insulin secretion. In contrast, this SNP polymorphism did not demonstrate any significant association with T2DM in the Caucasian population [45]. The inconsistent findings could be attributed to ethnic diversity or possible interactions among genetic and environmental factors.
Among environmental factors, age and the duration of bisphosphate use are widely acknowledged as potential risk factors associated with BRONJ [9,46]. In a large population-based study involving 58,069 participants, individuals who were 65 years or older exhibited a 3.0-fold (95% CI: 2.1–4.2) higher risk of BRONJ compared to those in the age group of 20–64 years [47]. In our study, hypertension was associated with a 2.1-fold higher risk of BRONJ, although statistical significance was not found in the multivariable analysis. There is no specific report providing detailed information on the relationship between these two factors, but this association could be explained by the higher prevalence of hypertension among older patients.
Another identified risk factor for BRONJ is the duration of bisphosphonate use. The onset of BRONJ differs depending on the potency of bisphosphonates as intravenous administration results in higher drug exposure compared to oral administration [48]. A prospective study investigating intravenous bisphosphonate (pamidronate and zoledronate) demonstrated that bisphosphonate use over 30 months increased the risk of developing ONJ [49]. However, a review of 103 patients who received oral bisphosphonates (alendronate and risedronate) reported a mean time to the onset of the disease as 4.6 years [50]. In our study, 74 (60.2%) patients were administered bisphosphonate orally, and 12 (9.8%) patients who initially received bisphosphonates intravenously switched to the oral route. Consistent with these findings, our study revealed that advanced age and long-term use of bisphosphonate were associated with a higher risk of BRONJ.
Some limitations need to be addressed. First, the sample size is relatively small. Second, potential confounding factors, such as smoking history and medication use, were not considered in the analysis, which could affect the occurrence of BRONJ. Third, the exclusion of two male patients to focus on gender-specific effects may limit the generalizability of the findings. Lastly, we could not measure the sole impact of the administration route because many patients who initially received bisphosphonates intravenously switched to the oral route. Additional comprehensive cohort studies are necessary to account for potential confounders and validate our findings.

4. Materials and Methods

4.1. Study Participants and Data Collection

We prospectively conducted this nested case–control study at the Ewha Womans University Mokdong Hospital between January 2014 and December 2018. Among the patients scheduled for dentoalveolar surgery, those who were currently or previously taking bisphosphonates were enrolled in this study. The inclusion criteria were patients aged 50 years or older with confirmed osteoporosis by a medical professional. Patients were excluded if they met any of the following criteria: (i) male sex, (ii) had cancer requiring antiresorptive drug treatment, or (iii) had insufficient clinical data or sample for DNA analysis. The diagnosis of BRONJ followed the guidelines outlined by the American Association of Oral and Maxillofacial Surgeons [48] and was made by maxillofacial surgeons. According to the guidelines, cases were considered to have developed BRONJ if all the following clinical manifestations were present: (i) the patient was receiving bisphosphonate treatment when BRONJ occurred, (ii) the exposed and necrotic bone in the maxillofacial area persisted for more than 8 weeks, and (iii) the patient had no history of head or neck radiation therapy. Controls consisted of patients under bisphosphonate treatment without any signs of BRONJ.
The study protocol was reviewed and approved by the Institutional Review Board (IRB) of Ewha Womans University Mokdong Hospital (IRB Number: 14-13-01). This study followed the ethical principles outlined in the 1694 Declaration of Helsinki and its subsequent amendments. We obtained approval from the ethics committee to ensure adherence to the required ethical standards for human research. All participants provided written informed consent before enrollment.
We reviewed the electronic medical records, gathering information such as age, gender, co-morbidities, the indication of bisphosphonates, types and doses of bisphosphonates, bisphosphonate therapy duration, and administration method.

4.2. Selection of Single-Nucleotide Polymorphisms (SNPs) and Genotyping Methods

We referenced previous research to select SNPs and determine the genotyping methods [31,51,52,53,54,55]. Based on this information, we identified five SNPs of PPARγ (one missense SNP (rs1801282), one synonymous SNP (rs3856806), one 3′-untranslated region (UTR) SNP (rs1152003), and two intron SNPs (rs1151999 and rs1175543)) and six SNPs of PPARGC1A (two missense SNPs (rs8192678 and rs3736265), two synonymous SNPs (rs2946385 and rs2970847), and two intron SNPs (rs10020457 and rs7665116). We obtained the details regarding these 11 SNPs, including chromosomal position, reference/alternate alleles, and functional information, from the SNP database of the National Center for Biotechnology Information. We collected Asian MAF and linkage disequilibrium (LD) data from Haploreg v4.2 [56].
Saliva samples were collected for genomic DNA using the tube format (OG300) of the Oragene® DNA Self-Collection Kit (DNAgenotek, Kanata, ON, Canada). The selected 11 SNPs were genotyped using the TaqMan assay. The RT-PCR was performed using the TaqMan® allele discrimination technique on an ABI 7300 instrument (Applied Biosystems, Carlsbad, CA, USA). The PCR was conducted in a 25 μL optical 8-cap strip, consisting of 11.25 μL of DNA samples and 13.75 μL of PCR mix. The PCR reagent mixture included 12.5 μL of the TaqMan Genotyping Master Mix and 1.25 μL of the 20× TaqMan SNP Genotyping Assay Mix (Applied Biosystems in Foster City, CA, USA). Following a denaturation step at 95 °C for 10 min, the PCR was carried out for 40 cycles at 92 °C for 15 s and 60 °C for 60 s.

4.3. Statistical Analysis

We used the chi-squared test and Fisher’s exact test to compare categorical variables between patients with and without BRONJ. We analyzed each SNP in both dominant and recessive models and selected the appropriate model considering the effect size and statistical significance. As PPARGC1A is a coactivator of PPARγ, we also investigated the combined genetic effects of PPARγ/PPARGC1A on BRONJ development by grouping genotypes with similar trends.
To determine the independent risk factors for BRONJ, we conducted a multivariable regression analysis, including variables with a p < 0.1 from the univariate analysis. We entered the variables by stepwise selection when the p-value was lower than 0.05 and removed them when it was higher than 0.1. The unadjusted OR and AOR with the 95% CI were calculated from univariate and multivariable analyses, respectively. The Hosmer–Lemeshow goodness-of-fit test was utilized to evaluate the adequacy of the prediction model’s fit. We evaluated the model’s discriminatory power using the AUROC analysis. A value of p < 0.05 indicated statistical significance. We used SPSS version 20 (IBM, Chicago, IL, USA) for all statistical analyses.

5. Conclusions

Our findings demonstrated an interactive association between the polymorphisms in PPARγ rs1151999/PPARGC1A rs2946385 with the development of BRONJ, along with age (≥70 years) and duration of bisphosphonate use (≥60 months). After validation in a larger cohort study, these results have the potential to facilitate early identification of high-risk patients and inform personalized osteoporosis treatment strategies for individuals undergoing bisphosphonate therapy.

Author Contributions

All the authors have made substantial contributions to the conception of the study. J.S.K., J.W.K., J.E.C. and H.S.G. contributed to designing the study. J.W.K., J.Y. and S.J.K. contributed to material preparation and data collection. J.S.K. and J.Y. performed data analysis and interpretation. J.S.K. and J.W.K. contributed to the drafting of the manuscript. J.E.C. and H.S.G. contributed to the critical revision of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI) funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI22C1377), and the National Research Foundation of Korea (NRF) grants funded by the Ministry of Education (grant number: 2022R1F1A1075439).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Ewha Womans University Mokdong Hospital (IRB number: 14-13-01).

Informed Consent Statement

Informed consent was obtained from all patients involved in the study.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare that they have no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Hong, J.W.; Nam, W.; Cha, I.H.; Chung, S.W.; Choi, H.S.; Kim, K.M.; Kim, K.J.; Rhee, Y.; Lim, S.K. Oral bisphosphonate-related osteonecrosis of the jaw: The first report in Asia. Osteoporos. Int. 2010, 21, 847–853. [Google Scholar] [CrossRef] [PubMed]
  2. Coleman, R.E. Risks and benefits of bisphosphonates. Br. J. Cancer 2008, 98, 1736–1740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Marx, R.E.; Sawatari, Y.; Fortin, M.; Broumand, V. Bisphosphonate-induced exposed bone (osteonecrosis/osteopetrosis) of the jaws: Risk factors, recognition, prevention, and treatment. J. Oral Maxillofac. Surg. 2005, 63, 1567–1575. [Google Scholar] [CrossRef] [PubMed]
  4. Khan, A.A.; Morrison, A.; Hanley, D.A.; Felsenberg, D.; McCauley, L.K.; O’Ryan, F.; Reid, I.R.; Ruggiero, S.L.; Taguchi, A.; Tetradis, S.; et al. Diagnosis and management of osteonecrosis of the jaw: A systematic review and international consensus. J. Bone Miner. Res. 2015, 30, 3–23. [Google Scholar] [CrossRef]
  5. Marx, R.E. Pamidronate (Aredia) and zoledronate (Zometa) induced avascular necrosis of the jaws: A growing epidemic. J. Oral Maxillofac. Surg. 2003, 61, 1115–1117. [Google Scholar] [CrossRef] [PubMed]
  6. Allen, M.R.; Kubek, D.J.; Burr, D.B.; Ruggiero, S.L.; Chu, T.-M.G. Compromised osseous healing of dental extraction sites in zoledronic acid-treated dogs. Osteoporos. Int. 2011, 22, 693–702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Sharma, D.; Ivanovski, S.; Slevin, M.; Hamlet, S.; Pop, T.S.; Brinzaniuc, K.; Petcu, E.B.; Miroiu, R.I. Bisphosphonate-related osteonecrosis of jaw (BRONJ): Diagnostic criteria and possible pathogenic mechanisms of an unexpected anti-angiogenic side effect. Vasc. Cell 2013, 5, 1. [Google Scholar] [CrossRef] [Green Version]
  8. Landesberg, R.; Cozin, M.; Cremers, S.; Woo, V.; Kousteni, S.; Sinha, S.; Garrett-Sinha, L.; Raghavan, S. Inhibition of oral mucosal cell wound healing by bisphosphonates. J. Oral Maxillofac. Surg. 2008, 66, 839–847. [Google Scholar] [CrossRef] [Green Version]
  9. Khan, A.A.; Sándor, G.K.; Dore, E.; Morrison, A.D.; Alsahli, M.; Amin, F.; Peters, E.; Hanley, D.A.; Chaudry, S.R.; Lentle, B.; et al. Bisphosphonate associated osteonecrosis of the jaw. J. Rheumatol. 2009, 36, 478–490. [Google Scholar] [CrossRef] [Green Version]
  10. Guo, Z.; Cui, W.; Que, L.; Li, C.; Tang, X.; Liu, J. Pharmacogenetics of medication-related osteonecrosis of the jaw: A systematic review and meta-analysis. Int. J. Oral Maxillofac. Surg. 2020, 49, 298–309. [Google Scholar] [CrossRef]
  11. Sandro Pereira da Silva, J.; Pullano, E.; Raje, N.S.; Troulis, M.J.; August, M. Genetic predisposition for medication-related osteonecrosis of the jaws: A systematic review. Int. J. Oral Maxillofac. Surg. 2019, 48, 1289–1299. [Google Scholar] [CrossRef]
  12. Yang, G.; Singh, S.; Chen, Y.; Hamadeh, I.S.; Langaee, T.; McDonough, C.W.; Holliday, L.S.; Lamba, J.K.; Moreb, J.S.; Katz, J.; et al. Pharmacogenomics of osteonecrosis of the jaw. Bone 2019, 124, 75–82. [Google Scholar] [CrossRef] [PubMed]
  13. Katz, J.; Gong, Y.; Salmasinia, D.; Hou, W.; Burkley, B.; Ferreira, P.; Casanova, O.; Langaee, T.Y.; Moreb, J.S. Genetic polymorphisms and other risk factors associated with bisphosphonate induced osteonecrosis of the jaw. Int. J. Oral Maxillofac. Surg. 2011, 40, 605–611. [Google Scholar] [CrossRef]
  14. La Ferla, F.; Paolicchi, E.; Crea, F.; Cei, S.; Graziani, F.; Gabriele, M.; Danesi, R. An aromatase polymorphism (g.132810C>T) predicts risk of bisphosphonate-related osteonecrosis of the jaw. Biomark. Med. 2012, 6, 201–209. [Google Scholar] [CrossRef] [PubMed]
  15. Arduino, P.G.; Menegatti, E.; Scoletta, M.; Battaglio, C.; Mozzati, M.; Chiecchio, A.; Berardi, D.; Vandone, A.M.; Donadio, M.; Gandolfo, S.; et al. Vascular endothelial growth factor genetic polymorphisms and haplotypes in female patients with bisphosphonate-related osteonecrosis of the jaws. J. Oral Pathol. Med. 2011, 40, 510–515. [Google Scholar] [CrossRef]
  16. Kastritis, E.; Melea, P.; Bagratuni, T.; Melakopoulos, I.; Gavriatopoulou, M.; Roussou, M.; Migkou, M.; Eleutherakis-Papaiakovou, E.; Terpos, E.; Dimopoulos, M.A. Genetic factors related with early onset of osteonecrosis of the jaw in patients with multiple myeloma under zoledronic acid therapy. Leuk. Lymphoma 2017, 58, 2304–2309. [Google Scholar] [CrossRef]
  17. Baroi, S.; Czernik, P.J.; Chougule, A.; Griffin, P.R.; Lecka-Czernik, B. PPARG in osteocytes controls sclerostin expression, bone mass, marrow adiposity and mediates TZD-induced bone loss. Bone 2021, 147, 115913. [Google Scholar] [CrossRef] [PubMed]
  18. Wei, W.; Wang, X.; Yang, M.; Smith, L.C.; Dechow, P.C.; Sonoda, J.; Evans, R.M.; Wan, Y. PGC1beta mediates PPARgamma activation of osteoclastogenesis and rosiglitazone-induced bone loss. Cell Metab. 2010, 11, 503–516. [Google Scholar] [CrossRef] [Green Version]
  19. Wagner, N.; Wagner, K.D. PPARs and Angiogenesis-Implications in Pathology. Int. J. Mol. Sci. 2020, 21, 5723. [Google Scholar] [CrossRef]
  20. Tontonoz, P.; Spiegelman, B.M. Fat and beyond: The diverse biology of PPARgamma. Annu. Rev. Biochem. 2008, 77, 289–312. [Google Scholar] [CrossRef]
  21. Yu, B.; Huo, L.; Liu, Y.; Deng, P.; Szymanski, J.; Li, J.; Luo, X.; Hong, C.; Lin, J.; Wang, C.Y. PGC-1α Controls Skeletal Stem Cell Fate and Bone-Fat Balance in Osteoporosis and Skeletal Aging by Inducing TAZ. Cell Stem Cell 2018, 23, 193–209. [Google Scholar] [CrossRef] [Green Version]
  22. Buccoliero, C.; Dicarlo, M.; Pignataro, P.; Gaccione, F.; Colucci, S.; Colaianni, G.; Grano, M. The Novel Role of PGC1α in Bone Metabolism. Int. J. Mol. Sci. 2021, 22, 4670. [Google Scholar] [CrossRef]
  23. Di Martino, M.T.; Arbitrio, M.; Guzzi, P.H.; Leone, E.; Baudi, F.; Piro, E.; Prantera, T.; Cucinotto, I.; Calimeri, T.; Rossi, M.; et al. A peroxisome proliferator-activated receptor gamma (PPARG) polymorphism is associated with zoledronic acid-related osteonecrosis of the jaw in multiple myeloma patients: Analysis by DMET microarray profiling. Br. J. Haematol. 2011, 154, 529–533. [Google Scholar] [CrossRef] [Green Version]
  24. Van Poznak, C.; Reynolds, E.L.; Estilo, C.L.; Hu, M.; Schneider, B.P.; Hertz, D.L.; Gersch, C.; Thibert, J.; Thomas, D.; Banerjee, M.; et al. Osteonecrosis of the jaw risk factors in bisphosphonate-treated patients with metastatic cancer. Oral Dis. 2022, 28, 193–201. [Google Scholar] [CrossRef] [PubMed]
  25. Soyal, S.; Krempler, F.; Oberkofler, H.; Patsch, W. PGC-1alpha: A potent transcriptional cofactor involved in the pathogenesis of type 2 diabetes. Diabetologia 2006, 49, 1477–1488. [Google Scholar] [CrossRef] [PubMed]
  26. Ruchat, S.-M.; Weisnagel, S.; Vohl, M.; Rankinen, T.; Bouchard, C.; Perusse, L. Evidence for interaction between PPARG Pro12Ala and PPARGC1A Gly482Ser polymorphisms in determining type 2 diabetes intermediate phenotypes in overweight subjects. Exp. Clin. Endocrinol. Diabetes 2009, 117, 455–459. [Google Scholar] [CrossRef]
  27. Ali, A.A.; Weinstein, R.S.; Stewart, S.A.; Parfitt, A.M.; Manolagas, S.C.; Jilka, R.L. Rosiglitazone causes bone loss in mice by suppressing osteoblast differentiation and bone formation. Endocrinology 2005, 146, 1226–1235. [Google Scholar] [CrossRef] [Green Version]
  28. Meunier, P.; Aaron, J.; Edouard, C.; Vignon, G. Osteoporosis and the replacement of cell populations of the marrow by adipose tissue: A quantitative study of 84 iliac bone biopsies. Clin. Orthop. Relat. Res. 1971, 80, 147–154. [Google Scholar] [CrossRef] [PubMed]
  29. Schwartz, A.V.; Sellmeyer, D.E.; Vittinghoff, E.; Palermo, L.; Lecka-Czernik, B.; Feingold, K.R.; Strotmeyer, E.S.; Resnick, H.E.; Carbone, L.; Beamer, B.A.; et al. Thiazolidinedione use and bone loss in older diabetic adults. J. Clin. Endocrinol. Metab. 2006, 91, 3349–3354. [Google Scholar] [CrossRef]
  30. Lecka-Czernik, B. Bone loss in diabetes: Use of antidiabetic thiazolidinediones and secondary osteoporosis. Curr. Osteoporos. Rep. 2010, 8, 178–184. [Google Scholar] [CrossRef] [Green Version]
  31. Harsløf, T.; Tofteng, C.L.; Husted, L.B.; Nyegaard, M.; Børglum, A.; Carstens, M.; Stenkjær, L.; Brixen, K.; Eiken, P.; Jensen, J.E.; et al. Polymorphisms of the peroxisome proliferator-activated receptor γ (PPARγ) gene are associated with osteoporosis. Osteoporos. Int. 2011, 22, 2655–2666. [Google Scholar] [CrossRef] [PubMed]
  32. Dragojevič, J.; Ostanek, B.; Mencej-Bedrač, S.; Komadina, R.; Preželj, J.; Marc, J. PPARG gene promoter polymorphism is associated with non-traumatic hip fracture risk in the elderly Slovenian population: A pilot study. Clin. Biochem. 2011, 44, 1085–1089. [Google Scholar] [CrossRef] [PubMed]
  33. Raponi, M.; Baralle, D. Alternative splicing: Good and bad effects of translationally silent substitutions. FEBS J. 2010, 277, 836–840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Pagani, F.; Baralle, F.E. Genomic variants in exons and introns: Identifying the splicing spoilers. Nat. Rev. Genet. 2004, 5, 389–396. [Google Scholar] [CrossRef]
  35. Kiel, D. Polymorphisms in the PPARr gene influence bone density in humans. J. Bone Miner. Res. 2005, 20, S234. [Google Scholar]
  36. Cirelli, T.; Nicchio, I.G.; Bussaneli, D.G.; Silva, B.R.; Nepomuceno, R.; Orrico, S.R.; Cirelli, J.A.; Theodoro, L.H.; Barros, S.P.; Scarel-Caminaga, R.M. Evidence Linking PPARG Genetic Variants with Periodontitis and Type 2 Diabetes Mellitus in a Brazilian Population. Int. J. Mol. Sci. 2023, 24, 6760. [Google Scholar] [CrossRef]
  37. Combarros, O.; Rodríguez-Rodríguez, E.; Mateo, I.; Vázquez-Higuera, J.L.; Infante, J.; Berciano, J.; Sánchez-Juan, P. APOE dependent-association of PPAR-γ genetic variants with Alzheimer’s disease risk. Neurobiol. Aging 2011, 32, 541.e1–547.e6. [Google Scholar] [CrossRef]
  38. Ding, Y.; Yang, H.; Wang, Y.; Chen, J.; Ji, Z.; Sun, H. Sirtuin 3 is required for osteogenic differentiation through maintenance of PGC-1α-SOD2-mediated regulation of mitochondrial function. Int. J. Biol. Sci. 2017, 13, 254–264. [Google Scholar] [CrossRef] [Green Version]
  39. Sánchez-de-Diego, C.; Artigas, N.; Pimenta-Lopes, C.; Valer, J.A.; Torrejon, B.; Gama-Pérez, P.; Villena, J.A.; Garcia-Roves, P.M.; Rosa, J.L.; Ventura, F. Glucose Restriction Promotes Osteocyte Specification by Activating a PGC-1α-Dependent Transcriptional Program. iScience 2019, 15, 79–94. [Google Scholar] [CrossRef]
  40. Colaianni, G.; Lippo, L.; Sanesi, L.; Brunetti, G.; Celi, M.; Cirulli, N.; Passeri, G.; Reseland, J.; Schipani, E.; Faienza, M.F.; et al. Deletion of the Transcription Factor PGC-1α in Mice Negatively Regulates Bone Mass. Calcif. Tissue Int. 2018, 103, 638–652. [Google Scholar] [CrossRef]
  41. Kadlec, A.O.; Chabowski, D.S.; Ait-Aissa, K.; Gutterman, D.D. Role of PGC-1α in vascular regulation: Implications for atherosclerosis. Arterioscler. Thromb. Vasc. Biol. 2016, 36, 1467–1474. [Google Scholar] [CrossRef] [Green Version]
  42. Bi, Y.; Gao, Y.; Ehirchiou, D.; Cao, C.; Kikuiri, T.; Le, A.; Shi, S.; Zhang, L. Bisphosphonates cause osteonecrosis of the jaw-like disease in mice. Am. J. Pathol. 2010, 177, 280–290. [Google Scholar] [CrossRef] [PubMed]
  43. Choi, H.; Lee, J.H.; Kim, H.J.; Park, W.; Lee, J.H.; Kim, J.H. Genetic association between VEGF polymorphisms and BRONJ in the Korean population. Oral Dis. 2015, 21, 866–871. [Google Scholar] [CrossRef] [PubMed]
  44. Sharma, R.; Matharoo, K.; Kapoor, R.; Bhanwer, A.J.S. Association of PGC-1α gene with type 2 diabetes in three unrelated endogamous groups of North-West India (Punjab): A case-control and meta-analysis study. Mol. Genet. Genom. 2018, 293, 317–329. [Google Scholar] [CrossRef]
  45. Ek, J.; Andersen, G.; Urhammer, S.A.; Gaede, P.H.; Drivsholm, T.; Borch-Johnsen, K.; Hansen, T.; Pedersen, O. Mutation analysis of peroxisome proliferator-activated receptor-gamma coactivator-1 (PGC-1) and relationships of identified amino acid polymorphisms to Type II diabetes mellitus. Diabetologia 2001, 44, 2220–2226. [Google Scholar] [CrossRef] [Green Version]
  46. Thumbigere-Math, V.; Tu, L.; Huckabay, S.; Dudek, A.Z.; Lunos, S.; Basi, D.L.; Hughes, P.J.; Leach, J.W.; Swenson, K.K.; Gopalakrishnan, R. A retrospective study evaluating frequency and risk factors of osteonecrosis of the jaw in 576 cancer patients receiving intravenous bisphosphonates. Am. J. Clin. Oncol. 2012, 35, 386–392. [Google Scholar] [CrossRef]
  47. Huang, Y.F.; Chang, C.T.; Muo, C.H.; Tsai, C.H.; Shen, Y.F.; Wu, C.Z. Impact of bisphosphonate-related osteonecrosis of the jaw on osteoporotic patients after dental extraction: A population-based cohort study. PLoS ONE 2015, 10, e0120756. [Google Scholar] [CrossRef] [PubMed]
  48. Ruggiero, S.L.; Dodson, T.B.; Fantasia, J.; Goodday, R.; Aghaloo, T.; Mehrotra, B.; O’Ryan, F. American Association of Oral and Maxillofacial Surgeons position paper on medication-related osteonecrosis of the jaw—2014 update. J. Oral Maxillofac. Surg. 2014, 72, 1938–1956. [Google Scholar] [CrossRef]
  49. Boonyapakorn, T.; Schirmer, I.; Reichart, P.A.; Sturm, I.; Massenkeil, G. Bisphosphonate-induced osteonecrosis of the jaws: Prospective study of 80 patients with multiple myeloma and other malignancies. Oral Oncol. 2008, 44, 857–869. [Google Scholar] [CrossRef]
  50. Palaska, P.K.; Cartsos, V.; Zavras, A.I. Bisphosphonates and time to osteonecrosis development. Oncologist 2009, 14, 1154–1166. [Google Scholar] [CrossRef]
  51. Chen, B.; Wang, Y.; Tang, W.; Chen, Y.; Liu, C.; Kang, M.; Xie, J. Association between PPARγ, PPARGC1A, and PPARGC1B genetic variants and susceptibility of gastric cancer in an Eastern Chinese population. BMC Med. Genom. 2022, 15, 274. [Google Scholar] [CrossRef]
  52. Gallicchio, L.; Kalesan, B.; Huang, H.Y.; Strickland, P.; Hoffman, S.C.; Helzlsouer, K.J. Genetic polymorphisms of peroxisome proliferator-activated receptors and the risk of cardiovascular morbidity and mortality in a community-based cohort in washington county, Maryland. PPAR Res. 2008, 2008, 276581. [Google Scholar] [CrossRef] [Green Version]
  53. Weydt, P.; Soyal, S.M.; Gellera, C.; Didonato, S.; Weidinger, C.; Oberkofler, H.; Landwehrmeyer, G.B.; Patsch, W. The gene coding for PGC-1alpha modifies age at onset in Huntington’s Disease. Mol. Neurodegener. 2009, 4, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lai, C.Q.; Tucker, K.L.; Parnell, L.D.; Adiconis, X.; García-Bailo, B.; Griffith, J.; Meydani, M.; Ordovás, J.M. PPARGC1A variation associated with DNA damage, diabetes, and cardiovascular diseases: The Boston Puerto Rican Health Study. Diabetes 2008, 57, 809–816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Geoffroy, P.A.; Etain, B.; Lajnef, M.; Zerdazi, E.H.; Brichant-Petitjean, C.; Heilbronner, U.; Hou, L.; Degenhardt, F.; Rietschel, M.; McMahon, F.J.; et al. Circadian genes and lithium response in bipolar disorders: Associations with PPARGC1A (PGC-1α) and RORA. Genes Brain Behav. 2016, 15, 660–668. [Google Scholar] [CrossRef] [Green Version]
  56. Ward, L.D.; Kellis, M. HaploReg v4: Systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016, 44, D877–D881. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The area under the receiver operating curve (ROC) for the development of bisphosphonate-related osteonecrosis of the jaw. The blue line presents the predicted probability of Model II (demographic factors and the combination of PPARγ rs1151999 GG and PPARGC1A rs2946385 GG genotypes), while the green line presents the reference.
Figure 1. The area under the receiver operating curve (ROC) for the development of bisphosphonate-related osteonecrosis of the jaw. The blue line presents the predicted probability of Model II (demographic factors and the combination of PPARγ rs1151999 GG and PPARGC1A rs2946385 GG genotypes), while the green line presents the reference.
Pharmaceuticals 16 01035 g001
Table 1. Demographic characteristics of study participants.
Table 1. Demographic characteristics of study participants.
CharacteristicsCase (n = 56)Control (n = 67)p-Value
Age (years) 0.003
 <709 (16.1)27 (40.9)
 ≥7047 (83.9)39 (59.1)
Co-morbidity
 Hypertension35 (62.5)28 (41.8)0.022
 Cardiovascular disease7 (12.5)8 (11.9)0.925
 Diabetes mellitus17 (30.4)16 (23.9)0.419
 Rheumatoid arthritis6 (10.7)2 (3.0)0.140 a
 Thyroid disease4 (7.1)2 (3.0)0.286
 Kidney disease2 (3.6)3 (4.5)1.000 a
 Liver disease0 (0)2 (3.0)0.500 a
 Cancer2 (3.6)6 (9.1)0.289 a
Types of bisphosphonates * 0.225 a
 Alendronate9 (34.6)12 (36.4)
 Ibandronate5 (19.2)12 (36.4)
 Pamidronate2 (7.7)0 (0)
 Risedronate10 (38.5)9 (27.3)
Route of administration 0.210
 Intravenous9 (20.5)18 (31.6)
 Oral35 (79.5)39 (68.4)
Treatment duration (months) 0.001
 <6026 (52.0)45 (83.3)
 ≥6024 (48.0)9 (16.7)
a Fisher’s exact test, * including alendronate, ibandronate, pamidronate, and risedronate.
Table 2. Association of PPARγ and PPARGC1A gene polymorphisms with bisphosphonate-related osteonecrosis of the jaw.
Table 2. Association of PPARγ and PPARGC1A gene polymorphisms with bisphosphonate-related osteonecrosis of the jaw.
Gene PolymorphismMAFGrouped GenotypesCase (n = 56)Control (n = 67)p-Value
PPARγ
rs11519990.390GG, GT52 (92.9)55 (82.1)0.077
G > T TT4 (7.1)12 (17.9)
rs38568060.107CC46 (82.1)53 (80.3)0.796
C > T CT, TT10 (17.9)13 (19.7)
rs11520030.402GG, GC47 (83.9)50 (74.6)0.208
G > C CC9 (16.1)17 (25.4)
rs18012820.041CC53 (94.6)61 (91.0)0.508 a
C > G CG, GG3 (5.4)6 (9.0)
rs11755430.402AA, AG50 (90.9)55 (82.1)0.162
A > G GG5 (9.1)12 (17.9)
PPARGC1A
rs81926780.415CC22 (39.3)19 (28.4)0.200
C > T CT, TT34 (60.7)48 (71.6)
rs29463850.291GG35 (62.5)29 (43.3)0.034
G > T GT, TT21 (37.5)38 (56.7)
rs100204570.102GG40 (71.4)59 (88.1)0.020
G > A GA, AA16 (28.6)8 (11.9)
rs76651160.260TT27 (48.2)41 (61.2)0.149
T > C TC, CC29 (51.8)26 (38.8)
rs29708470.240TT, TC23 (41.1)29 (43.3)0.805
T * > C CC33 (58.9)38 (56.7)
rs37362650.146GG38 (67.9)53 (79.1)0.157
G > A GA, AA18 (32.1)14 (20.9)
PPARγ/PPARGC1A
rs1151999/rs2946385N/AGG/GG16 (28.6)9 (13.4)0.038
Others40 (71.4)58 (86.6)
rs1151999/rs10020457N/AGT/GA12 (21.4)2 (3.0)0.001
Others44 (78.6)65 (97.0)
a Fisher’s exact test; * minor allele; MAF: minor allele frequency.
Table 3. Multivariable regression analyses to identify predictors of bisphosphonate-related osteonecrosis of the jaw.
Table 3. Multivariable regression analyses to identify predictors of bisphosphonate-related osteonecrosis of the jaw.
VariablesCrude OR
(95% CI)
Adjusted OR (95% CI)
Model IModel IIModel III
Age ≥ 70 years3.62 (1.52–8.59) **3.94 (1.50–10.30) **3.23 (1.18–8.84) *4.27 (1.57–11.64) **
Hypertension2.32 (1.12–4.80) * 2.13 (0.87–5.24)
Duration ≥ 60 months4.62 (1.87–11.42) ***3.32 (1.28–8.59) *4.27 (1.55–11.77) **3.07 (1.17–8.07) *
PPARγ rs1151999
GG, GT
2.83 (0.86–9.09)
PPARGC1A rs2946385
GG
2.18 (1.06–4.50) *
PPARGC1A rs10020457
GA, AA
2.95 (1.15–7.54) *
PPARγ/PPARGC1A rs1151999/rs2946385
GG/GG
2.41 (1.17–4.99) * 3.03 (1.01–9.11) *
PPARγ/PPARGC1A rs1151999/rs10020457
GT/GA
3.43 (1.29–9.08) * 5.12 (0.89–29.57)
Model I included all factors with p < 0.1 in univariate analysis. Model II and Model III included combined genotype of PPARγ/PPARGC1A rs1151999/rs2946385 and PPARγ/PPARGC1A rs1151999/rs10020457 instead of individual SNPs in Model I, respectively. OR: odds ratio, CI: confidence interval. * p < 0.05, ** p <0.01, *** p < 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, J.S.; Kim, J.W.; Yee, J.; Kim, S.J.; Chung, J.E.; Gwak, H.S. Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis. Pharmaceuticals 2023, 16, 1035. https://doi.org/10.3390/ph16071035

AMA Style

Kim JS, Kim JW, Yee J, Kim SJ, Chung JE, Gwak HS. Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis. Pharmaceuticals. 2023; 16(7):1035. https://doi.org/10.3390/ph16071035

Chicago/Turabian Style

Kim, Jung Sun, Jin Woo Kim, Jeong Yee, Sun Jong Kim, Jee Eun Chung, and Hye Sun Gwak. 2023. "Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis" Pharmaceuticals 16, no. 7: 1035. https://doi.org/10.3390/ph16071035

APA Style

Kim, J. S., Kim, J. W., Yee, J., Kim, S. J., Chung, J. E., & Gwak, H. S. (2023). Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis. Pharmaceuticals, 16(7), 1035. https://doi.org/10.3390/ph16071035

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