Next Article in Journal
Special Issue “Epigenetics in Neurodegenerative Diseases”
Previous Article in Journal
Enhanced Therapeutic Potential of Hybrid Exosomes Loaded with Paclitaxel for Cancer Therapy
Previous Article in Special Issue
Oral Pyogenic Granuloma: A Narrative Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of SIRT1 Gene SNPs and Clinical Characteristics in Medication-Related Osteonecrosis of the Jaw

1
Department of Internal Medicine and Oncology, Semmelweis University, 1083 Budapest, Hungary
2
Department of Oro-Maxillofacial Surgery and Stomatology, Semmelweis University, 1085 Budapest, Hungary
3
Hungarian Research Network SE-ENDOMOLPAT Research Group, 1085 Budapest, Hungary
4
Department of Plant Systematics, Ecology and Theoretical Biology, Eötvös Loránd University, 1117 Budapest, Hungary
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(7), 3646; https://doi.org/10.3390/ijms25073646
Submission received: 8 March 2024 / Revised: 22 March 2024 / Accepted: 22 March 2024 / Published: 25 March 2024

Abstract

:
Certain genetic factors, including single-nucleotide polymorphisms (SNPs) in the SIRT1 gene, have been linked to medication-related osteonecrosis of the jaw (MRONJ). This study examined four SNPs in the SIRT1 gene and implemented multivariate statistical analysis to analyze genetic and clinical factors in MRONJ patients. Genomic DNA was isolated from peripheral blood samples of 63 patients of European origin treated for MRONJ, and four SNP genotypes in the gene encoding the SIRT-1 protein were determined by Sanger sequencing. The allele frequencies measured in the MRONJ population were compared with allele frequencies measured in the European population in the National Center for Biotechnology Information Allele Frequency Aggregator (NCBI ALFA) database. Genetic and clinical factors were examined with multivariate statistical analysis. A C:A allele distribution ratio of 77.8:22.2 was measured in the rs932658 SNP. In the ALFA project, a C:A allele distribution ratio of 59.9:40.1 was detected in the European population, which was found to be a significant difference (p = 4.5 × 10−5). Multivariate statistical analysis revealed a positive correlation (0.275) between the genotype of SNP rs932658 and the number of stages improved during appropriate MRONJ therapy. It is concluded that allele A in SNP rs932658 in the SIRT1 gene acts as a protective factor in MRONJ.

1. Introduction

Medication-related osteonecrosis of the jaw (MRONJ) is a rare but serious side effect of mainly antiresorptive drugs, which significantly impairs the quality of life of affected patients [1]. MRONJ is defined as bone exposed or accessible through an extra or intraoral fistula in the maxillofacial region that has persisted for more than 8 weeks, current or previous treatment with antiresorptive agents alone or in combination with other drugs (e.g., antiangiogenic agents, immune modulators) and no history of metastatic disease or radiation exposure to the jaw [2]. Such antiresorptive agents include nitrogen-containing bisphosphonates and the receptor activator of nuclear factor kappa-Β ligand (RANKL) inhibitor drug denosumab. These antiresorptive drugs are effective in managing conditions, including osteoporosis and skeletal-related events (SREs), such as hypercalcemia of malignancy or pathological fractures due to bone metastasis of solid tumors (e.g., breast cancer, prostate cancer) or multiple myeloma [3,4]. MRONJ cases can be further classified into two distinct categories based on the assumed therapies. These are bisphosphonate-related osteonecrosis of the jaw (BRONJ) and non-bisphosphonate-related osteonecrosis of the jaw (non-BRONJ). Non-BRONJ encompasses instances induced by various pharmacological agents. These agents include denosumab, tyrosine kinase inhibitors (e.g., sunitinib, sorafenib or cabozantinib), everolimus (an mTOR-inhibitor), tocilizumab (an interleukin-6 receptor inhibitor) and anti-vascular endothelial growth factor (anti-VEGF) monoclonal antibodies, such as bevacizumab [5,6]. The prevalence of MRONJ ranges from 0 percent to 18 percent in cancer patients [7,8,9] and 0.02 percent to 0.3 percent in osteoporotic patients, depending on the antiresorptive agent used [10,11,12].
The exact pathomechanism of the disease is still unclear; however, several factors are thought to play a role in the development of MRONJ. Among others, the inhibition of bone remodeling, inflammation and infection are all possible key factors in the pathomechanism [13,14,15]. Over the last decade, a growing number of studies have suggested that genetic predisposition might also be a contributing factor to the development of MRONJ [16,17,18,19]. Several genes have been studied, including CYP2C8 or HERC4, and although whole exome and genome sequencing studies and analysis of candidate genes have not led to the clear identification of the gene that causes MRONJ, research has focused attention on several promising genes. In this context, SIRT1 stands out in view of its potential role in MRONJ development, although further studies are needed to determine its exact role in the process. A previous study suggested that the upregulation of the SIRT1 gene might be a protective factor in the development of MRONJ [17].
The SIRT1 gene encodes the Sirtuin-1 protein, which is a member of the sirtuins, a NAD+-dependent histone deacetylase protein family, involved in several crucial biological processes, such as inflammation, cell metabolism, aging or oxidative stress through the modulation of the activity of several transcription factors and enzymes [20,21]. SIRT1 is also thought to be a key regulator in bone homeostasis via epigenetically regulating several genes involved in bone metabolism [22].
Based on previous whole exome sequencing studies, four single-nucleotide polymorphisms (SNPs) were examined in or in proximity to the SIRT1 gene (rs932658, rs7896005, rs7894483 and rs3758391), which might have a causal role in the development of MRONJ. Figure 1 displays the schematic structure of the SIRT-1 gene with the SNPs examined. Clinical data were also collected and evaluated. The aim of this study was to further understand the possible association between MRONJ and the aforementioned SNPs in the SIRT1 gene in a single-center retrospective study.

