Abstract
In January 2020, the United Kingdom moved to a 1+1 schedule for the 13-valent pneumococcal conjugate vaccine (PCV13) with a single priming dose at 3 months and a 12month booster. We modelled the impact on invasive pneumococcal disease (IPD) out to 2030/31 of reductions in PCV13 coverage and population mixing associated with restrictions on non-essential health care visits and social distancing measures introduced in 2020/21 to reduce SARS-CoV-2 transmission.
Using an existing model of pneumococcal transmission in England and Wales we simulated the impact of a 40% reduction in coverage and a 40% reduction in mixing between and within age-groups during two lockdowns in spring 2020 and autumn/winter 2020/21. More and less extreme reductions in coverage and mixing were explored in a sensitivity analysis.
Predicted annual numbers of IPD cases under different coverage and mixing reduction scenarios with uncertainty intervals (UI) generated from minimum and maximum values of the model predictions using 500 parameter sets.
The model predicted that any increase in IPD cases resulting from a reduction in PCV13 coverage would be more than offset by a reduction in pneumococcal transmission due to social distancing measures and that overall reductions in IPD cases will persist for a few years after resumption of normal mixing. The net reduction in cumulative IPD cases over the five epidemiological years from July 2019 was predicted to be 13,494 (UI 12,211, 14,676) all ages. Similar results were obtained in the sensitivity analysis.
COVID-19 lockdowns are predicted to have had a profound effect on pneumococcal transmission resulting in a reduction in pneumococcal carriage prevalence and IPD incidence for up to five years after the end of the lockdown period. Carriage studies will be informative in confirming the predicted impact of the lockdown measures after they have been lifted.
Introduction
Following the first identification of local transmission of SARS-CoV-2 in England on February 28, 2020, [1] the United Kingdom (UK) government instituted a series of control measures to reduce spread of the virus. On March 18 schools along with restaurants and places of entertainment and leisure were closed, and on March 23 all non-essential travel was banned with members of the population only permitted to leave their home for exercise once a day, or in order to obtain food or medication [2]. The public was advised to only use the National Health Service when essential [3] and many general practices replaced face-to-face consultations with an online service [4]. This first national lockdown was eased in stages starting in May 2020 with children returning to schools in September 2020 [5]. As COVID-19 cases started to increase again in September 2020, a variety of locally-based measures were introduced in areas of increasing incidence in an attempt to reduce transmission. Despite these local measures cases continued to increase and a second national lockdown came into force in November 2020 which lasted for six months until May 2021 [5].
In the UK, paediatric vaccinations are given in general practice and there was initial concern that the measures taken to limit spread of SARS-CoV-2 would result in many children’s vaccinations being postponed with a resulting increase in vaccine-preventable infections [6]. Indeed early coverage data for measles/mumps/rubella vaccine for the first 3 weeks of the first lockdown in London suggested a drop of over 40% compared with the same period in 2019 [7]. Of particular concern is the impact of a drop in coverage of the 13 valent pneumococcal conjugate vaccine (PCV13) for which the national 2+1 schedule (2, 4, and 12 months) was changed on 1 January 2020 to a 1+1 schedule at 3 and 12 months [8]. The rationale for the schedule change was that after a decade of high coverage with PCV13 much of the protection of young children against vaccine-type invasive pneumococcal disease (IPD) is through herd immunity generated by the 12 month booster dose and that dropping a priming dose would have little effect on this. The supporting evidence for the schedule change came from an immunogenicity study showing similar post-booster antibody responses to PCV13 serotypes with a 2+1 and 1+1 schedule [9] and a modelling study which predicted little impact on IPD cases in children or adults if protection against vaccine-type carriage was similar post-booster with the two schedules [10].
Using our previously described pneumococcal transmission model [10] we investigated the impact on IPD if PCV13 coverage was substantially reduced during the lockdown periods in the UK. The model also allowed us to investigate the likely longer-term impact on pneumococcal transmission of social mixing restrictions put in place to limit the spread of the SARS-CoV-2 virus in England and Wales.
Model and assumptions
The same compartmental deterministic model and estimated parameter values as described by Choi et al. [10] were used in this study. Briefly, the model simulates the transmission dynamics of carried pneumococci between and within different age groups and predicts carriage changes by serotype grouping, either vaccine type (VT) or non-vaccine type (NVT), under different vaccination scenarios. Carriage changes are then translated into IPD cases using case-carrier ratios (CCRs) that describe the invasive potential of different carried serotypes. The model has a Susceptible-Infected-Susceptible (SIS) structure as it is assumed that a carriage episode does not generate any serogroup-specific protection. The model has 26 parameters which describe the propensity for vaccine-induced replacement of VTs with NVTs in carriage, the CCRs in different age groups, and the effectiveness of PCV vaccination against carriage of different vaccine serotypes. (see Supplementary Appendix for summary of model description and parameter estimation). The maximum likelihood values of these parameters and their uncertainty intervals are estimated by fitting the model to IPD data up to the end of the epidemiological year from July 2015 to June 2016 [10]. In implementing the effect of the 1+1 schedule, the base case in the model described by Choi et al. [10] assumed that one priming dose at 3 months of age provides half the protection against VT carriage as two doses while the effect of the 12 month booster remains unchanged. In the current analysis, the change to the 1+1 PCV13 schedule was implemented to reflect the actual start date of January 2020, not the assumed date of September 2018 as in the original model [10].
We assumed that, during the first two-month lockdown period starting from 23rd March 2020 and also during the second six-month lockdown starting from 5th November 2020, there was a 40% reduction in PCV13 coverage, as suggested by early declines in MMR vaccine coverage in London during the first lockdown [7]. For the sensitivity analysis we assumed two additional scenarios, one with no material decline in PCV13 coverage (0% decline) and a second with cessation of all PCV13 vaccination during the lockdown periods (100% decline). The 0% scenario was chosen to reflect a situation in which there was a rapid catch-up of missed PCV13 vaccinations after the first lockdown with no reduction during the second while the 100% reduction scenario was explored as an extreme scenario. We pessimistically assumed there would be no subsequent catch up vaccination for children under the 40% and 100% reduction scenarios as final coverage data for the affected birth cohorts will not be available until 2021/22.
To implement the impact of reduced social mixing during each of the lockdown periods we used the same mixing matrix as before [10] based on the POLYMOD study [11] and an additional contact survey among infants under one year [12] but reduced the number of contacts between and within each age group by 40% as a baseline scenario and explored 10% and 74% reduction scenarios for the sensitivity analysis. The 74% reduction scenario was based on early data from a study of the social contacts of 1,356 UK adults between 24 and 27 March [13] compared with the historical rates documented in the POLYMOD study [11]. However, this study [13] could not assess the impact on contacts made by children, particularly those of pre-school age who are the main drivers of pneumococcal transmission [10]. As the social distancing measures in children were less strict than adults, the estimated reduction in mixing would be lower in children and so we took the 74% scenario as the upper extreme scenario for the contact pattern reduction.
Results are shown as the median of the model outputs under the different scenarios. The uncertainty interval (UI) for each scenario represents the minimum and maximum values generated with 500 model parameters sets obtained from Choi et al. [10]. The impact on IPD of the temporary changes in PCV13 coverage and social mixing is presented by VT, which includes PCV13 serotypes minus serotypes 1 and 3, and NVT; serotype 3 is included with the NVT group in the model on the grounds that it has shown an increase post-PCV13 similar to that seem for NVTs while serotype 1 is excluded from the model due to secular changes which are independent of vaccine usage [10].
Patient and Public Involvement statement
Our study did not involve any direct participation by human subjects.
Results
Fig 1 shows the predicted number of IPD cases in England and Wales by epidemiological year out to 2030/31 in under 5 year olds and in 65+ year olds by VT and NVT group with four different scenarios. The first scenario presented in Figs 1 A&B is without any change in PCV13 coverage and contact patterns; the increase in IPD cases out to 2030/31 in 65+ year olds reflects the increasing mean age of individuals in this age group in the future together with the steep increase in IPD incidence with increasing age in the elderly population. The remaining three scenarios in Figs 1 C-H assume a 40% reduction in PCV13 coverage with respectively a 10% (C&D), 40% (E&F) or 74% (G&H) reduction in contact rates during the two lockdowns. Even with only a 10% reduction in contact rates in Figs 1 C&D, the predicted impact of reduced social mixing outweighed any increase in IPD due to reduction in PCV13 coverage. The predicted beneficial effect of the lockdown on IPD cases was even more marked with higher reductions in social mixing (Figs 1 E, F, G and H). Under all the scenarios, the model predicted that the impact of the two lockdowns on IPD would persist for up to five years after their cessation.
The cumulative number of additional IPD cases predicted in England and Wales over the first five epidemiological years starting from July 2019 is shown in Table 1 by age group and serotype grouping for a 40% reduction in contact rates and 40% reduction in PCV13 coverage during two lockdowns (base case scenario). Across all ages the model predicted a net saving of 13,494 IPD cases (UI 12,211, 14,676) over this period compared with the predicted number of IPD cases without any lockdowns. The results over the first five years by age group and different scenarios for the reduction in coverage and social mixing are presented in Table S1. Without any reduction in contact rates under the worst case of no PCV13 vaccination during the two lockdown periods and no catch up vaccination, the model predicted a cumulative additional 105 (UI 49, 197) VT IPD cases, of which 14 (UI 6, 26) were in children under 2 years of age over the first five epidemiological years since 2019/20. The overall predicted increase in IPD cases was smaller (73 (UI 33, 130)) due to a concomitant reduction in replacement with NVT IPD cases (Fig 2). While the predicted overall increase is small the perturbation in VT IPD cases as a result of the temporary interruption of vaccination is predicted to persist out to 2030/31.
Amin-Chowdbury et al. [14] reported that there were 2,440 IPD cases in England and Wales between February and June 2019 and 1,137 IPD cases during the same period in 2020 showing a 53% reduction in overall IPD cases in England and Wales. The model predicted that there would be 2,350 (UI 2,129, 2,531) IPD cases in the same period in 2019.
Table 2 shows the model predictions for the reduction in IPD cases in February to June 2020 with a 40% reduction in PCV13 coverage during the two lockdown periods and under various reduced mixing scenarios. The scenario with a 40% reduction in mixing pattern most closely matched the observed reduction in IPD during February to June 2020 [14]. Table S2 shows the predicted number of IPD cases during February to June 2020 with different assumptions about reductions in PCV13 coverage and in social mixing. Only the scenarios with a 40% reduction in mixing produced the observed reductions seen in IPD incidence with little impact of varying the PCV13 coverage reductions.
Discussion
In line with many other countries, the UK adopted a lockdown strategy to minimise social contacts and thereby reduce the spread of the SARS-CoV-2. The restriction on non-essential travel and GP visits raised concerns that vaccine coverage in children might be substantially reduced [6], echoing similar concerns raised at the global level [15]. With the recent change to a 1+1 PCV13 schedule in the UK, which is reliant on the booster dose to maintain herd immunity [9,10], a substantial reduction in PCV13 coverage could compromise control of IPD. Our modelling study suggests that any potential impact of the lockdown measures on PCV13 coverage and herd immunity will be more than offset by a reduction in pneumococcal transmission due to the reduction in population mixing such that a net reduction in IPD cases should occur. Our model also predicts that the net reduction in IPD cases will be sustained for up to five years after mixing has returned to normal levels in the population. This reduction in IPD would be the result of an overall reduction in pneumococcal carriage prevalence after the two lockdown periods, unlike the IPD reduction achieved with PCV. Carriage studies in England before and after the introduction of PCVs have shown little impact on overall carriage prevalence with the reduction in carriage of vaccine serotypes offset by an increase in non-vaccine serotypes [16]. The overall reduction in IPD cases achieved by PCVs is the consequence of the generally lower CCRs of the NVTs that replaced VTs in carriage.
Invasive pneumococcal disease is a consequence, albeit rare, of carriage of Streptococcus pneumoniae in the nasopharynx. Acquisition of the pneumococcus is by person-to-person transmission and occurs via respiratory droplets with transmission enhanced in settings with close mixing such as day care centres and military camps [17]. Visits to the GP for mild upper respiratory tract disease have also been shown to be a risk factor for pneumococcal acquisition [18]. It is therefore to be expected that the measures taken to reduce spread of SARS-CoV-2 would also reduce pneumococcal transmission. The impact on IPD in England and Wales of the first UK lockdown that lasted for two months from 23 March 2020 was recently reported [14]. Our model simulations of the impact of the first lockdown were similar and consistent with an overall reduction of 40% in contact patterns during this period.
Whether a reduction in coverage caused by the lockdowns would have a similar effect on other vaccine-preventable infections cannot be assessed with our model as it is specific to pneumococcal infection. However, the average number of weekly notifications of whooping cough and meningococcal septicaemia in England for one year from 23 March 2020 showed a 91.3% and 74.7% reduction respectively compared with the weekly average for the same period in 2019 [19]. As with IPD, nasopharyngeal colonisation is a pre-requisite for whooping cough and meningococcal disease. While the completeness of notifications by doctors may have suffered somewhat as a result of the COVID-19 epidemic, the magnitude of the decline in these two other infections is consistent with a substantial reduction in new colonisation episodes following the lockdown and suggests that, at least in the short term, this has outweighed any potential increase due to a reduction in coverage of pertussis or meningococcal vaccines.
Our model has limitations as previously described [10]. These include uncertainty over the future behaviour of non-vaccine serotypes and whether serotypes with different invasiveness potential will emerge, and also in the mixing patterns among the oldest age groups. However, given the magnitude of the predicted effect of the lockdowns, these uncertainties are unlikely to materially affect the results with respect to the effect of the lockdowns. It is also possible that the reduction in pneumococcal carriage predicted to occur as a result of the lockdowns may affect induction of natural immunity. Nasopharyngeal carriage can boost serotype-specific IgG antibody levels to the carried serotype and such naturally induced capsular antibodies have been shown to protect against carriage of that serotype [20]. If natural immunity against carriage was reduced as a result of decreased pneumococcal transmission this could offset the reduction in IPD predicted to persist after relaxation of the lockdown. Our model has an SIS structure so does not include a natural immunity compartment, so we were unable to explore the potential effect of such perturbations in carriage-induced immunity. We were also unable to incorporate the actual coverage reductions during the two lockdown periods as coverage data for completed 1+1 PCV13 courses is only available when a birth cohort reaches two years of age. Our baseline coverage reduction scenario may therefore have been extreme with catch-up vaccination rapidly offsetting any reduction during the two lockdown periods. However, as shown by the sensitivity analysis, the critical parameter in determining impact was reductions in mixing, with only a 10% reduction having a profound effect on IPD cases.
In conclusion, our pneumococcal model predicts that the effect of the COVID-19 lockdowns on pneumococcal carriage rates will result in a reduction in IPD cases over a five year period, thus obscuring the predicted small increase in VT IPD due to the adoption of a 1+1 PCV13 schedule [10]. Pneumococcal carriage studies conducted over the next few years could help elucidate the impact of the social distancing measures and of the adoption of the 1+1 PCV13 schedule on the future incidence of VT and NVT IPD.
Data Availability
The data used for this study are available on the previous published paper on https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.1002845.