Next Article in Journal
Treatment of Transthyretin Amyloid Cardiomyopathy: The Current Options, the Future, and the Challenges
Next Article in Special Issue
Effects of Type 2 Diabetes Mellitus on Osteoclast Differentiation, Activity, and Cortical Bone Formation in POSTmenopausal MRONJ Patients
Previous Article in Journal
Prediction of Complications and Surgery Duration in Primary Total Hip Arthroplasty Using Machine Learning: The Necessity of Modified Algorithms and Specific Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 11
Issue 8
10.3390/jcm11082145
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Identification of Potentially Pathogenic Variants Associated with Recurrence in Medication-Related Osteonecrosis of the Jaw (MRONJ) Patients Using Whole-Exome Sequencing

by
Songmi Kim
1,2,†,
Seyoung Mun
1,2,†,
Wonseok Shin
3,
Kyudong Han
1,2,* and
Moon-Young Kim
4,*
1
Center for Bio Medical Engineering Core Facility, Dankook University, Cheonan 31116, Korea
2
Department of Microbiology, Dankook University, Cheonan 31116, Korea
3
NGS Clinical Laboratory, Dankook University Hospital, Cheonan 31116, Korea
4
Department of Oral and Maxillofacial Surgery, College of Dentistry, Dankook University, Cheonan 31116, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
J. Clin. Med. 2022, 11(8), 2145; https://doi.org/10.3390/jcm11082145
Submission received: 16 March 2022 / Revised: 6 April 2022 / Accepted: 9 April 2022 / Published: 12 April 2022
(This article belongs to the Special Issue Breakthroughs in Oral and Maxillofacial Surgery)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: Bisphosphonates are antiresorptive and antiangiogenic drugs that prevent and treat bone loss and mineralization in women with postmenopausal osteoporosis and cancer patients. Medication-related osteonecrosis of the jaw (MRONJ) is commonly caused by tooth extraction and dental trauma. Although genetic and pathological studies about MRONJ have been conducted, the pathogenesis of MRONJ still remains unclear. Methods: We aimed to identify genetic variants associated with MRONJ, using whole-exome sequencing (WES). Ten MRONJ patients prescribed bisphosphonates were recruited for WES, and jawbone tissue and blood samples were collected from the patients. Results: The analysis of the WES data found a total of 1866 SNP and 40 InDel variants which are specific to MRONJ. The functional classification assay using Gene Ontology and pathway analysis discovered that genes bearing the MRONJ variants are significantly enriched for keratinization and calcium ion transport. Some of the variants are potential pathogenic variants (24 missense mutations and seven frameshift mutations) with MAF < 0.01. Conclusions: The variants are located in eight different genes (KRT18, MUC5AC, NBPF9, PABPC3, MST1L, ASPN, ATN1, and SLAIN1). Nine deleterious SNPs significantly associated with MRONJ were found in the KRT18 and PABPC3 genes. It suggests that KRT18 and PABPC3 could be MRONJ-related key genes.

1. Introduction

Medication-related osteonecrosis of the jaw (MRONJ) refers to an uncommon medical condition that can occur in patients receiving antiresorptive or antiangiogenic therapy to treat osteoporosis and diverse malignant conditions such as multiple myeloma, breast cancer, and prostate cancer [1]. When the first case of MRONJ was documented in 2003, only bisphosphonates were associated with the pathophysiology, and thus the condition was widely known as bisphosphonate-related osteonecrosis of the jaw (BRONJ) [2]. However, in light of the numerous recent reports of osteonecrosis cases associated with the use of various other antiresorptive and antiangiogenic agents such as denosumab, the current universal nomenclature for the disease has been modified to MRONJ. The antiangiogenic agents or targets of the VEGF pathway, such as the anti-VEGF monoclonal antibodies (Bevacizumab) and tyrosine kinase inhibitors (Sunitinib, Sorafenib), that are prescribed for the treatment of cancer-related conditions or organ rejection in renal transplant patients (Sirolimus) have also been indicated in the pathophysiology of the disease [3]. The key hypotheses of disease specificity for the jaws include bone remodeling inhibition, angiogenesis inhibition, inflammation, innate or acquired immune dysfunction, soft tissue toxicity, and genetic predisposition. Several hypotheses can explain the overall pathophysiology of MRONJ [4,5]. Antiresorptive and antiangiogenic reagents, including bisphosphonates and denosumab, inhibit bone remodeling by affecting osteoclast formation and differentiation, or they directly inhibit angiogenesis by reducing blood flow to the bone [5,6]. Genetic factors such as single nucleotide polymorphisms (SNPs) and mutations are significantly associated with MRONJ development by regulating bone remodeling, angiogenesis, and immune responses [7,8]. Prior attempts to determine the efficacy of a drug holiday in preventing MRONJ have yielded inconclusive results, and thus the evidence to support or refute the use of such a holiday remains limited [7]. The hallmark clinical feature of MRONJ is the presence of exposed necrotic bone in the oral cavity that persists for more than eight weeks without any signs of healing. This is frequently accompanied by concurrent infection signs such as pain, swelling, fistula formation, and in its severest form—pathological fractures (Figure 1). In some MRONJ patients, a specific symptom, known as the ‘numb chin syndrome’ may also manifest [8]. According to the current AAOMS guidelines, MRONJ can be classified into four distinct stages (stages 0 to III), for which the suggested methods of treatment differ. For early stages (stage I) of MRONJ patients, conservative treatment, including antibiotic and septic administration, and an antibacterial oral rinse, results in pain control and prevention of necrosis. Surgery may not be needed in the absence of disease progression, and if the patient has an adequate quality of life [4]. For patients suffering from advanced stage (stage II–III) MRONJ, invasive surgery, including jaw resection, can yield optimal mucosal healing; however, up to 25% of patients can suffer from severe postoperative morbidity following this type of extensive surgery. Therefore, to minimize morbidity by establishing novel prevention or early detection methods, recent studies have focused on determining the pathophysiology of MRONJ.
The incidence of MRONJ is relatively lower in the oral bisphosphonate group for osteoporosis than in the intravenous bisphosphonate group for cancer; hence, to date, many hypotheses concerning the pathogenesis of MRONJ have been raised by many authors. Previous studies reported that genetic polymorphisms in several candidate genes (TGFb1, MMP2, CYP2C8, VEGF, COL1A1, RANK, OPG, OPN, and PPARG) involved in the activation and differentiation of osteoclast and bone remodeling are significantly associated with MRONJ [1,9,10,11]. In addition, many researchers recently found that MRONJ is caused by multiple biological pathways involving immunoglobulin, cell adhesion, and cytoskeleton [12]. However, despite the numerous genetic studies, such as the genome-wide association studies (GWAS), whole-exome studies (WES), and the candidate gene studies (CGS), that have been performed to ascertain the condition’s exact pathogenic mechanism, the pathogenesis remains yet to be determined. Furthermore, although potential associations between the development of MRONJ and several genes have been documented in the previous literature, the GWAS, WES, and CGS have all failed to elucidate a single gene as the sole risk factor for MRONJ development.
Our study aimed to identify genetic variants associated with MRONJ by analyzing somatic and germline mutations through WES. By conducting an integrated analysis with the preceding research results, eight genes were determined to be closely related to MRONJ, and mutations in the genes were investigated in the present study. The result found that twelve and ten genetic mutations of KRT18 and PABPC3 genes, respectively, are predicted to modify the protein structure of the genes in MRONJ patients.

2. Materials and Methods

2.1. Ethics Statement

Experiments and analyses involving human subjects or human data in this study were approved by the Institutional Review Board of Dankook University Hospital (IRB No. 2018-01-009). All clinical investigations were performed in accordance with the Declaration of Helsinki. Written informed consent was obtained from all patients included in this study.

2.2. Patient Selections