2. Results

2.1. Study Population

The mean age of patients included in the study was 69.3 (range 45–91), 73% were female patients (n = 46), and 88.3% (n = 53) had a malignant disorder as primary disease, while 11.7% of patients (n = 7) suffered from osteoporosis as a reason for treatment. The majority of MRONJ cases (73.3%) were diagnosed as stage 2 disease (n = 44), 25% were deemed stage 3 (n = 15) and only one patient had stage 1 MRONJ. Descriptive statistics of the study population are shown in Table 1.

2.2. SIRT1 SNP Genotyping

There was a significant difference detected in the allele distribution of the rs932658 SNP by Chi-square test between the MRONJ group and the healthy European population measured in the National Center for Biotechnology Information Allele Frequency Aggregator (NCBI ALFA) database (p = 4.5 × 10−5) (Table 2). Representative Sanger’s sequence diagrams of SNP rs932658 are displayed in Figure 2.
No further significant differences were seen in the other SNPs. Minor allele frequencies are displayed in Figure 3.

2.3. Multivariate Statistical Analysis

By multivariate statistical analysis of clinical and genetic variables, several strong correlations were found, which indicate the reliability of this novel statistical method. These include a high correlation (0.7063) between the type of bisphosphonate used and the method of bisphosphonate administration (e.g., oral, intravenous) or a relatively high positive correlation (0.2789) between the presence of chemotherapy and the type of bisphosphonate used. These relationships are expressed on the vertical principal component (PCA) axis (Figure 4). As a new finding, a relatively high positive correlation between the genotype of rs932658 and the number of stages improved after the appropriate non-surgical and surgical treatment was detected (0.275). This suggests that genetic factors, specifically the rs932658 SNP in the SIRT1 gene, might play a role not only in the development of MRONJ but also in the clinical course of the disease. Patients with a more favorable genotype might be less likely to have MRONJ and have a higher rate of improvement once MRONJ has developed compared with less favorable genotype patients. As the diagram of Figure 4 shows, the other three genotypes are clustered; they are very highly correlated with one another (0.895 to 0.965). These appear uncorrelated with the other variables, thus explaining the contrast along the first PCA axis. A relatively high positive correlation (0.2417) was also detected between sex and healing. Previous dentoalveolar operation, corticosteroid therapy and diabetes mellitus did not show a correlation with other genetic or clinical variables. Smoking showed a positive correlation with the age of the patient (Supplement Table S2).
The PCoA of 60 patients is displayed in Figure 5 for axes 1 and 2. Apparently, the patients form two groups with a relatively large empty area separating them. To find the original variables best explaining this grouping, correlations between variables and the ordination scores for PCoA axes 1 and 2 were calculated. Then, variables with the highest correlations were selected. These were the three SNPs, rs7894483, rs7896005 and rs3758391, for axis 1 (r = 0.58, 0.61 and 0.59, respectively) and hormone therapy (r = 0.36) as well as the method of administration (r = −0.21) for axis 2. For visualization, the correlations were rescaled to best fit the ordination of patients, as usual in biplots, and arrows were drawn to point to their positions in Figure 5. As seen, the principal factors affecting the grouping of patients are the three genotypes, whereas none of the other clinical or genetic variables is influential.

3. Discussion

