Impact of the herpes zoster vaccination programme on hospitalised and general practice consulted herpes zoster in the 5 years after its introduction in England: a population-based study.

OBJECTIVES: To assess the impact of herpes zoster vaccination in the 5 years after introduction for 70- to 79-year-olds in England in September 2013.
DESIGN: Population based ecological impact assessment.
SETTING: Hospitals covering the whole English population for the period 2008 to 2018 and 293 general practices (GP) for the period 2005 to 2018, in England. PARTICIPANTS: Over the period the population contributed 117·5 million person-years for hospitalisation events and 6.96 million person-years for GP events in individuals aged 60 to 89.
INTERVENTIONS: Live attenuated herpes zoster vaccination (Zostavax), first used on 1st September 2013, in 70- and 79-year-olds with continued use in new 70 year-olds and with a staged catch-up of those aged 71 to 78 years in 2013. OUTCOME MEASURES: Herpes zoster and postherpetic neuralgia (PHN) consultation and hospitalisation rates in age-cohorts according to vaccine eligibility. Incidence rate ratios in age-cohorts eligible for vaccination compared with those non-eligible were calculated by Poisson regression. This was used to estimate prevented cases and, along with vaccine coverage, to estimate vaccine effectiveness.
RESULTS: Large and prolonged reductions in herpes zoster and PHN consultations and hospitalisations were observed in the 5 years post-implementation. For example, in 79 year-olds first eligible in 2013, the incidence rate ratio for consultations 5 years later was 0·65 (95% CI: 0·52 to 0·81). Over the whole period an estimated 40 500 fewer zoster consultations and 1840 fewer zoster hospitalisations occurred because of the vaccination programme. These reductions were consistent with effectiveness in the routine cohorts (vaccinated aged 70) of between 37% (for hospitalised zoster) and 75% (for PHN consultations) and, in catch up cohorts (vaccinated aged 78 to 79) of between 49% (for hospitalised PHN) and 66% (for PHN consultations).
CONCLUSION: Given the clear and sustained impact of herpes zoster vaccination over the 5-year period since introduction, optimising vaccination coverage is important to attain maximum benefit. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Use of a large general practice cohort and national hospital episode data has allowed precise estimation of the impact of the herpes zoster vaccination programme.

Five years of follow-up post vaccine implementation has allowed detailed assessment of the duration of protection against herpes zoster and postherpetic neuralgia.

The staged introduction of the vaccine with catch-up in different cohorts over time has allowed a robust analysis of impact accounting for time trends and age effects.

As an impact study it was only possible to indirectly assess vaccine effectiveness using impact estimates and coverage data.

Assessing duration of protection was approximate because up to 10% of individuals were vaccinated after the first year they were targeted.

Introduction

The herpes zoster (shingles) vaccine Zostavax (Zoster Vaccine Live; Merck & Co, Kenilworth, New Jersey, USA) was introduced in the UK on 1st September 2013 following the recommendation for its use in immunocompetent 70- to 79-year-olds by the UK’s Joint Committee on Vaccination and Immunisation (JCVI).1 2 The aim of the programme is to boost individuals’ pre-existing immunity to varicella zoster virus (VZV) and prevent herpes zoster which is caused by reactivation of latent VZV. Herpes zoster is typically characterised by a unilateral dermatomal vesicular rash which can lead to complications including postherpetic neuralgia (PHN) in which pain persists for more than 3 months after rash onset, especially in immunocompromised individuals and older adults.3 Zostavax is a live attenuated herpes zoster vaccine containing the same Oka VZV strain as the childhood varicella vaccine (Varivax) but at a higher dose. Single dose efficacy from clinical trials of 38% against shingles and 67% against PHN over a 3-year follow-up has been demonstrated in adults ≥70 years of age with more recent evidence indicating efficacy maintained at about 30% for herpes zoster and PHN up to 8 years after vaccination.4 5 The choice of target age group was based on a cost-effectiveness analysis which incorporated the age-specific incidence of herpes zoster and PHN,6 the decline in vaccine efficacy with age and the estimated duration of vaccine-induced protection.1

In the first year of the programme (2013 to 2014), a single dose of vaccine was offered to adults aged 70 (routine cohort) and 79 years (catch up) on 1st September 2013. In the following 4 years further routine and catch-up cohorts were offered the vaccine as outlined in table 1, with the most recent cohorts offered the vaccine as they turned 70 or 78 years of age rather than by age at September 1 of the year.

Table 1

Eligibility for shingles vaccination in England

Birth cohort	Age on 1 Sept 2013 (years)	First became eligible
Routine cohorts
2 Sept 1942 – 1 Sept 1943	70	2013–14
2 Sept 1943 – 1 Sept 1944	69	2014–15
2 Sept 1944 – 1 Sept 1945	68	2015–16
2 Sept 1945 – 1 Sept 1946	67	2016–17
2 Sept 1946 – 1 Sept 1947	66	2017–18
2 Sept 1947 – 1 Sept 1948	65	2017–18 (when turning 70)
Catch-up cohorts
2 Sept 1933 – 1 Sept 1934	79	2013–14
2 Sept 1934 – 1 Sept 1936	77 to 78	2014–15
2 Sept 1936 – 1 Sept 1937	76	2015–16
2 Sept 1937 – 1 Sept 1938	75	2016–17
2 Sept 1938 – 1 Sept 1939	74	2017–18
2 Sept 1939 – 1 Sept 1940	73	2017–18 (when turning 78)

