Abstract
Background To reduce the coronavirus disease burden in England, along with many other countries, the Government implemented a package of non-pharmaceutical interventions (NPIs) that have also impacted other transmissible infectious diseases such as norovirus. It is unclear what future norovirus disease incidence is likely to look like upon lifting these restrictions.
Methods Here we use a mathematical model of norovirus fitted to community incidence data in England to project forward expected incidence based on contact surveys that have been collected throughout 2020-2021.
Results We report that susceptibility to norovirus infection has likely increased between March 2020 to mid-2021. Depending upon assumptions of future contact patterns incidence of norovirus that is similar to pre-pandemic levels or an increase beyond what has been previously reported is likely to occur once restrictions are lifted. Should adult contact patterns return to 80% of pre-pandemic levels the incidence of norovirus will be similar to previous years. If contact patterns return to pre-pandemic levels there is a potential for the expected annual incidence to be up to 2-fold larger than in a typical year. The age-specific incidence is similar across all ages.
Conclusions Continued national surveillance for endemic diseases such as norovirus will be essential after NPIs are lifted to allow healthcare services to adequately prepare for a potential increase in cases and hospital pressures beyond what is typically experienced.
Background
In late January 2020 COVID-19 was first detected in the UK [1], and in response to increasing incidence and hospitalisations, the UK entered the first of three national lockdowns on the 23rd of March 2020. The lockdowns, and other social restrictions resulted in a reduction in contact patterns in the community, as shown through self-reported contact rates collected using online surveys [2] and other metrics [3]. Subsequently, reports to national surveillance for many endemic diseases, including norovirus, declined across England throughout the rest of 2020 [4]. This reduction in reports of infectious disease has been consistent across many pathogens and observed in other countries where non-pharmaceutical interventions (NPIs) for the control of COVID-19 have been implemented (eg. Australia [5, 6] and the USA [7, 8]). With the introduction of vaccination against COVID-19 in England a considerable reduction in COVID-19 hospitalisations and deaths have been observed and this intervention has superseded the blunt tool of NPIs. The UK roadmap out of lockdown details the gradual relaxation of NPIs in 2021 [9], which will likely result in contact patterns more consistent with those seen prior to 2020.
Norovirus is an endemic viral infection present globally across high- and low-income countries with – as yet - no commercially available vaccine [10]. Norovirus infection can be asymptomatic or can result in gastrointestinal disease symptoms such as stomach cramps, diarrhoea and vomiting. In more vulnerable groups, such as the elderly or immunocompromised, norovirus can result in more severe and/or prolonged disease, and affected individuals may require hospitalisation for treatment. Norovirus is the most common cause of gastrointestinal infections globally and has the highest burden of disease for intestinal infections in the UK [11]. Infections and outbreaks occur more frequently in the winter months, and in England are monitored through multiple national surveillance systems (eg. [12]) to ensure that unusual activity is detected, and that alerts to local health authorities are made for preventive actions. The probability of symptomatic disease is also influenced by host genetic factors [13], the infectious dose of norovirus [14], and previous exposure to infection. Norovirus consists of multiple genotypes and strains; since 2012 the GII.4/2012 strain has dominated in reports of symptomatic disease in England and globally. Prior to the GII.4/2012 strain, population strain replacement occurred every few years and was consistent across wide geographies. The majority of norovirus infection is considered to occur via community transmission, and consequently changes in community contact patterns will likely influence norovirus transmission. However, norovirus has multiple routes of transmission; for example foodborne transmission is thought to contribute to approximately 16-35% in the UK, based on microbial risk assessments and detailed case investigations [15].
Here we illustrate through modelling the impact of COVID-19 NPIs used throughout 2020-2021 on the dynamics and incidence of norovirus disease in the community, leading on from a recent report of a reduction in norovirus activity [4]. We provide indicators of surveillance to inform what early warning of a norovirus resurgence may look like.
Methods
Modelling norovirus infection prior to the COVID-19 pandemic (up until March 2020)
We built an age-structured SEIR-like model that follows a previously developed model for norovirus, focussing on infections with the GII.4/2012 strain [16] (Figure 1). We assume a population of 100,000 individuals with age-structure similar to that reported in England, and that 20% of the population remain resistant to norovirus infection due to nonsecretor status of human histoblood group antigen carbohydrates [17]. Upon infection individuals are assumed to enter a short stage of exposure (or pro-dromal infection), and then become fully infectious and symptomatic for on average of 2 days [18, 19] (full details of the model are given in SM). Individuals then enter the asymptomatic stage, where they remain moderately infectious [20] for on average 15 days when compared to the symptomatic stage but have no symptoms of disease [21]. Asymptomatic infection is assumed to correspond with norovirus shedding in stool that can be detected using RT-PCR. Upon recovery, individuals are immune to further symptomatic infection, but can develop asymptomatic infection. After an average of 5.1 years individuals are assumed to return to the susceptible class where infection will be symptomatic again [16]. To capture heterogeneity in transmission between ages, the transmission model is age-structured using eight age groups (0-4,5-14,15-24,25-34,35-44,45-54,55-64,65+). Symmetric contact rates between and within age groups prior to March 2020 were obtained for England from the European Polymod Study [22] and adapted for the age ranges used in here (using the R package socialmixr). All modelling and analyses were done in the software R (v 4.0.4).
A challenge of norovirus modelling is that there are broad observations in the epidemiology of norovirus which a model should ideally replicate, but these observations are from a wide spectrum of sources, with limitations and biases within each dataset. In the UK there have been two major studies in recent decades into infectious intestinal diseases (IID), which provide valuable information on community incidence. The first Study of Infectious Intestinal Disease (IID1) in England illustrated that norovirus is a common cause of IID in the community with considerable under-reporting to primary care [23]. A follow-on study (IID2) explored community incidence of gastrointestinal diseases (including norovirus), where incidence was highest in children and declined with age [24], and there was substantial under-reporting within the passive surveillance at the time. Serological data for norovirus is limited (due to challenges in establishing reliable correlates of protection [25] and enzyme immunoassays often being insensitive to strain-specific immunity), but where available, serological responses indicate seroconversion after acute illness [26] and suggest an increasing titre with exposure. In order to capture these general observations while needing to obtain reliable transmission rates within the community, we fit the probability of transmission given a contact to community age-stratified incidence data from IID2 [27] under different assumptions stated above, and also include a 20% under-reporting of symptomatic cases to account for a possible increase in incidence from the GII.4/2012 strain which has dominated symptomatic illness since this time (see SM for full details). For each scenario we provide estimates of incidence (ie. new cases of symptomatic infection) which were compared directly to incidence data. Additional model outputs were provided as a form of validation against what we would expect a realistic model of norovirus to output; estimates asymptomatic prevalence, a proxy for seroprevalence (assumed here to be the proportion of individuals that have recently recovered from symptomatic infection) and estimates of R0. While the fit of the model to the data is the main criteria for model selection, we consider these additional metrics when making decisions on the appropriateness of each set of model assumptions.
While the above model assumptions align well with much of the evidence, there are uncertainties with some assumptions, and alternative assumptions may affect disease dynamics. As part of sensitivity analysis, we additionally include other models which capture the main uncertainties and compare the fit of these models to the data and other summary metrics (see SM for full details).
Modelling norovirus infection during 2020-2021
From 23rd March 2020 onwards we assume that the age-specific contact rates change as reported in a longitudinal cross-sectional study (‘Comix’) of in-person contacts in England [27]. Weekly measures of age-specific contact rates were collected from 24th March 2020 and the study will continue to collect additional surveys until September 2021. The weekly contact rates are grouped into nine distinct periods of time representing 3 different lockdowns that saw different degrees of NPIs being introduced and restrictions easing again with varying contact rates [28]. The model is first run to endemic equilibrium using pre-pandemic values of age-specific contact rates. Between 23rd March 2020 until 19th July 2021 a contact matrix from Comix is used for each period of time (specified in the SM). For later projections we assume that contact rates return to pre-2020 values, an additional scenario that assumes adults only return to 80% of their pre-pandemic contact rates. We report the expected incidence of symptomatic norovirus infections in time, and the proportion of the population susceptible to symptomatic infection.
To support planning and national surveillance, we translate the incidence reported within the model to expected norovirus cases reported to the national SGSS. The SGSS is the national laboratory reporting system which records positive records of causative agents such as norovirus in England from frontline diagnostic laboratories [12]. As norovirus is not a notifiable causative agent in England under the Health Protection (Notification) regulations 2010 [29]. Data were extracted on the 28th June 2021 for all specimen types from 1st January 2014 onwards, and were deduplicated, providing an estimate of the weekly number of norovirus cases reported to SGSS. From the model the weekly incidence of symptomatic norovirus per 100,000 was scaled up to national incidence (using Office for National Statistics data on population size for England in 2019 of 56,290,000 individuals). Using the estimates within Tam et al. [11] we assume that for 287.6 (95% CI 239.1 to 346) symptomatic cases estimated nationally one case is reported within the SGSS. The data were compared to the model and the percentage difference reported, noting that the under-ascertainment reported in Tam et al. are from a laboratory reporting surveillance system in place prior to SGSS meaning that some changes in under-ascertainment are expected.
Results
The age-specific incidence of norovirus observed in England is best replicated assuming that primary infection is symptomatic, the duration of asymptomatic infectiousness is 15 to 20 days, and that asymptomatic infection is considerably less infectious than symptomatic infection (Figure 2A and Table 1). Simulations that correspond to these observations resulted in a shedding prevalence less than 1%, the proportion of the population recently recovered from infection ranging from 24.54-30.55% and an R0 between 1.81-2.01. The age-specific incidence of symptomatic infection is highest in children aged <5 years, and incidence declines with age. When we compare this modelled incidence to the weekly SGSS five-year average (from norovirus year 2014/15 to 2018/19), a further 27% reduction in reporting is needed for cases to be equivalent (Figure 1B). Other models with different assumptions had a poorer fit to the data and were not taken further in the analysis.
From 23rd March 2020 (“lockdown 1”) until 4th September 2020 (the “School reopening” stage) the rate of infection for norovirus is sufficiently low that new infections are rare (Figure 2A), and an increase in population susceptibility is observed (Figure 2B). This period corresponds with minimal cases reported to the laboratory surveillance SGSS system (Figure S1). The period between 4th September 2020 to 5th January 2021 corresponds to a higher average number of contacts within the population (up to 7.79 in average number of (unweighted) contacts) and an increase of incidence in the model. This increase in incidence was not observed in the SGSS system. The subsequent “lockdown 3” period corresponds to a reduction in contacts (from 6.61 to 3.47) and the rate of infection falls to low levels again, until schools are re-opened on the 9th March 2021 (Figure 2A). By this time point the model estimates that the proportion of the population susceptible to symptomatic infection has risen from 54% to 59% (a 9% increase) as a result of waning immunity (Figure 2B, age-specific values given in SM). Subsequently, model scenarios predict a rise in the rate of infection and a resurgence of norovirus in the community resulting in an annual incidence of cases up to 2-fold higher than simulations prior to 2020 (Figure 3A, full details in SM). However, this prediction is dependent upon assumptions on mixing patterns in the general population; if adult mixing is assumed to be 80% of pre-COVID levels a increase in incidence is not predicted. Instead, the rate of infection increases but at a lower level, resulting in a gradual increase in susceptibility in the population.
We also tested the impact of model assumptions regarding under-reporting and the duration of shedding (Figure S2). Assumptions about under-reporting affect the rate of infection; if incidence was under-reported in Harris [27] this will result in a higher rate of infection in the community and the potential for a more rapid variation in incidence. Assumptions on the duration of shedding have a lesser impact; simulation assuming a shorter duration of shedding suggest a higher incidence in the 2021/22 norovirus year.
Discussion
The reduced reporting of norovirus cases and outbreaks that was observed in 2020 as a consequence of COVID-19-associated NPIs has been described within a number of countries, including England [4], Australia [5] and the USA [7, 8]. Our modelling study quantifies the potential unintended consequences of the build-up of individuals susceptible to norovirus infection in England, with epidemic potential for the autumn/winter of 2021 and 2022. Norovirus is an example of an endemic infectious disease where changes in contact rates can affect both the short and long-term incidence of symptomatic disease. This resurgence is likely to further increase “winter pressures” within the NHS that may lead to more hospital admissions due to norovirus, hospital outbreaks of norovirus with bed-days lost as beds are kept empty for infection prevention and control, and treatment delays for non-norovirus patients [30, 31], therefore preparations for mitigating the effect of an increase in incidence are essential. This is especially true at the current time where there has been considerable disruption to healthcare services due to the COVID-19 pandemic and there are a record number of 5 million people on waiting lists to receive hospital care following the impact of COVID-19 in England [32]. While there are several caveats of the modelling to consider in the interpretation of this analysis, this is the first analysis of disease dynamics for norovirus affected by changes in NPIs.
More broadly, the unintended consequences of NPIs are likely to be widespread across other endemic diseases (eg. respiratory syncytial virus and influenza) where incidence pre-COVID-19 was limited by populations being largely immune to infection [33, 34]. A challenge in interpretation of the likely future trajectory of incidence is quantification of how NPI impact physical contacts in the community and consequently how this impacts transmission of infectious disease. Here we make use of the longitudinal cross-sectional study of contacts, which has been instrumental in quantifying incidence of COVID-19 [35], and show that these data can be very useful for quantifying dynamics of endemic diseases in addition to COVID-19. Continued collection of these data and incorporation into models of infectious disease will be instrumental in informing our understanding of disease transmission as the pandemic continues to unfold. Our analyses make the implicit assumption that contacts described within surveys correlate with contacts relevant for infectious disease transmission where the primary mode is direct contact; while this is apparent for infections such as ‘flu and SARS-CoV-2, there is less evidence for norovirus. The reported incidence of norovirus in the coming months, the changing behaviours of the community in contacts and infection prevention, and the actions public health officials will take to prevent a large epidemic, will be the test of this hypothesis.
In the analysis the incidence of norovirus reported was replicated well in the simulations, with the exception of the expected increase in cases at the end of 2020. This potentially suggests that the rate of infection assumed in the model is higher than is experienced in the community, despite the fact that the rate of infection was fitted to community data.
Alternatively, the effect of reduced to the SGSS due to changes in healthcare seeking behaviours and reduced laboratory testing of samples obtained may have impacted on the comparator data. This observation means that some caution should be applied to the findings presented here as the model predictions are sensitive to assumptions about contact rates and the probability of transmission, which remain uncertain. Further work will explore this in more detail through the use of additional contact surveys and exploration of other datasets that can improve our understanding of norovirus transmission in the community.
The majority of norovirus disease occurs within the community; there is a significant disease burden resulting in school closures and days of work lost (especially as infection can be experienced relatively frequently) [36]. Hygiene measures, such as hand washing with soap, enhanced environmental cleaning and staying at home when ill, have been part of usual recommended practice to reduce norovirus transmission, but the COVID-19 pandemic has indirectly enforced the adherence to these measures. It is possible that beyond 2021 these measures may remain as usual practice in the short to medium term, which may well mean that incidence of infectious diseases such as norovirus is lowered compared with historical norms. An additional challenge of norovirus largely being a community illness is that norovirus surveillance is affected by both under-ascertainment and under-reporting, resulting in considerable uncertainty in true estimates of burden. While we compared the model output to cases reported into the SGSS, surveillance is biased towards cases within hospitals, which will be a combination of community and hospital acquired infections. For this reason we have fitted to data from a community cohort study and used the SGSS data as a validation of the model estimates. The data from IID2 is over 10 years old, and may not fully represent incidence experienced in recent years, due to strain replacement in 2012 and hypothesised higher incidence of infection. Indeed, a recent community cohort study in the Netherlands report an incidence of symptomatic norovirus of 339.4 cases per 1000 person-years [37], which is higher than that reported in the IID2. We aimed to overcome the weaknesses in data by having multiple models with varying assumptions to cover a wide range of plausible present-day scenarios of norovirus dynamics. We identified a difference between model predictions and SGSS of 73%, which is well within a plausible range.
A further complexity in norovirus epidemiology is the possibility of strain replacement in the coming months or years; the GII.4/2012 strain has been in circulation in England since 2013 and could continue to dominate, but the emergence of other variants also have potential [38]. The likelihood and how a new variant may impact norovirus dynamics is currently uncertain and is likely to depend on the relative virulence of a new strain and the heterotypic immunity from infection with the GII.4/2012 strain. Further work will include consideration of strain-specific immunity profiles within the community to understand the impact of multi-strain dynamics and the potential impact of vaccination.
Conclusions
This modelling study suggests that norovirus incidence in the community is likely to remain approximately similar to or substantially increase beyond what has been experienced in years prior to 2020. The results are sensitive to future contact patterns in the community, and the adherence to preventive actions, which will affect the probability of transmission. The lower incidence of norovirus reported in surveillance aligns with model assumptions on reduced contact rates from March 2020 to June 2021, and is consistent with an increase in population susceptibility. The short to long term impact of this increased susceptibility places populations at risk of norovirus disease, but the scale of the impact remains to be seen.
Continued investment to maintain robust national surveillance systems will remain critical to enable measures to limit the impact of these resurgences and provide essential information to public health bodies to inform the implementation of preventive actions.
Data Availability
The code and data used to conduct these analyses are found at https://github.com/kath-o-reilly/norovirus_NPIs. For data access to the SGSS data used in the study applicants can apply to the PHE team by emailing noroOBK@phe.gov.uk
Declarations
Ethics approval and consent to participate
This study was an analysis of de-identified data from both public and
Consent for publication
not applicable
Availability of data and materials
The code and data used to conduct these analyses are found at https://github.com/kath-o-reilly/norovirus_NPIs. For data access to the SGSS data used in the study applicants can apply to the PHE team by emailing noroOBK{at}phe.gov.uk
Competing interests
The authors declare no competing interests.
Funding
This research was primarily funded by the Wellcome Trust “Noropatrol” [grant code203268/Z/16/Z]. Additional funding is acknowledged from Public Health England (PHE), which is an executive agency of the Department of Health and Social Care (DHSC), and National Institute for Health, USA (NIAID 1R01 AI148260). CoMix is funded by the EU Horizon 2020 Research and Innovations Programme - project EpiPose (Epidemic Intelligence to Minimize COVID-19’s Public Health, Societal and Economical Impact, No 101003688) and by the Medical Research Council (Understanding the dynamics and drivers of the COVID-2019 epidemic using real-time outbreak analytics MC_PC 19065).
Authors’ contributions
KO, FS, JB, and JE conceived the study. KO, CI, AG, AD, LL, KW carried out data analysis and KO carried out the modelling. Interpretation of the data findings were done in discussion with KO, FS, DA, CJ, AD, LL, LL, RG JB and JE, and all authors contributed to the final draft of manuscript, and approved the for publication.
Supplementary Methods
1. Full description of the mathematical model for norovirus
We built an age-structured SEIR-like model that follows a previously developed model for norovirus, focussing on GII.4 infections [1]. We assume a population of 100,000 individuals with age-structure similar to that reported in England, and that 20% of the population remain resistant to norovirus GII.4 infection due to nonsecretor status of human histoblood group antigen carbohydrates [2]. Births are assumed to either enter the genetically resistant class with probability g (G in Figure 1), or the susceptible class (S) with probability 1-g. Upon infection individuals are assumed to enter a short stage of exposure (or pro-dromal infection), and then become fully infectious and symptomatic for on average 2 days [3, 4]. The rate of infection is assumed to be a function of the proportion of individuals in the symptomatic and infectious class and the asymptomatic, where the relative contribution of the asymptomatic individuals is defined by δ. Individuals then enter the asymptomatic stage, where they remain moderately infectious for on average 15 days when compared to the symptomatic stage but have no symptoms of disease [5]. Asymptomatic infection is assumed to correspond with norovirus shedding in stool that can be detected using PCR. Upon recovery, individuals are immune to further symptomatic infection, but can develop asymptomatic infection. After an average of 5.1 years individuals are assumed to return to the susceptible class where re-infection will be symptomatic again. We make the simplistic assumption of there being one norovirus variant present in the population where immunity to further infection impacts only the probability of being symptomatic.
The system of ODEs consist of j(j = l, …, A) age categories which have the same number of infectious disease compartments. It is assumed that trans mission from age group i to age j occurs as specified using the time-varying force of infection (λ) below;
2. Dates for transition to ‘Polymod’ and ‘Comix’ contact matrices
The dates provided here are the corresponding dates used in Jarvis et al [27].
3. Exploration of alternative model assumptions for norovirus
The model structure developed by Lopman et al. [16] was our starting point for a model of norovirus transmission. Among the research group discussions of the model assumptions led us to trial alternative hypotheses by varying the model structure and comparing the model fit to the data. The hypothesis we explored are stated below with details of the justification. The log-likelihood fit of the model to the data is used to explore the evidence for using this model in further simulations; if the fit is substantially worse the model is not taken further. However, it should be acknowledged that this is not an exhaustive exploration of these hypotheses; the available data only enables us to vary one parameter at a time. With additional evidence these hypotheses should be further explored.
B. the duration of asymptomatic shedding is on average 20 days in duration
Studies where individuals with recent norovirus infection or in controlled experiments indicate that norovirus can be shed up to 3 weeks after the onset of illness [4]. The duration of shedding varies between individuals, and the quantity of virus shed also varies throughout the time-course of infection.
D. The first norovirus infection can be asymptomatic
While asymptomatic infection is well documented in older ages, it is less evident in younger ages where one might be confident that the infection is the first experienced. Within Han et al [6] a norovirus outbreak is reported within a paediatric ward where asymptomatic infection is reported but the extent of infections are not specified. Reports from paediatric wards are very useful in understanding the natural history of first infections, but children within paediatric wards may not be representative of the general population due to the medical reasons for them being present on the ward. To test this hypothesis we specify a model where a proportion of births (f) are assumed to be ‘born into’ the recovered class, resulting in the first infection being asymptomatic.
E. Asymptomatic infection is as infectious as symptomatic infection
The symptoms of norovirus infection (vomiting and diarrhoea) likely facilitate transmission to susceptible individuals through direct and fomite contact, and outbreak investigation has illustrated that asymptomatic individuals do not meaningfully contribute to transmission [7]. The default model assumes that for any given time period, asymptomatic infection is 5% as infectious as symptomatic infection. However, in a community setting, especially where parents may be caring for infected children, the evidence for asymptomatic transmission is less clear. To capture this alternative hypothesis that especially in community settings transmission of asymptomatically infected people may have a considerable contribution to transmission, we trial alternative values of relative infectiousness ranging from 0.05 to 1.0, and assess the fit of the model to the data.
4. Description of the model fitting procedure
The available data on age-specific incidence consists of 4 datapoints (in units of 1000-person-years), consequently only the probability of transmission given an infected contact is estimated. We assume that the likelihood of the model output given the data are normally distributed. The best fitting value of q is determined using a Metropolis-Hastings (MH) algorithm, using the following inputs;
Priors for q are assumed to be log-normally distributed with mean of 0.18 and standard deviation 0.5, resulting in wide 95% CI intervals of 0.064-0.490.
Proposed values for q (q’) are sampled from the current accepted value and a standard deviation of 0.003 (providing an acceptance rate varying from 0.5-0.7, depending upon simulations)
Iterations of the MH algorithm were run for at least 1,000 iterations, where convergence was assessed by there being no meaningful reduction in the log-likelihood with additional iterations run.
New values for q were accepted using probability
For each of the models that represented different assumptions about norovirus natural history, a comparison of the log-likelihood was performed. Formally, the model with the smallest log-likelihood would be the model to take forward to further simulations, and models with a difference in log-likelihood above 3.84 (p<0.05) should be excluded. However, several models (A0, B0 and A20, B20) have no biologically meaningful difference in reported incidence, and the additional outputs (seroprevalence, asymptomatic prevalence and R0) are not meaningfully different either, so we make the pragmatic decision to take all models forward to further simulations as we feel that the variability in projected incidence is an important consideration.
5. Comparison of SGSS data and model predictions for 2020
6. Table of estimated annual age-specific norovirus incidence (2019-2022)
These simulations assume that contact patterns will return to pre-COVID levels after July 2021.
7. Investigating the impact of reduced contact patterns from July 2021 compared to pre-pandemic
Acknowledgements
The authors are grateful for discussions within the wider Noropatrol consortium. The authors also acknowledge and are grateful for the data obtainable from public health diagnostic laboratories in England that submit laboratory reports of norovirus into SGSS.
Abbreviations
- IID
- Infectious Intestinal Disease
- NPI
- Non pharmaceutical interventions
- SGSS
- Second Generation Surveillance System