Ten MRONJ patients (stage II–IIII) from the oral and maxillofacial surgery department of Dankook University Hospital participated in this study. The patients visited the Dankook University Hospital for surgical intervention from March 2018 to December 2018 and received bisphosphonates (BPs) therapy, either risedronate (9 cases) or alendronate (1 case). Patients were considered to have MRONJ if all of the following three characteristics were present according to the position paper issued by the American Association of Oral and Maxillofacial Surgeons (AAOMS) in 2014; (i) current, previous treatment with an antiresorptive or antiangiogenic therapy, (ii) non-healing, exposed bone in the maxillofacial region that had persisted for more than eight weeks, and (iii) no history of radiation therapy to the jaw. We also analyzed the characteristics of the clinical findings (medications, smoking, sex, stage) in a panoramic radiographic study (Figure 1).

2.3. DNA Preparation

Peripheral venous blood (2 mL) and routine jawbone specimens were obtained from the ten MRONJ patients in the Department of Oral and Maxillofacial Surgery of Dankook University Hospital. Genomic DNA was extracted from the blood and the lesion jawbone tissues immediately after surgical treatment of ten MRONJ patients using a DNeasy Blood & Tissue Kit (Qiagen, Germany) in accordance with the manufacturer’s protocols, and stored at −80 °C until use.

2.4. Whole-Exome Sequencing (WES) and Individual Variant Calling