The findings of our study corroborated and validated the potential role of genetic factors, particularly SNPs associated with the SIRT1 gene, as well as other clinical factors in the development of MRONJ in a Central European population, which is different from those originally investigated. Our research group has already utilized principal component analysis to evaluate genetic and clinical data in MRONJ patients more than a decade ago [18]. In that study, a strong positive correlation was detected between CYP2C8 rs1934951 polymorphism and the anatomical localization of MRONJ. In an exome-wide association analysis of MRONJ patients with primary oncological disease, the strongest exome-level association was shown by the rs7896005 and rs375839 SNPs at the SIRT1/HERC4 locus, which would be associated with a lower risk of MRONJ [16]. A genome-wide association study of five thousand bisphosphonate-treated European MRONJ samples revealed a strong association of the rs2736308 SNP on chromosome 8 with an increased risk of MRONJ. In silico analyses suggested that this polymorphism indirectly regulates SIRT1 expression [23]. Consistent with the results in this work, the SIRT1 gene SNP rs932658 has also been demonstrated to be in association with MRONJ by regulating SIRT1 expression [17]. Other research groups have suggested potential roles of SNPs in the VEGF, CYP2C8, PPARG and other genes in the development of MRONJ [19,24,25,26,27].
In this study a significant difference was detected in allele distribution for rs932658 between MRONJ patients and the NCBI Allele Frequency Aggregator (ALFA) (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA507278, accessed on 1 January 2024) dataset subgroup measured in the European population, which suggests that individuals with certain genetic variations in the SIRT1 gene may have a lower chance of developing MRONJ. Moreover, the relatively high positive correlation found between the genotype of rs932658 and the number of stages improved during the therapy of MRONJ that was found may suggest that individuals with more favorable genetic variations not only have less chance for developing MRONJ but also have a higher likelihood of healing. This is the first observation of this kind and raises the notion that genetic factors might also influence the response to therapeutic interventions in MRONJ. A previous study suggested that the minor allele (allele A) of rs932658 is associated with higher gene expression of the SIRT1 gene, which might act as a protective factor in the development of MRONJ [17]. Our findings, in an independent Hungarian study cohort, confirm allele A in rs932658 as a potential protective factor in the development of MRONJ and predicts a greater degree of stage improvement in an already established MRONJ disease. This SNP, similar to other polymorphisms evaluated in this study (e.g., rs3758391), is located in the promoter region of the gene, which is universally understood to be crucial in the regulation of gene expression. One possible explanation of the functional importance of rs932658 might be its relative distance from the 5′-end of the SIRT1 gene. However, further studies are called for to understand the mechanism of the effect of the rs932658 SNP on SIRT1 gene expression.
Sirtuin proteins play an important role in bone biology. SIRT1 promotes bone formation by several molecular mechanisms. These include activating osteoblasts via stimulating the Wnt signaling pathway and decreasing the expression of sclerostin, an inhibitor of bone formation [28]. These findings established the start of several clinical trials investigating sirtuin-activating compounds for osteoporosis [28]. Despite significant progress over the past few years, many research questions have yet to be answered about sirtuin activators as potential therapeutic targets.
Consistent with the data in the literature [29], our study suggests that the prevalence of MRONJ is higher in women than in men. This is likely due to the higher prevalence of the primary diseases treated with antiresorptive agents in the female population (e.g., breast cancer, osteoporosis). However, the relatively high positive correlation between the sex of the patient and healing suggests that the patient’s sex might not be irrelevant in the clinical course of MRONJ. In our study group, female patients also had a higher tendency to heal from MRONJ than male patients. Between the sex of the patient and the genotype of rs932658, there was no positive correlation detected, so there might be factors other than genetic factors determining the correlation between the sex of the patient and healing. Other previously described risk factors like corticosteroid therapy or diabetes mellitus did not show a correlation with other clinical or genetic variables in our study.
Principal coordinates analysis revealed distinct patient groupings, primarily influenced by the genotypes of the three SNPs, i.e., rs7894483, rs7896005 and rs3758391, respectively (r = 0.58, 0.61, 0.59). The patients in the distinct groups might act differently in some clinical respects, so further studies may be needed to identify the relevance of this apparent patient grouping.
In discussing the outcomes of our investigation, it is imperative to acknowledge the advantages that underpin the significance of the findings presented in this article. This study was conducted on a relatively large and independent cohort of Central European MRONJ patients. This enabled us to detect genetic associations with greater reliability. Furthermore, a novel multivariate statistical analysis was deployed to assess the possible correlations between clinical and genetic variables. The results of the multivariate statistical analysis further validated the findings made by genotyping. This innovative methodological approach not only adds knowledge to current discourse but also might set a precedent for further investigations in this field.
It is also crucial to acknowledge the limitations of this study. The sample size of sixty-three patients, while providing valuable insights, may not fully capture the diversity of genetic variations that could contribute to MRONJ. Additionally, the study was based on a slightly incomplete dataset, and further research with larger and more diverse cohorts is required to validate and extend these findings.
The findings presented here might be useful in the future due to multiple reasons. Firstly, the possible role of SIRT1 in MRONJ and, more broadly, in bone biology was further validated by our findings. The association of MRONJ with the identified SNP implicates the important role of sirtuins in the pathomechanism of this adverse drug reaction, resulting in a severe decrease in the quality of life of affected patients. These results suggest that the modulation of sirtuins by pharmacological agents (e.g., STACs) might present a promising new field of therapeutics in MRONJ. Moreover, with the improvement of knowledge on the genetic background of MRONJ, new risk stratification tools based on the genetic predisposition of the patient can be developed, further aiding personalized medicine. Furthermore, these SNPs can be determined relatively easily and cost-effectively and might help patients avoid MRONJ, which could ease the load on caregivers and help patients and doctors concentrate on the management of the often malignant primary disease.
There is still a definite need for further research regarding the wide genetic landscape of MRONJ and, specifically, the role of the SIRT1 gene in the disease. Since genetic factors may vary slightly in the investigated population, further validation of the association between rs932658 and the risk of MRONJ should be conducted in populations with different genetic backgrounds.
In summary, this research provides insights into the genetic and clinical factors underlying MRONJ, thus promoting a better understanding of its development and progression. The introduction of novel multivariate statistical analysis might also contribute to increasing our knowledge of MRONJ. There are some MRONJ risk assessment questionnaires published in the literature [30,31], whereas no prevalent, internationally accepted and validated a priori risk estimation algorithm for MRONJ has been available. In addition, there is absolutely no procedure that takes any genetic data into account in the assessment. The potential integration of genetic data into clinical management offers the possibility of personalized treatment in the future. A future goal would be to screen out high-risk patients for whom antiresorptive therapy should or should not be considered. Contribution to that goal with the results of this study is intended.

4. Materials and Methods

