Article Text

Extended report
Hormone replacement therapy and mid-term implant survival following knee or hip arthroplasty for osteoarthritis: a population-based cohort study
  1. D Prieto-Alhambra1,2,3,4,5,
  2. M K Javaid1,
  3. A Judge1,2,
  4. J Maskell2,
  5. C Cooper1,2,
  6. N K Arden1
  7. on behalf of the COASt Study Group
  1. 1Oxford NIHR Musculoskeletal Biomedical Research Unit, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Oxford, UK
  2. 2MRC Lifecourse Epidemiology Unit, University of Southampton, Southampton General Hospital, Southampton, UK
  3. 3GREMPAL Research Group, IDIAP Jordi Gol Primary Care Research Institute, Universitat Autònoma de Barcelona, Barcelona, Spain
  4. 4URFOA, Institut Municipal d'Investigacions Mèdiques (IMIM), Parc de Salut Mar, Barcelona, Spain
  5. 5RETICEF (Red Temática de Investigación Cooperativa en Envejecimiento y Fragilidad), Instituto Carlos III, Madrid, Spain
  1. Correspondence to Professor Nigel Arden, Oxford NIHR Musculoskeletal Biomedical Research Unit, University of Oxford, Windmill Road, Headington, Oxford, OX3 7LD, UK; nigel.arden{at}


Objectives Osteolysis and subsequent prosthesis loosening is the most common cause for revision following total knee arthroplasty (TKA) or total hip arthroplasty (THA). Hormone replacement therapy (HRT) could reduce osteolysis through its antiresorptive effects. We studied whether HRT use is associated with reduced revision rates in a community-based cohort of women undergoing TKA or THA for osteoarthritis.

Methods Female participants in the General Practice Research Database undergoing a primary TKA or THA from 1986 to 2006 were included. We excluded patients aged <40 years at the date of primary, and those with a history of previous hip fracture or rheumatoid arthritis. Women with at least 6 months of HRT were identified as HRT users. We further explored the associations among HRT use of ≥12 months, adherence (medication possession ratio) and cumulative use and revision risk. Cox models were fitted to model implant survival in years. Propensity score matching was used to control for confounding.

Results We matched 2700 HRT users to 8100 non-users, observed for a median (IQR) of 3.3 (1.5–6.1) years after TKA/THA. HR for HRT ≥6 months was 0.62 (95% CI 0.41 to 0.94), whereas HR for ≥12 months was 0.48 (0.29 to 0.78). Higher adherence and therapy duration were associated with further reductions in revision rates. Preoperative HRT appeared unrelated to implant survival.

Conclusions HRT use is associated with an almost 40% reduction in revision rates after a TKA/THA. These findings require replication in external cohorts and experimental studies.

  • Osteoarthritis
  • Epidemiology
  • Knee Osteoarthritis
  • Orthopedic Surgery

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.


Lower limb osteoarthritis (OA) is increasingly common worldwide, and accounted for more than 90% of the 153 000 total knee arthroplasty (TKA) and total hip arthroplasty (THA) procedures carried out in the UK in 2010.1 There is no cure for OA and TKA and THA are the most effective treatments for severe knee and hip OA.2 ,3 The rates of primary arthroplasties are increasing4 ,5 due to an increasingly elderly population in addition to an increasing prevalence of obesity.

Together with patient reported outcomes,3 implant survival constitutes the most important element in the evaluation of TKA/THA surgery: revision surgery offers poorer clinical results6 ,7 and is more costly than primary arthroplasty.8 According to the 8th annual report of the National Joint Registry for England and Wales, revision rates were estimated at 2.3% for THA and 2.2% for TKA after 3 years of follow-up.1 The main causes for failure after the first year are osteolysis and subsequent aseptic loosening, which account for approximately 75% and 40% of revision surgeries after THA and TKA respectively in this same report.1 These occur as a consequence of osteolysis around the prosthesis leading to its migration and, eventually, to function impairment and symptomatic failure.9

Strategies aimed at reducing periprosthetic osteolysis and the consequent bone loss and migration would seem a logical way to reduce arthroplasty failure and hence the need for revision. Antiresorptive agents, due to their antiosteoclastic activity, have been tested for this purpose: several randomised controlled trials10–13 but not all14 have suggested that bisphosphonates could be effective in improving surrogate outcomes. Further, a recent observational study carried out by our group has demonstrated that bisphosphonate use significantly improves implant survival.15

