# Patent

## Title
Pharmacogenetic test anti-resorptive therapy-associated osteonecrosis of the jaw

## Abstract

The present invention relates to methods and compositions for testing individuals to determine whether they are at increased risk of developing anti-resorptive therapy-associated osteonecrosis of the jaw.

## Background
PRIORITY CLAIM
This application claims priority to U.S. Provisional Application Ser. No. 61/471,532, filed Apr. 4, 2011, the contents of which is incorporated by reference in its entirety herein.
SEQUENCE LISTING
The instant application contains a Sequence Listing which has been submitted in ASCII format via EFS-Web and is hereby incorporated by reference in its entirety. Said ASCII copy, created on Nov. 19, 2012, is named 07005045.txt and is 59,951 bytes in size.


GRANT INFORMATION
This invention was made with government support under Grant No. DE018143-01 awarded by the National Institute of Dental and Craniofacial Research, National Institutes of Health. The government has certain rights in the invention.


1. INTRODUCTION
The present invention relates to methods and compositions for determining whether a subject is at increased risk for developing anti-resorptive therapy-related osteonecrosis of the jaw.
2. BACKGROUND OF THE INVENTION
Bisphosphonates (BPs) are widely prescribed anti-osteoclastic medications. The intravenously administered BPs pamidronate and zoledronic acid are used in oncology to control bone metastasis and hypercalcemia. Oral BPs are used to control or prevent bone loss in osteoporosis, including osteoporosis associated with menopause. An estimated 3 million American women are currently being treated with oral bisphosphonates [1]. The monoclonal antibody Denosumab is also used to treat these conditions.
BPs are synthetic analogs of pyrophosphate that readily localize to bones due to their affinity for hydroxyapatite, and reduce osteoclastic activity. They are not readily metabolized, and thus, have long-lasting effects that might extend for several years. BPs are especially attracted to, and localize in, areas of the bone undergoing inflammation or resorption. They are subsequently phagocytozed and internalized by osteoclasts. These internalized bisphosphonates, in turn, trigger apoptosis (cell death) of the osteoclasts, thus inhibiting osteoclast-mediated bone resorption [2]. Osteoclasts seem to be affected by BPs both in terms of number and function. Animal studies have also demonstrated some antiangiogenic properties, which may partially explain the development of osteonecrosis due to limited healing ability of the bone because of reduced vasculature [3].
BPs, especially zoledronic acid, have been associated with a serious adverse effect, osteonecrosis of the jaw. According to the American Association of Oral and Maxillofacial Surgeons (AAOMS), BP-related osteonecrosis of the jaw (BRONJ) is defined as exposed bone in the maxillofacial region for more than eight weeks in patients treated with a bisphosphonate that have no prior history of radiation therapy to the jaws [4]. The non-healing exposed necrotic lesions may involve the mandible or the maxilla or both, and can be painful, persistent, and resistant to treatment. The incidence of BRONJ varies in different studies. BRONJ affects as many as 5-10% of zoledronic acid users and far fewer users of oral bisphosphonates.
Osteonecrosis of the jaw has also been reported in association with denosumab treatment [29]. It would be desirable to identify individuals who are at risk for developing osteonecrosis of the jaw, so that subjects at greater risk could be considered for alternatives to anti-resorptive therapy or, if suitable alternatives are not available, could be monitored more closely.
3. SUMMARY OF THE INVENTION
The present invention relates to methods and compositions for testing individuals to determine whether they are at increased risk of developing anti-resorptive therapy-associated osteonecrosis of the jaw (“ARONJ”). It is based, at least in part, on the results of a genome wide association analysis which revealed that certain Single Nucleotide Polymorphisms (“SNPs”) are significantly associated with osteonecrosis of the jaw among bisphosphonate users, including SNPs in the RNA-binding motif, single-stranded-interacting protein 3 (“RBMS3”) gene (for example, the SNP rs17024608) as well as SNPs in other genes, including but not limited to those for insulin-like growth factor I receptor (“IGF1R”), insulin-like growth factor binding protein 7 (“IGFBP7”), dihydropyrimidine dehydrogenase (“DPYD”), ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (“ABCC4”), and glutathione S-transferase mu 2 (“GSTM2”) and other SNPs as listed in TABLES 1-5.



4. BRIEF DESCRIPTION OF THE FIGURES
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
FIG. 1A-B: Population structure of Caucasian subjects from the combined ONJ/Hap Map 3/European collections. The red circles represent ONJ subjects, azure dots CEU Hap Map 3, blue dots CEU Hap Map III, yellow dots North European POPRES, green dots central European POPRES, gray dots Spanish POPRES, light yellow dots South European POPRES. (A) PC axis 2 versus PC axis 1 (ordinate versus abscissa). (B) PC axis 1 versus PC axis 3 (ordinate versus abscissa).
FIG. 2. QQ plot for Caucasian ONJ Study Group.
FIG. 3. Manhattan Plot for Caucasian ONJ Study Group.
FIG. 4: Manhattan Plot of the region surrounding rs17024608. Red dots represent the imputed markers; Gray dots represent the original markers.
FIG. 5A-B. Results of SNP analysis within 20 kb of IGF gene family candidate genes. (A) QQ for candidate SNPs (logistic test, drug-exposed study group) and (B) SNP quality graph of the top hit.
FIG. 6A-B. Results of SNP analysis of ADME genes. (A) QQ for ADME SNPs (logistic test, drug-exposed study group) and (B) SNP quality graph of the top hit.
FIG. 7A-B. SNP analysis of IGF-related SNPs in the extended study group. (A) QQ plot for candidate SNPs (logistic test, extended study group) and (B) SNP quality graph of the top hit.
FIG. 8A-B. SNP analysis of ADME-related SNPs in the extended study group. (A) QQ plot for ADME SNPs (logistic test, extended study group) and (B) SNP quality graph of the top hit.
FIG. 9. QQ plots for IGF-candidate SNPs (logistic test, final study group).
FIG. 10A-B. (A) QQ plots for ADME SNPs (logistic test, final study group) and (B) SNP quality graph of the top hit.
FIG. 11A-Q. Nucleotide sequence (SEQ ID NO: 86)of the intron between positions 29941247 to 29977576 of the human RBMS3 gene on chromosome 3 (numbered consistently with their the numbering of the RBMS3 gene as represented in NCBI Ace. No. NC—000003 Region 29322803 . . . 30051886 GPC—000000027). Sequence flanking the A/G substitution in rs17024608 is bolded and underlined, and the variant base (A) involved in the A to G substitution is shown by a capital letter. For clarity, residue 618421 in the figure is residue 29941247 of human chromosome 3.



5. DETAILED DESCRIPTION OF THE INVENTION
For purposes of clarity of description and not by way of limitation, the detailed description of the invention is divided into the following subsections:

    
    
        (i) biomarker genes and SNPs;
        (ii) methods of treatment/diagnosis; and
        (iii) kits.
    
    


The term “anti-resorptive therapy” as used herein refers to therapy in which an agent (“anti-resorptive agent”) is administered that inhibits resorption of bone, for example in the treatment of osteoporosis or in oncology to control bone metastasis and/or hypercalcemia. Non-limiting examples of anti-resorptive agents include, but are not limited to, bisphosphonates and monoclonal antibodies specific for RANK ligand, such as, but not limited to, denosumab.
“Anti-resorptive therapy-associated osteonecrosis of the jaw”, or “ARONJ”, includes BRONJ as well as ONJ associated with other anti-resorptive agents.
The term “bisphosphonates” as used herein refers to a class of drugs comprising a bisphosphonate structure that may be used to treat osteoporosis or hypercalcemia. Non-limiting examples of bisphosphonates include orally-administered drugs such as alendronate (Fosamax), etidronate (Didronel), ibadronate (Boniva), and risedronate (Actonel) and intravenously administered drugs such as pamidronate (Aredia) and zoledronic acid (Zometa).
The term “osteonecrosis of the jaw” (“ONJ”) means a clinical condition characterized by, in a subject, exposed bone in the maxillofacial region that persists for more than eight weeks. Typically one or more region of the affected bone area is necrotic.
The term “allelic variation” refers to the presence, in a population, of different forms of the same gene characterized by differences in their nucleotide sequences (sequences in genomic DNA). The variation may be in the form of one or more substitution, insertion, or deletion of a nucleotide. Different alleles may be functionally the same, or may be functionally different. In one subset of allelic variations a single nucleotide is different between alleles, and is referred to as a Single Nucleotide Polymorphism (“SNP”). Allelic variation in a known sequence may be identified by standard sequencing techniques. A “variation” or “variant,” as those terms are used herein, is relative to the ancestral gene found in the majority of the population. Unless specified otherwise, the presence of a SNP means that at the single nucleotide position for which alleles have been identified, the nucleotide present is the variant nucleotide, not the nucleotide found in the majority of the population (however, TABLE 4 lists SNPs and then specifies which allele is associated with ARONJ, and either the allele found in the majority of the population or a variant may be specified). The variation (variant) is comprised of a substituted nucleotide or nucleotides or an insertion or deletion of a nucleotide or nucleotides. Herein, generally the ancestral nucleotide is listed first and the variation (variant) nucleotide is listed second (for example, in A/G A is the ancestral nucleotide and G is the variation (variant) nucleotide). If there is an insertion, the ancestral nucleotide is represented by a hyphen (e.g., -/G). If there is a deletion, the variation (variant) nucleotide is represented by a hyphen (e.g., G/-). Numerous allelic variations (variants), captured in SNPs, of gene s are known in the art and catalogued (for example, in the National Center for Biotechnology Information “Entrez SNP”). Allelic variations that are not SNPs include deletions or insertions or substitutions of multiple consecutive nucleotides.
In non-limiting embodiments of the invention, the presence of an allelic variation, for example a SNP, may be determined using a technique such as, but not limited to, primer extension or polymerase chain reaction, using primer(s) designed based on sequence in the proximity of the variation, followed by sequencing. For example, and not by way of limitation, the presence of a SNP indicative of increased risk of ARONJ may be determined by a method comprising using at least one primer sequence complementary to a sequence flanking the location of the SNP (for example, within 80 nucleotides, or within 50 nucleotides, or within 30 nucleotides, or within 20 nucleotides, or within 10 nucleotides, of the SNP) in a primer extension reaction or polymerase chain reaction to generate a test fragment that contains the location of the SNP and determining the nucleotide present at the location of the single nucleotide polymorphism, for example by sequencing all or a portion of the test fragment.
A subject may be a human or a non-human subject. In a specific non-limiting embodiment, the subject suffers from osteoporosis. In a specific non-limiting embodiment, the subject is a postmenopausal woman. In a specific non-limiting embodiment, the subject has a cancer and has or is at risk for bone metastasis. In a specific non-limiting embodiment the subject has hypercalcemia.
To assess whether the subject carries a ARONJ biomarker as described herein, a sample of nucleic acid from the subject may be used. The nucleic acid may be genomic DNA or RNA reflective of the allelic variation or a cDNA copy thereof. For example, a sample comprising a cell from the subject may be collected. For example, the sample may be a tissue or body fluid, including but not limited to saliva, blood or its components, skin, hair follicle, urine, etc. The sample may be obtained by scraping the inside of the subject's mouth or cervix (eg in the context of a Pap smear). In a non-limiting example, as part of the detection process, nucleic acid may be at least partially purified from the sample.
5.1 Biomarker Genes and SNPS
In various non-limiting embodiments, the following genes and SNPs have been related to ARONJ, such that these genes and SNPs may be used as biomarkers for increased risk of ARONJ. Allelic variation and SNP variants may be used as indicators that a subject is at increased risk for ARONJ. The genes listed below, SNPs associated with those genes, and the further SNPs listed below are collectively referred to herein as “ARONJ biomarkers”.
5.1.1 RBMS3
RBMS3 is an RNA-binding protein that belongs to the c-myc gene single-strand binding protein family. Its transcripts are alternatively spliced to form different mRNAs. In humans it is located at about positions 29322803-30051886 on chromosome 3 (NCBI Reference Sequence NC—00003.11).
In one non-limiting embodiment of the present invention, an allelic variation in the RBMS3 gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of RBMS3 from about positions 29900000 to 29990000 of the human gene on chromosome 3 (see NC—000003.11), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in the intron of RBMS3 from about positions 29941247 to 29977576 of the human gene on chromosome 3, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ. The sequence of this intron is set forth in FIG. 11.
In another non-limiting embodiment of the present invention, an allelic variation in a region of the intron of RBMS3 from about positions 29941247 to 29977576 of the human gene on chromosome 3, said region being between about 29954000 to 29955000 (between about nucleotides numbered 631198 and 632198 in FIG. 11), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variant which is a SNP of RBMS3 selected from those set forth in TABLE 1 is an indicator that a subject is at increased risk for developing ARONJ. For clarity, the SNP is shown in bold text and underlined; for example, in rs17024608 the SNP is a substitution of ancestral nucleotide A by G, as shown by “A/G”.