For NGS library construction, each genomic DNA (1 μg) was sheared to 180–280 bp fragments using a DNA shearing system (Covaris, Woburn, MA, USA). Sequencing libraries were generated using a SureSelect Human All Exon Kit (Agilent Technologies, Santa Clara, CA, USA) following the manufacturer’s instructions. DNA that was ligated with adapters were enriched during PCR reaction. The WES libraries were captured and purified using streptavidin-coated magnetic beads and the AMPure XP system (Beckman Coulter, Beverly, MA, USA). WES was performed on the Illumina NovaSeq 6000 with a 2 × 150 bp paired end (PE) read.
Sequencing reads were mapped to the human reference genome (GRCh37/hg19) using the Burrows–Wheeler Alignment tool (BWA 0.7.8) with default parameters [13]. PCR duplicates were removed using Picard (http://picard.sourceforge.net/index.shtml, accessed on 8 November 2019). Variant calling was performed using the Genome Analysis Tool Kit (GATK v3.8) [14]. The functional annotation of variants was conducted using ANNOVAR (2015Dec14) [15]. To obtain genetic variants highly related to MRONJ, we compared two datasets with the variant candidates (Korean variant archive; KOVA and GSK project). KOVA was used as the control dataset to exclude the common Korean variants. The GSK project investigating the MRONJ-related variants found in the saliva of MRONJ patients was used to collect variants shared between MRONJ diagnosed patients.

2.5. Functional Gene Classification and Pathway Analysis

To investigate the biological relevance of the candidate genes, we performed gene ontology (GO) enrichment analysis to classify them into three classes, including biological process (BP), cellular component (CC), and molecular function (MF), by using Metascape software (https://metascape.org/gp/index.html, accessed on 19 February 2020) [16]. Then, the enrichment of candidate genes was conducted to identify pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genom e.jp/kegg/pathw ay.html, accessed on 19 February 2020). The significant function of genes was filtered by statistical cut-off using the adjusted p-value (FDR; False Discovery Rate).

2.6. Pathogenicity of Variants

To investigate the occurrence and the frequency of the variants in the general population, three public databases were used as follows: ExAc (n = 60,706, http://exac.broadinstitute.org/, accessed on 14 April 2021), 1000 Genomes project phase 3 database (1000 G; n = 2504, http://www.internationalgenome.org/, accessed on 14 April 2021), and NHLBI Exome Sequencing Project (n = 6503, http://evs.gs.washington.edu/EVS/, accessed on 14 April 2021). Using the comparison of the three control databases, we investigated the frequency of the genetic variations identified in MRONJ patients and selected unique variants with minor allele frequencies (MAF) < 0.01. In addition, the pathogenicity of the variants was analyzed to predict functional effects sorting intolerant from tolerant (SIFT) (http://sift.jcvi.org, accessed on 14 April 2021), using Polymorphism Phenotyping v2 (PolyPhen-2) (http://genetics.bwh.harvard.edu/pph2/, accessed on 14 April 2021), a Mutation Assessor, and a Mutation Taster.
Furthermore, we predicted protein structures using the AlphaFold protein structure database (https://alphafold.ebi.ac.uk/, accessed on 10 August 2021), Phyre2 (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index, accessed on 10 August 2021). In addition, we used UCSF chimera (https://www.cgl.ucsf.edu/chimera/, accessed on 10 August 2021) and DynaMut2 (http://biosig.unimelb.edu.au/dynamut2/, accessed on 10 August 2021) to analyze the protein stability and protein conformational changes [17,18]. The filtration and prioritization framework used for data analysis is described in Supplementary Figure S1.

3. Results

3.1. Clinical Findings

The clinical characteristics of the ten MRONJ patients are shown in Table 1. One man and nine women were enrolled in this study, and none of them were smokers. The mean age was 77.8 ± 5.64 years. Of the ten MRONJ patients, the mandible bone was involved in seven (70%), and the rest were maxillary bones. One of the ten MRONJ patients was spontaneous, eight had teeth extractions, and the other had root canal therapy prior to the onset of MRONJ. Only one patient was taking alendronate, and the others were taking risedronate (Table 1). Ten patients developed the disease after long-term administration for at least four years.

3.2. Whole-Exome Sequencing and Quality Controls

WES data were generated with an average 11.4 Gb of sequences per sample on the Illumina NovaSeq6000 platform. Sequencing quality for Q20 and Q30 was an average of 96.75% and 91.99%, respectively (Supplementary Table S1). Of the WES reads, 97.34% were mapped to the human reference genome. The blood and lesion samples (i.e., the lesion jawbone tissues) were sequenced to an average of 104.16× depth on target, and an average of 99.57% of exonic regions were covered with at least 10× depth on target (Supplementary Table S1).

3.3. Identification of Single Nucleotide Polymorphism (SNPs) Related to MRONJ

From WES data, an average of 433,602 germline SNP variants and 367,806 somatic SNP variants were identified in ten patients (Supplementary Table S2). To identify variants associated with MRONJ, we used a strict filtering method that collects common variants from all ten patients. We identified a total of 39,720 germline SNP variants and 32,225 somatic SNP variants shared by all ten patients. Then, to rule out the common variants in Korean, we compared the total variants with the KOVA dataset containing 293,049 variants in 1055 healthy Korean individuals [19]. The 219,772 MRONJ-associated variants discovered from the GSK project with 142 saliva samples of MRONJ patients [20] were used for variant recalibration and separation of somatic mutations from germline mutations (Figure 2). Through the analysis, we identified a total of 4171 SNPs in 2321 genes associated with MRONJ, including 3876 SNPs detected in all three datasets (saliva, blood, and lesion), 216 SNPs detected in both blood and saliva datasets, and 79 SNPs detected in both saliva and lesion datasets. Additionally, we collected 1151 lesion-specific SNPs to identify somatic variants (Table 2, Figure 2). We found 1020 nonsynonymous SNP variants, including the missense (n = 946), nonsense (n = 8), and unknown mutations (n = 66) in the exonic and splicing region (Table 3, Supplementary Table S3).
In the MRONJ common group, we found 942 nonsynonymous variants, including missense, nonsense, and unknown variants, and ten frameshift InDels in the MRONJ common group (saliva, blood, and lesion). The top ten genes that accumulated a large number of nonsynonymous variants were PDE4DIP, HLA-C, HLA-B, GPRIN2, KCNJ12, KCNJ18, HLA-A, and MST1L genes (Table 3, Supplementary Table S3). Twenty-three of 942 nonsynonymous variants were distributed in PDE4DIP, which encodes anchoring the primary c-AMP-hydrolyzing enzyme and the phosphodiesterase 4D (PDE4D) protein, to the Gogi/centrosome region of the cell. Nine nonsynonymous variants were identified in KCNJ18 or KCNJ12, which regulate potassium concentration in the skeletal muscle. KCNJ12 is an inwardly rectifying K+ channel activated by phosphatidylinositol 4,5-bisphosphate that controls the concentration of calcium ions in the cytosol of skeletal muscle [21]. KCNJ18 encodes the skeletal muscle inward-rectifier K+ channel Kir2.6, which is regulated by thyroid hormones (T3) [22,23]. Three variants (R39Q, R40H, I249V) in the KCNJ18 gene were previously reported, and they are related to thyrotoxic periodic paralysis (TPP) [22,24].
We identified nine lesion-specific missense variants in eight genes, including PRAMEF1 (c.386 C>T, p.T129M and c.412T>C, p.C138R), IGFN1 (c. 5047G>A, p.G1683R), OR2T12 (c. 205A>G, p.M69V), PABPC3 (c.878T>G, p.V293G), PLIN4 (c. 787G>A, p.G263S), LONRF2 (c.330C>G, p.D110E), MUC12 (c.2564G>A, p.R855H), FOXD4 (c.614C>T, p.P205L) (Table 3; Supplementary Table S3). PRAMEF1, a member of PRAME, is related to reproductive tissues during development. PABPC3, a cytoplasmic PABP, is associated with mRNA translation and stability that participates in nonsense-mediated decay. MUC12, a member of the mucin family, encodes an integral membrane glycoprotein that contributes to host cell defense infected by the bacterium and is associated with epithelial cell protection and regulation of cell growth. Ladeberg et al. (2008) reported that bisphosphonate inhibits the proliferation of keratinocytes in the oral mucosa [25]. Although no significance has been identified in the pathogenic effects and the phenotypic difference of these genes in MRONJ, they may affect biological functions such as cell morphology, proliferation, and differentiation, broken by protein structural alteration in protein–protein interactions (PPIs).

3.4. Identification of Insertion/Deletions (InDels) Variants Related to MRONJ

As for InDel variants, we identified an average of 72,754 germline InDels and 72,715 somatic InDels from WES data of ten patients (Supplementary Table S2). A total of 4901 germline InDels and 3892 somatic InDels are common in all ten MRONJ patients (Figure 2). Finally, a total of 1641 InDels in 1142 genes were only discovered in MRONJ patients, except for common Korean variants as we described above, which include 265 InDels detected in all three datasets (saliva, blood, and lesion), 241 InDels detected in both saliva and lesion dataset, and 46 variants detected in both saliva and blood dataset. Additionally, we included 1089 lesion-specific InDels, which were somatic variants. Of these, we collected 132 exonic InDel variants to identify variants that affect MRONJ, including common InDels in the saliva and blood (n = 4), common InDels in the saliva and lesion (n = 42), common InDels in the saliva, blood, and lesion (n = 40), and lesion-specific InDels (n = 46). They contain 32 frameshift InDels, 53 non-frameshift InDels, 46 unknown mutations, and one nonsense mutation (Table 3, Supplementary Table S4).
In MRONJ common group, ten frameshift InDels are identified in nine genes, including MST1L (c.894_897del, p.G298fs and c.119_123del, p.Q40fs), HRNR (c.1delA, p.M1fs), NUDT18 (c.111delG, p.R37fs), GPATCH4 (c.1083_1084insGT, p.F362fs), ATG3 (c.920dupT, p.L307fs), CYFIP2 (c.280dupC, p.I93fs), NOP16 (c.586_587insAC, p.R196fs), SSPO (c.9204dupC, p.S3068fs), and OR2T35 (c.609_615del, p.C203fs) (Table 3). HRNR is a member of the S-100 fused protein family, involved in protein phosphorylation regulation, cytoskeleton component, calcium homeostasis, and cell proliferation [26]. ATG3, a ubiquitin-like-conjugating enzyme, is involved in autophagy-related to osteoblast differentiation and mineralization, which plays an important role in bone homeostasis [27]. The representative genes harbored severe frameshift InDels loading that can impair these roles in genes that are vital for bone formation and control.
Eleven lesion-specific frameshift InDels were also identified in eight genes, which were SARM1 (c.549_550del, p.G183fs; c.544_545insGC, p.S182fs), SLAIN1 (c.219_220insGG, p.A73fs; c.229dupC, p.Q76fs), ZFPM (c.1330_1331del, p.E444fs; c.1335delT, p.P445fs), CTU2 (c.1097_1097del, p.R366fs), GOLGA6L2 (c.2080delA, p.R694fs), LOC283710 (c.75_76del, p.P25fs), PRDM15 (c.262dupC, p.R88fs), and SCARF2 (c.2247dupG, p.P750fs). SARM1 is highly associated with the neuron degeneration that mediates neuronal cell death. A loss-of-function (LOF) mutation (C844Y) in PRDM15, which regulates NOTCH and WNT signaling pathways, was identified in patients with brain malformations [28]. SCARF2, the scavenger receptor type F family, containing epithermal growth factor (EGF)-like domains and putative N-glycosylation sites contributes to lipid metabolism and intracellular signaling pathways in pathophysiological conditions. Two genetic variants (c.773G>A in exon 4 and c.1328_1329delTG in exon 8) of SCARF2, influencing cytoskeletal organization, were identified in patients with Van Den Ende-Gupta syndrome [29,30]. Interestingly, SCARF2 is expressed in the mandibular, maxillary, and urogenital ridge tissues in mice at ten embryonic days and human endothelial cells [30].
We also identified 33 lesion-specific non-frameshift InDels in 26 genes, which were GOLGA6L2 (n = 3), NCOR2 (n = 3), ASPN (n = 2), ATN1 (n = 2), MPRIP (n = 2), GOLGA8R (n = 1), and ZFPM1 (n = 1) (Table 3; Supplementary Table S4). NOCR2, a member of thyroid hormone and retinoic acid receptor-associated corepressors, regulates androgen receptor signaling and androgen-induced cell proliferation. ASPN is related to the regulation of periodontal ligament differentiation and mineralization that maintains homeostasis of the tooth-supporting system [31]. Asporin protein binds to collagen fiber, resulting in induction of mineralization and inhibition of TGF-b in the cartilage differentiation [32]. ATN1 is related to Dentatorubral pallidoluysian atrophy (DRPLA) disease, a rare neurodegenerative disorder, which occurs with the expansion from 7–35 copies to 49–93 copies of a trinucleotide repeat (AAG/CAA).

3.5. Enrichment and Pathway Analysis

To understand the functional relevance of the isolated genetic variants focusing on MRONJ, we performed Gene Ontology (GO) and Kegg pathway analyses on 665 genes containing 1020 nonsynonymous SNP variants. As shown in Figure 3, the most significant gene set enrichment in the biological processes (BP) was ‘detection of stimulus’, ‘keratinization’, and ‘calcium ion transport’. Interestingly, the ‘protein kinase C signaling pathway’ was also elected in the biological processes (Figure 3; Supplementary Table S5). As referred to in the previous study, the bisphosphonates inhibit protein kinase C signaling, resulting in the accumulation of calcium and degradation of osteoclast [33]. The most significant enriched function was ‘ECM structural constituent’ for the molecular function (MF), followed by ‘olfactory receptor activity’, and ‘calcium ion binding’ (Figure 3; Supplementary Table S5). For the cellular component (CC), the function of genes associated with ‘keratin filament’ was most significantly enriched (Figure 3; Supplementary Table S5). For the Kegg pathway, the top-ranked significant pathways were ‘olfactory transduction’, ‘autoimmune thyroid disease’, and ‘ECM-receptor interaction’ (Supplementary Table S5). One of the primary mechanisms, the extracellular matrix (ECM), comprises over 90% of type I collagen, the most abundant collagen, in human bone. It is a quite meaningful result because Type I collagen is synthesized by osteoblasts, which drives bone tissue remodeling during wound healing.
Furthermore, we examined the interaction between the classified functions of MRONJ-specific variants by function network analysis using ClueGO (p-value < 0.05; Benjamini-Hochberg). The integrative function network analysis showed that the eight parent functions, including ‘integral component of luminal side of endoplasmic reticulum membrane’, ‘O-glycan processing’, ‘sensory perception of the chemical stimulus’, ‘keratin filament’, ‘calcium ion transport’, ‘protein kinase A signaling’, ‘phenylpropanoid metabolic process’, and ‘post-synaptic specialization organization’ were mainly clustered with 665 genes (Supplementary Figure S2).

3.6. Pathogenicity Analysis

For the private or rare variants (MAF < 0.01), we also utilized three databases, including the 1000 genomes project database, NHLBI exome sequencing project (ESP), and ExAC database. We filtered in the candidate variants with allele frequencies lower than 0.01 or no allele frequency in all three databases. A total of 24 candidate SNPs and seven candidate InDels were selected in four genes (KRT18, MUC5AC, NBPF9, and PABPC3) and five genes (MST1L, ASPN, ATN1, PABPC3, and SLAIN1), respectively (Table 4).
In addition, we analyzed 24 rare SNP variants for the annotation of potentially deleterious effects using SIFT and polyphen2. Of these, three SNP variants (p.G38C, p.G43R in KRT18, and p.K231E in PABPC3) were predicted to be ‘deleterious’ in the SIFT, and ‘possibly damaging or probably damaging’ in Polyphen-2 (Table 4). Six SNP variants (p.R27W, p.P28Q, p.A32S, p.S34R in KRT18, and p.A181T, p.L207F in PABPC3) were predicted to be either deleterious or benign in the SIFT and Polyphen-2 (Table 4). KRT18 encoding the type 1 intermediate filament chain regulates protein synthesis, organelle positioning, and protection from apoptosis and necrosis. PABPC3 is a member of RNA binding proteins involved in the cytoplasmic regulatory process of mRNA metabolisms, such as RNA stability and translation initiation. RBMS3 (the c-myc gene single-strand binding protein), a paralog for PABPC3, regulates myc/ras cooperative transforming activity, apoptosis, DNA replication, and cytoplasmic RNA metabolism. GWAS studies showed that two SNPs (rs17024608 and rs10510628) in the RBMS3 were associated with a bisphosphonates-induced reduction in collagen formation, osteonecrosis disruption, and bone turnover [34,35].
To predict protein structural alteration by amino acid (AA) changes, protein conformational changes in WT and mutants of KRT18 and PABPC3 were predicted using the AlphaFold protein structure database, phyre2 web tool. Nine rare SNP variants (p.R27W, p.P28Q, p.A32S, p.S34R, p.G38C, and p.G43R in KRT18, and p.A181T, p.L207F, and p.K231E in PABPC3) were structurally changed (Figure 4A). In addition, we calculated the distance between the atoms of two residues to predict the effect of the mutations on the distances between each pair of residues using the UCSF chimera tool. The distances between wild-type and mutations were changed in KRT18 and PABPC3 (Figure 4B). Notably, we identified pairs of charged residues that were significantly skewed from 8.68 Å to 4.67 Å by the R27W mutation (Figure 4B).
Furthermore, we used DynaMut2 to analyze the entropy energy change between wild-type and mutant structures (in kcal/mol/K). All six substitutions (p.R27W, p.P28Q, p.A32S, p.S34R, p.G38C, and p.G43R), located in the head domain, which is a primary phosphorylation site of KRT18, have weak effects (−0.5 < ΔΔG < 0.5 kcal/mol) on KRT18 stability (Figure 5). The phosphorylation of KRT18 at serine 33 and 34 regulate filament stability, keratin reorganization, binding to 14-3-3 protein during mitosis, and keratin protein turnover during apoptosis [36,37,38]. The glycosylation at serine 30 and 31 protects epithelial tissue from injury [37]. The loss and mutations of KRT18 contribute to the phenotypes of several human diseases, including skin disorders, chronic renal disease, and liver cirrhosis [38].
In addition, three substitutions (A181T, L207F, and K231E) are located in the RNA recognition motifs (RRMs) of PABPC3 that bind the poly(A)tail of single-stranded RNAs (Figure 4A). Unfortunately, the entropy change of PABPC3 could not be analyzed because the intact protein structure of PABPC3 could not be obtained from the protein data bank (PDB). However, a p.K231E substitution located in the RNP1 motif ((K/R)-G-(F/Y)-(G/A)-F-V-X-(F/Y)), a highly conserved region, may interact with the translation initiation factor and regulatory proteins, affecting translation initiation and nonsense-mediated decay.

4. Discussion

Bisphosphonates, known as a primary drug for the onset of MRONJ, is the antiresorptive agent that prevents loss of bone density. They are used to treat osteoporosis, osteogenesis imperfecta, primary hyperparathyroidism, and cancer. The bisphosphonates bind preferentially to calcium ions in the bones and accumulate in high concentrations. However, the bisphosphonates cause mucosal ulceration and exposure to the necrotic bone, leading to osteonecrosis of the jaw. Previous studies reported that the bisphosphonates induce senescence and apoptosis of the cells by blocking the cholesterol biosynthetic pathways (mevalonate pathway) that indirectly inhibit the blood supply of the bone [20,39,40]. Bisphosphonates, which are closely related to anti-angiogenesis, affect the structure of soft tissues in the jaw, causing soft tissue cell damage. In this study, we identified the candidate genes and sorted out the variants closely associated with MRONJ in lesion and blood samples through WES. Although previous studies have constantly attempted to screen the candidate genes and variants that show a significant correlation with MRONJ using GWAS, WES, and CGS, the single genes and variants associated with MRONJ development have not been identified; therefore, the need for in-depth research on the onset, development, diagnosis, and treatment of MRONJ is drawing attention by securing various genetic data according to various races, drugs, living environments, and clinical aspects.
In function enrichment analysis for MRONJ-specific variants, the protein kinase C (PKC) signaling pathway, that can either read out lipid signals alone or combine the ability to read out simultaneous calcium ions, is highly enriched in the biological processes. The calcium ion is a second messenger that regulates differentiation and activation of osteoclasts according to calcium concentration [41]. The calcium ion signaling regulates the activation of PKC and NFATc1 in osteoclasts, which is involved in cell proliferation and differentiation [33]. Previous studies reported that treatments of bisphosphonates, such as alendronate and risedronate, were accumulated in the bone matrix. They trigger osteoclast apoptosis, resulting in cytoskeletal regulation [42,43]. Human parathyroid hormone (PTH) is highly associated with the regulation of calcium metabolism; thus, when calcium ions increase in the serum, the osteoclast activates, which degrades bone. These results suggest that the calcium ion is closely related to MRONJ development, regulating cell proliferation and differentiation. Through our study, we were able to observe variations in the following regulatory factors: FLT4, MC1R, ULK4, ADGRV1, TTN, MYOM1, and AIP, closely related to PKC, which respond immediately to calcium signals in osteoclast activity.
Through the systematic pathogenicity analysis, we identified 24 candidate SNPs and seven candidate InDels in four genes (KRT18, MUC5AC, NBPF9, and PABPC3) and five genes (MSTIL, ASPN, ATN1, PABPC3, and SLAIN1), respectively. KRTs are intermediate filament proteins found in skin and epithelial tissues that play an important role in constructing the cytoskeleton of the cells. KRTs maintain the structure of skin cells and other epithelial tissues by regulating electrolyte transport, post-translational modification, and prevention of degradations; however, disruption of KRTs leads to skin and mucosal diseases. As shown in the results of previous studies, the genetic effect of KRT4 and KRT13 mutations, a member of the keratin family, is mainly implicated in the pathogenesis of white sponge nevus (WSN) disorder, a rare autosomal dominant disorder in oral mucosa [44,45].
KRT18 is also involved in the cytoskeletal signaling pathway, estrogen signaling pathway, and staphylococcus aureus infection pathway [46,47,48,49]. Estrogen generated from the mevalonate pathway contributes to lipid metabolism, producing cholesterol, fatty acids, and steroid hormones, inducing osteoblast proliferation, osteoclast apoptosis, and stimulating matrix synthesis associated with MRONJ [50,51]. Bisphosphonates could influence the production of estrogen, a type of steroid hormone, by inhibiting the mevalonate pathway, and reduced estrogen exacerbates age-related bone loss in women [52]. In addition, a previous study reported that the microbial films of predominant bacteria, including Fusobacterium, Actinomyces, Staphylococcus, Streptococcus, Selenomonas, and Trepenemes, were present in osteonecrosis of MRONJ patients [53,54,55,56]. Disruption of KRT18 could be made to fail to interact with CDH1 (calcium-dependent cell adhesion proteins) and EPCAM (epithelial cell adhesion molecule), which are involved in the proliferation of epithelial cells and defense against mucosal infection, respectively, which affects filament reorganization; therefore, disruption of KRT18 may influence the cytoskeleton of the skin and mucosal, leading to oral mucosa diseases, including MRONJ.
Here, we accomplished the first rare-variant association study for MRONJ in lesions isolated from human jawbone tissue and blood samples. Our genetic expectation demonstrated that the KRT18 and PABPC3 could be potential candidate genes related to MRONJ and developing the MRONJ phenotype. This information could help clinicians to exclude patients who may develop MRONJ during dental procedures and take special preventive measures. Furthermore, it would be an innovative action to prevent drug side effects caused by genotypes through drug sensitivity diagnosis in patients, before the prescription of osteoporosis treatment drugs.
As a limitation of the current research, first, the number of samples in this study was not large enough to detect statistically significant genetic variants. If more robust data is accumulated from patient cohorts with similar clinical environments and conditions, it is believed that it would provide a strong basis from which to derive genetic implications. Second, examining the functional association with germline and somatic variants focusing on MRONJ might provide a crucial message to identify alterations that can be issued in the transcription or translation stage.

5. Conclusions

We performed whole-exome sequencing to identify MRONJ-associated variants. We identified 954 nonsynonymous SNP variants and 40 InDels accompanied with 691 genes, shared by MRONJ patients in this study. Of these, we found significant genetic variations between MRONJ patients and control sets, and we examined their function through analysis of GO classification and KEGG pathways. The genes bearing the variants were highly related to biological functions of the cytoskeleton, calcium-ion binding, and keratinization. Notably, nine variants located in KRT18 and PABPC3 genes were considered pathogenic variations, influencing immoral protein changes. Our finding suggests that KRT18 and PABPC3 may be implicated in MRONJ pathophysiology, and that they may be used as a therapeutic target in the treatment of osteonecrosis. It is clear that further studies using more cases and controls are needed to confirm the molecular biological association between KRT18 and PABPC3 in the development and occurrence of MRONJ, but we have demonstrated new proven genetic changes to understand the pathogenesis of MRONJ.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm11082145/s1, Figure S1: Outline of experimental design and analysis workflow for this study; Figure S2: Molecular enrichment networks analysis of genes associated with MRONJ; Table S1: Summary of sequencing statistics; Table S2: Number of variants identified from ten MRONJ patients; Table S3: List of SNP candidates; Table S4: List of InDel candidates; Table S5: Enrichment of GO terms and KEGG pathways of MRONJ-associated genes.

Author Contributions

Conceptualization, S.K. and M.-Y.K.; methodology, S.M.; validation, S.K., S.M. and W.S.; formal analysis, S.M.; investigation, W.S.; resources, M.-Y.K.; data curation, S.M.; writing—original draft preparation, S.K.; writing—review and editing, K.H.; visualization, S.K. and S.M.; supervision, K.H.; project administration, M.-Y.K.; funding acquisition, M.-Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

M.-Y.K. was supported by a National Research Foundation of Korea (NRF) grant funded by the Ministry of Science and ICT (2021R1C1C1008843), Priority Research Center Program provided by the Ministry of Education (2019R1A6A1A11034536) and Medical Research Center Program NRF (2021R1A5A2022318).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Dankook University Hospital (DKUH IRB. 2018-01-009, approved on 7 March 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

The WES data generated in this study are available from the corresponding author upon request.

Acknowledgments

The authors gratefully acknowledge the Center for Bio-Medical Engineering Core Facility at Dankook University for providing critical reagents and equipment.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AA, Amino acid; AAOMS, the American Association of Oral and Maxillofacial Surgeons; BP, Biological process; BPs, Bisphosphonates; CC, Cellular component; CGS, Candidate gene studies; ECM, Extracellular matrix; FDR, False discovery rate; GO, gene ontology; GWAS, The Genome-Wide Association Studies; MAF; Minor allele frequencies; MRONJ, Medication-related osteonecrosis of the jaws; MF, Molecular function; WES, Whole-exome studies.

References

  1. 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] [PubMed]
  2. Wang, J.; Goodger, N.M.; Pogrel, M.A. Osteonecrosis of the jaws associated with cancer chemotherapy. J. Oral Maxillofac. Surg. 2003, 61, 1104–1107. [Google Scholar] [CrossRef]
  3. Chang, J.; Hakam, A.E.; McCauley, L.K. Current Understanding of the Pathophysiology of Osteonecrosis of the Jaw. Curr. Osteoporos. Rep. 2018, 16, 584–595. [Google Scholar] [CrossRef] [PubMed]
  4. Ruggiero, S.L.; Dodson, T.B.; Aghaloo, T.; Carlson, E.R.; Ward, B.B.; Kademani, D. American Association of Oral and Maxillofacial Surgeons’ Position Paper on Medication-Related Osteonecrosis of the Jaw—2022 Update. J. Oral Maxillofac. Surg. 2022; in press. [Google Scholar] [CrossRef]
  5. He, L.; Sun, X.; Liu, Z.; Qiu, Y.; Niu, Y. Pathogenesis and multidisciplinary management of medication-related osteonecrosis of the jaw. Int. J. Oral Sci. 2020, 12, 30. [Google Scholar] [CrossRef]
  6. Beaumont, S.; Harrison, S.; Er, J. Review of Myeloma Therapies and Their Potential for Oral and Maxillofacial Side Effects. Cancers 2021, 13, 4479. [Google Scholar] [CrossRef]
  7. Otto, S.; Pautke, C.; Arens, D.; Poxleitner, P.; Eberli, U.; Nehrbass, D.; Zeiter, S.; Stoddart, M. A Drug Holiday Reduces the Frequency and Severity of Medication-Related Osteonecrosis of the Jaw in a Minipig Model. J. Bone Miner. Res. 2020, 35, 2179–2192. [Google Scholar] [CrossRef]
  8. Fortunato, L.; Amato, M.; Simeone, M.; Bennardo, F.; Barone, S.; Giudice, A. Numb chin syndrome: A reflection of malignancy or a harbinger of MRONJ? A multicenter experience. J. Stomatol. Oral Maxillofac. Surg. 2018, 119, 389–394. [Google Scholar] [CrossRef]
  9. 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]
  10. 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]
  11. Lee, K.H.; Kim, S.H.; Kim, C.H.; Min, B.J.; Kim, G.J.; Lim, Y.; Kim, H.S.; Ahn, K.M.; Kim, J.H. Identifying genetic variants underlying medication-induced osteonecrosis of the jaw in cancer and osteoporosis: A case control study. J. Transl. Med. 2019, 17, 381. [Google Scholar] [CrossRef]
  12. Kim, K.M.; Rhee, Y.; Kwon, Y.D.; Kwon, T.G.; Lee, J.K.; Kim, D.Y. Medication Related Osteonecrosis of the Jaw: 2015 Position Statement of the Korean Society for Bone and Mineral Research and the Korean Association of Oral and Maxillofacial Surgeons. J. Bone Metab. 2015, 22, 151–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The genome analysis toolkit: A mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic. Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  16. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef] [PubMed]
  17. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Zidek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef]
  18. Rodrigues, C.H.M.; Pires, D.E.V.; Ascher, D.B. DynaMut2: Assessing changes in stability and flexibility upon single and multiple point missense mutations. Protein Sci. 2021, 30, 60–69. [Google Scholar] [CrossRef]
  19. Lee, S.; Seo, J.; Park, J.; Nam, J.Y.; Choi, A.; Ignatius, J.S.; Bjornson, R.D.; Chae, J.H.; Jang, I.J.; Lee, S.; et al. Korean Variant Archive (KOVA): A reference database of genetic variations in the Korean population. Sci. Rep. 2017, 7, 4287. [Google Scholar] [CrossRef] [Green Version]
  20. Kim, J.H.; Ko, Y.J.; Kim, J.Y.; Oh, Y.; Hwang, J.; Han, S.; Kim, S.; Lee, J.H.; Han, D.H. Genetic investigation of bisphosphonate-related osteonecrosis of jaw (BRONJ) via whole exome sequencing and bioinformatics. PLoS ONE 2015, 10, e0118084. [Google Scholar] [CrossRef]
  21. Cheng, J.; Peng, W.; Cao, X.; Huang, Y.; Lan, X.; Lei, C.; Chen, H. Differential Expression of KCNJ12 Gene and Association Analysis of Its Missense Mutation with Growth Traits in Chinese Cattle. Animals 2019, 9, 273. [Google Scholar] [CrossRef] [Green Version]
  22. Soufi, M.; Ruppert, V.; Rinne, S.; Mueller, T.; Kurt, B.; Pilz, G.; Maieron, A.; Dodel, R.; Decher, N.; Schaefer, J.R. Increased KCNJ18 promoter activity as a mechanism in atypical normokalemic periodic paralysis. Neurol. Genet. 2018, 4, e274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Maggi, L.; Bonanno, S.; Altamura, C.; Desaphy, J.F. Ion Channel Gene Mutations Causing Skeletal Muscle Disorders: Pathomechanisms and Opportunities for Therapy. Cells 2021, 10, 1521. [Google Scholar] [CrossRef] [PubMed]
  24. Paninka, R.M.; Mazzotti, D.R.; Kizys, M.M.; Vidi, A.C.; Rodrigues, H.; Silva, S.P.; Kunii, I.S.; Furuzawa, G.K.; Arcisio-Miranda, M.; Dias-da-Silva, M.R. Whole genome and exome sequencing realignment supports the assignment of KCNJ12, KCNJ17, and KCNJ18 paralogous genes in thyrotoxic periodic paralysis locus: Functional characterization of two polymorphic Kir2.6 isoforms. Mol. Genet. Genom. 2016, 291, 1535–1544. [Google Scholar] [CrossRef]
  25. 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]
  26. Donato, R. S100: A multigenic family of calcium-modulated proteins of the EF-hand type with intracellular and extracellular functional roles. Int. J. Biochem. Cell Biol. 2001, 33, 637–668. [Google Scholar] [CrossRef]
  27. Guo, Y.F.; Su, T.; Yang, M.; Li, C.J.; Guo, Q.; Xiao, Y.; Huang, Y.; Liu, Y.; Luo, X.H. The role of autophagy in bone homeostasis. J. Cell. Physiol. 2021, 236, 4152–4173. [Google Scholar] [CrossRef]
  28. Mzoughi, S.; Fong, J.Y.; Papadopoli, D.; Koh, C.M.; Hulea, L.; Pigini, P.; Di Tullio, F.; Andreacchio, G.; Hoppe, M.M.; Wollmann, H.; et al. PRDM15 is a key regulator of metabolism critical to sustain B-cell lymphomagenesis. Nat. Commun. 2020, 11, 3520. [Google Scholar] [CrossRef]
  29. Ishii, J.; Adachi, H.; Aoki, J.; Koizumi, H.; Tomita, S.; Suzuki, T.; Tsujimoto, M.; Inoue, K.; Arai, H. SREC-II, a new member of the scavenger receptor type F family, trans-interacts with SREC-I through its extracellular domain. J. Biol. Chem. 2002, 277, 39696–39702. [Google Scholar] [CrossRef] [Green Version]
  30. Anastasio, N.; Ben-Omran, T.; Teebi, A.; Ha, K.C.; Lalonde, E.; Ali, R.; Almureikhi, M.; Der Kaloustian, V.M.; Liu, J.; Rosenblatt, D.S.; et al. Mutations in SCARF2 are responsible for Van Den Ende-Gupta syndrome. Am. J. Hum. Genet. 2010, 87, 553–559. [Google Scholar] [CrossRef] [Green Version]
  31. Kalamajski, S.; Aspberg, A.; Lindblom, K.; Heinegard, D.; Oldberg, A. Asporin competes with decorin for collagen binding, binds calcium and promotes osteoblast collagen mineralization. Biochem. J. 2009, 423, 53–59. [Google Scholar] [CrossRef]
  32. Aki, T.; Hashimoto, K.; Ogasawara, M.; Itoi, E. A whole-genome transcriptome analysis of articular chondrocytes in secondary osteoarthritis of the hip. PLoS ONE 2018, 13, e0199734. [Google Scholar] [CrossRef]
  33. Tatsuda, Y.; Iguchi, K.; Usui, S.; Suzui, M.; Hirano, K. Protein kinase C is inhibited by bisphosphonates in prostate cancer PC-3 cells. Eur. J. Pharmacol. 2010, 627, 348–353. [Google Scholar] [CrossRef] [PubMed]
  34. Nicoletti, P.; Cartsos, V.M.; Palaska, P.K.; Shen, Y.; Floratos, A.; Zavras, A.I. Genomewide pharmacogenetics of bisphosphonate-induced osteonecrosis of the jaw: The role of RBMS3. Oncologist 2012, 17, 279–287. [Google Scholar] [CrossRef] [Green Version]
  35. Kiel, D.P.; Demissie, S.; Dupuis, J.; Lunetta, K.L.; Murabito, J.M.; Karasik, D. Genome-wide association with bone mass and geometry in the Framingham Heart Study. BMC Med. Genet. 2007, 8 (Suppl. 1), S14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Busch, T.; Armacki, M.; Eiseler, T.; Joodi, G.; Temme, C.; Jansen, J.; von Wichert, G.; Omary, M.B.; Spatz, J.; Seufferlein, T. Keratin 8 phosphorylation regulates keratin reorganization and migration of epithelial tumor cells. J. Cell Sci. 2012, 125, 2148–2159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ku, N.O.; Azhar, S.; Omary, M.B. Keratin 8 phosphorylation by p38 kinase regulates cellular keratin filament reorganization: Modulation by a keratin 1-like disease causing mutation. J. Biol. Chem. 2002, 277, 10775–10782. [Google Scholar] [CrossRef] [Green Version]
  38. Toivola, D.M.; Ku, N.O.; Resurreccion, E.Z.; Nelson, D.R.; Wright, T.L.; Omary, M.B. Keratin 8 and 18 hyperphosphorylation is a marker of progression of human liver disease. Hepatology 2004, 40, 459–466. [Google Scholar] [CrossRef]
  39. Reszka, A.A.; Halasy-Nagy, J.; Rodan, G.A. Nitrogen-bisphosphonates block retinoblastoma phosphorylation and cell growth by inhibiting the cholesterol biosynthetic pathway in a keratinocyte model for esophageal irritation. Mol. Pharmacol. 2001, 59, 193–202. [Google Scholar] [CrossRef] [Green Version]
  40. 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–8. [Google Scholar] [CrossRef] [Green Version]
  41. George, E.L.; Lin, Y.L.; Saunders, M.M. Bisphosphonate-related osteonecrosis of the jaw: A mechanobiology perspective. Bone Rep. 2018, 8, 104–109. [Google Scholar] [CrossRef] [PubMed]
  42. Allen, M.R.; Burr, D.B. Mandible matrix necrosis in beagle dogs after 3 years of daily oral bisphosphonate treatment. J. Oral Maxillofac. Surg. 2008, 66, 987–994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Roelofs, A.J.; Coxon, F.P.; Ebetino, F.H.; Lundy, M.W.; Henneman, Z.J.; Nancollas, G.H.; Sun, S.; Blazewska, K.M.; Bala, J.L.; Kashemirov, B.A.; et al. Fluorescent risedronate analogues reveal bisphosphonate uptake by bone marrow monocytes and localization around osteocytes in vivo. J. Bone Miner. Res. 2010, 25, 606–616. [Google Scholar] [CrossRef] [PubMed]
  44. Kimura, M.; Nagao, T.; Machida, J.; Warnakulasuriya, S. Mutation of keratin 4 gene causing white sponge nevus in a Japanese family. Int. J. Oral Maxillofac. Surg. 2013, 42, 615–618. [Google Scholar] [CrossRef]
  45. Cai, W.; Chen, Z.; Jiang, B.; Yu, F.; Xu, P.; Wang, M.; Wan, R.; Liu, J.; Xue, Z.; Yang, J.; et al. Keratin 13 mutations associated with oral white sponge nevus in two Chinese families. Meta Gene 2014, 2, 374–383. [Google Scholar] [CrossRef]
  46. Tao, G.Z.; Looi, K.S.; Toivola, D.M.; Strnad, P.; Zhou, Q.; Liao, J.; Wei, Y.; Habtezion, A.; Omary, M.B. Keratins modulate the shape and function of hepatocyte mitochondria: A mechanism for protection from apoptosis. J. Cell Sci. 2009, 122, 3851–3855. [Google Scholar] [CrossRef] [Green Version]
  47. Toivola, D.M.; Boor, P.; Alam, C.; Strnad, P. Keratins in health and disease. Curr. Opin. Cell Biol. 2015, 32, 73–81. [Google Scholar] [CrossRef]
  48. Gilbert, S.; Loranger, A.; Daigle, N.; Marceau, N. Simple epithelium keratins 8 and 18 provide resistance to Fas-mediated apoptosis. The protection occurs through a receptor-targeting modulation. J. Cell Biol. 2001, 154, 763–773. [Google Scholar] [CrossRef] [Green Version]
  49. Kovacic, N.; Grcevic, D.; Katavic, V.; Lukic, I.K.; Marusic, A. Targeting Fas in osteoresorptive disorders. Expert Opin. Targets 2010, 14, 1121–1134. [Google Scholar] [CrossRef] [Green Version]
  50. Galea, G.L.; Meakin, L.B.; Sugiyama, T.; Zebda, N.; Sunters, A.; Taipaleenmaki, H.; Stein, G.S.; van Wijnen, A.J.; Lanyon, L.E.; Price, J.S. Estrogen receptor alpha mediates proliferation of osteoblastic cells stimulated by estrogen and mechanical strain, but their acute down-regulation of the Wnt antagonist Sost is mediated by estrogen receptor beta. J. Biol. Chem. 2013, 288, 9035–9048. [Google Scholar] [CrossRef] [Green Version]
  51. Hayashi, M.; Nakashima, T.; Yoshimura, N.; Okamoto, K.; Tanaka, S.; Takayanagi, H. Autoregulation of Osteocyte Sema3A Orchestrates Estrogen Action and Counteracts Bone Aging. Cell Metab. 2019, 29, 627–637.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Camuzard, O.; Santucci-Darmanin, S.; Breuil, V.; Cros, C.; Gritsaenko, T.; Pagnotta, S.; Cailleteau, L.; Battaglia, S.; Panaia-Ferrari, P.; Heymann, D.; et al. Sex-specific autophagy modulation in osteoblastic lineage: A critical function to counteract bone loss in female. Oncotarget 2016, 7, 66416–66428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hansen, T.; Kunkel, M.; Weber, A.; James Kirkpatrick, C. Osteonecrosis of the jaws in patients treated with bisphosphonates—histomorphologic analysis in comparison with infected osteoradionecrosis. J. Oral Pathol. Med. 2006, 35, 155–160. [Google Scholar] [CrossRef]
  54. Kos, M.; Junka, A.; Smutnicka, D.; Szymczyk, P.; Gluza, K.; Bartoszewicz, M. Bisphosphonates enhance bacterial adhesion and biofilm formation on bone hydroxyapatite. J. Craniomaxillofac. Surg. 2015, 43, 863–869. [Google Scholar] [CrossRef]
  55. Sedghizadeh, P.P.; Stanley, K.; Caligiuri, M.; Hofkes, S.; Lowry, B.; Shuler, C.F. Oral bisphosphonate use and the prevalence of osteonecrosis of the jaw: An institutional inquiry. J. Am. Dent. Assoc. 2009, 140, 61–66. [Google Scholar] [CrossRef] [Green Version]
  56. Thumbigere-Math, V.; Michalowicz, B.S.; Hodges, J.S.; Tsai, M.L.; Swenson, K.K.; Rockwell, L.; Gopalakrishnan, R. Periodontal disease as a risk factor for bisphosphonate-related osteonecrosis of the jaw. J. Periodontol. 2014, 85, 226–233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jcm 11 02145 g001 550