Hormone replacement therapy (HRT) has antiresorptive effects, and was widely used in the 1980s and 1990s, offering an excellent opportunity for research on the potential effects of this kind of drug on long-term implant survival wherever follow-up data for more than 5 years are available.

We took advantage of the existence of a validated source of such data in the UK within the General Practice Research Database (GPRD),14 and aimed to test whether HRT is associated with reduced rates of revision surgery following lower limb total joint replacement.


Study design

This is a population-based retrospective cohort study.

Study population and setting

We screened the GPRD to identify all women with a medical diagnosis code for primary THA or TKA from 1986 to the end of 2006. The GPRD comprises of computerised records of all clinical and referral events in both primary and secondary care in addition to comprehensive demographic information, medication prescription data, clinical events, specialist referrals, hospital admissions and their major outcomes in a sample of 6.5 million patients from 433 contributing practices, chosen to be representative of the wider UK population. The GPRD is administered by the Medicines and Healthcare products Regulatory Agency. Only practices that pass quality control are used as part of the GPRD database. Deletion or encoding of personal and clinic identifiers ensures the confidentiality of information in the GPRD. Data are stored using OXMIS and Read codes for diseases that are cross-referenced to the International Classification of Diseases (ICD-9) and we used these to identify primary TKA and THA. The used list of codes is available online as supplementary table 1, and shown to be valid to identify patients undergoing TKA/THA in GPRD data for the study period.4

Women were included in this study if they were aged at least 40 years at the time of primary surgery, and those with a medical diagnosis code for rheumatoid arthritis were excluded. We also excluded those with a history of hip fracture prior to primary hip arthroplasty as we could not accurately ascertain whether participants were undergoing THA for OA or for the hip fracture or its late complications. Using these criteria, we identified 24 733 patients.

Ascertainment of outcome

The main outcome of this study was time from primary THA/TKA to implant failure. Subjects with a revision surgery procedure were identified using ‘Read/OXMIS codes’.

Participants were followed up until revision date (outcome), loss to follow-up (death or transfer out of the primary care practice), end of year 10 following THA/TKA or end of study (31/12/2006), whichever came first.

Main exposure: HRT use

Participants who had been prescribed HRT by their general practitioner either for at least 6 months before implant revision with a high adherence (medication possession ratio (MPR) >80%, calculated as the number of daily doses received divided by the number of days of follow-up) or if at least six prescriptions had been filled in the first 6 months after treatment initiation were identified as HRT users. Conversely, HRT non-users were defined as those who did not fulfil any of these two criteria, as well as those who had their first prescription of HRT after implant revision.


The following potential confounders were considered in the analyses: age, gender, body mass index (BMI), joint replaced (hip/knee), year of joint replacement operation, recorded diagnosis of OA (yes/no), fracture presurgery (yes/no), calcium and vitamin D supplements, use of bisphosphonates, use of selective oestrogen receptor modulators, oral glucocorticosteroid therapy, smoking status and alcohol intake recorded closest to the date of the primary surgery, General Practitioner (GP) practice deprivation score (as defined by the Index of Multiple Deprivation), region of UK, comorbid conditions registered by the physician from the following list (asthma, malabsorptive syndromes, inflammatory bowel disease, hypertension, hyperlipidaemia, ischaemic heart disease, stroke, chronic obstructive pulmonary disease, chronic kidney failure, neoplasms, diabetes), and use of drugs which can affect fracture risk (proton pump inhibitors, antiarrhythmics, anticonvulsants, antidepressants, anti-Parkinson drugs, statins, thiazide diuretics, anxiolytics).

Statistical analyses

To address the issue of confounding by indication, we used propensity score matching methods.16 The propensity score represents the probability that a patient is an HRT user, and was estimated for the whole study population using multivariate logistic regression models.17 In such models, HRT use was a binary outcome, and all the covariates as listed above were introduced as potential confounders. We used multiple imputation (ICE procedure in Stata18) for important missing covariates (BMI, smoking and alcohol intake) for the propensity score models. We imputed the values of these three variables in 10 different datasets using the following covariates as potentially explanatory variables: fracture history, socio-economic status, year of TKA/THA, age at surgery, use of concomitant medications, cardiovascular disease, OA, type 2 diabetes, malignancies, chronic kidney disease, joint replace (knee/hip), number of comorbid conditions and region.