Individuals remain eligible until their 80th birthday.

Evaluation of the vaccine programme 3 years following its introduction demonstrated a clear impact in these target ages on general practice (GP) consultations for herpes zoster and PHN with herpes zoster incidence reducing from about nine to six per 1000 person-years.7 In addition, vaccine effectiveness in the 3-year period has been demonstrated as 64% (95% CI: 60% to 68%) against GP consultations for herpes zoster and 81% (95%CI: 61% to 91%) for PHN.8 In this paper, we extend the previous 3-year evaluation to a 5-year period post introduction, which allows a fuller assessment of waning. We also assess impact on hospitalisations for herpes zoster and PHN allowing a more complete evaluation of the reduction in disease burden.

Methods

Data sources

General practice data on consultations and vaccine uptake

The Royal College of General Practitioners (RCGP) Research and Surveillance Centre (RSC) is a sentinel primary care network representing over 1% of the English population that is geographically representative of the general population.9 In March 2019 we obtained data on herpes zoster, PHN and herpes zoster vaccinations for the period from 1st October 2005 to 30th September 2018 for patients aged 60 to 89 years from all 293 practices in the network. We also obtained denominator data for patients registered each month, stratified by age at 1st September 2013, year/month, gender and GP practice. Events were identified using Read codes (coded thesaurus of clinical terms used for recording patient findings and procedures—see online supplementary file 1).10 For analysis we aggregated these data using the same strata as for the denominator data. Herpes zoster consultations and PHN consultations were dropped if they occur within 12 months of a previous consultation.

Hospital admissions data

For hospital admissions we used Hospital Episode Statistics (HES) Admitted Patient Care, which contains individual level hospital admissions for all at National Health Service hospitals in England.11 In August 2019, we obtained all Finished Consultant Episodes for the period from 1st October 2008 to 30th September 2018, for patients aged 60 to 89 years with an International Classification of Diseases (ICD) code B022 herpes zoster or with ICD code G530 PHN in any of the patients 20 diagnosis fields on discharge. We required a minimum gap of 12 months between episodes. We calculated age in September before the year they became ill and aggregated data by age and month/year of admission. To obtain denominators we obtained England Office for National Statistics population statistics by year of age for the period 2008 to 2017 and matched these to the appropriate numerator.12 We then generated the age on 1st September 2013 so that the vaccine eligible cohorts could be allocated.

Statistical analysis

We defined age cohorts according to participants age on 1st September 2013, which allowed us to identify those targeted in each year of the programme as given in table 1. Those becoming eligible by turning 70 or 78 in 2017 had low uptake in 2017 to 2018 and we therefore excluded these cohorts from the impact analysis from 1st October 2017 by treating them as a separate group for this period. We calculated cumulative vaccine uptake from the RCGP data by summing the vaccine uptake within each month from September 2013 to August 2018 for all cohorts.

To assess impact of the vaccine programme, we modelled the incidence of GP consultations as well as of hospital admissions according to age (in years), time period (2005 to 2018 for RCGP and 2008 to 2018 for HES) and vaccine eligibility. Vaccine eligibility is the factor that measures impact and we determined this for each incidence strata according to age and period and assigned this as not eligible, first eligible in that year (defined as 1st October to 30th September), and first eligible in 1, 2, 3 and 4 years previously. The reason for starting in October before measuring impact was to allow time for immunity to develop and to take into account low uptake in the September of each year. Individual vaccination status was therefore not used in this impact analysis, just eligibility by birth cohort. For the RCGP data we also assessed gender and region (London, Midlands, north England and south England) but did not include these in the final models as they did not change impact estimates.

We modelled the data using multivariable Poisson regression with an offset for person time. This model measures impact as the relative incidence in the vaccine eligible cohorts compared with unvaccinated cohorts. We investigated trends by age and period as linear, quadratic and cubic effects as well as an interaction effect between age and period, with the simplest model retained if the simplification did not result in more than a three percentage point change in impact estimates, although we retained a linear period effect in all models. To check robustness of impact estimates we fitted a random effects model with GP as the random effect as well as a model with age and period included as factors rather than polynomials.

To predict the expected incidence and number of events in the absence of vaccination we fitted the final model using data from the cohorts who were not eligible for vaccination and used this to extrapolate to the vaccine eligible cohorts. We also used this model to generate figures of observed compared with modelled incidence as well as a table of observed and expected numbers of events along with expected incidence and incidence reduction. We did sample size calculations prior to obtaining data to ensure reductions of 25% could be detected (80% power and 5% significance). We estimated vaccine effectiveness (VE) consistent with impact and coverage by dividing the impact (percentage reduction) by the average cumulative coverage in the period/age cohort assessed. For example, a VE of 60% is consistent with a 30% reduction and 50% coverage. Precision of VE is based on the 95% CIs for impact. We did all analyses in Stata 14.

