Introduction

Osteonecrosis (ON), also known as avascular necrosis and aseptic necrosis, is defined as bone cell death following a compromise of blood flow to the bone. ON is most common in the femoral head (i.e., hip) but can occur at any skeletal site (e.g., knee, shoulder, and ankle) [1, 2]. The majority of ON cases are secondary to trauma [1]. Non-traumatic ON can also occur, but the underlying pathology is unclear [1, 3]. In published literature, non-traumatic ON has been associated with a number of risk factors including corticosteroid use, alcohol consumption, immunosuppressive therapy, autoimmune diseases such as systemic lupus erythematosus and rheumatoid arthritis, hematologic/thrombotic disorders, malignancies and metabolic disorders such as diabetes mellitus, and renal failure [1, 35]. Patients who experience non-traumatic ON usually have more than one risk factor, which indicates the pathogenesis of non-traumatic ON is probably multifactorial [2].

The majority of studies to date have assessed risk factors for ON in specific diseases with corticosteroid use, e.g., systemic lupus erythematosus and organ transplantation [47]. Few studies have been conducted in a general population [3, 8]. The purpose of this study was to examine the incidence of ON, patient characteristics, and selected potential risk factors for ON in two general population health record databases in the UK: the General Practice Research Database (GPRD) and The Health Improvement Network (THIN) database.

Methods

Study population and databases

The GPRD database contains computerized information entered by approximately 450 general practitioners in the UK. Data on approximately 3.4 million active patients (total of approximately 13 million) are systematically recorded, anonymized, and sent to GPRD where the data are collated and organized for research purposes. Symptoms and diagnoses are coded using the Oxford Medical Information System (OXMIS) and the READ clinical classification system. The THIN database contains similar information entered by general practitioners in the UK and contains information on over six million patients from 358 general practice offices, including data on approximately 2.8 million active patients.

Only data from medical practices that passed quality control checks are included in the GPRD database [9]. Both the GPRD and THIN databases have been validated in several independent studies [1014]. Both databases use the READ classification to code specific diagnoses; a drug dictionary based on the MULTILEX classification is used to code drugs. Information collected in both of the databases includes patient demographics and records of primary care visits as well as diagnoses from specialist referrals, hospital admissions, and the results of laboratory, radiographic, and diagnostic tests. Prescriptions issued by general practitioners are also recorded. Practices selected from THIN did not contribute to the GPRD during the study period, thereby avoiding duplication of ON cases.

Each database was screened for all permanently registered adults (aged 18 years or older) from 1989 to 2003. ON was defined as a patient with a record of at least one of the READ codes listed in Table 1. For each identified case, the first record of ON during the period of data collection was considered the index date. Within each database, each case was matched to up to six controls with no record of ON. The matching criteria included age (± 5 years), sex, and medical practice (registered at the same practice at the index date of the case). The index date of each control patient was assigned the same date as the corresponding matched case. Cases and controls were required to have a minimum of 3 months (i.e., 91 days) enrollment prior to the index date.

Table 1 List of READ/OXMIS codes used for identifying osteonecrosis cases

The overall study design was a case–control study that combined information from each of the two databases (GPRD and THIN). Cases with a diagnosis of ON were further assessed by examining the free text fields with key search terms for each subject. After identifying all diagnoses of ON, the incidence of ON was computed over time, and analyses were carried out to explore potential risk factors for ON.

Statistical methods and analysis

Incidences were calculated using midyear population counts. Possible risk factors, selected a priori, were considered for inclusion based on a review of the potential risk factors previously cited in the published literature [1, 47, 15]. Risk factors of interest included (1) drug exposure (including systemic corticosteroids, immunosuppressants, anti-infectives, bisphosphonates, statins, and hormone replacement therapy (HRT; women only)) in the 2 years prior to the index date and (2) medical history (including hospitalization, referral or specialist visit, bone fracture, any cancer including hematological cancer, inflammatory bowel disease (IBD), gout, solid organ or bone transplantation, asthma, renal failure or dialysis, congenital or acquired hip dislocation, diabetes mellitus, osteoporosis, connective tissue disease, and osteoarthritis) in the 5 years prior to the index date. Each of these potential risk factors was separately entered into a regression model. Additionally, alcohol consumption was considered (depending on the proportion of subjects with data for this variable). Baseline demographic characteristics for cases and controls were compared.