TABLE 1





SEQ
SNP




ID
ref
Chr 3



No.
no.
position
Variation (bold, underlined)







 1
rs17024608
29954690
GATAGAATAGAACTATTTGATGTGGA/GCCATGAG





AATTTAAAAGTATCTGCA


 


 2
rs3821577
29910544
AAGGCAGACGTATGGTGCCATGATGAA/GATTGGGTCCCA





TGGAACACCCAATG


 


 3
rs9820707
29954837
TGAGAAATGACAGAAAAACTATATCAT/CCAGCAGAGAA





AACATCTTCTTGCTG


 


 4
rs9875937
29954551
ATTGAGGGCAGGCATAAGACCTCTTGG/TCTATCATCACT





AGGTGGCTCATAAA


 


 5
rs9876178
29954887
TGACCTAGGAAGGCATACAAATAGTTC/AAGTTGGCGTCA





CACTGTTTATGCAT


 


 6
rs13319154
29954187
ACTTTCAATTTTCTGAAGAAATGTATA/TGGGCAACTACAT





TAATATTCTTCAG


 


 7
rs13326291
29954399
CTCATGGAAAAGATGGCAAGATTGTTT/CAGTTTGTGCCA





GAATTTTCCTTAAT


 


 8
rs75830538
29954957
GTGGTACTTGTTCCCTTTTATATCACT/CTAACATAAGTAT





AATTTCATTGCAC


 


 9
rs114730671
29954785
TCTATCTTACAATAATTACTCAATATG/TTTCAGCAGACTC





TGAAATCTGCACT


 


10
rs116600197
29954066
TCAGCTGAAACAACTGTGTTATTCAGC/TGTAGTGTCTTAA





TTGGCAGTTACAT


 


11
rs116863073
29954716
ACCATGAGAATTTAAAAGTATCTGCAT/CAAATTGATGAT





TTGGCATCAGATAA


 


12
rs118013282
29954445
CTTAATAAATGTCTTCAAAGCTCTCAC/TAGGAGCCCAAG





CATGTACACAGTTG









In another non-limiting embodiment of the present invention, an allelic variant which is a SNP of RBMS3 as occurs in rs17024608, shown in TABLE 1, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variant which is a substitution of A by G at position 29954690 of chromosome 3 is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variant which is a substitution of A by G in SEQ ID NO: 1: (GATAGAATAGAACTATTTGATGTGGA/GCCATGAGAATTTAAAAGTATCTGCA) at the position indicated by A/G, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variant which is a SNP of RBMS3 selected from those set forth in TABLE 2 is an indicator that a subject is at increased risk for developing ARONJ.










TABLE 2





SEQ
SNP




ID
ref
Chr 3



NO.
no.
position
Variation







13
rs10510628
29853403
GCTCTGCCGTTTCTTAGGAAGTTGTGG/AAAGATATTGGT





CTTTTTGTGAATAT


 


14
rs4599260
29853269
CCATTAGGATAAAGAT/CGAGTGACCTCAAAAA


 


15
rs10514681
29853882
TAGTGGGAAAGTTTTAAAGAGACCGTC/TATTGAGTGCTT





TGATATGTTTGTTC


 


16
rs35393422
29853420
GAAGTTGTGGAAGATATTGGTCTTTTT/-





GTGAATATGTATGACACTATTCATT


 


17
rs79049188
29853041
CTCTCATGGAGTCTACATTCTAAGGTC/TTCATAGGAAAC





ACATGTAACTTTAC


 


18
rs115136555
29853407
TGCCGTTTCTTAGGAAGTTGTGGAAGA/GTATTGGTCTTTT





TGTGAATATGTAT









5.1.2 IGF1R
The IGF1R gene in humans is located at about positions 99145510 to 99555008 on chromosome 15 (NCBI Reference Sequence NG—009492).
In one non-limiting embodiment of the present invention, an allelic variation in the IGF1R gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
5.1.3 IGFBP7
The IGFBP7 gene in humans is located at about positions 57897244 to 57976539 on chromosome 4 (NCBI Reference Sequence NC—000004.11).
In one non-limiting embodiment of the present invention, an allelic variation in the IGFBP7 gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of IGFBP7 from about positions 57920000 to 57950000 of the human gene on chromosome 4 (see NC—000004.11), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of IGFBP7 from about positions 57930000 to 57940000 of the human gene on chromosome 4 (see NC—000004.11), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of IGFBP7 from about positions 57940000 to 57950000 of the human gene on chromosome 4 (see NC—000004.11), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation which is a SNP of IGFBP7 as occurs in rs11934877 is an indicator that a subject is at increased risk for developing ARONJ. The SNP of rs11934877 is located at position 57941026 of chromosome 4, where the ancestral nucleotide T is replaced by C, for example as in SEQ ID NO:19







TAATCTGTGTTAAAACAATATAGCATT/CATCTGCTTTGAATGCACTA


 


GGCACC.






In another non-limiting embodiment of the present invention, an allelic variation which is a SNP of IGFBP7 as occurs in rs17761305 is an indicator that a subject is at increased risk for developing ARONJ. The SNP of rs17761305 is located at position 57934091 of chromosome 4, where the ancestral nucleotide C is replaced by T, for example as in SEQ ID NO 20:







CCCCTGGAGAATAATTGATAGGGTAGC/TGAAAAATGTGGATATCATA


 


AAATAT.






5.1.4 DPYD
The DPYD gene in humans is located at about positions (−)98513111 to (−) 97416801 on the negative strand of chromosome 1 (NCBI Reference Sequence NG—008807.1)
In one non-limiting embodiment of the present invention, an allelic variation in the DPYD gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of DPYD gene from about positions (−)97800000 to (−)97500000 of the human gene on chromosome 1 (see NG—008807.1), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of DPYD gene from about positions (−)97700000 to (−)97650000 of the human gene on chromosome 1 (see NG—008807.1), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation which is a SNP of DPYD as occurs in rs10875055 is an indicator that a subject is at increased risk for developing ARONJ. The SNP of rs10875055 is located at position (−)97683997 of chromosome 1, where the ancestral nucleotide C is replaced by T, for example as in SEQ ID NO:21:







TTCATCTCACTAATAAGAGCTACCCAC/TCCGCCTTTATACAGAGGTT


 


CTCAGA






5.1.5 ABCC4
The ABCC4 gene in humans is located at about positions 95672083-95953687 on chromosome 13 (NCBI Reference Sequence NC—000013.10).
In one non-limiting embodiment of the present invention, an allelic variation in the ABCC4 gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of ABCC4 gene from about positions 95730000 to 95740000 of the human gene on chromosome 13 (see NC—000013.10), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of ABCC4 gene from about positions 95732000 to 95738000 of the human gene on chromosome 13 (see NC—000013.10), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation which is a SNP of ABCC4 as occurs in rs1189437 is an indicator that a subject is at increased risk for developing ARONJ. The SNP of rs1189437 is located at position 95735604 of chromosome 13, where the ancestral nucleotide A is replaced by C, for example as in SEQ ID NO:22:







GAGTGTAATCCTAACAACAACTCATGA/CAAGTATTTTTGAAAAGAA


 


TACTTGA.






5.1.6 GSTM2
The GSTM2 gene in humans is located at about positions 110210644 to 110226619 on chromosome 1 (NCBI Reference Sequence NC—000001.10).
In one non-limiting embodiment of the present invention, an allelic variation in the GSTM2 gene, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation in a region of GSTM2 gene from about positions 110210000 to 110215000 of the human gene on chromosome 3 (see NC—000001.10), where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In another non-limiting embodiment of the present invention, an allelic variation which is a SNP of GSTM2 as occurs in rs673151 is an indicator that a subject is at increased risk for developing ARONJ. The SNP of rs673151 is located at position 110213458 of chromosome 1, where the ancestral nucleotide G is replaced by A, for example as in SEQ ID NO:23:









GAAAGATGAGGAGATATTCAGAGGATG/AAGTGGAAGAAAGGAGGG



 



GGAAAAAG.






5.1.7 Further SNP Biomarkers
In further non-limiting embodiments of the present invention, the presence of an allelic variation which is a SNP selected from those set forth in TABLE 3 is an indicator that a subject is at increased risk for developing ARONJ.










TABLE 3





SEQ
SNP
Chromosome



ID
ref
location/



NO.
no.
gene
Variation







24
rs8012823
72343297chr14
CTGGAGGGCTTCCAAT/CGGACTGATCTCTTTG




DPF3 gene



 


25
rs11802277
Chr. 1
GAGAGTGAGAACTGAGTGGGCTGGGGAGCT





G/TATGGACTTTACCCTGCCATTTCCTAACCC


 