Patient and public involvement

Patients or the public were not involved in the design or planning of the study.

Results

Vaccine uptake

Uptake based on the 293 practices contributing to the RCGP RSC shows the routine cohorts in the first 4 years of the programme achieving about 60% to 65% uptake in the first year they are targeted and then having further gradual increases to over 75% (figure 1). The figure does show that more recent routine cohorts have achieved lower coverage, for example, by the August at the end of the first year targeted the first routine cohort (age 70 in September 2013 and vaccinated in 2013/2014) had 63% uptake compared with 46% for the fourth routine cohort (age 66 in September 2013 and vaccinated in 2017/2018). The catch-up cohorts achieve 60% to 65% uptake which plateaus as they are no longer eligible once they reach age 80 (figure 2). Uptake in these practices was similar to that reported nationally which was 61·8%, 59·0%, 54·9%, 48·3% and 49·1%, respectively, for the first 5 years of the programme in the routine cohorts (aged 70, 69, 68, 67 and 66 years in September 2013) and 59·6%, 57·8%, 58·5%, 55·5%, 49·4% and 50·8%, respectively, for the catch-up cohorts (79, 78, 77, 76, 75 and 74 years in September 2013).13–17

Figure 1

Cumulative coverage for each routine age cohort throughout 5-year vaccination programme. Cohorts are grouped by age on 1st September 2013.

Figure 2

Cumulative coverage for each catch-up age cohort throughout 5-year vaccination programme. Cohorts are grouped by age on 1st September 2013.

Vaccine impact on GP consultations

From the 293 RCGP practices the number of patients registered aged 60 to 89 years at October 1st each year increased from 447 567 in 2005 to 621 993 in 2018 with an average of 542 039 over the 14 years. These patients contributed a total 6·96 million person-years of data and are described in table 2 to show herpes zoster and PHN events and rates by vaccine eligibility, year and age. The incidence data were then modelled to estimate vaccine impact. Following the fitting process the final model included the variable to measure impact (vaccine eligibility), a log-linear time trend (a 0·2% increase per year, independent of age) and a quadratic age effect, with incidence increasing from 5·9 per 1000 person-years at 60 to a plateau of 10·2 by age 85. No significant interactions were identified. The model with factors for age and period gave similar results, as did the random effects model.

Table 2

Description of the aggregated RCGP and hospitalised cohort person-years and herpes zoster events according to vaccine eligibility, year and age

		GP consultations (RCGP RSC)						Hospitalisations (England)
Factor	Level	Eligible for vaccination			Not eligible for vaccination			Eligible for vaccination			Not eligible for vaccination
		Total person-years* (1000s)	Herpes zoster events (rate per 1000 person- years)	PHN events (rate per 1000 person-years)	Total person-years (1000s)	Herpes zoster events (rate per 1000 person- years)	PHN events (rate per 1000 person-years)	Total person-years* (1000s)	Herpes zoster events (rate per 10 000 person- years)	PHN events (rate per 10 000 person-years)	Total person-years (1000s)	Herpes zoster events (rate per 10 000 person- years)	PHN events (rate per 10 000 person-years)
Year	2005–06				446	3514 (7.9)	510 (1.1)
	2006–07				463	3738 (8.1)	548 (1.2)
	2007–08				485	3907 (8.1)	562 (1.2)
	2008–09				501	4010 (8.0)	573 (1.1)				11 104	3658 (3.3)	592 (0.5)
	2009–10				516	4187 (8.1)	598 (1.2)				11 285	3987 (3.5)	739 (0.7)
	2010–11				529	4266 (8.1)	595 (1.1)				11 440	4042 (3.5)	773 (0.7)
	2011–12				542	4386 (8.1)	625 (1.2)				11 593	4142 (3.6)	821 (0.7)
	2012–13				556	4586 (8.3)	652 (1.2)				11 753	4280 (3.6)	878 (0.7)
	2013–14	38	252 (6.6)	39 (1.0)	532	4357 (8.2)	621 (1.2)	816	227 (2.8)	66 (0.8)	11 118	4223 (3.8)	950 (0.9)
	2014–15	94	555 (5.9)	77 (0.8)	489	3946 (8.1)	584 (1.2)	1969	666 (3.4)	182 (0.9)	10 132	3919 (3.9)	859 (0.8)
	2015–16	133	828 (6.2)	117 (0.9)	462	3536 (7.7)	511 (1.1)	2777	949 (3.4)	224 (0.8)	9502	3688 (3.9)	802 (0.8)
	2016–17	173	1074 (6.2)	140 (0.8)	432	3410 (7.9)	478 (1.1)	3666	1417 (3.9)	339 (0.9)	8804	3514 (4.0)	791 (0.9)
	2017–18	219	1320 (6.0)	211 (1.0)	351	2688 (7.7)	352 (1.0)	4412	1904 (4.3)	481 (1.1)	7123	2830 (4.0)	584 (0.8)