4.1. Patients

Sixty-three consecutive patients of Hungarian origin suffering from MRONJ were enrolled in this study. All MRONJ patients were treated at the Department of Oro-Maxillofacial Surgery and Stomatology, Semmelweis University. After detailed medical history taking, physical and laboratory examination, peripheral blood samples were taken from the patients. For three of them, the collection of previous medical data was unsuccessful; therefore, these patients were excluded from multivariate data analysis. All patients were informed and asked for consent before enrolling them in this study. The 1964 Declaration of Helsinki and its subsequent amendments were followed in this study. The study had the approval of the institutional ethics committee of Semmelweis University of Budapest (Semmelweis University Regional and Institutional Committee of Science and Research Ethics, protocol code: 10862-1/2016/EKU; IV/4613-4/2020/EKU). Access to human data was restricted to members of the research team. A password-protected medium was used to store the data. The authors confirm that all the methods were carried out in accordance with the applicable guidelines and regulations.

4.2. Genotyping

Four SNPs in the SIRT1 gene were genotyped in the 63 subjects. Peripheral blood samples were collected, and genomic DNA was extracted from each subject using High Pure PCR Template Purification kit (Roche Diagnostics, GmbH, Mannheim, Germany). The PCR mixture (20 μL) contained 1 μL genomic DNA (50 ng/μL), 2 × 0.50 μL primer, 10 μL of JumpStart Taq Readymix (SIGMA-ALDRICH, Co., 3050 Spruce Street, St. Louis, MO 63103 USA) and 8 μL ultrapure PCR water. Cycling conditions comprised an initial cycle at 94 °C for 5 min followed by 40 cycles of 94 °C for 30 s, 60 °C for 30 s and by a final step at 72 °C for 1 min. Sanger dideoxy sequencing was performed to determine the genotype of the SNPs (Eurofins Genomics Europe GmbH, Ebersberg, Germany).

4.3. Statistical Analysis

The distribution of SIRT1 SNP genotypes and alleles in the NCBI Allele Frequency Aggregator (ALFA) reference database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA507278, accessed on 1 January 2024) and MRONJ populations was analyzed by Chi-square tests. Statistics with a p-value < 0.05 were considered significant. Chi-square tests were performed using IBM SPSS 28.0 software (SPSS Inc., Chicago, IL, USA).
For further analysis, the data were summarized in a matrix with 60 rows (patients) and 26 columns (variables) (Supplement Table S1). Conventional multivariate techniques cannot be applied to this dataset for two reasons. The variables are not homogeneous regarding the measurement scale. A few are ratio-scale variables (age, duration of treatment and number of recurrences) for which arithmetic differences between values are meaningful. For the ordinal variables (stage, stage improvement), only the sequence of possible states can be interpreted. The majority of variables are nominal; that is, the only information that can be obtained by comparing two values is their equality or non-equality. Examples are sex (male or female), primary disease (lung cancer, osteoporosis, etc.) or smoking (yes or no). Also, the SNP genotypes are nominal, each value corresponding to a combination of alleles. For example, for SIRT1 rs932658, 1 means CC, 2 means AC and 3 means AA in position. Furthermore, the data are incomplete; 138 data points (8.5%) are unknown due to various factors (lack of measurement, etc.). In such cases, to reveal data structure with patients or variables as units of comparison, special methods are called for. Objects may be compared using the Gower dissimilarity index [32], whereas correlations among variables of the mixed type may be calculated by the d-correlation formula recently suggested by Podani et al. (2023) [33]. Reducing dimensionality, thus revealing underlying data structure, is then made available through principal coordinates analysis (PCoA) of Gower dissimilarities and principal components analysis (PCA) from the correlation matrix. The results of both ordination procedures are displayed by scatter diagrams with axes as dimensions, most commonly selecting the first two. The relative importance of these dimensions is expressed by percentage contributions of eigenvalues to the total variance. The larger the sum of the eigenvalues for the axes shown, the more faithful the two-dimensional representation with respect to the entire data space. Calculations were performed by the DCORR application [33] and the SYN-TAX package [34].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25073646/s1.

Author Contributions

Conceptualization, P.L. and J.P.K.; methodology, P.L., J.P.K., H.P. and A.I.; software, J.P.; formal analysis, B.B. (Bence Bojtor), R.A., B.L., A.K. and B.T.; investigation, M.V., S.S., B.B. (Bence Bojtor) and B.B. (Bernadett Balla); resources, M.V., Z.P., J.P.K. and P.L.; data curation, B.B. (Bence Bojtor) and J.P.; writing—original draft preparation, B.B. (Bence Bojtor); writing—review and editing, B.B. (Bence Bojtor), B.B. (Bernadett Balla), P.L., J.P.K. and J.P.; visualization, B.B. (Bence Bojtor) and J.P.; supervision, P.L. and J.P.K.; funding acquisition, P.L. and I.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the Ministry of Innovation and Technology (2020-4.1.1.-TKP2020-MOLORKIV) as well as the Hungarian Research Network (SE-ELKH ENDOMOLPAT).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Semmelweis University (protocol code: 10862-1/2016/EKU; IV/4613-4/2020/EKU; 3 March 2016, 3 July 2020).

Informed Consent Statement

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

Data Availability Statement

Original data generated and analyzed during this study are included in this published article or in the data repositories listed in References.