26
rs10754178
196416252chr 1
TTCTTTCCATCTTCAA/GTGAAATATTTTGCCT




KCNT2



 


27
rs13096022
7425350 chr3
TCAGGA/GTTAAC




GRM7



 


28
rs6900513
66333105chr6
GATAACCTCCAAAGAC/TTTGGTTGTAATTTTC




EYS



 


29
rs1873291
66344855chr6
TTAAAAATCAATTTGC/TCTTCATCACAGACAG




EYS



 


30
rs10781262
77650696chr9
TTAATAAAGGGTAAGATTGGGCTATCA/GTATTTGAATT





AGCAGAATCACTCTA


 


31
rs1471646
chr 1
GACTCATCTGACTTAGAAATGGGTGG[A/G]TGAAAAG




RP11-382E9.1
AATCTTCACTCACTATGT


 


32
rs4870310
chr 6
ancestral gene nucleotide is T; variant is C




RP11-
CTCGTTGTTTTTCTGGCACTACAAGAT/CGTTCCAGGTTC




15G8.1
ATCTTATATATTCT


 


33
rs10875148
98859508chr1
ancestral gene nucleotide is G; variant is A





TATTTTCTAACTCTTTCTGGTATAATG/AGGAACAGTCAA





GATCTGAACAAGAG


 


34
rs4562759
83835639chr10
ancestral gene nucleotide is T; variant is C




NRG3
AGAGAGTAGAATTAAGTGGTTTTTCAT/CTGGAATTATG





GAGGGAGAATGAAAT


 


35
rs4878512
27630418chr9
ancestral gene nucleotide is A; variant is G




RP11-
AAGAATGTAAAGCATTCTCTAGTCCTA/GTTTCCTTTCCT




53518.1
TGTTTGTTTTCTCT


 


36
rs12613966
AC113618.1
ancestral gene nucleotide is C; variant is T





AAATGTCACCTTTGAGTAGTGAAGTTC/TTGGATGATTTT





TATTTTCTTATTTT


 


37
rs11189381
RP11-459F3.3
ancestral gene nucleotide is T, variant is C





GCTGGCTTCTTCTTAATCAGAACTGCT/CCTTAGCTTCAA





GAGAGGCTGGAAAA


 


38
rs6861122
AC022120.1
ancestral gene nucleotide is C; variant is T





TTAGTAATATGCCTTTAAGGTACCTCC/TATCTAAACTGA





ACATGCTCATTA


 


39
rs4431170
MARCH1
ancestral gene nucleotide is A; variant is G





TGAAACCAAACTATAAACTACTTCTTA/GTCTTTGTGAGA





GAATTCCAGGGCAC









In another non-limiting embodiment of the present invention, an allelic variant which is a SNP of the human EYS gene in a region on chromosome 6 between 66325000 and 66350000, where the allelic variation may be a SNP, is an indicator that a subject is at increased risk for developing ARONJ.
In further non-limiting embodiments of the present invention, an allelic variation which is a SNP selected from those set forth in TABLE 4 is an indicator that a subject is at increased risk for developing ARONJ, where the nucleotide in the column marked “allele” is the allele indicative of the increased risk (among the SNPs listed in TABLE 4, sometimes the allele associated with ARONJ is the variant allele, and sometimes it is the allele found in the majority of the population tested). The SNPs listed in TABLE 4, except for rs17024608 which is discussed above, are denoted SNP ID #40-147. Commercial probes or sequences corresponding to the SNPs listed in TABLE 4 are publicly available, either for the SNP named (e.g. “rs17024608”) or at the chromosomal position specified (e.g., “12-85466404” denotes a SNP found at position 85466404 on human chromosome 12). For convenience, and not by way of limitation, the sequences of certain regions containing SNPs are set forth in TABLE 5 below as SEQ ID NOS: 40-85.


















TABLE 4





CHR
SNP
LOCATION
ALLELE
TESTED
SA|OR
SE
L95
U95
STAT
P. value
SNP ID#


























3
rs17024608
29954690
G
1744
5.825
0.3276
3.065
11.07
5.379
7.47E−08



22
rs5768434
46977516
T
1752
12.61
0.4784
4.937
32.2
5.298
1.17E−07
40


12
rs11064477
6944626
A
1689
21.66
0.6126
6.52
71.98
5.02
5.16E−07
41


12
12-85466404
85466404
A
1696
9.296
0.4769
3.65
23.67
4.675
2.95E−06
42


8
8-58133986
58133986
T
1746
7.326
0.427
3.173
16.92
4.664
3.10E−06
43


1
rs1886629
194421521
C
1755
3.698
0.2878
2.104
6.501
4.544
5.53E−06
44


2
rs7588295
166115757
G
1716
8.681
0.4783
3.399
22.17
4.518
6.24E−06
45


4
rs4431170
165504024
G
1748
5.176
0.3665
2.523
10.62
4.485
7.28E−06
46


6
rs7740004
120897902
A
1757
5.952
0.3992
2.722
13.02
4.469
7.87E−06
47


10
rs11189381
99553188
C
1728
6.816
0.4303
2.933
15.84
4.461
8.17E−06
48


15
rs12903202
56094085
G
1690
4.036
0.3145
2.179
7.477
4.436
9.15E−06
49


18
rs17751934
47455812
T
1773
5.009
0.3632
2.458
10.21
4.436
9.16E−06
50


11
11-23990403
23990403
C
1708
12.17
0.5655
4.016
36.86
4.419
9.94E−06
51


3
rs7613687
60124976
T
1733
13.44
0.6107
4.061
44.49
4.255
2.09E−05
52


3
rs7634338
8864809
C
1708
12.55
0.6001
3.87
40.68
4.215
2.50E−05
53


3
rs7612449
181586884
T
1718
11.79
0.5708
3.853
36.1
4.323
1.54E−05
54


4
rs10029016
8496903
A
1685
11.73
0.5612
3.905
35.24
4.388
1.15E−05
55


12
rs10505722
1102632
G
1716
11.39
0.5651
3.763
34.48
4.305
1.67E−05
56


14
rs2332834
71952831
C
1690
11.3
0.6162
3.378
37.82
3.935
8.31E−05
57


7
rs10271074
78589212
G
1713
11.08
0.557
3.719
33.02
4.318
1.57E−05
58


3
3-85058749
85058749
A
1729
10.99
0.6046
3.36
35.95
3.965
7.35E−05
59


5
5-159725177
159725177
C
1744
10.91
0.5466
3.736
31.84
4.371
1.24E−05
60


1
rs4951362
202580842
C
1719
10.66
0.5452
3.662
31.04
4.341
1.42E−05
61


2
rs62153910
96382119
T
1766
10.46
0.5396
3.632
30.11
4.35
1.36E−05
62


1
1-55273566
55273566
G
1740
10.37
0.5975
3.216
33.46
3.915
9.05E−05
63


2
2-15831165
15831165
A
1726
10.33
0.5364
3.609
29.55
4.352
1.35E−05
64


2
2-53263505
53263505
C
1722
10.03
0.5434
3.456
29.09
4.242
2.22E−05
65


2
rs4062819
95777963
A
1730
9.872
0.5353
3.457
28.19
4.277
1.89E−05
66


7
7-13058926
13058926
T
1696
9.791
0.528
3.479
27.56
4.321
1.55E−05
67


15
rs7176436
56752223
G
1731
9.767
0.5365
3.413
27.95
4.248
2.16E−05
68


2
2-95978053
95978053
A
1731
9.685
0.5337
3.403
27.57
4.254
2.10E−05
69


1
rs17127107
65099496
G
1685
9.492
0.5126
3.476
25.92
4.391
1.13E−05
70


7
7-78589823
78589823
A
1717
9.426
0.5521
3.194
27.82
4.063
4.84E−05
71


2
rs62156621
79426168
A
1727
9.31
0.524
3.334
26
4.258
2.06E−05
72


2
rs921245
137980811
C
1720
9.295
0.5267
3.311
26.09
4.233
2.30E−05
73


12
12-85423495
85423495
C
1701
9.284
0.5088
3.425
25.17
4.379
1.19E−05
74


12
12-85114303
85114303
C
1742
9.104
0.5056
3.379
24.53
4.368
1.25E−05
75


9
9-29422005
29422005
T
1739
9.06
0.5455
3.111
26.39
4.04
5.34E−05
76


1
rs1736563
169371174
A
1691
8.941
0.5421
3.09
25.87
4.041
5.32E−05
77


5
5-65503171
65503171
A
1765
8.934
0.5453
3.068
26.02
4.016
5.93E−05
78


7
7-78597753
78597753
G
1721
8.788
0.5503
2.989
25.84
3.949
7.83E−05
79


10
10-94242917
94242917
G
1743
8.785
0.5164
3.193
24.17
4.208
2.58E−05
80


10
10-94239487
94239487
G
1743
8.779
0.5164
3.191
24.15
4.207
2.59E−05
81


10
10-94264341
94264341
A
1739
8.758
0.5164
3.183
24.1
4.202
2.65E−05
82


2
2-15818709
15818709
A
1728
8.724
0.5391
3.033
25.1
4.018
5.87E−05
83


12
12-85440122
85440122
G
1702
8.663
0.5067
3.209
23.39
4.261
2.04E−05
84


12
rs4639998
85442168
T
1702
8.663
0.5067
3.209
23.39
4.261
2.04E−05
85


4
rs7669796
160863851
G
1716
8.662
0.5458
2.972
25.25
3.955
7.64E−05
86


12
rs17653326
85528244
A
1748
8.661
0.5037
3.227
23.24
4.286
1.82E−05
87


9
rs17057133
73509497
C
1692
8.616
0.5078
3.184
23.31
4.241
2.23E−05
88


16
rs3135009
12556820
T
1721
8.608
0.5224
3.092
23.96
4.121
3.78E−05
89


2
2-15819709
15819709
C
1729
8.597
0.539
2.989
24.73
3.991
6.57E−05
90


2
rs2544530
15820699
G
1728
8.569
0.5392
2.978
24.65
3.984
6.77E−05
91


8
rs16915413
94070482
A
1730
8.569
0.5422
2.961
24.8
3.962
9.30E+01
92


10
10-59772956
59772956
T
1735
8.567
0.5454
2.942
24.95
3.938
8.21E−05
93


9
9-17649682
17649682
G
1743
8.399
0.5149
3.062
23.04
4.133
3.58E−05
94


9
9-17652481
17652481
T
1742
8.394
0.5149
3.06
23.03
4.132
3.60E−05
95


3
3-85899151
85899151
T
1739
8.389
0.5348
2.941
23.93
3.977
6.98E−05
96


12
12-39900069
39900069
A
1696
8.387
0.5327
2.952
23.83
3.992
6.55E−05
97


16
rs2856790
12557384
C
1726
8.375
0.5203
3.02
23.22
4.084
4.42E−05
98