Age	60–69†				3415	23 381 (6.8)	2402 (0.7)				57 270	9897 (1.7)	1849 (0.3)
	70–73	346	1940 (5.6)	234 (0.7)	715	6154 (8.6)	862 (1.2)	7164	1874 (2.6)	434 (0.6)	10 951	3303 (3.0)	755 (0.7)
	74–77	26	139 (5.4)	14 (0.5)	843	7696 (9.1)	1299 (1.5)	519	153 (3.0)	38 (0.7)	14 016	6287 (4.5)	1453 (1.0)
	78–82	271	1850 (6.8)	321 (1.2)	608	5946 (9.8)	1135 (1.9)	5669	2925 (5.2)	774 (1.4)	9240	5888 (6.4)	1352 (1.5)
	83–89†	14	100 (7.0)	15 (1.1)	722	7354 (10.2)	1511 (2.1)	288	211 (7.3)	46 (1.6)	12 378	12 908 (10.4)	2380 (1.9)

*Vaccine eligible person-years includes all time from the point that each age-cohort vaccination first became eligible to receive vaccine.

†Age group not eligible for vaccination.

GP, general practice; PHN, postherpetic neuralgia; RCGP, Royal College of General Practitioners; RSC, Research and Surveillance Centre.

Figures 3 and 4 show the observed and predicted incidence from the model in the absence of vaccination for herpes zoster and PHN by year (2005/2006 to 2017/2018) for those cohorts aged between 66 years and 81 years as of 1st September 2013. The figures clearly show the impact in the first year each cohort is targeted for vaccination and its persistence in the following years.

Observed GP-consulted herpes zoster incidence and model-predicted incidence in the absence of vaccination by year for each age cohort. Cohorts are grouped by age at 1 September 2013. Each year runs from October 1 of the first year to September 30 of the second. Blue points indicate observed zoster incidence and red points indicate observed incidence in age cohorts eligible for vaccination. Red lines show the model-predicted incidence in the absence of vaccination. Note that ages 71, 72, 80 and 81 are cohorts not eligible in the period. GP, general practice.

Observed GP-consulted postherpetic neuralgia incidence and model-predicted incidence in the absence of vaccination by year for each age cohort. Cohorts are grouped by age at 1 September 2013. Each year runs from October 1 of the first year to September 30 of the second. Blue points indicate observed zoster incidence and red points indicate observed incidence in age cohorts eligible for vaccination. Red lines show the model-predicted incidence in the absence of vaccination. Note that ages 71, 72, 80 and 81 are cohorts not eligible in the period. GP, general practice.

Cumulative uptake observed and predicted cases and relative incidence estimates from the fitted models are summarised for herpes zoster (table 3) and PHN (table 4). The tables also show the vaccine effectiveness required to generate the observed reductions, which is approximately 50% to 60% for herpes zoster. The incidence reduction, if applied to the England population translates to approximately 40 500 fewer herpes zoster (including PHN) episodes over the study period. Reductions were greater for PHN (table 4) with a 47% reduction in PHN incidence in the routine cohorts across the first 5 years and a 38% reduction in the catch-up cohorts. This would be equivalent to a reduction in PHN episodes of 0·6 to 0·7 per 1000 person-years consistent with vaccine effectiveness of approximately 75% for routine and 66% for catch-up cohorts. The incidence reduction, if applied to the England population, translates to approximately 8700 fewer PHN episodes over the 5-year study period.

Table 3

Impact of routine and catch-up herpes zoster vaccination on GP-diagnosed herpes zoster by time since cohorts were first eligible for vaccination

Period	Age on 1 Sept 2013 (age when first eligible for vaccination)	Average cumulative uptake*	Expected events†	Observed events	Incidence rate ratio‡(95% CI)	Expected incidence per 1000 person-years	Incidence reduction per 1000 person-years(95% CI)	Vaccine effectiveness§(95% CI)
Routine cohorts
First year after vaccine eligibility	66–70 years	47%	1098	851	0.77 (0.72 to 0.83)	8.24	1.87 (1.42 to 2.30)	48% (37% to 59%)
	(70–71 years)
Second year after vaccine eligibility	67–70 years	64%	838	509	0.61 (0.55 to 0.66)	8.45	3.33 (2.85 to 3.76)	61% (53% to 69%)
	(70–71 years)
Third year after vaccine eligibility	68–70 years	71%	616	378	0.61 (0.55 to 0.68)	8.66	3.36 (2.79 to 3.88)	54% (45% to 63%)
	(70–71 years)
Fourth year after vaccine eligibility	69–70 years	75%	413	217	0.52 (0.46 to 0.6)	8.86	4.22 (3.55 to 4.80)	63% (53% to 72%)
	(70–71 years)
Fifth year after vaccine eligibility	70 years	77%	199	124	0.62 (0.52 to 0.74)	9.05	3.42 (2.33 to 4.34)	49% (33% to 62%)
	(70–71 years)
All years of the programme	66–70 years	62%	3164	2079	0.66 (0.63 to 0.69)	8.50	2.92 (2.66 to 3.18)	55% (50% to 60%)
	(70–71 years)
Catch-up cohorts
Second year after vaccine eligibility	75–79 years	62%	757	481	0.63 (0.58 to 0.69)	9.82	3.60 (3.00 to 4.14)	59% (49% to 68%)
	(78–80 years)
