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.