12
12-85389759
85389759
G
1768
8.369
0.5024
3.126
22.4
4.228
2.35E−05
99


3
3-29881391
29881391
C
1745
8.209
0.4884
3.152
21.38
4.31
1.63E−05
100


3
3-29881550
29881550
G
1745
8.207
0.4884
3.151
21.37
4.31
1.63E−05
101


3
rs6443519
179112920
A
1723
8.199
0.5176
2.973
22.61
4.065
4.81E−05
102


2
2-79440423
79440423
G
1736
8.146
0.5165
2.96
22.42
4.061
4.89E−05
103


9
9-26597745
26597745
A
1692
8.107
0.4926
3.087
21.29
4.248
2.16E−05
104


10
10-55210159
55210159
A
1734
8.027
0.5261
2.862
22.51
3.959
7.54E−05
105


2
rs7568908
79462538
T
1707
7.903
0.5157
2.876
21.72
4.008
6.11E−05
106


7
7-54640120
54640120
T
1753
7.771
0.4985
2.925
20.64
4.113
3.90E−05
107


7
7-54633747
54633747
G
1752
7.766
0.4984
2.924
20.63
4.112
3.92E−05
108


7
7-54635902
54635902
A
1752
7.766
0.4984
2.924
20.63
4.112
3.92E−05
109


7
7-54638869
54638869
G
1752
7.766
0.4984
2.924
20.63
4.112
3.92E−05
110


7
7-54638903
54638903
G
1752
7.766
0.4984
2.924
20.63
4.112
3.92E−05
111


2
rs2544538
15825910
A
1737
7.707
0.5184
2.79
21.29
3.939
8.19E−05
112


8
8-58130916
58130916
A
1732
7.601
0.4593
3.09
18.7
4.416
1.01E−05
113


2
2-15828311
15828311
A
1738
7.563
0.5184
2.738
20.89
3.903
9.51E−05
114


19
19-55065301
55065301
C
1732
7.557
0.5053
2.807
20.34
4.003
6.26E−05
115


2
2-21318018
21318018
A
1728
7.553
0.4896
2.893
19.72
4.13
3.63E−05
116


2
2-15838147
15838147
G
1716
7.549
0.5177
2.737
20.82
3.905
9.43E−05
117


20
20-55029650
55029650
C
1702
7.534
0.4877
2.897
19.6
4.141
3.46E−05
118


20
20-55029675
55029675
C
1702
7.534
0.4877
2.897
19.6
4.141
3.46E−05
119


2
rs57446887
21316772
C
1725
7.528
0.4896
2.884
19.65
4.123
3.74E−05
120


2
rs10203465
81260623
C
1748
7.519
0.4849
2.907
19.45
4.161
3.17E−05
121


2
2-15840134
15840134
C
1709
7.507
0.5177
2.722
20.71
3.894
9.86E−05
122


14
14-98055995
98055995
A
1742
7.502
0.4786
2.936
19.17
4.211
2.55E−05
123


14
14-98056427
98056427
C
1742
7.502
0.4786
2.936
19.17
4.211
2.55E−05
124


1
rs2805873
57713578
A
1713
7.499
0.4683
2.995
18.78
4.302
1.69E−05
125


2
2-15833678
15833678
G
1724
7.499
0.5177
2.718
20.69
3.891
9.96E−05
126


7
7-54624309
54624309
G
1754
7.487
0.4964
2.83
19.81
4.056
5.00E−05
127


7
rs11982678
54625449
A
1754
7.487
0.4964
2.83
19.81
4.056
5.00E−05
128


2
2-21321893
21321893
A
1729
7.413
0.4883
2.847
19.3
4.102
4.09E−05
129


2
rs7602629
21323199
A
1729
7.413
0.4883
2.847
19.3
4.102
4.09E−05
130


14
14-98042176
98042176
G
1743
7.407
0.4779
2.903
18.9
4.19
2.79E−05
131


3
rs12636997
29874808
G
1712
7.351
0.5042
2.737
19.75
3.957
7.60E−05
132


14
rs12431810
98044037
A
1744
7.33
0.4777
2.874
18.7
4.17
3.05E−05
133


8
8-91414191
91414191
A
1707
7.312
0.4836
2.834
18.86
4.114
3.88E−05
134


2
rs312047
21328107
A
1686
7.236
0.4866
2.788
18.78
4.067
4.76E−05
135


20
rs2297429
61699516
T
1726
7.208
0.506
2.673
19.43
3.903
9.49E−05
136


3
3-29883350
29883350
A
1751
7.182
0.4669
2.876
17.93
4.222
2.42E−05
137


2
2-81311658
81311658
C
1750
7.18
0.4836
2.783
18.52
4.077
4.57E−05
138


3
rs12638932
29882677
T
1749
7.171
0.4669
2.872
17.91
4.219
2.45E−05
139


3
3-29881826
29881826
C
1748
7.166
0.4669
2.87
17.89
4.218
2.47E−05
140


2
rs1986414
201199695
A
1711
7.164
0.4573
2.924
17.56
4.306
1.66E−05
141


3
3-29884382
29884382
C
1758
7.126
0.4674
2.851
17.81
4.202
2.65E−05
142


2
rs2219366
81318138
A
1748
7.118
0.4834
2.76
18.36
4.06
4.90E−05
143


19
19-55069295
55069295
C
1732
7.111
0.5022
2.658
19.03
3.906
9.37E−05
144


9
9-17634306
17634306
T
1755
7.102
0.5027
2.651
19.02
3.9
9.63E−05
145


10
rs1886951
2425384
A
1686
7.085
0.456
2.899
17.32
4.294
1.75E−05
146


19
19-55072907
55072907
T
1735
7.054
0.502
2.637
18.87
3.892
9.95E−05
147






















TABLE 5










SEQ






AL-
SNP
ID



CHR
SNP
LOCATION
LELE
ID#
NO:
Sequence





















22
rs5768434
46977516
T
 40
40
TTTGTATCTAATGTTTATATTcttta[C/T]gcttgaaactatttgaagcattcta


 


12
rs11064477
6944626
A
 41
41
AGCCTTCCCAGCTTGGGTCCTGTTTC[A/G]GAGCCCAGGCCTTGCTTCCCCTTAG


 


 1
rs1886629
194421521
C
 44
42
GAATCTCCAAGTTAATTTTTTCATTG[C/T]TTACCCTTTATTACTGCATGAAGTC


 


 2
rs7588295
166115757
G
 45
43
ACTTTTTTTTTTCAATTCTGAAATCA[A/G]CTTCTTAGATCTACTTTTCCTTACC


 


 4
rs4431170
165504024
G
 46
44
TGAAACCAAACTATAAACTACTTCTT[A/G]TCTTTGTGAGAGAATTCCAGGGCAC


 


 6
rs7740004
120897902
A
 47
45
TACAACTCTCAAGAAGCAATCTTTTT[A/T]TGGCAGGAAGCTGCATGAGTGATAA


 


10
rs11189381
99553188
C
 48
46
GCTGGCTTCTTCTTAATCAGAACTGC[C/T]CTTAGCTTCAAGAGAGGCTGGAAAA


 


15
rs12903202
56094085
G
 49
47
CACTGATGCTAAGGCAAGAGTTATCG[A/G]TAGCCCAAGCTCTGGCTAGTATCCA


 


18
rs17751934
47455812
T
 50
48
ATAGAATAAATAAGGTAACAGGTAAT[C/T]TAAACAAAGAGAAAAATATCTATTT


 


 3
rs7613687
60124976
T
 52
49
GGGAAATGAATCCCTGGCAGGATGGA[C/T]ATTAACCAAACACCTAACACATCAG


 


 3
rs7634338
8864809
C
 53
50
GGGAATATATGGTATTAATGAAAGCA[C/T]TGCTACACCTAGGTACCCCATACCT


 


 3
rs7612449
181586884
T
 54
51
TTAAGTTATTGCTATCCTATCCTTCT[C/T]CAAAGAGCTTTTCTTCTGTTTTTTA


 


 4
rs10029016
8496903
A
 55
52
gtggaaggatggttaccagagggtgg[A/G]aaggggagtggagggttgggggaga


 


12
rs10505722
1102632
G
 56
53
AATGATAAGCTTAAGAAAATACAAGC[A/C/G/T]TAAGTATGGCTCTTAATACGCTAGA


 


14
rs2332834
71952831
C
 57
54
GTCTCCTTTTCTTCTTGTCCCCTGCT[C/T]TCAGCTCTGAATGAGAAAAGTTTTG


 


 7
rs10271074
78589212
G
 58
55
TATAAAGTCATCCTACTTTCTCTTTC[A/G]TACTTTCAAAAGTTGGTATTCAGTA


 


 1
rs4951362
202580842
C
 61
56
GGCGAGAAAGACATATGTGGATGGAA[C/T]ATTTCAGAACAACTTGTATTTCCAA


 


 2
rs62153910
96382119
T
 62
57
CTACACTTGGGGTCCCAGAAGAGCTG[C/T]CCCCTCCTGCACATTTCCAATGCAG


 


15
rs7176436
56752223
G
 68
58
TTATGAAAAGGCTGTGAAGCTGAAGA[A/G]AAACTAAGAAATGGATATTGCTGCA


 


 1
rs17127107
65099496
G
 70
59
GTGGTTCATCTCATAGCAGACTTGCT[C/G]TAGAAGGTGAAACTCCCGGATTTTA


 


 2
rs62156621
79426168
A
 72
60
CACTCACACTTGTCCAGAGTTTCCTC[A/T]TGCTTTCTACTCTTCTGAAATCCTG


 


 2
rs921245
137980811
C
 73
61
TAGCTTCCACAACATTCCCAGGCTAC[A/G]AGAGCTTACAGTCCATTAGCACTGA


 


 1
rs1736563
169371174
A
 77
62
TTAAATCACAAATGCAGTCTCAATCC[A/G]GAAAATAGATCCCATCATATGTGAT


 


12
rs4639998
85442168
T
 85
63
tctttctcatatttatgtgcttcctt[C/T]aggagctctagtaaggcaggtctgg


 


 4
rs7669796
160863851
G
 86
64
gccaaggggtgcctgcaggcccatgc[G/T]gagccaccctcagacctccttcagc


 


12
rs17653326
85528244
A
 87
65
AGGAGGATATATAACCCTGGCTTGAA[A/G]AAGATGGAGATAGCTACAAGAGATG


 


 9
rs17057133
73509497
C
 88
66
GTTGGACGGATACTCATCTCGTGTAA[C/T]GGTCATAGAAAGATTCTGAGTGCTC


 


16
rs3135009
12556820
T
 89
67
CTCATGTACGCAGGGTGTTTCCCTAG[C/T]TGACATGTCTGAGGACGTCTTTCGT


 


 2
rs2544530
15820699
G
 91