Propensity scores were used to match each HRT user to three comparable non-users using a 0.02 SD calliper.19 This is a standard method for minimising confounding by indication which provides participants with balanced baseline characteristics in both treatment arms and eliminates drug users with no comparable controls as well as unexposed participants with measurable contraindications.20

HRT users and matched non-users were included in Cox regression survival models to estimate the effect of HRT use on implant survival. Proportional hazards assumptions were checked using the formal Schöenfeld's residuals test.

In these models, HRT users who started therapy postoperatively were imputed the date when they were defined as drug users (6 or 12 months after first prescription) as index date, and corresponding matched non-users were assigned the same index date. HRT users who started HRT treatment before surgery and their matched non-users had joint replacement date defined as their index date. This method has been proposed as a solution to overcome immortal time bias in pharmaco-epidemiological studies.21

In order to investigate whether the timing of therapy initiation has an impact on implant survival, we differentiated HRT users with a first prescription before primary arthroplasty from those who started HRT after surgery. We then compared these with HRT non-users using similar Cox regression models.

Sensitivity analyses

A series of sensitivity analyses were carried out in order to test whether the existing association between HRT use and time to revision could be causal. First, HRT use was redefined: only patients with at least 12 prescriptions were here considered HRT users with controls rematched to the new index date (date of first prescription of HRT plus 12 months). Second, we categorised HRT users into low adherence (MPR <0.4), intermediate adherence (MPR 0.4 to <0.8) and high adherence (MPR ≥0.8), and estimated HRs for each of these categories when compared with HRT non-users after adjusting for propensity scores. Finally, in order to study whether a higher cumulative use of HRT provided better protection, we classified HRT users into: <365 daily defined doses (DDD), 365–1065 DDDs and >1065 DDDs. Propensity-adjusted HRs were calculated for each of these using the methods described above.

One final sensitivity analysis was run to measure the potential impact that possible unobserved confounders could have on our findings: the Rosenbaum boundaries sensitivity analysis22 ,23 measures the required imbalance of a confounder to explain the observed associations (if any) in propensity-matched cohort analyses.

All statistical analyses were carried out in Stata IC (V.12).


Out of the eligible 24 733 women identified, 3644 were defined as HRT users, and the remaining 21 089 were non-users. Out of these, propensity-matching methods were used to select 2700 HRT users and three comparable controls (8100 women). These were followed up for a median (inter-quartile range) of 3.3 (1.5–6.1) years after primary arthroplasty (figure 1), making a total of 43 481 person-years (32 845 HRT non-users and 10 636 HRT users) of observation. Baseline characteristics for both arms are shown in table 1. Matched HRT users and non-users only differed significantly in the use of medications with a potential deleterious effect on fracture risk (see Covariates section above), which was higher in HRT users. The final propensity score logistic equation yielded a c statistic of 0.81 (95% CI 0.80 to 0.81), indicating a good ability to predict HRT use.

Table 1

Baseline characteristics

Figure 1

Population flow-chart.

Out of the 10 800 study participants, 10 632(98.4%) were censored before revision (8609 subjects at end of study, 803 at the end of year 10 postsurgery and 1220 lost to follow-up) and 168 (1.6%) were revised. Overall cumulative revision rate at 3 years were 0.97% (0.74% to 12.5%) for THA and 0.76% (0.53% to 10.5%) for TKA. HRT use for at least 6 months was associated with a significant reduction in risk of failure: failure incidence was 2.61/1000 person-years-at-risk (pyar) (1.79 to 3.61) in HRT users, and 4.25 (3.81 to 5.02) pyar among HRT non-users, with corresponding HR of 0.62 (0.41 to 0.94; p=0.023) (figure 2). Adjustment for use of drugs with a potential deleterious effect on bone metabolism did not change these results: HR 0.61 (0.40 to 0.92; p=0.019). Similarly, use of HRT for a year or more was related to a further reduction in failure risk: HRT 0.48 (0.29 to 0.78; p=0.003). When HRT users were stratified into those starting therapy before or after primary arthroplasty, no significant protection was observed for those starting preoperatively (HR 1.06 (0.66 to 1.70; p=0.80)), but a strong protective effect was present for those with a first prescription after surgery (HR 0.24 (0.10 to 0.55); p=0.001) (table 2).

Table 2

Cox models on the effect of hormone replacement therapy (HRT) use on implant survival

Figure 2

Kaplan–Meier estimates of probability of revision surgery according to hormone replacement therapy (HRT) use.