Conflicts of Interest

The authors declare no conflicts 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. Miksad, R.A.; Lai, K.C.; Dodson, T.B.; Woo, S.B.; Treister, N.S.; Akinyemi, O.; Bihrle, M.; Maytal, G.; August, M.; Gazelle, G.S.; et al. Quality of life implications of bisphosphonate-associated osteonecrosis of the jaw. Oncologist 2011, 16, 121–132. [Google Scholar] [CrossRef]
  2. Ruggiero, S.L.; Dodson, T.B.; Fantasia, J.; Goodday, R.; Aghaloo, T.; Mehrotra, B.; O’Ryan, F.; American Association of Oral and Maxillofacial Surgeons. 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]
  3. Cadieux, B.; Coleman, R.; Jafarinasabian, P.; Lipton, A.; Orlowski, R.Z.; Saad, F.; Scagliotti, G.V.; Shimizu, K.; Stopeck, A. Experience with denosumab (XGEVA(R)) for prevention of skeletal-related events in the 10 years after approval. J. Bone Oncol. 2022, 33, 100416. [Google Scholar] [CrossRef] [PubMed]
  4. Ross, J.R.; Saunders, Y.; Edmonds, P.M.; Patel, S.; Broadley, K.E.; Johnston, S.R. Systematic review of role of bisphosphonates on skeletal morbidity in metastatic cancer. BMJ 2003, 327, 469. [Google Scholar] [CrossRef]
  5. De Cicco, D.; Boschetti, C.E.; Santagata, M.; Colella, G.; Stagliano, S.; Gaggl, A.; Bottini, G.B.; Vitagliano, R.; D’Amato, S. Medication-Related Osteonecrosis of the Jaws: A Comparison of SICMF-SIPMO and AAOMS Guidelines. Diagnostics 2023, 13, 2137. [Google Scholar] [CrossRef]
  6. Kaczoruk-Wieremczuk, M.; Adamska, P.; Adamski, L.J.; Wychowanski, P.; Jereczek-Fossa, B.A.; Starzynska, A. Oral Surgery Procedures in a Patient with Hajdu-Cheney Syndrome Treated with Denosumab—A Rare Case Report. Int. J. Environ. Res. Public Health 2021, 18, 9099. [Google Scholar] [CrossRef] [PubMed]
  7. Ng, T.L.; Tu, M.M.; Ibrahim, M.F.K.; Basulaiman, B.; McGee, S.F.; Srikanthan, A.; Fernandes, R.; Vandermeer, L.; Stober, C.; Sienkiewicz, M.; et al. Long-term impact of bone-modifying agents for the treatment of bone metastases: A systematic review. Support. Care Cancer 2021, 29, 925–943. [Google Scholar] [CrossRef]
  8. Gnant, M.; Mlineritsch, B.; Stoeger, H.; Luschin-Ebengreuth, G.; Knauer, M.; Moik, M.; Jakesz, R.; Seifert, M.; Taucher, S.; Bjelic-Radisic, V.; et al. Zoledronic acid combined with adjuvant endocrine therapy of tamoxifen versus anastrozol plus ovarian function suppression in premenopausal early breast cancer: Final analysis of the Austrian Breast and Colorectal Cancer Study Group Trial 12. Ann. Oncol. 2015, 26, 313–320. [Google Scholar] [CrossRef] [PubMed]
  9. Raje, N.; Terpos, E.; Willenbacher, W.; Shimizu, K.; Garcia-Sanz, R.; Durie, B.; Legiec, W.; Krejci, M.; Laribi, K.; Zhu, L.; et al. Denosumab versus zoledronic acid in bone disease treatment of newly diagnosed multiple myeloma: An international, double-blind, double-dummy, randomised, controlled, phase 3 study. Lancet Oncol. 2018, 19, 370–381. [Google Scholar] [CrossRef]
  10. Saag, K.G.; Petersen, J.; Brandi, M.L.; Karaplis, A.C.; Lorentzon, M.; Thomas, T.; Maddox, J.; Fan, M.; Meisner, P.D.; Grauer, A. Romosozumab or Alendronate for Fracture Prevention in Women with Osteoporosis. N. Engl. J. Med. 2017, 377, 1417–1427. [Google Scholar] [CrossRef]
  11. Hallmer, F.; Andersson, G.; Gotrick, B.; Warfvinge, G.; Anderud, J.; Bjornland, T. Prevalence, initiating factor, and treatment outcome of medication-related osteonecrosis of the jaw-a 4-year prospective study. Oral. Surg. Oral. Med. Oral. Pathol. Oral. Radiol. 2018, 126, 477–485. [Google Scholar] [CrossRef]
  12. Bone, H.G.; Wagman, R.B.; Brandi, M.L.; Brown, J.P.; Chapurlat, R.; Cummings, S.R.; Czerwinski, E.; Fahrleitner-Pammer, A.; Kendler, D.L.; Lippuner, K.; et al. 10 years of denosumab treatment in postmenopausal women with osteoporosis: Results from the phase 3 randomised FREEDOM trial and open-label extension. Lancet Diabetes Endocrinol. 2017, 5, 513–523. [Google Scholar] [CrossRef]
  13. Soundia, A.; Hadaya, D.; Esfandi, N.; de Molon, R.S.; Bezouglaia, O.; Dry, S.M.; Pirih, F.Q.; Aghaloo, T.; Tetradis, S. Osteonecrosis of the jaws (ONJ) in mice after extraction of teeth with periradicular disease. Bone 2016, 90, 133–141. [Google Scholar] [CrossRef]
  14. Soma, T.; Iwasaki, R.; Sato, Y.; Kobayashi, T.; Nakamura, S.; Kaneko, Y.; Ito, E.; Okada, H.; Watanabe, H.; Miyamoto, K.; et al. Tooth extraction in mice administered zoledronate increases inflammatory cytokine levels and promotes osteonecrosis of the jaw. J. Bone Miner. Metab. 2021, 39, 372–384. [Google Scholar] [CrossRef]
  15. Mawardi, H.; Giro, G.; Kajiya, M.; Ohta, K.; Almazrooa, S.; Alshwaimi, E.; Woo, S.B.; Nishimura, I.; Kawai, T. A role of oral bacteria in bisphosphonate-induced osteonecrosis of the jaw. J. Dent. Res. 2011, 90, 1339–1345. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, G.; Hamadeh, I.S.; Katz, J.; Riva, A.; Lakatos, P.; Balla, B.; Kosa, J.; Vaszilko, M.; Pelliccioni, G.A.; Davis, N.; et al. SIRT1/HERC4 Locus Associated With Bisphosphonate-Induced Osteonecrosis of the Jaw: An Exome-Wide Association Analysis. J. Bone Miner. Res. 2018, 33, 91–98. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, G.; Collins, J.M.; Rafiee, R.; Singh, S.; Langaee, T.; McDonough, C.W.; Holliday, L.S.; Wang, D.; Lamba, J.K.; Kim, Y.S.; et al. SIRT1 Gene SNP rs932658 Is Associated With Medication-Related Osteonecrosis of the Jaw. J. Bone Miner. Res. 2021, 36, 347–356. [Google Scholar] [CrossRef] [PubMed]
  18. Balla, B.; Vaszilko, M.; Kosa, J.P.; Podani, J.; Takacs, I.; Tobias, B.; Nagy, Z.; Lazary, A.; Lakatos, P. New approach to analyze genetic and clinical data in bisphosphonate-induced osteonecrosis of the jaw. Oral. Dis. 2012, 18, 580–585. [Google Scholar] [CrossRef]
  19. Sarasquete, M.E.; Garcia-Sanz, R.; Marin, L.; Alcoceba, M.; Chillon, M.C.; Balanzategui, A.; Santamaria, C.; Rosinol, L.; de la Rubia, J.; Hernandez, M.T.; et al. Bisphosphonate-related osteonecrosis of the jaw is associated with polymorphisms of the cytochrome P450 CYP2C8 in multiple myeloma: A genome-wide single nucleotide polymorphism analysis. Blood 2008, 112, 2709–2712. [Google Scholar] [CrossRef]
  20. Yang, Y.; Liu, Y.; Wang, Y.; Chao, Y.; Zhang, J.; Jia, Y.; Tie, J.; Hu, D. Regulation of SIRT1 and Its Roles in Inflammation. Front. Immunol. 2022, 13, 831168. [Google Scholar] [CrossRef] [PubMed]
  21. Chen, C.; Zhou, M.; Ge, Y.; Wang, X. SIRT1 and aging related signaling pathways. Mech. Ageing Dev. 2020, 187, 111215. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, Y.; Zhou, F.; Liu, H.; Li, J.; Che, H.; Shen, J.; Luo, E. SIRT1, a promising regulator of bone homeostasis. Life Sci. 2021, 269, 119041. [Google Scholar] [CrossRef]
  23. Yang, G.; Singh, S.; McDonough, C.W.; Lamba, J.K.; Hamadeh, I.; Holliday, L.S.; Wang, D.; Katz, J.; Lakatos, P.A.; Balla, B.; et al. Genome-wide Association Study Identified Chromosome 8 Locus Associated with Medication-Related Osteonecrosis of the Jaw. Clin. Pharmacol. Ther. 2021, 110, 1558–1569. [Google Scholar] [CrossRef]
  24. 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]
  25. 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] [PubMed]
  26. Kim, S.; Mun, S.; Shin, W.; Han, K.; Kim, M.Y. Identification of Potentially Pathogenic Variants Associated with Recurrence in Medication-Related Osteonecrosis of the Jaw (MRONJ) Patients Using Whole-Exome Sequencing. J. Clin. Med. 2022, 11, 2145. [Google Scholar] [CrossRef] [PubMed]
  27. Stockmann, P.; Nkenke, E.; Englbrecht, M.; Schlittenbauer, T.; Wehrhan, F.; Rauh, C.; Beckmann, M.W.; Fasching, P.A.; Kreusch, T.; Mackensen, A.; et al. Major histocompatibility complex class II polymorphisms are associated with the development of anti-resorptive agent-induced osteonecrosis of the jaw. J. Craniomaxillofac Surg. 2013, 41, 71–75. [Google Scholar] [CrossRef]
  28. Li, Q.; Cheng, J.C.; Jiang, Q.; Lee, W.Y. Role of sirtuins in bone biology: Potential implications for novel therapeutic strategies for osteoporosis. Aging Cell 2021, 20, e13301. [Google Scholar] [CrossRef]
  29. McGowan, K.; McGowan, T.; Ivanovski, S. Risk factors for medication-related osteonecrosis of the jaws: A systematic review. Oral. Dis. 2018, 24, 527–536. [Google Scholar] [CrossRef]
  30. SDCEP. Oral health management of patients at risk of medication-related osteonecrosis of the jaw. Br. Dent. J. 2017, 222, 930. [Google Scholar] [CrossRef]
  31. Drudge-Coates, L.; Van den Wyngaert, T.; Schiodt, M.; van Muilekom, H.A.M.; Demonty, G.; Otto, S. Preventing, identifying, and managing medication-related osteonecrosis of the jaw: A practical guide for nurses and other allied healthcare professionals. Support. Care Cancer 2020, 28, 4019–4029. [Google Scholar] [CrossRef] [PubMed]
  32. Gower, J.C. A general coefficient of similarity and some of its properties. Biometrics 1971, 27, 857–871. [Google Scholar] [CrossRef]
  33. Podani, J.; Schmera, D.; Bagella, S. Correlating variables with different scale types: A new framework based on matrix comparisons. Methods Ecol. Evol. 2023, 14, 1049–1060. [Google Scholar] [CrossRef]
  34. Podani, J. SYN-TAX 2000. Computer programs for data analysis in ecology and systematics. In User’s Manual Scientia, Budapest; Scientia Publishing: Budapest, Hungary, 2001. [Google Scholar]