Crude odds ratios (ORs) and 95% confidence intervals (CIs) were calculated for each risk factor in a univariate analysis using conditional logistic regression, comparing cases and controls. After excluding risk factors that had an insignificant OR or did not reach an overall 1% prevalence, a final, multivariable logistic regression model was derived.

Results

A total of 792 cases and 4,660 controls were included in the analysis, with 99% of cases having at least five matched controls. Fifty-three percent of the cases and 53.1% of the controls were female, with a mean age of 57.5 years among cases and 57.6 years among controls. Mean observation time was 8.9 person-years for cases and 9.4 person-years for controls. The most common site of ON was the hip, representing 75.9% of the cases (Table 2).

Table 2 Baseline characteristics of cases and controls

The age-adjusted annual incidence rates of ON by sex and the osteonecrosis incidence rates by sex and age cohort are shown in Figs. 1 and 2. Overall, the recorded incidence of ON increased over time from approximately 1.4/100,000 in 1989 to approximately 3/100,000 in 2003. The incidence of ON in women was higher than men from 1989 to 1998; however, from 1999 through 2003, there was a shift to men having a higher incidence than women (Fig. 1). Overall and for women, the incidence of ON increased with age. The incidence of ON in men remained constant from age 40 to 79 (around 2/100,000), increasing to 3/100,000 at age 80 years and older. From ages 18–59, men had a higher incidence than women; however, women 60 years and older had a higher incidence than men (Fig. 2).

Fig. 1
figure 1

Age-adjusted annual incidence rates by sex (GPRD and THIN research databases)

Fig. 2
figure 2

Osteonecrosis incidence rates by sex and age cohort (1989–2003). Incidence rates are weighted average of the annual sex- and age-cohort-specific incidence rates (GPRD and THIN research databases)

Table 3 shows descriptive statistics for each of the potential risk factors of interest. Drug exposure was captured over the prior 2-year period and classified a priori, as follows: None; Exposed (2+ prescriptions within 120 days in the previous 2 years); or Intermittent (all other possible exposure scenarios). In the study population, anti-infectives were the most commonly prescribed therapy (22.9% among cases and 15.3% among controls). Relevant medical history was captured for the previous 5 years. The most commonly reported disease condition was osteoarthritis in 21.7% of cases and 7.8% of controls. A large proportion of subjects were missing data for alcohol consumption (46.3% of cases and 51.2% of controls; Table 3), and it was, therefore, decided to exclude this variable from multivariable modeling (Tables 4 and 5).

Table 3 Potential risk factors of interest
Table 4 Multivariable logistic regression modeling: selected potential risk factors of osteonecrosis at any site
Table 5 Multivariable logistic regression modeling: selected potential risk factors for osteonecrosis of the hip (subset N = 4134)

Statistically elevated crude ORs were observed for bisphosphonates, systemic corticosteroids, immunosuppressants (intermittent only), anti-infectives, and HRT (exposed only; Table 4). Marked differences between cases and controls were noted for hospitalization, referral or specialist visit, and bone fracture in the prior 5-year period, and, despite a relatively small proportion of exposed subjects, the crude ORs were also significantly increased for cancer, IBD, gout, solid organ or bone transplantation, asthma, renal failure or dialysis, osteoporosis, connective tissue disease, and osteoarthritis (Table 4). Statistically significant risk factors for ON from the final multivariable logistic regression model were systemic corticosteroid use (intermittent and exposed), hospitalization, referral or specialist visit, bone fracture, any cancer, osteoporosis, connective tissue disease, and osteoarthritis (Table 4).

An additional analysis was performed in the subset of cases with hip ON and their matched controls because these represented a potentially more homogeneous population and also included the majority (75.9%) of the identified ON cases overall (Table 2). A total of 601 cases and 3,533 controls were included in the hip ON subset analysis. Approximately 54% of cases and controls in the hip ON subset were female with a mean age of 58.3 years. Statistically significant risk factors for hip ON from the adjusted multivariable logistic regression model were the same as the overall ON population except for the inclusion of immunosuppressant use (intermittent) and the exclusion of osteoporosis (Table 5).

Of recent interest is the use of bisphosphonates and a postulated association with osteonecrosis of the jaw (ONJ) [1619]. In our case–control study, only 4.4% of ON cases were bisphosphonate users within the previous 2 years (Table 3). Across all cases, only three had the jaw mentioned as the site of ON, and none of them had been exposed to bisphosphonates (Table 2).