In addition, there was a significant association between HRT adherence and implant survival (p for trend 0.007): HR 0.70 (0.44 to 1.14; p=0.15) for HRT users with an MPR of <0.4; HR 0.66 (0.29 to 1.50; p=0.33) for those with an MPR of 0.4 to <0.8; and HR 0.22 (0.05 to 0.89; p=0.033) for users with an MPR ≥0.8 (table 3, figure 3A).

Table 3

Sensitivity analyses: effect of adherence and cumulative use of HRT on implant survival

Figure 3

Kaplan–Meier estimates of probability of revision surgery according to adherence to hormone replacement therapy (HRT) (A) and HRT cumulative use (B).

Similarly, cumulative use of HRT was also directly related to prosthesis survivorship (p for trend 0.008): HR 1.23 (0.63 to 2.42; p=0.54) for users of <365 DDDs; HR 0.40 (0.19 to 0.86; p=0.019) for users of 365 to 1065 DDDs; and HR 0.53 (0.29 to 0.99; p=0.046) for HRT users with >1065 DDDs (table 3, figure 3B).

Finally, a Rosenbaum boundaries sensitivity analysis suggested a γ of 1.30 for the observed association between 6-month HRT use and implant survival. This means that for any unmeasured confounder to explain the observed association, that confounder would need to produce a 30% increase in the odds of receiving HRT. This same analysis, when applied to the 12-month HRT analysis, suggested a γ of 1.65, suggesting that a confounder related to a 65% imbalance in the probability of receiving HRT would explain the observed association.


Key findings

We have shown that HRT use for at least 6 months after primary arthroplasty is associated with a significant reduction in prosthesis failure of up to almost 40% after lower limb total joint arthroplasty. Hormone therapy for at least 1 year appeared related to further reduction in failure rates to about 50% of that in HRT non-users. In addition, higher adherence to HRT, as well as higher duration of the treatment, appeared directly related to improved effects on implant survival.

Preoperative HRT use did not appear related to implant survival in these data.


Aseptic loosening, leading to symptomatic failure, is the most common indication of late revision surgery after a THA/TKA.24 The mechanisms underlying loosening are still obscure, but it is widely accepted that periprosthetic bone loss secondary to chronic inflammation and osteoclastic activity is the main pathway. Osteolysis within this process can be at least partially prevented with antiresorptive agents,25 and bisphosphonates seem to be currently the group of drugs with a wider body of evidence to support them: most of the clinical trials carried out up to date10 ,11 pointed towards a beneficial effect on surrogates for loosening, such as periprosthetic bone loss or implant migration. In addition, a recent cohort study has suggested that bisphosphonate therapy could also extend implant survival and almost halve the risk of failure.15

Regarding HRT, animal studies have suggested that oestrogen deficiency exerts a negative influence on bone tissue around knee implants, while HRT minimises periprosthetic bone loss26 and improves bone ingrowth around the implant,27 comparable with alendronate therapy28 and more effectively than calcitonin.29 Some studies have shown that both in animals30 and in humans,31 treatment with HRT improves osseointegration after dental implantation surgery. These are the first data, to our knowledge, that have been published describing the effect of HRT on prothesis survival in women undergoing a TKA/THA.

The observed lack of association between preoperative HRT use and implant survival can be due to several reasons. We speculate that reducing bone resorption before surgery might have deleterious effects on the integration of the prosthesis immediately after surgery. In contrast, postoperative initiation of HRT appeared strongly associated with an improved implant survival in our data. These findings are novel and require validation in future studies.

As the rate of revision at 7 years is under 5% (NJR), it is likely that any antiresorptive drugs would be used in a targeted manner to patients identified as being at high risk of failure.32 ,33 Given the risk of venous thromboembolic events with HRT,34 HRT should be stopped 6 weeks before and after surgery. This is supported by our observation that the protective effect of HRT was only evident when started after and not before the primary arthroplasty. These findings are encouraging for the research of antiresorptive agents that might improve outcomes following total joint replacement surgery in the future.

Strengths and limitations