Third year after vaccine eligibility	77–79 years	63%	578	397	0.68 (0.62 to 0.76)	9.95	3.14 (2.42 to 3.79)	50% (38% to 60%)
	(78–80 years)
Fourth year after vaccine eligibility	78–79 years	63%	409	273	0.67 (0.59 to 0.75)	10.07	3.36 (2.51 to 4.13)	53% (40% to 65%)
	(78–80 years)
Fifth year after vaccine eligibility	79 years	59%	124	81	0.65 (0.52 to 0.81)	10.18	3.54 (1.91 to 4.84)	58% (32% to 80%)
	(79–80 years)
All years of the programme	74–79 years	57%	2811	1950	0.69 (0.66 to 0.73)	9.85	3.02 (2.68 to 3.34)	53% (47% to 59%)
	(78–80 years)

*Calculated by taking the mean of the cumulative uptake values in each month from October to September of the relevant years and cohorts.

†Expected if vaccine not introduced; based on model results for unvaccinated cohorts.

‡Estimated from the Poisson regression model with a log-linear time trend, quadratic age effect and the factor for vaccine eligibility.

§Effectiveness required to generate the observed reduction given the impact and coverage.

GP, general practice.

Table 4

Impact of routine and catch-up herpes zoster vaccination on GP-diagnosed postherpetic neuralgia by time since cohorts were first eligible for vaccination

Period	Age on 1 Sept 2013 (age when first eligible for vaccination)	Average cumulative uptake*	Expected events†	Observed events	Incidence rate ratio‡ (95% CI)	Expected incidence per 1000 person-years	Incidence reduction per 1000 person-years(95% CI)	Vaccine effectiveness§(95% CI)
Routine cohorts
First year after vaccine eligibility	66–70 years	47%	151	103	0.67 (0.55 to 0.82)	1.13	0.37 (0.20 to 0.51)	70% (38% to 95%)
	(70–71 years)
Second year after vaccine eligibility	67–70 years	64%	121	52	0.43 (0.32 to 0.56)	1.22	0.70 (0.53 to 0.82)	89% (68% to 100%)
	(70–71 years)
Third year after vaccine eligibility	68–70 years	71%	93	46	0.49 (0.37 to 0.66)	1.30	0.66 (0.45 to 0.83)	71% (48% to 89%)
	(70–71 years)
Fourth year after vaccine eligibility	69–70 years	75%	65	35	0.53 (0.38 to 0.75)	1.39	0.65 (0.35 to 0.86)	62% (34% to 82%)
	(70–71 years)
Fifth year after vaccine eligibility	70 years	77%	32	12	0.37 (0.21 to 0.64)	1.47	0.93 (0.52 to 1.17)	82% (46% to 100%)
	(70–71 years)
All years of the programme	66–70 years	62%	461	248	0.53 (0.46 to 0.61)	1.24	0.58 (0.49 to 0.66)	75% (63% to 86%)
	(70–71 years)
Catch-up cohorts
First year after vaccine eligibility	74–79 years	48%	174	133	0.76 (0.63 to 0.90)	1.78	0.44 (0.18 to 0.65)	51% (21% to 77%)
	(78–80 years)
Second year after vaccine eligibility	75–79 years	62%	143	87	0.60 (0.48 to 0.74)	1.86	0.75 (0.48 to 0.96)	64% (41% to 83%)
	(78–80 years)
Third year after vaccine eligibility	77–79 years	63%	112	48	0.42 (0.32 to 0.56)	1.93	1.12 (0.84 to 1.32)	91% (69% to 100%)
	(78–80 years)
Fourth year after vaccine eligibility	78–79 years	63%	81	55	0.67 (0.51 to 0.88)	2.00	0.66 (0.25 to 0.98)	53% (20% to 78%)
	(78–80 years)
Fifth year after vaccine eligibility	79 years	59%	25	13	0.51 (0.30 to 0.88)	2.06	1.01 (0.24 to 1.45)	82% (20% to 100%)
	(79–80 years)
All years of the programme	74–79 years	57%	536	336	0.62 (0.55 to 0.69)	1.88	0.71 (0.59 to 0.84)	66% (54% to 78%)
	(78–80 years)

*Calculated by taking the mean of the cumulative uptake values in each month from October to September of the relevant years and cohorts.

†Expected if vaccine not introduced; based on model results for unvaccinated cohorts.

‡Estimated from the Poisson regression model with a quadratic time trend, quadratic age effect and the factor for vaccine eligibility.

§Effectiveness required to generate the observed reduction given the impact and coverage.

GP, general practice.

Vaccine impact on hospitalisations

A total of 117·5 million person-years of data for those diagnosed with herpes zoster and PHN aged 60 to 89 from October 2008 to September 2018 were included in the analysis of impact. Table 2 shows person-years and herpes zoster/PHN events contributed by cohorts eligible and not eligible for vaccination by year and age group. Modelling the hospitalisation incidence data indicated that the age followed a cubic trend and time-period a quadratic trend. The cubic age term allowed for a small decline in PHN incidence in the oldest ages (>85) following the usual increase with age. The time trend showed increases over the period (table 2). As with GP consultations no significant interactions were identified and the models with factors for age and time effects and with GP random effects gave similar results.