Figure 1. Schematic diagram showing the SIRT1 gene with the SNPs tested. One single-nucleotide polymorphism (SNP) is located in the intronic region (rs7896005). One SNP is located in the 5′ upstream region (rs7894483). Two SNPs are located in the promoter region of the gene (rs932658, rs3758391). The distribution of alleles in the SNP highlighted in red is significantly different in the study population compared to the average population.
Figure 1. Schematic diagram showing the SIRT1 gene with the SNPs tested. One single-nucleotide polymorphism (SNP) is located in the intronic region (rs7896005). One SNP is located in the 5′ upstream region (rs7894483). Two SNPs are located in the promoter region of the gene (rs932658, rs3758391). The distribution of alleles in the SNP highlighted in red is significantly different in the study population compared to the average population.
Ijms 25 03646 g001
Figure 2. Representative Sanger sequencing diagrams of the rs932658 SNP. All sequencing diagrams were verified by its reverse base pair. (A) Sanger sequencing diagram of a homozygous C genotype. (B) Sanger sequencing diagram of a heterozygous CA genotype. (C) Sanger sequencing diagram of a homozygous A genotype.
Figure 2. Representative Sanger sequencing diagrams of the rs932658 SNP. All sequencing diagrams were verified by its reverse base pair. (A) Sanger sequencing diagram of a homozygous C genotype. (B) Sanger sequencing diagram of a heterozygous CA genotype. (C) Sanger sequencing diagram of a homozygous A genotype.
Ijms 25 03646 g002
Figure 3. Bar chart showing the minor allele frequencies of the SNPs tested. The lighter color represents the Medication-Related Osteonecrosis of the Jaw (MRONJ) patients studied; the darker color represents the average population’s value based on the NCBI Allele Frequency Aggregator (ALFA) database. (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA507278, accessed on 1 January 2024) A significant difference between MRONJ and the control group was detected in rs932658. *** p = 4.5 × 10−5. ns = not significant.
Figure 3. Bar chart showing the minor allele frequencies of the SNPs tested. The lighter color represents the Medication-Related Osteonecrosis of the Jaw (MRONJ) patients studied; the darker color represents the average population’s value based on the NCBI Allele Frequency Aggregator (ALFA) database. (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA507278, accessed on 1 January 2024) A significant difference between MRONJ and the control group was detected in rs932658. *** p = 4.5 × 10−5. ns = not significant.
Ijms 25 03646 g003
Figure 4. Principal components ordination of 26 variables used in the present study. Eigenvalues associated with axes 1 and 2 account for 12.8% and 8.5% of total variance, respectively. Necrosis elsewhere = osteonecrosis at a different site compared with the anatomical location of the primarily diagnosed MRONJ; Perprim/persec = per primam or per secundam wound healing after surgical treatment of MRONJ; Healing = complete healing after treatment (yes or no); Max/mand = maxillary or mandibular localization of MRONJ; Area = frontal, premolar or molar localization of MRONJ; Tooth extraction = prior tooth extraction in medical history (yes or no); Hormone therapy = prior hormone therapy in medical history (yes or no); Chemotherapy = prior chemotherapy in medical history (yes or no).
Figure 4. Principal components ordination of 26 variables used in the present study. Eigenvalues associated with axes 1 and 2 account for 12.8% and 8.5% of total variance, respectively. Necrosis elsewhere = osteonecrosis at a different site compared with the anatomical location of the primarily diagnosed MRONJ; Perprim/persec = per primam or per secundam wound healing after surgical treatment of MRONJ; Healing = complete healing after treatment (yes or no); Max/mand = maxillary or mandibular localization of MRONJ; Area = frontal, premolar or molar localization of MRONJ; Tooth extraction = prior tooth extraction in medical history (yes or no); Hormone therapy = prior hormone therapy in medical history (yes or no); Chemotherapy = prior chemotherapy in medical history (yes or no).
Ijms 25 03646 g004
Figure 5. Principal coordinates ordination of 60 patients (full symbols). Eigenvalues associated with axes 1 and 2 account for 15% and 11% of total variance, respectively. Blue arrows and labels refer to variables with the highest correlations with either axis. Note that variable–axis correlation values were scaled down for best fit to the ordination, so that their relative positions and relative lengths matter only.
Figure 5. Principal coordinates ordination of 60 patients (full symbols). Eigenvalues associated with axes 1 and 2 account for 15% and 11% of total variance, respectively. Blue arrows and labels refer to variables with the highest correlations with either axis. Note that variable–axis correlation values were scaled down for best fit to the ordination, so that their relative positions and relative lengths matter only.
Ijms 25 03646 g005
Table 1. Summary of the clinical characteristics of the study population.
Table 1. Summary of the clinical characteristics of the study population.
CharacteristicsMRONJ Patients n = 63
Age (years) (mean ± SD) *69.33 ± 11.06
Sex
Female46
Male17
Duration of treatment (months) (mean ± SD) *46.42 ± 39.77
Disorder *
Breast cancer26
Prostate cancer11
Myeloma multiplex11
Osteoporosis7
Lung cancer2
Melanoma malignum1
Colon cancer1
Cervical cancer1
Chemotherapy *
Yes38
No22
Hormone deprivation therapy
Yes32
No28
Antiresorptive agent *, **
Alendronate (p.o.)1
Ibandronate (p.o.)4
Risedronate (p.o.)1
Zoledronate (iv./inj.)29
Denosumab (inj.)29
MRONJ localization *, ***
Maxilla19
Mandibula42
MRONJ stage *, ****
11
244
315
* There were no available clinical data for three patients. ** Four patients had both bisphosphonate and denosumab. *** One patient had both maxillar and mandibular involvement. **** Based on the American Association of Oral and Maxillofacial Surgeons’(AAOMS) Position Paper on Medication-Related Osteonecrosis of the Jaws—2022 update.
Table 2. Summary of the genotypes and allele frequencies. NCBI ALFA: National Center for Biotechnology Information Allele Frequency Aggregator for European population. p-value is for the statistical difference between the study and European populations.
Table 2. Summary of the genotypes and allele frequencies. NCBI ALFA: National Center for Biotechnology Information Allele Frequency Aggregator for European population. p-value is for the statistical difference between the study and European populations.
SNPGenotype 1Genotype 2Genotype 3Allele FrequencyALFAp-Value
rs932658CC: 65.1% (n = 41)CA: 25.4% (n = 16)AA: 9.5% (n = 6)C:A 77.8:22.2C:A 59.9:40.1p = 4.5 × 10−5
rs7896005 *GG: 47.4% (n = 27)GA: 40.3% (n = 23)AA: 12.3% (n = 7)G:A 67.5:32.5G:A 65.5:34.5ns. (p = 0.64)
rs7894483 **AA: 9.9% (n = 6)AT: 39.3% (n = 24)TT: 50.8% (n = 31)A:T 29.5:70.5A:T 32.9:67.1ns. (p = 0.43)
rs3758391 **CC: 50.8% (n = 31)CT: 37.7% (n = 23)TT: 11.5% (n = 7)C:T 69.7:30.3C:T 67.7:32.3ns. (p = 0.64)
* Genotyping of six samples was unsuccessful. ** Genotyping of two samples was unsuccessful. ns. = not significant.
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