Our study has several strengths and limitations. The main limitation of this study is its observational nature. We have used propensity score matching to control for confounding by indication, which has been proposed as the gold standard analytical strategy to reduce such confounding.17 ,35 Unsolved issues are residual confounding (due to unobserved variables, such as bone mineral density, ethnicity or design of implant or type of fixation) and the lack of a placebo drug to which HRT can be compared. The GPRD has been widely used for epidemiological research, and has previously been shown to be valid for prescription data and discrete outcomes such as fracture and surgery. However, we cannot exclude that in a small number of cases, the side of primary or revision surgery may be contralateral to the index primary. This would lead to a misclassification of the outcome, a random bias and therefore would drive any association of HRT with revision towards the null hypothesis.

Although propensity scores are one of the gold standard methods to account for confounding by indication, there is always the risk of unobserved (residual) confounding. Based on our Rosenbaum sensitivity analysis, one such unmeasured confounder would need to be associated with a 30%–65% increase in revision rates to explain the observed association between HRT use (6 and 12 months, respectively) and implant survival. The most likely unobserved predictor of revision in our data is implant type and fixation: there are no data available in GPRD either about the type of prosthesis used or on the use of cement. According to data from the National Joint Registry,36 3-year revision rates can vary from 0.98% to 2.2% for THA and from 1.4% to 2.1% for TKA based on these two elements. Therefore, only an unexpected strong association between HRT use and implant used, for which there is no previous evidence or biological plausibility, would explain our findings.

Finally, the revision rates observed in our data (about 2% for a median 3.3 years of follow-up) are similar to those reported by the National Joint Registry (1.7% for THAs and 1.8% for TKAs for 3 years of follow-up),36 suggesting completeness of the data. However, it must be admitted that the study period (1986–2006) could explain some differences in revision rates in our data compared with the most current National Joint Registry report (clinical activity in 2012). For all these reasons, our findings require replication in external cohorts and, if possible, in randomised controlled trial data, before these results can be implemented in clinical practice.

The main strengths of this study are the number of participants studied and the duration of follow-up. Furthermore, it is a population-based study, which increases the generalisability of the results. Data from GPRD include all population groups, including older and high complexity patients, and are fully representative of clinical practice in the UK, as clinical information is collected through actual primary care practice.


In our observational cohort study, HRT use after arthroplasty is associated with an almost 40% reduction in implant failure after a TKA/THA for OA. Higher adherence and longer treatment duration further improve implant survival. These findings require replication in experimental studies.


We would like to acknowledge all the patients and professionals who participate in data collection within the GPRD/CPRD Database, as well as the COAST Study participants.


Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Files in this Data Supplement:

  • Press release

    Press release

    Files in this Data Supplement:


  • Handling editor Tore K Kvien

  • Collaborators The Coast Study Group—The COASt Study Group consists of the following members: Mark Mullee, James Rafferty, Andrew Carr, Andrew Price, Kassim Javaid, David Beard, Douglas Altman and Nicholas Clarke.

  • Contributors All authors were involved in drafting the article or revising it critically for important intellectual content, and all authors approved the final version to be published. DP-A had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study conception and design: DP-A, NKR, AJ, MKJ and CC. Acquisition of data: DP-A, CC, JM and NKA. Analysis and interpretation of data: DP-A, JM, MKJ, AJ, CC and NKA.

  • Funding This article presents independent research commissioned by the National Institute for Health Research (NIHR) under its Programme Grants for Applied Research funding scheme (RP-PG-0407-10064). The views expressed in this article are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. Support was also received from the NIHR Biomedical Research Unit into Musculoskeletal Disease, Nuffield Orthopaedic Centre and University of Oxford. Partial funding by: Institut Catala de la Salut-IDIAP Jordi Gol (grant for a research fellowship at the IMIM). Unrestricted educational grants from Merck, Sharpe and Dohme, Novartis and Southampton Rheumatology Trust.

  • Competing interests All authors have completed the Unified Competing Interests form at (available on request from the corresponding author) and declare: this article presents independent research commissioned by the National Institute for Health Research (NIHR) under its Programme Grants for Applied Research funding scheme (RP-PG-0407-10064) and unrestricted educational grants from Merck, Sharpe and Dohme, Novartis and Southampton Rheumatology Trust; Daniel Prieto Alhambra received partial funding by the Institut Catala de la Salut-IDIAP Jordi Gol; all authors declare no financial relationships with any organisations that might have an interest in the submitted work in the previous 3 years; no other relationships or activities that could appear to have influenced the submitted work.

  • Ethics approval All research undertaken using data obtained from CPRD is approved by, as appropriate, an ethics committee, a scientific committee and the National Information Governance Board Ethics and Confidentiality Committee.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No data are available for sharing.