Quantification of risk factors for herpes zoster: population based case-control studyBMJ 2014; 348 doi: https://doi.org/10.1136/bmj.g2911 (Published 13 May 2014) Cite this as: BMJ 2014;348:g2911
- Harriet J Forbes, research fellow in epidemiology,
- Krishnan Bhaskaran, lecturer in statistical epidemiology,
- Sara L Thomas, clinical senior lecturer,
- Liam Smeeth, professor of clinical epidemiology,
- Tim Clayton, senior lecturer in medical statistics,
- Sinéad M Langan, NIHR clinician scientist and clinical senior lecturer
- 1Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK
- Correspondence to: H Forbes
- Accepted 13 April 2014
Objectives To quantify the effects of possible risk factors for herpes zoster at different ages.
Design Case-control study.
Setting UK Clinical Practice Research Datalink primary care data.
Participants 144 959 adults diagnosed with zoster between 2000 and 2011; 549 336 age, sex, and practice matched controls.
Main outcome measures Conditional logistic regression was used to generate adjusted odds ratios to estimate the strength of association of each potential risk factor with zoster and assess effect modification by age.
Results The median age of the cases and controls was 62 years. Factors associated with increased risk of zoster included rheumatoid arthritis (3111 (2.1%) v 8029 (1.5%); adjusted odds ratio 1.46, 99% confidence interval 1.38 to 1.55), inflammatory bowel disease (1851 (1.3%) v 5118 (0.9%); 1.36, 1.26 to 1.46), chronic obstructive pulmonary disease (6815 (4.7%) v 20 201 (3.7%); 1.32, 1.27 to 1.37), asthma (10 243 (7.1%) v 31 865 (5.8%); 1.21, 1.17 to 1.25), chronic kidney disease (8724 (6.0%) v 29 437 (5.4%); 1.14, 1.09 to 1.18), and depression (6830 (4.7%) v 22 052 (4.0%); 1.15, 1.10 to 1.20). Type 1, but not type 2, diabetes showed some association with zoster (adjusted odds ratio 1.27, 1.07 to 1.50). The relative effects of many assessed risk factors were larger in younger patients. Patients with severely immunosuppressive conditions were at greatest risk of zoster—for example, patients with lymphoma (adjusted odds ratio 3.90, 3.21 to 4.74) and myeloma (2.16, 1.84 to 2.53), who are not eligible for zoster vaccination.
Conclusions A range of conditions were associated with increased risk of zoster. In general, the increased risk was proportionally greater in younger age groups. Current vaccines are contraindicated in people at the greatest risk of zoster, highlighting the need for alternative risk reduction strategies in these groups.
Herpes zoster (commonly known as shingles) is caused by the reactivation of latent varicella zoster virus when specific cell mediated immunity becomes compromised. It is a common disease in older people,1 with a lifetime risk of up to 30% rising to 50% among those living to 85 years.2 Zoster typically presents as a painful unilateral vesicular dermatomal rash that causes acute morbidity lasting two to four weeks.3 A severe complication—postherpetic neuralgia—although uncommon in patients aged under 50, develops in 12% of zoster patients aged 50 years or over1 4; it causes intense pain that can last from months to years and is associated with considerable impairment of quality of life.
A live vaccine that reduces the risk of zoster and postherpetic neuralgia has recently been developed.4 The vaccine is licensed in peoople aged over 50 years. In the United States and Australia, it is recommended in people aged 60 years or over. In the United Kingdom, it is routinely available for patients aged 70, with a catch-up campaign for 72-79 year olds (on the basis of cost effectiveness studies5); the current catch-up cohort is those aged 79. Table 1⇓ shows some examples of recommedations on vaccination in major countries. As the zoster vaccine is expensive, targeting vaccination towards groups at high risk of zoster is necessary. Age is the most important risk factor for zoster and postherpetic neuralgia, so it drives vaccination policies. The incidence of zoster rises from 3.5 per 1000 person years among 50-54 year olds to 7.1 per 1000 person years among 75-79 year olds in the UK.6 However, whether people with other risk factors, particuarly in younger age groups, might also benefit from vaccination is not clear.
Several clinical conditions not listed as contraindications for the zoster vaccine have been associated in some studies with an increased risk of zoster; these include autoimmune conditions such as rheumatoid arthritis,7 systemic lupus erythematosus,8 and inflammatory bowel disease,9 and chronic conditions such as diabetes mellitus,10 11 chronic obstructive pulmonary disease,11 12 chronic kidney disease,13 14 15 asthma,16 and depression.2 17 Large, adequately powered studies investigating the association of these conditions with zoster are lacking.
This study therefore aimed to quantify the effects of various proposed risk factors for zoster and explore whether their effects differ by age group, to help to identify any groups of patients at high risk who are not currently targeted for vaccination. We aimed to add to the existing literature by using a very large data source and a consistent methodological approach to help comparison between various risk factors.
We did a matched case-control study to quantify the effects of a range of individual level factors on the risk of zoster in a general population.
This study used data from the UK Clinical Practice Research Datalink, a primary care database of anonymised patients’ records containing complete prescribing and diagnostic information and feedback from hospital referrals. Clinical diagnoses are coded in the Clinical Practice Research Datalink with Read codes.18 It is one of the largest sources of continuous patients’ records in the world, containing data on approximately 7% of the UK population, and is broadly representative of patients’ and practices’ characteristics in the UK.19 Sixty per cent of patients in the Clinical Practice Research Datalink have data available in the Hospital Episode Statistics database, which is a linked computerised database of hospital attendances in England from 1997.
Selection of cases
The base study population consisted of all patients aged 18 years or over, under follow-up between 1 January 2000 and 31 December 2011, with no evidence of previous zoster (no codes for zoster or postherpetic neuralgia before the start of follow-up). We identified zoster cases by using the Clinical Practice Research Datalink and Hospital Episodes Statistics. Cases in the Clinical Practice Research Datalink were those with a zoster Read code and at least 12 months’ follow-up before a first diagnosis of zoster, to exclude past cases of zoster recorded retrospectively after registration at a general practice.20 Incident zoster in Hospital Episodes Statistics was identified by ICD-10 (international classification of diseases, 10th revision) codes (B02, B02.0, B02.1, B02.31, B02.7, B02.8, B02.9, G53.0) that appeared in the primary diagnosis field; we excluded patients with these codes recorded in secondary diagnosis fields, as this could reflect either incident zoster or sequelae of past zoster. For cases identified in Hospital Episodes Statistics, we took the date of hospital admission as the index date. If zoster was identified in both Hospital Episodes Statistics and the Clinical Practice Research Datalink, we took the earliest recorded zoster diagnosis as the index date.
Selection of controls
We identified all potential controls for each case and used incidence density sampling to select up to four controls per case at random.21 Controls were registered with the practice at the index date of the case and for at least 12 months before and were matched to cases by practice, sex, and age (within 1 year). We selected controls before their exposure status for other risk factors was known. By matching on practice, we controlled for practice level socioeconomic status. Controls took the index date of their matched case. Controls had no history of zoster or postherpetic neuralgia at the index date, but they could go on to have zoster and therefore also be included as a case.22 Controls with no contact with their practice (that is, no consultation record of any kind, including repeat prescriptions and face to face or telephone consultations) any time from six months before to 12 months after the index date were assumed inactive in the database and excluded.
Information on risk factors was based on records entered in the Clinical Practice Research Datalink any time before the index date. For each risk factor, we developed a definition and a list of relevant Read codes. A clinician reviewed all code lists, and these are available on request. Patients without Read codes for a clinical condition were assumed not to have the condition.
Key risk factors of interest
Our key risk factors of interest were rheumatoid arthritis, systemic lupus erythematosus, inflammatory bowel disease, chronic obstructive pulmonary disease, asthma, chronic kidney disease, depression, and diabetes (overall and by type) (see appendix section A for a full descriptions of how these risk factors were identified).
Our models included severely immunosuppressive conditions determined by the Advisory Committee on Immunization Practices to be contraindications to vaccination—namely, recent history (less than two years before index date) of leukaemia, lymphoma, or bone marrow transplant, or any history of HIV, haematopoietic stem cell transplantation, myeloma, or “other unspecified cellular immune deficiencies” (for example, pancytopenia).23 We also included use of immunosuppressive treatment; we extracted all relevant prescriptions before the index date and calculated duration of prescription (using data on quantity of tablets prescribed and numeric daily dose). We defined oral corticosteroid exposure as a 14 day course of high dose (≥20 mg/day) oral corticosteroids in the month before the index date. We included exposure to other immunosuppressive treatments in the month before the index date as an additional covariate (see appendix section A for further details).
Other characteristics of patients included smoking (ex-smoker, current smoker, non-smoker), alcohol consumption (non-drinker, current drinker, ex-drinker), and body mass index (underweight (<18.5), normal weight (18.5-24.9), overweight (≥25), obese (≥30)) (see appendix section A for further details). We captured exposure to inhaled corticosteroids within three months before the index date (see appendix section A for further details), as some evidence suggests that they may increase the risk of zoster.24 We compared follow-up time (in years) for cases and controls, to check that the opportunity to record the exposures of interest in the two groups was equal.25
To explore whether the association between our risk factors of interest and zoster could be affected by ascertainment bias (that is, visiting a general practice for the risk factor may increase the likelihood of receiving a diagnosis of zoster), we investigated the association between zoster and treated epilepsy (an epilepsy Read code before index date and an epilepsy treatment (British National Formulary chapter 4.8.1) 12 months before the index date). Epilepsy is a chronic condition requiring high level healthcare use26; it is not thought to be associated with zoster.
We described the characteristics of the study population by case-control status. The analysis accounted for the matched design by using conditional logistic regression, so all odds ratios accounted for the matching factors of age (within 1 year), sex, practice, and calendar time. Initially, we calculated univariable odds ratios to estimate the strength of association of each potential risk factor with zoster; 99% confidence intervals were generated to account for multiple testing and reduce the chance of picking up small and clinically unimportant associations. We then did multivariable analyses on patients with complete data for all variables (see also below for approach to missing data).
We investigated a potentially effect modifying role of age at diagnosis (index) date by calculating stratum specific odds ratios for each variable in the multivariable model by age groups (>50 years, 50-59 years, 60-69 years, and ≥70 years). We added interaction terms between age and other risk factors one at a time into the multivariable regression model and used likelihood ratio tests to investigate interactions.
To quantify the effect of our assessed risk factors, we estimated the age specific absolute rate of zoster for each factor by multiplying the age specific effect estimates by the age specific rate of zoster in the general population in 2010. We estimated the age specific general population rates within the Clinical Practice Research Datalink by dividing the number of incident cases of zoster in 2010 by the total person time contribution of the population at risk in that year, which included all patients in the Clinical Practice Research Datalink with at least one year of follow-up and no history of zoster.
We then evaluated to what extent combinations of risk factors could lead to a high risk of zoster, specifically within age groups not targeted for vaccination. We calculated the predicted risk of zoster for each person by summing the coefficients for their risk factors from the final model and multiplying their predicted odds ratio by the age specific rate of zoster in the general population. We were thus able to determine the number of cases aged under 70 but with a predicted risk of zoster equal to or higher than that of the general population aged at least 70 years (vaccination is currently offered to 70 and 79 year olds in the UK). We repeated this analysis with a 60 years age cut-off, to reflect vaccination policy from different countries. We restricted these calculations to people without severely immunosuppressive conditions.
Sensitivity analyses included use of a wider definition of exposure to oral corticosteroids and other immunosuppressive treatment (anyone exposed three months before the index date). We used multiple imputation by chained equations to account for missing data (11% of patients had missing data for alcohol or smoking).27 The imputation model included all covariates from the main outcome model, together with matching variables of age and sex, and extra comorbidities as additional markers of alcohol or smoking related diseases. We created five imputed datasets and combined them for analysis. We further adjusted for patient level socioeconomic status in the subgroup of patients registered at English practices for whom these data were available (see appendix sections B-D for further details). We applied two different definitions of an “active” control. Firstly, we considered controls to be inactive if they had no contact with their general practice during any time from one year before to two years after the index date. Secondly, we considered controls to be inactive if they had not consulted their general practice in the previous three years. Finally, we calculated the mean yearly consultation rate before the index date (by dividing the total number of face to face or telephone consultations during follow-up by the total years of follow-up before the index date) among patients with our risk factors of interest and among patients with epilepsy.
We identified a total of 145 397 incident cases of zoster. After matching, we removed inactive controls (n=27 928; 4.8% of all 577 264 matched controls). We excluded a small number of cases, as they had no eligible controls (n=438; 0.3%), largely owing to advanced age, leaving 144 959 cases.
Table 2⇓ shows descriptive details of the cases and controls; 59.4% of cases and 61.0% of controls were female (the higher proportion among controls was due to the variable number of controls per case). The median age at diagnosis of zoster was 62 (interquartile range 48-73) years for cases and 62 (49-74) years for controls. Approximately 45% of zoster cases occurred in patients under 60 years and 65% in those under 70 years. Registered follow-up time before the index date was equivalent in cases and controls (8.6 (4.3-12.1) years).
Table 3⇓ shows the proportion of cases and controls with our key risk factors of interest and other covariates. The strongest risk factors for zoster were the severely immunosuppressive conditions, although their overall prevalence was low—2.8% of cases and 1.2% of controls (table 3⇓). Odds ratios varied from 1.78 to 13.46, although the latter estimate (for haematopoietic stem cell transplantation) had a particularly wide 99% confidence interval, as it was based on only 29 exposed patients.
Among our key risk factors of interest, systemic lupus erythematosus had the strongest relative risk of zoster (adjusted odds ratio 1.72, 99% confidence interval 1.45 to 2.04) (table 3⇑), but the condition was very rare (0.3% of cases and 0.1% of controls). The risk of zoster was increased by more than 30% among patients with rheumatoid arthritis (adjusted odds ratio 1.46, 1.38 to 1.55) and chronic obstructive pulmonary disease (1.32, 1.27 to 1.37). Asthma, chronic kidney disease, and depression were associated with a greater than 10% increased risk of zoster. Although no association existed between zoster and diabetes overall (adjusted odds ratio 1.02, 0.99 to 1.05), we found strong evidence that patients classified as having type 1 diabetes were at increased risk of zoster (1.27, 1.07 to 1.50).
To assess whether some of the overall effect of these conditions was mediated by their treatments, we additionally adjusted for immunosuppressive treatments and inhaled corticosteroids (table 3⇑, final column). The associations between zoster and rheumatoid arthritis, systemic lupus erythematosus, inflammatory bowel disease, chronic obstructive pulmonary disorder, and asthma were attenuated after adjustment but did not disappear. Among conditions not commonly treated with these drugs (chronic kidney disease, depression, and diabetes), the association did not vary after this additional adjustment.
Treated epilepsy, not expected to be associated with zoster, was marginally more common among cases than controls (1420 (0.98%) v 5142 (0.94%)). We found little evidence of association with zoster in crude (odds ratio 1.05, 0.97 to 1.14) or adjusted analyses (1.06, 0.97 to 1.15), suggesting that other associations observed in this study are unlikely to be explained by ascertainment bias.
We explored how the effect of our risk factors varied by age (tables 4⇓ and 5⇓). We found some evidence of age interaction. In general, the relative effects of these risk factors on zoster decreased with increasing age. For example, among patients with chronic kidney disease the odds ratio point estimate was 1.53, 1.21, and 1.11 among patients aged 18-49, 50-59, and 60-69 years respectively.
To quantify the effect of our key risk factors of interest, particularly among age groups not targeted for vaccination, we estimated absolute rates of zoster (table 6⇓). The rates remained low among patients aged 18-49 years (for example, for those with rheumatoid arthritis: 3.51 (99% confidence interval 2.40 to 5.13) per 1000 person years), despite the high relative risks in this age group. The absolute rates of zoster among patients aged 60-69 years were relatively high, particularly for those with chronic obstructive pulmonary disease and rheumatoid arthritis.
We calculated the predicted risk of zoster among 143 620 zoster cases (having excluded 1339 patients with severely immunosuppressive conditions) on the basis of their risk factor profile. Among these cases, 40 409 (28.1%) had at least one risk factor of interest. Of 97 789 zoster cases aged under 70 years, 3619 (3.7%) had a predicted risk as high as that for the general population aged 70 or over; 643 (0.96%) of 67 000 patients under 60 years had a predicted risk as high as those aged 60-69 years.
Sensitivity analyses imputing missing data for smoking and alcohol (see appendix section B), using different definitions for immunosuppressive treatment (see appendix section C), adjusting for patient level socioeconomic status (see appendix section D), and applying different definitions of an active control (see appendix section E) did not change the study findings. The consultation rate was similar among patients with epilepsy and those with our risk factors of interest (see appendix section F).
In this large matched case-control study, a range of conditions were associated with increased risk of zoster, including rheumatoid arthritis, systemic lupus erythematosus, inflammatory bowel disease, chronic obstructive pulmonary disease, asthma, chronic kidney disease, type 1 diabetes, and depression. The relative effects of many of these risk factors were larger among younger patients. This raises the question of whether vaccination of certain younger, high risk groups may be beneficial; cost effectiveness studies, also considering the risk of postherpetic neuralgia, would be needed to answer this question. Among a wide range of risk factors, people with contraindications to currently available vaccines were at the highest risk of zoster, highlighting the need to identify strategies to reduce the risk of zoster among these groups.
Strengths and limitations of study
This large study included all incident cases of zoster in the Clinical Practice Research Datalink over an 11 year period, providing high power to detect small effects, except for rare potential risk factors. We assessed a wide range of factors, allowing adjustment for many of the known as well as potential confounders, in addition to exploring interactions with age, which is of key public health relevance in terms of vaccination policy.
The study may be subject to ascertainment bias; regular general practice visits for a chronic condition such as diabetes may increase the likelihood of receiving a zoster diagnosis. However, we believe that most zoster cases would present to their general practice, owing to the extensive rash and considerable pain associated with zoster and free general practice attendance in the UK.28 A survey about immunisation practices in the United States among people aged 60 and over found that 95% of those who knew they had zoster sought care.29 Furthermore, we found little evidence to suggest that patients with epilepsy (who were visiting their general practice as frequently as patients with our risk factors of interest) were at increased risk of zoster.
As with any observational study, residual confounding may be present. For example, exposure to varicella contacts, hypothesised to naturally boost varicella zoster virus specific immunity in people with latent varicella zoster virus infection and thus protect against zoster,30 was not available.
Some misclassification of zoster is possible. In UK primary care, diagnosis of zoster is clinically based with no laboratory testing. Validation of zoster diagnoses in primary care in the Netherlands found that 90.8% of cases had raised antibodies indicating recent zoster infection, suggesting that diagnosis of zoster has a high positive predictive value.31 The frequency with which general practitioners wrongly reject the diagnosis of zoster has not, however, been investigated; this may be more common among younger patients in whom zoster is uncommon. However, the most common misclassification of zoster is as recurrent herpes simplex, which very rarely presents in a dermatomal distribution.32 Zoster is typically a straightforward clinical diagnosis, so the extent of misclassification is likely to be minimal. Any misclassification is likely to be non-differential with respect to exposure status, thus leading to an underestimation of associations.
Finally, excluding “inactive” controls to remove patients who would not attend their general practice if they did develop zoster, may have resulted in healthier controls not being included and led to an underestimation of the odds ratios.
Comparison with other studies
The strength of association between the autoimmune conditions in this study and risk of zoster is broadly in line with previous studies.7 8 9 33 The increased risk associated with these conditions may partly reflect exposure to immunosuppressive drugs, as well as the conditions themselves. When we adjusted for immunosuppressive treatments, the increased risk of zoster was attenuated, suggesting that these drugs contribute to the overall increased risk of zoster in these patients.
The effect of chronic obstructive pulmonary disease is also consistent with previous studies.11 12 After adjustment for oral and inhaled corticosteroids commonly used to treat chronic obstructive pulmonary disease, the association reduced in these patients, again suggesting that these treatments add to their overall risk of zoster. Larger associations between zoster and depression have been reported.11 17 34 No consensus exists on how to capture depression in electronic records; previous studies used different definitions, which may contribute to the varied findings. For example, a matched case-control study using commercial claims and encounters data in the United States found a higher risk of zoster among patients aged over 65 with depression (odds ratio 1.52, 95% confidence interval 1.46 to 1.58) but only included very severe cases (depressive disorder, chronic depressive personality disorder, or dysthymic disorder).11 The association of zoster with chronic kidney disease was also smaller than previous studies have suggested.15
The risks of zoster are reportedly increased in patients with diabetes of unspecified type.10 11 34 35 We found some evidence that type 1, but not type 2, diabetes was associated with an increased risk of zoster (although differentiating between type was not based on diagnostic codes, so misclassification is possible). Patients with type 1 diabetes may be at higher risk of zoster owing to the immune dysregulation resulting from autoimmune pathogenesis.36 The lack of effect found with type 2 diabetes must be interpreted cautiously, as negative confounding by ethnic group is possible; the type 2 diabetes population may have included South Asian patients (known to be at lower risk of zoster2 37 and higher risk of diabetes). Ethnicity is poorly recorded for older patients not recently registered in the Clinical Practice Research Datalink, so we could not capture this characteristic.
The observed increased risks of zoster suggest that all these conditions are associated with some cell mediated immunosuppression. Disease induced immunosuppression has been proposed in chronic obstructive pulmonary disease and asthma patients,38 39 40 41 42 although long term exposure to corticosteroids in these patients is also likely to contribute to immunosuppression. Patients with chronic kidney disease are known to have immune dysfunction,43 particularly when the disease reaches end stage. For depression, stress may increase the risk of zoster as it affects several neuroendocrine functions that can affect cell mediated immunity,44 and depression itself has been associated with poorer zoster immunity.45 46
Conclusions and policy implications
Severe immunosuppression is known to be associated with an increased risk of zoster.11 47 48 49 What this study has highlighted, however, is that the strongest clinical risk factors for zoster are contraindications to its vaccine; the people arguably in most need of protection against zoster cannot currently benefit from vaccination. Alternative risk reduction strategies in these patients would help those at greatest risk of this disease and its complications.
Another important consideration regarding vaccination policy is whether the risk factors for zoster are also risk factors for zoster’s main morbidity, postherpetic neuralgia. To date, very few studies have looked at this research question, and most have been underpowered. The role of severe immunosuppression as a risk factor for postherpetic neuralgia is inconclusive; three studies found that “general” immunosuppression was a risk factor for postherpetic neuralgia,35 50 51 whereas another study found no association.52 Autoimmune conditions have scarcely been assessed; one large cohort study in 34 280 zoster patients identified in Taiwanese health insurance records found that those with systemic lupus erythematosus were at greater risk of postherpetic neuralgia (rate ratio 2.27, 95% confidence interval 1.75 to 2.94).35 Three studies found point estimates for the association between diabetes and postherpetic neuralgia greater than one in multivariate analyses, although the evidence was insufficient to confirm the association.51 53 54 Some evidence suggests that depression is associated with postherpetic neuralgia.55 Larger studies assessing these and other risk factors for postherpetic neuralgia are needed to improve our understanding of the risk factors for postherpetic neuralgia.
This study also raises the question of whether younger age groups at high risk of zoster may benefit from vaccination. Most cases of zoster occur in age groups too young for vaccination; this is because vaccination policy considers the risk of both zoster and postherpetic neuralgia, which is very uncommon in patients under 50. Deciding whether vaccination of groups of patients currently too young for vaccination would be cost effective will therefore require additional research on the risk factors for postherpetic neuralgia, specifically among patients aged under 50.
Our study has identified a range of conditions that are associated with an increased risk of zoster, raising the question of whether targeted zoster vaccination of specific high risk groups at younger ages is warranted. However, this study has also shown that patients at the greatest risk of zoster are not currently eligible to receive the vaccine.
What is already known on this topic
A new herpes zoster vaccination campaign introduced in the United Kingdom in September 2013 is focused on older people only
Recent literature has suggested that a range of clinical conditions are associated with an increased risk of zoster, raising the possibility that some younger people may be at high risk
However, large, highly powered studies investigating the association of these clinical conditions with zoster are lacking
What this study adds
Conditions associated with increased risk of zoster included rheumatoid arthritis, systemic lupus erythematosus, inflammatory bowel disease, chronic obstructive pulmonary disease, asthma, chronic kidney disease, type 1 diabetes, and depression
The increased risks were generally greater among younger age groups
The strongest risk factors are contraindications to vaccination, emphasising the need for alternative risk reduction strategies among these groups
Cite this as: BMJ 2014;348:g2911
Contributors: All authors were involved in the design of the study. HJF did the statistical analysis and wrote the first draft. All authors contributed to further drafts and approved the final manuscript. HJF is the guarantor.
Funding: This article presents independent research supported by a National Institute for Health Research (NIHR) clinician scientist fellowship (to SML; grant number NIHR/CS/010/014). KB is supported by an NIHR postdoctoral fellowship (PDF-2011-04-007). LS is supported by a Wellcome Trust senior research fellowship in clinical science. SLT is supported by an NIHR career development fellowship (CDF-2010-03-32). All authors carried out this research independently of the funding bodies. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the UK Department of Health.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: Approval was obtained from the Clinical Practice Research Datalink Independent Scientific Advisory Committee (11_028A2RA) and the LSHTM Research Ethics Committee (A241 5930).
Transparency declaration: The lead author (the manuscript’s guarantor) affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.
Data sharing: No additional data available.
This is an Open Access article distributed in accordance with the Creative Commons Attribution (CC BY-3.0) which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited: http://creativecommons.org/licenses/by/3.0/.