Figure 1. Panoramic radiographic images of MRONJ patients in this study. (A) Extraction of the mandibular molars. (B) Three months after the extraction. Red arrows point to the necrotic area.
Figure 1. Panoramic radiographic images of MRONJ patients in this study. (A) Extraction of the mandibular molars. (B) Three months after the extraction. Red arrows point to the necrotic area.
Jcm 11 02145 g001
Jcm 11 02145 g002 550
Figure 2. Comparisons of the variants among four datasets. The number of variants identified from bidirectional approaches (blood and lesion tissue) compared with the previous data set from the KOVA and GSK-saliva. Overlapping areas in the Venn diagram represent common variants between the comparison groups. (A) SNPs variants; (B) InDels variants.
Figure 2. Comparisons of the variants among four datasets. The number of variants identified from bidirectional approaches (blood and lesion tissue) compared with the previous data set from the KOVA and GSK-saliva. Overlapping areas in the Venn diagram represent common variants between the comparison groups. (A) SNPs variants; (B) InDels variants.
Jcm 11 02145 g002
Jcm 11 02145 g003 550
Figure 3. GO functional analysis of candidate genes. GO enrichment analysis of WES genes was retrieved using Metascape software. Significantly (p < 0.05) enriched GO terms involved in biological process (pink), cellular component (blue), and molecular function (yellow) and branches are presented. The adjusted statistically significant values were negative 10-base log-transformed. The X-axis represents the enrichment scores of the terms −log10(p-value), and the Y-axis represents the enriched GO terms in the pathway. P-value was adjusted by Benjamini–Hochberg FDR. GO, gene ontology.
Figure 3. GO functional analysis of candidate genes. GO enrichment analysis of WES genes was retrieved using Metascape software. Significantly (p < 0.05) enriched GO terms involved in biological process (pink), cellular component (blue), and molecular function (yellow) and branches are presented. The adjusted statistically significant values were negative 10-base log-transformed. The X-axis represents the enrichment scores of the terms −log10(p-value), and the Y-axis represents the enriched GO terms in the pathway. P-value was adjusted by Benjamini–Hochberg FDR. GO, gene ontology.
Jcm 11 02145 g003
Jcm 11 02145 g004 550
Figure 4. Comparison of the wild-type protein structure with the mutant. (A) Gene structures. Blue characters indicate rare variants that have MAF < 0.05. (B) Structural comparison of the distance between wild-type (left side) and mutants (right side) in KRT18 (a–f) and PABPC3 (g–i) protein. The adjacent amino acid is noted by black words, and the predicted distance between two residues in the protein structure is noted by orange words.
Figure 4. Comparison of the wild-type protein structure with the mutant. (A) Gene structures. Blue characters indicate rare variants that have MAF < 0.05. (B) Structural comparison of the distance between wild-type (left side) and mutants (right side) in KRT18 (a–f) and PABPC3 (g–i) protein. The adjacent amino acid is noted by black words, and the predicted distance between two residues in the protein structure is noted by orange words.
Jcm 11 02145 g004
Jcm 11 02145 g005 550
Figure 5. Analysis of stability and dynamic changes between WT and mutants in KRT18. Stability and dynamic changes between WT (A–C) and five mutations (D–F) of KRT18 were analyzed. The figure shows the change in Gibbs free energy (in Kcal/mol) for mutations and non-covalent interactions established by the mutated residue. The orange and blue dash lines indicate polarity and van der Waals interactions, respectively.
Figure 5. Analysis of stability and dynamic changes between WT and mutants in KRT18. Stability and dynamic changes between WT (A–C) and five mutations (D–F) of KRT18 were analyzed. The figure shows the change in Gibbs free energy (in Kcal/mol) for mutations and non-covalent interactions established by the mutated residue. The orange and blue dash lines indicate polarity and van der Waals interactions, respectively.
Jcm 11 02145 g005
Table
Table 1. Clinical characteristics of patients with MRONJ.
Table 1. Clinical characteristics of patients with MRONJ.
Patient No.GenderAgeBPs
Administered		Duration of
Administration (Years)		Affected
Area		Type of Dental
Intervention		StageSmoking
1Male70Risedronate10Maxilla and MandibleTooth extractionIIINo
2Female77Risedronate10MandibleTooth extractionIINo
3Female86Risedronate15MandibleTooth extractionIIINo
4Female81Risedronate20MaxillaTooth extractionIIINo
5Female80Alendronate5MaxillaSpontaneousIIINo
6Female77Risedronate4MandibleTooth extractionIIINo
7Female67Risedronate5MaxillaTooth extractionIIINo
8Female82Risedronate10MandibleTooth extractionIIINo
9Female76Risedronate20MandibleTooth extractionIINo
10Female80Risedronate10MandibleRoot canal therapyIIINo
Table
Table 2. The number of variant candidates in MRONJ patients.
Table 2. The number of variant candidates in MRONJ patients.
VariantSNPsInDels
Saliva;
Blood		Saliva;
Lesion		Saliva;
Blood;
Lesion		Lesion-SpecificSaliva;
Blood		Saliva;
Lesion		Saliva;
Blood;
Lesion		Lesion-Specific
Intragenic15752319961138202238750
 Exonic73311866184414046
 Intronic6716108255631138180624
 Exonic; Splicing--1--1--
 Splicing--8-1692
 5’UTR5414416112536
 3’UTR121982114442
 5’UTR;3’UTR--------