68
TTCCAAACCCGGCTTTCCCTCATTTG[A/G]TCCTCAAACACCCCATCGGGGGGCC


 


 8
rs16915413
94070482
A
 92
69
TAGCCATTGACAAACCCTGTAGACAT[A/G]AGAATTTAATATGTGATAAAGATAA


 


16
rs2856790
12557384
C
 98
70
CACTTGTCCTGTGGACTCATGCCATG[C/T]CACCGTAGTGCTGAGTGACGCTTAA


 


 3
rs6443519
179112920
A
102
71
ttcatgctaaaaactctgaataaact[A/C]ggttttgatggaatgtatctcaaaa


 


 2
rs7568908
79462538
T
106
72
TAGCAATGCAAGCATGTTTGCCTTCA[G/T]ATAGCTAAATGACTGCgtgattgct


 


 2
rs57446887
21316772
C
120
73
TTATTCCTATGCTTGTTGGCTGCATG[C/T]ATGTCTTCTTTTGAAAAGTGTATCT


 


 2
rs10203465
81260623
C
121
74
tggtaataattacaagaaaaatgtct[C/G]tacatgttcagtataaatgcaacct


 


 1
rs2805873
57713578
A
125
75
TGGCTTCCTGGGTGCAACATCCGGAC[C/T]GATGTCTTTTATTGTTGTTATTGtt


 


 7
rs11982678
54625449
A
128
76
tatacacacacatataATTAATGCTA[A/G]AAGGCTATACACAGGAACACTATTG


 


 2
rs7602629
21323199
A
130
77
gaggcctcaggaagcttttactgatg[A/G]ctgtaggcaaagtgagagcaggcac


 


 3
rs12636997
29874808
G
132
78
AAAAAATTTCATCTTATATGTAGTAC[A/G]TAGTAATCTATAAATATAAAATACA


 


14
rs12431810
98044037
A
133
79
ATGTTGAGATAGTAGCACAGGAGGCC[A/C]GGGATTTATTTGGGTTATACACAGA


 


 2
rs312047
21328107
A
135
80
AAGAGAGAATCAAAAGGCAGGTCCTC[A/G]CAGCAGGCGCTGGGACATCTGTATC


 


20
rs2297429
61699516
T
136
81
ACAAGAAAGCAAGAGCTGCCAGGGCC[C/T]CTTCCAGCAGGGAGGCTGACCCTGC


 


 3
rs12638932
29882677
T
139
82
AAGCTTCTGactgttaagggtgatgg[A/T]tatgttcattatcttgacaatggtg


 


 2
rs1986414
201199695
A
141
83
CCTGTGCAAATGAAATGCTCATCCCC[A/G]CAAAGAAGGAATATGGGGCTGGCAG


 


 2
rs2219366
81318138
A
143
84
CACCAAACACACATACATACACATTT[A/G]TAGCATTTTGGAGCTAGAAAAGCT


 


10
rs1886951
2425384
A
146
85
CAGCAGCCACTCCTGGCAGAACCCCT[C/T]CTCCCATGCCAGCCACCCTCTTGAG