Bojtor, B.; Vaszilko, M.; Armos, R.; Tobias, B.; Podani, J.; Szentpeteri, S.; Balla, B.; Lengyel, B.; Piko, H.; Illes, A.; et al. Analysis of SIRT1 Gene SNPs and Clinical Characteristics in Medication-Related Osteonecrosis of the Jaw. Int. J. Mol. Sci. 2024, 25, 3646. https://doi.org/10.3390/ijms25073646

AMA Style

Bojtor B, Vaszilko M, Armos R, Tobias B, Podani J, Szentpeteri S, Balla B, Lengyel B, Piko H, Illes A, et al. Analysis of SIRT1 Gene SNPs and Clinical Characteristics in Medication-Related Osteonecrosis of the Jaw. International Journal of Molecular Sciences. 2024; 25(7):3646. https://doi.org/10.3390/ijms25073646

Chicago/Turabian Style

Bojtor, Bence, Mihaly Vaszilko, Richard Armos, Balint Tobias, Janos Podani, Szofia Szentpeteri, Bernadett Balla, Balazs Lengyel, Henriett Piko, Anett Illes, and et al. 2024. "Analysis of SIRT1 Gene SNPs and Clinical Characteristics in Medication-Related Osteonecrosis of the Jaw" International Journal of Molecular Sciences 25, no. 7: 3646. https://doi.org/10.3390/ijms25073646

APA Style

Bojtor, B., Vaszilko, M., Armos, R., Tobias, B., Podani, J., Szentpeteri, S., Balla, B., Lengyel, B., Piko, H., Illes, A., Kiss, A., Putz, Z., Takacs, I., Kosa, J. P., & Lakatos, P. (2024). Analysis of SIRT1 Gene SNPs and Clinical Characteristics in Medication-Related Osteonecrosis of the Jaw. International Journal of Molecular Sciences, 25(7), 3646. https://doi.org/10.3390/ijms25073646

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