Figures 5 and 6 show the observed and predicted incidence from the models in the absence of vaccination for herpes zoster and PHN by year (2008/2009 to 2017/2018) for those cohorts aged between 66 years and 81 years as of September 2013. The figures clearly show the impact in the first year each cohort is targeted for vaccination and in all the following years. Observed and predicted cases and relative incidence estimates from the fitted models are summarised for herpes zoster (table 5) and PHN (table 6). Reduction for herpes zoster incidence were 23% (range: 19 to 27) across the routine cohorts and 29% (range: 26 to 32) across the catch-up cohorts in the 5 years of the programme. These relative reductions translated into incidence reductions of 0·78 per 10 000 for routine and 2·1 per 10 000 for catch-up cohorts and a total of 1840 prevented herpes zoster (including PHN) hospitalisations for the period 2013/2014 to 2017/2018.

Table 5

Impact of routine and catch-up herpes zoster vaccination on hospitalised herpes zoster by time since cohorts were first eligible for vaccination

Period	Age on 1 Sept 2013 (age when first eligible for vaccination)	Average cumulative uptake*	Expected events†	Observed events	Incidence rate ratio‡(95% CI)	Expected incidence per 10 000 person-years	Incidence reduction per 10 000 person-years(95% CI)	Vaccine effectiveness§(95% CI)
Routine cohorts
First year after vaccine eligibility	66–70 years	47%	823	686	0.83 (0.77 to 0.90)	2.94	0.50 (0.30 to 0.68)	36% (22% to 49%)
	(70–71 years)
Second year after vaccine eligibility	67–70 years	64%	658	519	0.78 (0.72 to 0.86)	3.27	0.71 (0.47 to 0.92)	34% (22% to 44%)
	(70–71 years)
Third year after vaccine eligibility	68–70 years	71%	538	395	0.73 (0.66 to 0.81)	3.66	0.99 (0.71 to 1.25)	38% (27% to 48%)
	(70–71 years)
Fourth year after vaccine eligibility	69–70 years	75%	394	311	0.78 (0.70 to 0.87)	4.10	0.90 (0.51 to 1.25)	29% (17% to 40%)
	(70–71 years)
Fifth year after vaccine eligibility	70 years	77%	205	116	0.56 (0.46 to 0.67)	4.61	2.04 (1.52 to 2.47)	57% (43% to 69%)
	(70–71 years)
All years of the programme	66–70 years	62%	2619	2027	0.77 (0.73 to 0.81)	3.41	0.78 (0.65 to 0.91)	37% (31% to 43%)
	(70–71 years)
Catch-up cohorts
First year after vaccine eligibility	74–79 years	48%	1277	976	0.76 (0.71 to 0.81)	6.26	1.49 (1.16 to 1.79)	50% (39% to 60%)
	(78–80 years)
Second year after vaccine eligibility	75–79 years	62%	1136	770	0.67 (0.63 to 0.73)	7.03	2.29 (1.93 to 2.63)	52% (44% to 60%)
	(78–80 years)
Third year after vaccine eligibility	77–79 years	63%	961	656	0.68 (0.62 to 0.73)	7.92	2.57 (2.12 to 2.98)	51% (42% to 59%)
	(78–80 years)
Fourth year after vaccine eligibility	78–79 years	63%	756	554	0.72 (0.66 to 0.79)	9.00	2.50 (1.90 to 3.05)	44% (34% to 54%)
	(78–80 years)
Fifth year after vaccine eligibility	79 years	59%	251	180	0.70 (0.61 to 0.82)	10.17	3.02 (1.86 to 4.01)	50% (31% to 66%)
	(79–80 years)
All years of the programme	74–79 years	57%	4382	3136	0.71 (0.68 to 0.74)	7.36	2.10 (1.88 to 2.32)	50% (45% to 55%)
	(78–80 years)

*Calculated by taking the mean of the cumulative uptake values in each month from October to September of the relevant years and cohorts.

†Expected if vaccine not introduced; based on model results for unvaccinated cohorts.

‡Estimated from the Poisson regression model with a quadratic time trend, cubic age effect and the factor for vaccine eligibility.

§Effectiveness required to generate the observed reduction given the impact and coverage.

Table 6

Impact of routine and catch-up herpes zoster vaccination on hospitalised postherpetic neuralgia by time since cohorts were first eligible for vaccination

Period	Age on 1 Sept 2013 (age when first eligible for vaccination)	Average cumulative uptake*	Expected events†	Observed events	Incidence rate ratio‡ (95% CI)	Expected incidence per 10 000 person-years	Incidence reduction per 10 000 person-years (95% CI)	Vaccine effectiveness§ (95% CI)
Routine cohorts
First year after vaccine eligibility	66–70 years	47%	204	154	0.75 (0.64 to 0.88)	0.73	0.18 (0.08 to 0.27)	53% (25% to 77%)
	(70–71 years)
Second year after vaccine eligibility	67–70 years	64%	170	133	0.78 (0.65 to 0.93)	0.84	0.19 (0.06 to 0.29)	34% (11% to 54%)
	(70–71 years)