Table 6 reports the type of bisphosphonate exposure for cases and controls in this study. Etidronate was the most common compound reported; this was the only oral bisphosphonate marketed for the treatment of osteoporosis in the UK in the early 1990s. Further, the distribution by type of bisphosphonate is overall consistent with market share in the UK during the study period. No cases or controls with intravenous bisphosphonate use were identified in this study. Exposure to bisphosphonates was not associated with an increased risk of ON in the adjusted model of all skeletal sites combined (Table 4) or in the adjusted model for the hip subset (Table 5).

Table 6 Types of bisphosphonates used by cases and controls within the previous 2-year study period

Discussion

From 1989 to 2003, in this study population, the observed incidence of ON ranged from approximately 1.4–3.0/100,000 within the combined GPRD/THIN dataset. The reason for the increased incidence over time is not known but could be due in part to the increasing use of more advanced radiographic techniques, especially MRI, that are more sensitive in detecting ON. Consistent with the literature [1, 5], we found that the hip was the most commonly diagnosed site for ON (75.9%) and that corticosteroid use was associated with a 3-fold increased risk for ON. Significant risk factors for ON at all skeletal sites combined did not differ substantially from those for ON of the hip. While we did not assess trauma specifically, bone fracture in the prior 5 years was associated with a 5.8-fold increased risk of ON at all skeletal sites both combined and at the hip.

As observed in other studies, a history of connective tissue disease or cancer were significant risk factors for ON. This may be confounded by the frequent use of corticosteroids in these populations [46, 20]. In addition, overall disease severity/morbidity may also contribute to a higher rate of ON in these populations [1, 4]. There were two risk factors that showed a risk reduction (70% with statin use and 60% with diabetes mellitus); however, neither was statistically significant and neither met the criteria for inclusion in the multivariable model.

Our study population was 53% female. This contrasts with previous findings that ON is more common in men in the general population (with the exception of systemic lupus erythematosus populations) [1]. In addition, the age of our study population ranged between 42 and 73 years (mean = 57.6 years; median = 59.0 years), which is older than previously reported in the literature [1, 21].

Although a history of osteoporosis in the prior 5 years was a significant risk factor in this study, bisphosphonate use was not. Only three cases had the jaw mentioned as the site of ON, and none of these had been exposed to bisphosphonates in the previous 2 years. In this study, there were no cases of ON with intravenous bisphosphonate use, which has been reported with ONJ in the treatment of multiple myeloma and metastatic carcinoma in the literature [1619]. It should also be noted that the study period was prior to the recent literature and recent awareness of ONJ. Given that prior bone fracture was the strongest risk factor observed in this study and that bisphosphonates are indicated for the prevention and treatment of osteoporosis that is often first identified after a fracture occurs, confounding by indication may explain the observation of bisphosphonate use and ON in the univariate analysis (elevated crude OR).

There are several limitations to this study. As with the use of any medical records database, misclassification bias is possible. The case definition was developed to include all available READ codes in order to minimize the likelihood that true cases of ON were missed (i.e., sensitive) and that cases were not falsely classified (i.e., specific). Some cases of ON may not have been recorded or diagnosed; the diagnosis of non-traumatic ON is difficult because the disease is silent until pain presents [1]. In general, cases of ONJ identified by dental professionals may not be consistently recorded in the medical records databases. Further, ONJ had no separate READ code during the study period. To address this limitation, possible cases were assessed from a review of the text fields for ON cases with any mention of “jaw.” Another limitation is that prescriptions written by specialists may not have been recorded by the general practitioner. The study design was based on an a priori selection of risk factors that have been previously cited in the literature [1, 47, 15] with particular focus on those that were highly correlated; therefore, this study may have excluded other potentially important risk factors.

In conclusion, using data from the UK GPRD and THIN databases, we found that significant predictors of ON at any skeletal site included use of systemic corticosteroids in the previous 2 years, hospitalization, referral or specialist visit, bone fracture, any cancer, osteoporosis, connective tissue disease, and osteoarthritis within the past 5 years. Bisphosphonate use was not a significant predictor of ON. This study aimed to provide a broader perspective on the descriptive epidemiology of ON risk factors than previous published studies. Studies utilizing more recent data may further elucidate the understanding of key predictors of ON.