5.2 Methods of Treatment/Diagnosis
The present invention provides for a method of determining whether or not a human subject is at increased risk for developing ARONJ comprising determining whether nucleic acid, e.g. genomic DNA, of the subject carries an allelic variation or SNP that is a ARONJ biomarker, as described above, for example, by determining the nucleotide sequence of at least a portion of the DNA of the subject and determining whether a variant SNP which is a ARONJ biomarker as set forth above is present, where the presence of said ARONJ biomarker indicates that the subject is at increased risk for developing ARONJ. If a subject is at increased risk of developing ARONJ, the subject may be cautioned/warned to abstain from use of an anti-resorptive agent, such as a BP or Denosumab, avoid dental surgical procedures, or be frequently examined for early detection of ONJ (and then cessation of BP treatment or other therapeutic measures).
In non-limiting embodiments, the present invention provides for a method of determining whether a human subject is at increased risk for developing ARONJ comprising determining whether the genomic DNA of the subject carries an allelic variation of a gene selected from the group consisting of the human RBMS3 gene, the human IGF1R gene, the human IGFBP7 gene, the human DPYD gene, the human ABCC4 gene, or the human GSTM2 gene, where the presence of said allelic variation indicates that the subject is at increased risk for developing ARONJ. In particular non-limiting embodiments, the allelic variations may be in regions or subregions of those genes as described above, in SNPs represented in TABLES 1, 2 and 3, or in SNPs rs11934877, rs17761305, rs10875055, rs1189437, or rs673151.
In non-limiting embodiments, the present invention provides for a method of determining whether a human subject is at increased risk for developing ARONJ comprising determining whether the genomic DNA of the subject carries a SNPs set forth in TABLE 4 (i.e., any one of SNP ID #40-147 or combinations thereof), where the presence of the indicated allele of the SNP indicates that the subject is at risk for developing ARONJ.
The present invention further provides for a method of treating a human subject suffering from osteoporosis, comprising: (i) determining whether the human subject is at increased risk for developing ARONJ comprising determining whether the genomic DNA of the subject carries an allelic variation that is a ARONJ biomarker, as described above, for example, by determining the nucleotide sequence of at least a portion of the DNA of the subject and determining whether a variant SNP which is a ARONJ biomarker as set forth above is present, where the presence of said ARONJ biomarker indicates that the subject is at increased risk for developing ARONJ; and (ii) if the subject is at increased risk for developing ARONJ, recommending that the subject not be treated with an anti-resorptive agent or be treated with an anti-resorptive agent with relatively lower incidence of ONJ and optionally recommending an alternative treatment for osteoporosis, such as, but not limited to, calcium supplementation, exercise, and/or (iii) if the subject is at increased risk for developing ARONJ, recommending that the subject not undergo bone invasive dental procedures; and/or (iv) if the biomarker studies suggest that the subject is not at increased risk of ARONJ, initiating or continuing anti-resorptive therapy.
In a non-limiting embodiment, the present invention provides for a method of treating a subject suffering from osteoporosis, comprising obtaining the sequence of a portion of nucleic acid collected from the subject to determine whether the subject carries one or more single nucleotide polymorphism indicative of an increased risk of anti-resorptive therapy-associated osteonecrosis of the jaw selected from the group consisting of the single nucleotide polymorphisms set forth in TABLES 1, 2, 3 4 and combinations thereof, and if the one or more single nucleotide polymorphism is absent, treating the subject with an anti-resorptive agent.
The present invention further provides for a method of treating a human subject suffering from hypercalcemia, comprising: (i) determining whether the human subject is at increased risk for developing ARONJ comprising determining whether the genomic DNA of the subject carries an allelic variant that is a ARONJ biomarker, as described above, for example, by determining the nucleotide sequence of at least a portion of the DNA of the subject and determining whether a variant SNP which is a ARONJ biomarker as set forth above is present, where the presence of said ARONJ biomarker indicates that the subject is at increased risk for developing ARONJ; and (ii) if the subject is at increased risk for developing ARONJ, recommending that the subject be treated with a BP that carries a relatively lower risk of ARONJ or that the subject not be treated with BP and optionally recommending an alternative treatment for hypercalcemia, such as, but not limited to, gallium nitrate, plicamycin (formerly mithramycin), calcitonin, hemodialysis or peritoneal dialysis and/or (iii) if the subject is at increased risk for developing ARONJ, recommending that the subject not undergo bone invasive dental procedures and/or (iv) if the biomarker studies suggest that the subject is not at increased risk of ARONJ, initiating or continuing BP therapy.
In certain non-limiting embodiments, the present invention provides for a method of determining whether a human subject is at increased risk for developing bisphosphonate related osteonecrosis of the jaw comprising determining whether the genomic DNA of the subject carries one or more allelic variation which is a SNP selected from the group consisting of the SNPs set forth in TABLES 1, 2, 3 4 and combinations thereof, where the presence of said allelic variation indicates that the subject is at increased risk for developing ARONJ. In certain non-limiting embodiments, the SNPs include one, two, three or four or at least two, at least three, or at least four of the following: rs17024608 wherein A is substituted by G; rs17761305 wherein C is substituted by T; rs11934877 wherein T is substitued by C; rs10875055 wherein C is substitued by T; rs1189437 wherein A is substitued by C; rs673151, wherein G is substitued by A; and combinations thereof.
5.3 Kits
The present invention provides for kits that may be used to practice the above methods for determining whether the genomic DNA of a subject carries an allelic variation that is a ARONJ biomarker.
In certain non-limiting embodiments of the invention, a kit may comprise one or more primer nucleic acid having a sequence that is complementary to a nucleotide sequence containing or in proximity to the location of a SNP that is a ARONJ biomarker, as described above. For example, and not by way of limitation, the primer may be extended across the sequence having a SNP.
A “primer” as that term is used herein is a polynucleotide that is at least 8 nucleotides, at least 10 nucleotides, at least 15 nucleotides, at least 20 nucleotides, or at least 25 nucleotides in length (and may be, for example but not limitation, up to 20 nucleotides, up to 30 nucleotides, up to 40 nucleotides, up to 50 nucleotides, up to 100 nucleotides, up to 200 nucleotides, up to 500 nucleotides, up to 1000 nucleotides, in length) and, under reaction conditions, forms a hybrid structure with its target sequence, due to complementarity of at least one sequence in the probe or primer with a sequence in the target sequence. The target sequence, in non-limiting embodiments, may be at least 8 nucleotides, at least 10 nucleotides, at least 15 nucleotides, at least 20 nucleotides, or at least 25 nucleotides in length (and may be, for example but not limitation, up to 20 nucleotides, up to 30 nucleotides, up to 40 nucleotides, up to 50 nucleotides, up to 100 nucleotides, up to 200 nucleotides, up to 500 nucleotides, or up to 1000 nucleotides, in length. In non-limiting embodiments, the primer may be identical to the complement of its target sequence, may be at least 99 percent identical to the complement of its target sequence, may be at least 98 percent identical to the complement of its target sequence, or may be at least 95 percent identical to the complement of its target sequence, and may optionally be fused to a second nucleic acid or other molecule that is non-specific to the subject's nucleic acid but that is used in the detection assay (for example, for purification of extended or amplified primer).
The target sequence may span the location of the SNP or other allelic variation or may be in proximity to it (for example, but not by way of limitation, within up to 20 nucleotides, or up to 50 nucleotides, or up to 100 nucleotides, or up to 200 nucleotides, or up to 500 nucleotides, or up to 1000 nucleotides in genomic DNA or, if the SNP is located in a transcribed region, in RNA or cDNA).
In a specific non-limiting embodiment of the invention, the target sequence of a primer may be within 50 or within 100 nucleotides on either side of (i) the nucleotide that is the location of the SNP or (ii) the other allelic variation, in a genomic DNA, RNA or cDNA sequence.
One primer, as described above, may be used to generate a test fragment by primer extension. Two such primers may be used to generate a test fragment by polymerase chain reaction. The placement of the primer(s) is/are such that the test fragment comprises the location of the SNP or other allelic variation. The resulting test fragment may be sequenced to determine whether a SNP or other allelic variation of the ARONJ biomarker is present.
In non-limiting embodiments, a primer used according to the invention has a target sequence in (or is complementary to at least a portion of or is at least 90 or at least 95 or at least 99 percent homologous to (as determined by standard software such as BLAST or FASTA)) one of SEQ ID Nos 1-85.
In non-limiting embodiments of the invention, primer extension may be used to extend a primer as described above to generate a test fragment that comprises a nucleotide that is the location of a SNP that is a ARONJ biomarker. In non limiting embodiments, the sequence of the test fragment may be determined so as to determine what nucleotide is present at the site of the SNP.
In non-limiting embodiments of the invention, polymerase chain reaction may be used to amplify a test fragment between two primers (as described above), where the test fragment comprises a nucleotide that is the location of a SNP that is a ARONJ biomarker.
In non-limiting embodiments of the invention. a kit may comprise at least one, or at least two, or at least three, or at least four, or at least five, or at least six, or at least seven, or at least eight, or at least nine, or at least ten, and optionally up to five, optionally up to ten, optionally up to twenty, or optionally up to fifty, primer(s) selected from the group of:
primer(s) that may be used to generate test fragments that comprise a nucleotide that is the location of a SNP listed in TABLE 1, or in TABLE 2, or in TABLE 3; or in TABLE 4, and
primer(s) that may be used to generate test fragments that comprise a nucleotide that is the location of a SNP or allelic variation in the human RBMS3 gene, human IGF1R gene, human IGFBP7 gene, human DPYD gene, human ABCC4 gene, or human GSTM2 gene, including the regions and subregions described above, and including the SNPs associated particularly with each of these genes.
In a specific non-limiting embodiment, the kit comprises at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of the SNP of rs17024608, namely, the substitution of A by G in SEQ ID NO:1=GATAGAATAGAACTATTTGATGTGGA/GCCATGAG AATTTAAAAGTATCTGCA (bolded and underlined), where the presence of a G rather than A is indicative on increased risk of ARONJ, together with a package insert that describes the association between this SNP and ARONJ. Said kit may optionally comprise at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of a SNP of a human gene selected from the group consisting of IGF1R, UGFBP7, DPYD, ABCC4, GSTM2, SNP ID #40-147, or combinations thereof. For example, said kit may further comprise at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of the SNP of rs17761305 corresponding to position 57934091 of chromosome 4, as represented by the substitution of C by T in SEQ ID NO:20. As another example, said kit may further comprise at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of the SNP of rs11934877 corresponding to position 57941026 of chromosome 4, as represented by the substitution of T by C in SEQ ID NO:19.
In certain non-limiting embodiments of the invention a kit as described above comprises primers having, as target sequences, ARONJ biomarkers as set forth above, where said primers for ARONJ biomarkers are at least 20 percent, or at least 30 percent, or at least 40 percent, or at least 50 percent, or at least 60 percent, or at least 70 percent, or at least 80 percent, or at least 90 percent, of the primers present in the kit. In such embodiments, the kit is directed toward detecting ARONJ-associated markers and does not include a majority of primers that are not ARONJ-associated, although primers to serve as controls, for example, may be included (in non-limiting embodiments, the percentage of non-ARONJ associated primers may be up to 50 percent, or up to 40 percent, or up to 30 percent, or up to 20 percent, or up to 10 percent, or up to 5 percent).
In certain non-limiting embodiments, a kit for detecting ARONJ biomarkers may comprise at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of the SNP of rs17024608, namely, the substitution of A by G in SEQ ID NO:1=GATAGAATAGAACTATTTGATGTGGA/GCCATGAG AATTTAAAAGTATCTGCA (bolded and underlined), where the presence of a G rather than A is indicative on increased risk, together with one or at least one primer for another ARONJ biomarker set forth above.
In certain non-limiting embodiments, a kit for detecting ARONJ biomarkers may comprise at least one primer that may be used to generate a test fragment that comprises the nucleotide that is the location of the SNP of rs17024608, namely, the substitution of A by G in SEQ ID NO:1=GATAGAATAGAACTATTTGATGTGGA/GCCATGAG AATTTAAAAGTATCTGCA (bolded and underlined), where the presence of a G rather than A is indicative on increased risk, together with two or at least two primers for other ARONJ biomarkers set forth above, for example, but not limited, to one, two, three or four or at least two, at least three, or at least four of the following: rs17761305 wherein C is substituted by T; rs11934877 wherein T is substitued by C; rs10875055 wherein C is substitued by T; rs1189437 wherein A is substitued by C; rs673151, wherein G is substitued by A; and combinations thereof.
6. EXAMPLE
6.1 Materials and Methods
This research involves an observational, hospital-based, epidemiologic case control study. The research protocol was reviewed and approved by the institutional review boards (IRB) of the participating institutions. Human subject participation required the signing of a written informed consent, as approved by each institution's IRB. The study base of this case control study consisted of individuals who had received bisphosphonates, and who had received care in the clinics of the Massachusetts General Hospital (MGH), the Brigham & Women's Hospital (BWH), the Harvard School of Dental Medicine (HSDM) and its affiliated clinics, and the Nova University Dental School in Florida. In addition to providing their signed consent, individual participants were required to have the ability to answer to a questionnaire, and to provide a saliva sample. No individual was excluded from the study on the basis of gender, religion, political or sexual orientation, or minority group membership.
The identification of cases and controls occurred at the level of the recruiting clinic. Initially, electronic medical records and clinical notes were searched to identify bisphosphonate users. Among the bisphosphonate users, confirmed ONJ cases and unaffected controls were identified and invited to participate in the study by means of an introductory letter. Letters were mailed at the home address of record for all subjects. Three weeks after the initial mailings, a second wave of follow-up letters were sent, followed with telephone calls. Because of multiple co-morbidities, participants were offered the option to participate over the mail or to visit the clinic for an in-person session. Further, research visits were arranged to coincide with scheduled visits in the Oncology wards, when possible. To avoid misclassification of the disease among the controls (avoiding classifying persons with osteonecrosis as unaffected healthy controls) intra-oral examinations were performed on all controls that expressed a willingness to participate. Both cases and controls used the exact same research instrument and same saliva collection method. The research instrument contained questions on the following fields: demographics, including gender, race, ethnicity and availability of medical insurance; recent radiation to the head and neck (a positive answer would exclude them from the study); exposure to risk factors such as tobacco and alcohol use; co-morbidities, including cancer, osteoporosis and various autoimmune diseases; use of certain medications such a steroids, statins or thalidomide; details of their bisphosphonate use; having interventional dental procedures prior to ONJ, such as implant placement or dental extractions; and details on the osteonecrosis of the jaw, including symptomatology and recurrence. The saliva collection method utilized the Oragene DNA collection kit (DNA Genotek, Canada).
Following recruitment, the saliva kits were mailed in one batch to the genotyping facility. DNA was extracted following the manufacturer's recommended protocol. High throughput genotyping was performed using the Human Omni Express 12v1.0 Beadchip (Illumina, San Diego) according to the manufacturer's protocol. The Human Omni Express 12v1.0 Beadchip captures 731,442 markers, representing more than 91% of human variation for major alleles with frequencies above 5% in Caucasians.
All genetic data was imported at the Columbia University Medical Center Division of Bio-informatics computer cluster for statistical analysis. After converting the Final Report file format (standard export format from Illumina's GenomeStudio) into ped and map files, all downstream analyses was carried out in PLINK software [5].
To test the quality of the genotyping and to decide the call rate thresholds (both per sample and per marker), the missingness rate was checked by individual and by SNP (locus), respectively. All samples had call rates greater than 95%. 39,456 SNPs with MAF (minor allele frequency) less than 0.01 were excluded as alleles with such low frequency they would have no chance of approaching significance in this study. PLINK was used to test for Hardy Weinberg Equilibrium and SNPs were excluded that deviated from HWE at a p<0.0000001. Cryptic relatedness was tested for by estimating the identity-by-descent (IBD) for all possible pairs of individuals. To estimate the effect of population structure, the smartPCA program from the EIGENSTRAT package (version 3.0) [6] was used to conduct Principal Components Analysis (PCA) in order to expose population structure of the ONJ study group. This process was repeated when looking for additional genetically—matched population controls in publically available GWAS datasets. Subjects were selected from POPRES [7], Wellcome Trust Case Control Consortium [8], Illumina iControlDB [9], and the international Serious Adverse Events Consortium (iSAEC) [10]. All subjects except the ones from iControlDB were genotyped using Illumina 1M or 1M-duo chips, and the subjects from iControlDB were genotyped using Illumina 500K chip. SNPs from known regions of long-range linkage disequilibrium (LD) [11] were removed before conducting PCA.
The association of single SNPs were tested primarily using logistic regression with the PCA eigen values as covariates under an additive model. The Cochran-Mantel-Haenszel stratified test was also utilized. Both tests take into account the population structure to minimize inflation of test scores. Standard case-control association analyses set the significance p-value at the p<5×10-8 level, with the exception of a candidate gene analysis that used less stringent correction. The candidate gene sub-analysis focused on certain genes that were considered to be of interest, including the Insulin-like Growth Factor (IGF) family, and several ADME genes [12].
Subsequently the genotypes of 30 Caucasian ONJ cases and 1,743 controls from the “extended study group” (described in more detail in the Results section) were imputed using IMPUTE2 [13] (version February 2009), with data from the 1000 Genomes Project (112 individuals, release number March 2010) and HapMap III (June 2010, all ethnicity) as the reference panels. Only the imputed genotypes with posterior probability (reference) of greater than 0.9 were retained. All known SNPs with poor quality were pruned before the imputation to avoid false positives. The genome was divided into 5000 bp length segments and was imputed using ethic mixed panels to increment the quality of the imputation for rare variants. Stringent QC was carried out on the imputed genotypes. Copy number variations (CNVs) were subsequently inferred from SNP chip data using PennCNV software (April 2009 version) [14]. To ensure the accuracy of CNV calling, stringent sample and CNV filtering procedures were applied. All samples were included that had a LRR standard deviation <0.5, maximum number of total CNV calls <50, BAF median >0.55 or <0.45, BAF drift >0.01 or WF>0.05 or <−0.05 (default parameters). Additionally, to ensure high-confidence CNVs, individual CNVs with PennCNV-generated confidence score <10; those with calls based on fewer than 10 SNPs/CNV probes; and those with span within 1 Mb from centromeres or telomeres, were excluded.
6.2 Results
Recruitment. A total of 67 individuals were recruited in the period 2008-09. Of those, 32 were female cases with a mean age of 62.8 years; 15 were female controls with a mean age of 64.8 years; 5 were male controls with a mean age of 63.6 years and 15 were male cases with a mean age of 64.8 years. The majority of the cases (28/47) and controls (13/20) had received zoledronic acid, with an average duration of 22.5 months. The mean months on zoledronic acid was higher in cases than in controls, but the difference was not statistically significant. Similarly, there was no significant difference between cases and controls in mean months on zoledronic acid for the 14 subjects that reported a positive history of osteoporosis. Of the 67 individuals that participated in the study, we were able to extract DNA from 53 samples; 35 patients with osteonecrosis of the jaw and 18 treatment-tolerant controls.
Population structure and selection of genetically matched population controls. In order to identify the ethnicity of the members of the ONJ study group, the genotype data of the original study group, also referred here as “ONJ study group,” was combined with that of 987 HapMap III subjects, which include subjects from 11 populations. Six individuals were found to not cluster with the Caucasian (CEU and TSI) HapMap III samples, including three apparent African Americans, one possible Mexican and two subjects with mixed ethnicity. To increase sample size, additional genomic data was introduced for various classes of controls. More specifically, the genotypic results of the 47 Caucasian subjects (cases and controls) in the original study were merged with that of a selected European study group to separate sub-populations among Europeans. The members of the European study group were selected from publicly available datasets genotyped with the Illumina 1M or 1M-duo chips (POPRES, WTCCC and iSAEC study groups), representing sub-populations among Europeans. This analysis showed that the Caucasian subjects clustered with individuals of northwestern, southern and eastern European descent (FIG. 1A-B). Based on the eigen scores of the first two principal components, 1,122 genetically-matched population controls were selected to form the “combined study group”.
To further increase the number of genetically matched controls, especially subjects of Eastern European origin, the combined study group was merged with 2,978 Caucasian samples from the iControlDB dataset and PCA was used to cluster all samples together. From this merged study group, the closest controls for each case were selected based on the eigen scores of the first six principal components to form the “extended study group”. Finally, the Caucasian ONJ study group was merged with 101 treatment-tolerant cancer subjects from a dataset available via dbGAP to form a “treatment-matched study group” [15] that contains our 30 ONJ cases and a total of 118 treatment-tolerant controls; of the 118 controls, 101 come from the dbGap phs000210.v1 and 17 come from our original ONJ study. Considering recent evidence about the rarity of ONJ, the Extended study group was combined with 27 additional Eastern European subjects taken from phs000210.v1 study group to form the final study group.
Discovery Phase. The Caucasian ONJ study group contains data from 30 Caucasian cases and 17 Caucasian treatment exposed controls. In total, 631,507 SNPs passed quality control. Logistic regression was used to quantify the Odds Ratio and the 95% Confidence Interval of each SNP using the eigen score from the six significant components as covariate to control for population structure. Given the small sample size, single marker association was tested by the Fisher Exact test. TABLE 6, below, summarizes the findings from the genome wide association analysis. This process as repeated comparing the data from the 30 Caucasian cases with various sets of controls, as explained above. Statistical tests of association in all genome-wide association sub-analyses failed to reach statistical significance (genome-wide level of p<5×10-8). Several markers reached borderline significance at the p<5×10-7level, such as rs8012823 at DPF3; rs11802277 at AL365331.2; rs1075417818.1; rs11189381 at RP11-459F3.3; rs6861122 at AC022120.1; rs4431170 at MARCH1; rs13096022 at GRM7, and several SNPs appeared consistently in the lists of the 10 top associated genes across the various analyses, suggesting possible involvement in the etiology of ONJ. See Appendix for the 10 most associated SNPs of each analysis, along with ORs, p-values, QQ plots and Manhattan plots for each of the different case-control comparisons (caucasian ONJ study group; combined study group; extended study group; treatment matched study group; and, final study group).
Imputation Analysis. Prior to statistical evaluation a quality control routine was followed that included first using the -png-miss option to fill in the missing genotype and then calculating the Average Posterior Probability (APP) of imputed genotypes for each SNP. APP gives a broadly indication of the accuracy of imputed genotypes of a particular SNP. This parameter ranges from 0 to 1, where 1 means the complete certainty of the called genotype. SNPs with APP <0.9 were discarded. The difference of missingness was tested between cases and controls and the signal intensity was assessed manually using GenomeStudio to access the quality of genotype calls on subjects where the SNPs were genotyped.
Using imputation, 3,542,142 SNPs were analyzed and their association with the risk of ONJ was tested using logistic regression on the extended study group, with six eigenvalues as covariates to control for population structure. Two significant SNPs at the p<5×10-8 level, rs17425952 and rs233723, did not pass quality control, and were therefore discarded as non-significant. One significant SNP was found located in an intron of gene RBMS3, rs17024608, to be associated with an Odds Ratio of 5.4 (p=7.5×10-8). Rs17024608 was present in the genotype data of the combined study group, with a p-value of 7.3×10-6 (logistic regression). It was not present in the genotype data of the extended study group because it was pruned for not passing the missing rate threshold (>2% of missing data). Rs3821577, the best proxy of rs17024608 on the SNP chips (r2=0.14), showed association P-value of 1.0×10-5 (logistic regression). In the treatment matched study group, rs17024608 showed an association P-value=1.4×10-5 (OR=5.6, logistic regression). TABLES 7 and 8 summarize the genetic characteristics of rs17024608 and its statistical parameters across the various datasets. FIG. 4 presents the Manhattan plot of the region surrounding the SNP.
CNV Association analysis. Burden and common copy number variants association analyses were performed. Associations were tested using two tails permuted (×10,000) Fisher's exact test analysis using the PLINK software, by considering duplications and deletions separately. Singleton oversized CNVs larger than 700 kb were investigated to find evidence for individual predisposition to ONJ. All CNVs were excluded that had coverage smaller than 20 genetic markers/CNVs. All analyses were performed on the Caucasian subjects.
After the stringent QC, one subject was excluded from the association analysis. Fifty two individuals (33 cases and 19 controls) passed stringent quality-control criteria for CNV calling; 431 CNVs were called, of which 71 were duplications and 360 were deletions. Cases and controls did not differ significantly in their rate of CNVs for both deletions and duplications. After multi-test correction, none of the common CNVs had a significant association. However, two unique oversized (greater than 700 kb) duplications were found in cases, and none in controls (TABLE 9). The duplications were found on chromosomes 2 (925,407 bp) and 22 (730,236 bp) respectively.
Candidate Genes Analyses. Genes of the IGF gene family and genes related to drug-Absorption, Distribution, Metabolism and Excretion (ADME) were considered of interest at the inception of this study. With regard to the IGF gene family, 1,083 SNPs located within 20 kb of the candidate genes were available in the treatment-matched study group. FIG. 5A shows the QQ-plot for this set of SNPs. The most significantly associated SNPs were rs11934877 (OR 4.128; 95% CI: 1.918-8.885, p-value=0.0002), intronic to IGFBP7 (FIG. 5B. 2,564 SNPs close to genes related to drug-Absorption, Distribution, Metabolism and Excretion (ADME) [16] were also extracted. FIG. 6A shows the QQ-plot from those SNPs. The most significantly associated SNP was rs10875055 (OR=4.324; 95% CI: 1.999-9.353, p value=0.0002001), which is intronic to DPYD (FIG. 6B). The same two panels of markers were retrieved for the extended study group logistic regression results. Data on 867 candidate IGF SNPs and 1,247 ADME SNPs were available for analysis. FIGS. 7A and 8A show QQ plots respectively for the TOP candidate and ADME SNP panels. Top associated candidate SNPs were rs11934877 (OR=2.95; 95% CI: 1.66-5.25, p-value=0.00022) intronic close to gene IGFBP7, rs17761305 (OR=2.88; 95% CI: 1.64-5.07, p value=0.00023) intergenic SNP close to IGF1R (FIG. 7B). Top associated ADME SNPs were rs1189437 (OR=4.64; 95% CI: 2.06-10.46, p value=0.00021), intronic within ABCC4 (FIG. 8B). The same two panels of markers were analyzed in the Final study group (FC) via logistic regression. Data on 967 IGF candidate SNPs and 1,382 ADME SNPs were used. The top associated SNP was rs673151, intronic in GSTM2 gene, with a OR=3.57; 95% CI: 1.77-7.21, p-value of 0.00035. FIGS. 9 and 10A show QQ plots respectively for the IGF candidate- and ADME SNP panels.







TABLE 6







Genome Wide Association Study Results (ten most significant SNPs) for Caucasian ONJ Study Group




















ancestral



SNP
P
OR (%95)
chrom
Coordinate
type
allele
closest gene


















rs8012823
5.66E−07
0.09
(0.03-0.24)
14
73273544
INTRONIC
T
DPF3


rs11802277
7.38E−07
33
(4.2-257.1)
1
1.18E+08
DOWNSTREAM
G
AL365331.2


rs6900513
1.17E−05
0.12
(0.044-0.32)
6
66333105
INTRONIC
C
EYS


rs10781262
1.80E−05
0.12
(0.046-0.33)
9
77650696
INTERGENIC
A
C9orf41


rs1471646
1.93E−05
11.05
(3.04-40.08)
1
1.99E+08
WITHIN_NON_CODING_GENE
C
RP11-382E9.1


rs4870310
2.29E−05
0.11
(0.040-0.31)
6
1.55E+08
DOWNSTREAM
T
RP11-15G8.1


rs10875148
2.42E−05
23.57
(3.02-183.9)
1
98859508
INTERGENIC
G
AL160056.1


rs4562759
2.65E−05
8.7
(2.95-25.63)
10
83835639
INTRONIC
T
NRG3


rs4878512
2.71E−05
9.16
(2.87-29.26)
9
27630418
INTERGENIC
A
RP11-53518.1


rs1873291
3.18E−05
0.1423
(0.055-0.36)
6
66344855
INTRONIC
C
EYS
















TABLE 7







MAF of rs17024608 among different groups of subjects










caucasian selection (# samples)
MAF







cases (30)
0.2833



controls (1743)
0.08226



exposed controls (118)
0.06356



EE exposed and general population controls (122)
0.02



EE general population controls
0.015

















TABLE 8







OR (CI ± 95) of the top hit in all analyses.










COHORT
Model
OR 95% CI
P value





EXTENDED (30 vs 1743)
ADD
5.371 (2.8-10.3)
4.25E−07


DRUG EXPOSED (30 vs
ADD
7.102 (2.731-18.47)
5.80E−05


112)


IMPUTATION (30 vs 1743) *
ADD
5.825 (3.065-11.07)
7.47E−08



DOM
6.483 (2.966-14.17)
2.79E−06



REC
23.79 (5.518-102.5)
2.13E−05





* On the imputation dataset the association of rs17024608 under dominant and recessive models was tested.













TABLE 9







CVN Details on Large Duplications (found solely in cases)
















FIID
CHR
BP1
BP2
TYPE
SCORE
SITES
length
Start SNP
End SNP



















ojn5111
2
132144891
133070297
DUP
53.9
102
925,407
rs850234
rs16837705


ojn1304
22
19063495
19793730
DUP
360.7
205
730,236
rs6003971
rs2845421









6.3 Discussion
Osteonecrosis of the jaw is a serious adverse effect of bisphosphonates, especially among cancer patients on zoledronic acid. For this vulnerable group, osteonecrosis of the jaw adds yet another burden to their already compromised health, negatively affecting their quality of life. A test able to screen subjects for ONJ susceptibility prior to initiating bisphosphonates would have a great clinical utility as it would reduce the incidence of osteonecrosis. The present study identified SNPs in the genes IGF1R, IGFBP7, DPYD, ABCC4, and GSTM2, one or more of which may be used in such a test.
IGFs, especially IGF1 with its tyrosine kinase domain, are growth factors with potent signal transduction capabilities. Insulin like growth factors are molecules with important roles in normal growth and development. IGF1-deficient children fail to achieve appropriate height and pharmacologic therapies now exist to correct such deficiencies [17]. IGF1 and IGF2 are able to influence the replication and differentiation of bone cells through activation of their receptors, especially IGF1R, which plays a role in the cell cycle [18-19]. However, IGF2R seems to have a pro-apoptotic effect since it binds IGF2 and thus reduces available ligand levels for IGF1R [20]. IGF-binding proteins (IGFBPs), produced by bone cells, compete with the receptors in binding the ligands and thus affect the bioavailability of IGF1 and IGF2. IGFBP-4 binds IGFs and blocks their action, whereas IGFBP-5 promotes the stimulatory effects of IGFs [21].
DPYD polymorphisms have been associated with fluorouracil toxicity, especially with bone marrow and gastrointestinal toxicity, mucositis and leucopenia. DPYD has also been linked to autism spectrum disorders, Barett esophagus and adenocarcinoma [22-24]. A PubMed search on “DPYD gene and bone” or “DPYD gene and necrosis” returned no results; this genotype may be potentially involved in the soft tissue aspect of osteonecrosis' pathogenesis.
ABCC4, an ATP-binding cassette transporter gene, codes for Multidrug Resistance Protein, (MRP4/ABCC4) a transporter that actively effluxes endogenous and xenobiotic substrates out of cells. Inherited variation in ABC transporters has been associated with the occurrence of serious adverse effects. For example, ABCC4 has been linked to cyclophosphamide-induced adverse drug reactions in breast cancer patients, especially leucopenia/neutropenia. Currently there is no published information on MRP4 in bone phenotypes. [25-26].
GSTM2 codes for the phase II detoxifying enzyme glutathione-s-transferase, an enzyme that protects cells against toxic insults and enhances cell survival. While a direct effect on bone necrosis has not been described in the published literature, Owur and Kong (Biochem Pharmacol, 2002) have raised the hypothesis that increased concentrations of certain xenobiotics leads to gst-mediated apoptosis, with extreme increases in concentrations leading to rapid cell necrosis. [27] Knowing that bisphosphonates are attracted to areas of active bone loss or trauma, it is then plausible that localized spikes in bisphosphonate concentrations would lead to gst-mediated cell toxicity.
The high throughput analysis of the present study was able to identify only one strong signal in RBMS3, and several weak signals in various genes. RBMS3 is a binding protein that belongs to the c-myc family of genes. The protein is located in the cytoplasm and it has two RNA binding domains. It has been shown to bind Prx1, a homeobox transcriptional factor that upregulates collagen 1 [28]. The effect of the specific genotype in the etiology of osteonecrosis is currently unknown; however, it is plausible that RBMS3 rs17024608 may be involved in reduced collagen formation and the disruption of tissue repair.
7. REFERENCES


    
    
        1. Cartsos V M, Seamanduras A, Koo S, Zavras A I. Implications of bisphosphonate use in dentistry. Analecta Periodontologica 2009; 20:1-25
        2. Rodan G A, Fleisch H A: Bisphosphonates: mechanisms of action. J Clin Invest 1996; 97:2692-6.
        3. Fournier P, Boissier S, Filleur S, Guglielmi 5, Cabon F, Colombel M, Clezardin P. Bisphosphonates inhibit angiogenesis in vitro and testosterone-stimulated vascular regrowth in the ventral prostate in castrated rats. Cancer Res 2002; 62:6538-6544.
        4. American Association of Oral and Maxillofacial Surgeons (AAOMS). Position paper on bisphosphonate-related osteonecrosis of the jaw. Can be found electronically at: http://www.aaoms.org/docs/position_papers/bronj_update.pdf (accessed last on Feb. 10, 2011)
        5. Purcell, S., et al., PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 2007. 81(3): p. 559-75.
        6. Price, A. L., et al., Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet, 2006. 38(8): p. 904-9.
        7. Nelson, M. R., et al., The Population Reference Sample, POPRES: a resource for population, disease, and pharmacological genetics research. Am J Hum Genet, 2008. 83(3): p. 347-58.
        8. 3. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature, 2007. 447(7145): p. 661-78.
        9. http://www.illumina.com/science/icontroldb.ilmn
        10. Y., et al., Genome-Wide Association Study of Serious Blistering Skin Rash Caused by Drugs. Pharmacogenomics Journal, 2010 (In press.)
        11. Novembre, J., et al., Genes mirror geography within Europe. Nature, 2008. 456(7218): p. 98-U5.
        12. Ahmadi, K. R., et al., A single-nucleotide polymorphism tagging set for human drug metabolism and transport. Nature Genetics, 2005. 37(1): p. 84-89.
        13. Marchini, J., et al., A new multipoint method for genome-wide association studies by imputation of genotypes. Nature Genetics, 2007. 39(7): p. 906-913.
        14. Wang, K., et al., PennCNV: An integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Research, 2007. 17(11): p. 1665-1674.
        15. Ingle J N, Schaid D J, Goss P E, Liu M, Mushiroda T, Chapman J A, Kubo M, Jenkins G D, Batzler A, Shepherd L, Pater J, Wang L, Ellis M J, Stearns V, Rohrer D C, Goetz M P, Pritchard K I, Flockhart D A, Nakamura Y, Weinshilboum R M. Genome-wide associations and functional genomic studies of musculoskeletal adverse events in women receiving aromatase inhibitors. J Clin Oncol. 2010 Nov. 1; 28(31):4674-82. Accessed with permission via dbGAP phs000210.v1.p1
        16. Ahmadi, K. R., et al., A single-nucleotide polymorphism tagging set for human drug metabolism and transport. Nature Genetics, 2005. 37(1): p. 84-89.
        17. Patient Information Sheet Mecasermin [rDNA origin] (marketed as Increlex). Can be found electronically at the FDA web site: http://www.fda.gov/cder/drug/InfoSheets/patient/mecaserminPIS.htm
        18. Gordeladze J O, Reseland J E, Drevon C A. Pharmacological interference with transcriptional control of osteoblasts: a possible role for leptin and fatty acids in maintaining bone strength and body lean mass. Curr Pharm Des 2001; 7(4):275-90
        19. Andreassen T T, Oxlund H. The effects of growth hormone on cortical and cancellous bone. J Musculoscel Neuronal Interact 2001; 1:49-58
        20. Pollak M N, Schernhammer E S, Hankinson S E. Insulin-like growth factors and neoplasia. Nat Rev Cancer 2004; 4(7):505-18.
        21. Canalis E. Skeletal growth factors. In; Marcus R, Feldman D., Kelsey J. Osteoporosis. Academic Press, San Diego, Calif., pp 261-133
        22. Coate L, Cuffe S, Horgan A, Hung R J, Christiani D, Liu G. Germline genetic variation, cancer outcome, and pharmacogenetics. J Clin Oncol 2010; 28(26):4029-37
        23. Botelho N K, Schneiders F I, Lord S J, Freeman A K, Tyagi S, Nancarrow D J, Hayward N K, Whiteman D C, Lord R V. Gene expression alterations in formalin-fixed, paraffin-embedded Barrett esophagus and esophageal adenocarcinoma tissues. Cancer Biol Ther 2010; 10(2):172-9.
        24. Schwab M, Zanger U M, Marx C, et al. Role of genetic and nongenetic factors for fluorouracil treatment-related severe toxicity: a prospective clinical trial by the German 5-FU Toxicity Study Group. J Clin Oncol 2008; 26(13):2131-8.
        25, Tagami M, Kusuhara S, Imai H, Uemura A, Honda S, Tsukahara Y, Negi A. MRP4 knockdown enhances migration, suppresses apoptosis, and produces aggregated morphology in human retinal vascular endothelial cells. Biochem Biophys Res Commun 2010; 400(4):593-8.
        26. Low S K, Kiyotani K, Mushiroda T, Daigo Y, Nakamura Y, Zembutsu H. Association study of genetic polymorphism in ABCC4 with cyclophosphamide-induced adverse drug reactions in breast cancer patients. J Hum Genet. 2009 October; 54(10):564-71
        27. Owuor E D, Kong A N. Antioxidants and oxidants regulated signal transduction pathways. Biochem Pharmacol. 2002; 64(5-6):765-70.
        28. Fritz D, Stefanovic B. RNA binding protein BRMS3 is expressed in activated hepatic stellate cells and liver fibrosis and increases expression of transcription factor Prx1. J Mol Biol 2007; 371(3):585-595.
        29. Aghaloo et al., T. Oral Maxilofac Surg. 2011; 68(5):959-963.
    
    


Various publications are cited herein, the contents of which are hereby incorporated by reference in their entireties.

## Claims

1. A method of treating a human subject suffering from osteoporosis, comprising determining, from a nucleic acid sample from the subject, whether the subject carries single nucleotide polymorphism rs17024608 wherein A is substituted by G, where if the polymorphism is present, the subject is at increased risk for developing bisphosphonate-related osteonecrosis of the jaw relative to a subject that does not carry the polymorphism; and either (i) where the subject carries the polymorphism, treating the subject with an osteoporosis therapy that is an alternative to bisphosphonate treatment and not treating the subject with a bisphosphonate agent or (ii) where the subject does not carry the polymorphism, treating the subject with an osteoporosis therapy that is a bisphosphonate agent.

2. The method of claim 1, where the osteoporosis therapy that is an alternative to bisphosphonate treatment is calcium supplementation.

3. The method of claim 1 where the bisphosphonate agent is selected from the group consisting of alendronate, etidronate, ibadronate, risedronate, pamidronate and zoledronic acid.

4. The method of claim 1 where the bisphosphonate agent is selected from the group consisting of pamidronate and zoledronic acid.

5. The method of claim 1 where the bisphosphonate agent is zoledronic acid.