Third year after vaccine eligibility	68–70 years	71%	143	92	0.64 (0.52 to 0.79)	0.97	0.35 (0.21 to 0.47)	51% (30% to 68%)
	(70–71 years)
Fourth year after vaccine eligibility	69–70 years	75%	108	63	0.58 (0.45 to 0.75)	1.12	0.47 (0.28 to 0.62)	56% (34% to 73%)
	(70–71 years)
Fifth year after vaccine eligibility	70 years	77%	57	30	0.52 (0.36 to 0.75)	1.28	0.61 (0.32 to 0.82)	62% (32% to 83%)
	(70–71 years)
All years of the programme	66–70 years	62%	682	472	0.69 (0.62 to 0.77)	0.89	0.27 (0.21 to 0.33)	50% (38% to 61%)
	(70–71 years)
Catch-up cohorts
First year after vaccine eligibility	74–79 years	48%	343	274	0.80 (0.70 to 0.90)	1.68	0.34 (0.16 to 0.5)	42% (20% to 62%)
	(78–80 years)
Second year after vaccine eligibility	75–79 years	62%	303	191	0.63 (0.54 to 0.73)	1.88	0.7 (0.51 to 0.86)	60% (44% to 74%)
	(78–80 years)
Third year after vaccine eligibility	77–79 years	63%	252	176	0.69 (0.59 to 0.81)	2.08	0.63 (0.39 to 0.84)	48% (30% to 64%)
	(78–80 years)
Fourth year after vaccine eligibility	78–79 years	63%	192	141	0.73 (0.61 to 0.87)	2.28	0.62 (0.3 to 0.89)	43% (21% to 62%)
	(78–80 years)
Fifth year after vaccine eligibility	79 years	59%	60	38	0.62 (0.45 to 0.86)	2.45	0.93 (0.34 to 1.35)	64% (24% to 93%)
	(79–80 years)
All years of the programme	74–79 years	57%	1150	820	0.72 (0.66 to 0.78)	1.93	0.55 (0.43 to 0.66)	49% (38% to 59%)
	(78–80 years)

*Calculated by taking the mean of the cumulative uptake values in each month from October to September of the relevant years and cohorts.

†Expected if vaccine not introduced; based on model results for unvaccinated cohorts.

‡Estimated from the Poisson regression model with a quadratic time trend, cubic age effect, and the factor for vaccine eligibility.

§Effectiveness required to generate the observed reduction given the impact and coverage.

Observed hospitalised herpes zoster incidence and model-predicted incidence in the absence of vaccination by year for each age cohort. Cohorts are grouped by age at 1 September 2013. Each year runs from October 1 of the first year to September 30 of the second. Blue points indicate observed zoster incidence and red points indicate observed incidence in age cohorts eligible for vaccination. Red lines show the model-predicted incidence in the absence of vaccination. Note that ages 71, 72, 80 and 81 are cohorts not eligible in the period.

Observed hospitalised postherpetic neuralgia incidence and model-predicted incidence in the absence of vaccination by year for each age cohort. Cohorts are grouped by age at 1 September 2013. Each year runs from October 1 of the first year to September 30 of the second. Blue points indicate observed zoster incidence and red points indicate observed incidence in age cohorts eligible for vaccination. Red lines show the model-predicted incidence in the absence of vaccination. Note that ages 71, 72, 80 and 81 are cohorts not eligible in the period.

Reductions were greater for PHN (table 6) with a 31% reduction in PHN incidence in the routine cohorts across the first 5 years and a 28% reduction in the catch-up cohorts. This would be equivalent to reduction in PHN episodes of 0·27 and 0·55 per 10 000 person-years and gives an estimated number of PHN hospitalisations prevented of 540. Tables 5 and 6 also show the vaccine effectiveness of approximately 37% to 50% for herpes zoster and 50% for PHN, a little lower than that against GP consultations.

Discussion

In this paper, we provide important evidence of the sustained population impact of the herpes zoster vaccination programme on GP consultations for herpes zoster and postherpetic neuralgia, 5 years following its implementation in England. Despite some challenges in programme implementation and accurately identifying eligible individuals for vaccination, our evaluation demonstrates the continued impact observed in the first cohorts targeted for vaccination in 2013/2014, indicating that the effectiveness of the vaccine programme has been maintained. This impact is equivalent to approximately 40 500 fewer GP consultations for herpes zoster and herpes zoster-related disease including 8700 fewer PHN consultations in the first 5 years of the national programme. These levels of disease prevention are particularly relevant in the context of a slowly increasing incidence of herpes zoster over time, reflecting changes in population demographics and is similar to trends observed in other countries.18