Intergenic331528141641215213
Upstream; Downstream--11---4
Upstream212419-4131
Downstream-1913--17
ncRNA exonic12712815-11427
ncRNA splicing---1----
ncRNA intronic1237175412657
ncRNA exonic; splicing--------
NA--163-----
Total variants2167938761151462412651089
(Total Genes)14254212559841205230687
Table
Table 3. Effects of variants in the exonic regions.
Table 3. Effects of variants in the exonic regions.
Variant EffectsSingle Nucleotide Polymorphisms (SNPs)Insertion/Deletions (InDels)
Saliva;
Blood		Saliva;
Lesion		Saliva;
Blood;
Lesion		Lesion-SpecificSaliva;
Blood		Saliva;
Lesion		Saliva;
Blood;
Lesion		Lesion-Specific
Missense43128829----
Synonymous29109255----
Unknown19524121222
Stop gain--8-1---
Nonframeshift insertion-----6313
Nonframeshift deletion-----6520
Frameshift insertion----2655
Frameshift deletion---- 356
Total variants73311866184424046
(Total Genes)(49)(21)(1261)(12)(3)(37)(38)(34)
Table
Table 4. Annotation of rare variants with deleterious effects.
Table 4. Annotation of rare variants with deleterious effects.
TypePositionIDREFALTGeneNucleotide
Change		Aminoacid
Change		Exonic
Function		SIFTPolyphen2
HVAR/HDIV
SNPsChr1:
145368473		rs1043749GCNBPF9./../.unknown./../.
Chr1:
145368518		rs61813437CTNBPF9./../.unknown./../.
Chr11:
1213275		rs71251383GAMUC5AC./../.unknown./../.
Chr12:
53343033		rs78514003GAKRT18c.G76Ap.A26Tmissense SNV0.549, T0.003, B
/0.005, B
Chr12:
53343036		rs77825282CTKRT18c.C79Tp.R27Wmissense SNV0.032, D0.001, B
/0.001, B
Chr12:
53343040		rs74379840CAKRT18c.C83Ap.P28Qmissense SNV0.098, T0.716, P
/0.982, D
Chr12:
53343051		rs74953757GTKRT18c.G94Tp.A32Smissense SNV0.171, T0.688, P
/0.91, P
Chr12:
53343059		rs78343594CAKRT18c.C102Ap.S34Rmissense SNV0.001, D0.081, B
/0.087, B
Chr12:
53343069		rs77999286GTKRT18c.G112Tp.G38Cmissense SNV0.009, D0.923, D
/0.988, D
Chr12:
53343084		rs75441140GCKRT18c.G127Cp.G43Rmissense SNV0.001, D0.554, P
/0.949, P
Chr12:
53343099		.GAKRT18c.G142Ap.V48Mmissense SNV0.147, T0.11, B
/0.642, P
Chr12:
53343105		rs78479490CTKRT18c.C148Tp.R50Cmissense SNV0.299, T0.004, B
/0.003, B
Chr12:
53343124		rs76183244GAKRT18c.G167Ap.G56Dmissense SNV0.054, T0.015, B
/0.017, B
Chr12:
53343148		.CGKRT18c.C191Gp.A64Gmissense SNV0.718, T0.009, B
/0.028, B
Chr12:
53343158		rs77364359AGKRT18c.A201Gp.I67Mmissense SNV0.493, T0.0, B
/0.0, B
Chr13:
25670877		rs112107735GAPABPC3c.G541Ap.A181Tmissense SNV0.042, D0.253, B
/0.627, P
Chr13:
25670907		rs76264750CAPABPC3c.C571Ap.P191Tmissense SNV1.0, T 0.0, B
/0.0, B
Chr13:
25670919		rs76861216AGPABPC3c.A583Gp.I195Vmissense SNV1.0, T0.023, B
/0.005, B
Chr13:
25670953		rs74040928GAPABPC3c.G617Ap.R206Hmissense SNV0.089, T0.043, B
/0.114, B
Chr13:
25670955		rs79397892CTPABPC3c.C619Tp.L207Fmissense SNV0.192, T0.978, D
/0.999, D
Chr13:
25670988		rs74564616TGPABPC3c.T652Gp.L218Vmissense SNV0.287, T0.041, B
/0.027, B
Chr13:
25671027		rs78826513AGPABPC3c.A691Gp.K231Emissense SNV0.001, D0.953, D
/0.995, D
Chr13:
25671172		rs79593984CAPABPC3c.C836Ap.T279Kmissense SNV1.0, T0.0, B
/0.0, B
Chr13:
25671214		rs201081849TGPABPC3c.T878Gp.V293Gmissense SNV1.0, T0.0, B
/0.0, B
IndelsChr1:
17085999		rs59375146CCCCGCMST1L./../.frameshift deletion./../.
Chr9:
95237024		.CTCATCACTCATCATCAASPN./../.Non-frameshift deletion./../.
Chr9:
95237024		.CTCATCACTCAASPN./../.Non-frameshift deletion./../.
Chr12:
7045891		rs797045323ACAGCAGCAG
CAGCAGAGCAG
CAGCAG		ACAGCAGCAGCAGATN1./../.Non-frameshift insertion./../.
Chr12:
7045891		rs797045323ACAGCAGCAG
CAGCAGCAGCA
GCAGCAG		ACAGCAGCAG
CAGCAGCAGCA
GCAGCAGCAGCAG		ATN1./../.Non-frameshift insertion./../.
Chr13:
25671149		rs368285293ACGGAPABPC3./../.Non-frameshift deletion./../.
Chr13:
78272267		rs201380414TTGGSLAIN1./../.frameshift insertion./../.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/jcm11082145

AMA Style

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. Journal of Clinical Medicine. 2022; 11(8):2145. https://doi.org/10.3390/jcm11082145

Chicago/Turabian Style

Kim, Songmi, Seyoung Mun, Wonseok Shin, Kyudong Han, and Moon-Young Kim. 2022. "Identification of Potentially Pathogenic Variants Associated with Recurrence in Medication-Related Osteonecrosis of the Jaw (MRONJ) Patients Using Whole-Exome Sequencing" Journal of Clinical Medicine 11, no. 8: 2145. https://doi.org/10.3390/jcm11082145

APA Style

Kim, S., Mun, S., Shin, W., Han, K., & Kim, M. -Y. (2022). Identification of Potentially Pathogenic Variants Associated with Recurrence in Medication-Related Osteonecrosis of the Jaw (MRONJ) Patients Using Whole-Exome Sequencing. Journal of Clinical Medicine, 11(8), 2145. https://doi.org/10.3390/jcm11082145

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