For the first time, we additionally demonstrate the impact of the vaccine programme on hospital admissions for herpes zoster and PHN with the impact of hospitalised herpes zoster at a rate similar to that observed for GP consultations for herpes zoster. In total, an estimated 1840 fewer herpes zoster hospitalisations including 540 PHN consultations were prevented in the first 5 years of the programme in England. The reduction in GP consultations and hospitalisations will have large cost savings. Using a cost of £75·60 for a non-PHN GP case, £340·00 for a PHN GP case and £2816 for a hospitalised herpes zoster or PHN case gives a total cost saving of over 10·5 million pounds (£75·60×31 800 + £340×8700 + £1840×2816 = £10 543 520) over the first 5 years of the programme. These costs per case are derived from those used in recent cost-effectiveness analysis (Alicia Rosello, personal communication) and the age distribution of the cases in our analysis. The average annual reduction in hospital costs £1·6 million per year ((£340×8700 + £1840×2816)÷5) compares to an estimated cost prior to vaccination for herpes zoster hospitalisation of £13 million per year across all ages.19 These estimates of reduced burden and cost are likely to be conservative because they do not account for any lack of sensitivity in the coded episodes, which is likely for PHN, and also do not allow for the fact that some individuals will have had repeat episodes with a year of one another as these were excluded to ensure episodes included were independent.

In our evaluation, the vaccine effectiveness against herpes zoster to generate the observed incidence reductions in GP consultations in the first 5 years following vaccination is approximately 50% to 60%. This compares with long-term efficacy estimated from the long-term persistence study of approximately 38% and suggests that the protection from the vaccine does not wane as rapidly in clinical practice compared with the trial settings.4 5 In an 8-year follow-up cohort study in the US between 2007 and 2014, VE was 49·1% (95% CI: 47·5 to 50·6) across all follow-up and VE waned from 67·5% during the first year after vaccination to 31·8% (95% CI: 15·1 to 45·2) by year 8.20 Our higher estimates may reflect differences in threshold for health seeking behaviour and the fact that milder cases may not present to general practice, although a US cohort study found lower effectiveness in outpatient data with a decline from 38% (95% CI: 37 to 40) within a year of vaccination to 21% (95% CI: 11 to 31) after seven or more years.21 That study found higher estimates and less waning for hospitalised zoster (VE 77% (95% CI: 39 to 67) within a year and 55% (95% CI: 39 to 67) at four or more years). Vaccine effectiveness against PHN was notably higher at 75% for routine cohorts and 66% for catch-up cohorts in the first 5 years of the programme in England. This is more closely aligned to the efficacy generated in the clinical trials with long-term efficacy of 67%, as well as post licensure US studies of PHN with VE at about 50% between 4 and 7 years after vaccination.21 22 This effectiveness is consistent with evidence that the major benefits from this vaccine are through preventing the complications from severe herpes zoster, rather than preventing herpes zoster itself. Our VE estimates are also similar to those estimated through a formal evaluation in the UK using a cohort design based on an alternative primary care data source. In that study VE was 64% (95% CI: 60% to 68%) against incident herpes zoster and 81% (95% CI: 61% to 91%) against PHN, with very similar VE estimates in the routine and catch-up cohorts and with some evidence of waning for the first 3 years of the programme.8 Similar results have also been reported in other primary care effectiveness studies in the UK with VE seen not to vary by age but being lower in those with type 2 diabetes or prior Zoster.23 24A key strength of this study is the use of a large population-based sentinel primary care data source with robust historical data on consultations which are geographically representative of the country and provide consistent recording over time. In addition, the use of hospital admissions data across the whole of England add further weight to this evaluation. Despite this, it is important to note however, that in our study, time since first eligibility for vaccination is only an approximate measure of time since vaccination, particularly in the routine cohorts who continue to be eligible. Nevertheless, the continued low levels in the catch-up cohorts are particularly reassuring.

While, the findings from this 5-year evaluation are reassuring it is of concern that there has been a declining trend in vaccine uptake. Efforts to improve uptake have included simplifying the eligibility criteria for vaccination from April 2017 (ie, eligible on 70th and 78th birthday); although there are challenges in accurately monitoring uptake since this change, it is notable that as yet, this does not appear to have translated into a measurable improvement in uptake. To encourage improved uptake RCGP RSC has added feedback to member practices about their vaccine administration compared with the rest of the network as well as about rates of shingles and PHN; these data are also reported as a national observatory and show impact through the lower proportion of herpes zoster and PHN cases vaccinated than the whole eligible population.25 26

In 2019, the UK JCVI made specific recommendations for the use of a new recombinant subunit vaccine, Shingrix.27 This vaccine, licensed as a two-dose schedule, administered at 2- to 6-month intervals, has been shown to have efficacy of >90% in older adults with limited waning to 4 years and has been introduced into the US programme in 2018.28 As a recombinant vaccine, this vaccine is being recommended in the UK, for use in immunocompromised individuals known to be at high risk of herpes zoster and currently unable to benefit from the existing herpes zoster vaccination programme. In addition, Shingrix has been recommended for use in the routine programme for adults aged 60 to 70 years based on a recent cost-effectiveness analysis. However due to the limited availability to meet high global demand, it is currently not known when Shingrix will be available for use in the UK programme.

The findings from this evaluation of the current herpes zoster vaccination programme in England however should provide considerable reassurance to health professionals and patients and serves as a reminder of the significant health benefits conferred to a target population that is both increasing in size and that is particularly affected from the severe consequences of herpes zoster and herpes zoster-related complications.