Next Article in Journal
SARS-CoV-2 Post Vaccinated Adverse Effects and Efficacy in the Egyptian Population
Previous Article in Journal
Protective Immunization of Atlantic Salmon (Salmo salar L.) against Salmon Lice (Lepeophtheirus salmonis) Infestation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network Study

1
Mathematics and Statistics, Centre for Disease Modelling, York University, Toronto, ON M3J 1P3, Canada
2
Department of Mathematics, Virginia Tech, Blacksburg, VA 24061, USA
3
Department of Mathematics, Purdue University, West Lafayette, IN 46202, USA
4
National Science Foundation, Alexandria, VA 22314, USA
5
Department of Mathematics, California State University, Northridge, CA 91330, USA
6
Department of Mathematics, University of Szeged, 6720 Szeged, Hungary
7
Epidemiology and Biostatistics, School of Population and Global Health, McGill University, Montreal, QC H3A 0G4, Canada
8
National Microbiology Laboratory, Public Health Risk Sciences Division, Public Health Agency of Canada, St. Hyacinthe, QC J2S 2M2, Canada
*
Author to whom correspondence should be addressed.
Vaccines 2022, 10(1), 17; https://doi.org/10.3390/vaccines10010017
Submission received: 25 September 2021 / Revised: 23 October 2021 / Accepted: 24 November 2021 / Published: 23 December 2021

Abstract

:
COVID-19 seroprevalence changes over time, with infection, vaccination, and waning immunity. Seroprevalence estimates are needed to determine when increased COVID-19 vaccination coverage is needed, and when booster doses should be considered, to reduce the spread and disease severity of COVID-19 infection. We use an age-structured model including infection, vaccination and waning immunity to estimate the distribution of immunity to COVID-19 in the Canadian population. This is the first mathematical model to do so. We estimate that 60–80% of the Canadian population has some immunity to COVID-19 by late Summer 2021, depending on specific characteristics of the vaccine and the waning rate of immunity. Models results indicate that increased vaccination uptake in age groups 12–29, and booster doses in age group 50+ are needed to reduce the severity COVID-19 Fall 2021 resurgence.

1. Introduction

The COVID-19 pandemic continues to affect the lives of Canadians. Despite increasing vaccination uptake, of first and second doses, and great decreases in COVID-19 cases of late, questions remain as to the future of the COVID-19 pandemic in this country, and any future need for COVID-19 vaccines to tackle COVID-19 resurgence.
A required step in understanding possibilities of resurgence or future vaccine needs lies in the determination and quantification of immunity in the Canadian population. Seroprevalence studies in different population cohorts have been conducted (see [1,2] for examples, and [3] for more information, many others are underway—see [4] for details) which can inform immunity distribution calculations. A recent statistical study by the COVID-19 Immunity Task Force (CITF), which incorporates different population seroprevalence studies into their analysis, estimated that the Canadian population has some immunity against COVID-19: due to infection, 5.4% (95% CrI: 0.6 to 15.8) as of 31 May 2021; due to infection and vaccination, 44.9% (95% CrI: 44.2 to 45.8) as of same date [5] (Credible interval (CrI)). The CITF [5] is the only working group that we know of that has attempted to quantify seroprevalence in the Canadian population.
It is of interest to provide other estimates of COVID-19 immunity in the Canadian population. Mathematical models of COVID-19 infection and vaccination can be used to estimate immunity distributions. In a previous study, we developed a mathematical model of COVID-19 infection and vaccination [6]. The model tracks infection and immunity status by age. An important difference of our model over other models of COVID-19 is that it incorporates differential outcomes of immunity, determined by infection severity, which is ultimately related to the prevalence of comorbidities in the Canadian population by age. A second important difference is that our model includes the effects of waning immunity, whereby immunity protection gained from infection or vaccination can decay over time. The model therefore has the capacity to provide time varying estimates of COVID-19 seroprevalence, which in turn can be used to inform public health decision-makers on vaccination policy i.e., targeting primary coverage to specific age groups, or booster doses.
In the current study we use our mathematical model to determine distributions of immunity in the Canadian population, by age, from infection, and from vaccination. The model is fit to daily COVID-19 incidence data up to 27 June 2021 [7], and incorporates actual (up to 27 June 2021) and projected (to September 2021) coverage of the first and second doses of COVID-19 vaccines [7,8]. We then use the model to quantify distributions of immunity from January 2020 to March 2022, given different assumed characteristics of the vaccines against various variants of concern (i.e., protection from infection, protection from severe disease), and different rates of waning immunity. In summary, we find that 60–80% of the Canadian population has some immunity to COVID-19 by late Summer 2021, depending on specific characteristics of the vaccine and the waning rate of immunity. Model results also indicate that increased vaccination uptake in age groups 12–29, and booster doses in age group 50+ are needed to reduce the severity COVID-19 Fall 2021 resurgence.

2. Methods

We implemented a model of COVID-19 infection with age structure (i.e., groups 0–4, 5–9, ⋯, 75+ years). A flow diagram of the model is shown in Figure 1 for one age group. The model is based on a Susceptible-Exposed-Infected-Vaccinated-Susceptible model structure (SEIVS). We use S i , E j k , I j , and V i to denote the number of susceptible, exposed, infected and vaccinated individuals in each age group, where i ( 1 i 4 ) denotes immune status, j ( 2 j 4 ) denotes symptom severity, k ( 1 k 3 ) represents stages in the exposed class (to obtain Gamma-distributed exposed sojourns), and = 1 , 2 denotes the number of doses of vaccine that individuals have received.
A detailed description of the mathematical model can be found in [6] and in the Appendix. Briefly, we assume that mild ( I 2 ), moderate ( I 3 ), or severe ( I 4 ) disease can be experienced upon infection, and that the probability of mild, moderate, or severe disease is determined by the comorbidity status in each age group [9]. We assume that all I 4 infecteds will be reported. We also assume that some fraction of I 3 will get tested and will be reported. Finally, we account for a small number of reported cases in I 3 and I 2 that will be tested because of contact tracing.
We assume that immunity gained after infection correlates with the severity of infection, such that higher levels of immunity are gained in individuals that have experienced more severe disease [10,11,12]. Immunity gained from vaccination is also implemented in the model, using actual Canadian vaccination roll-out data and projections [7,8]. Three different types of vaccine characteristics are considered which reflect protective capacities against infection and/or disease (given different variants of concern (VOC)) of the vaccines used in the Canadian population [13,14,15,16,17,18,19]. Characteristics of the vaccines are listed in Table 1.
In addition to modelling the gain of immunity, we also consider immunity decay in the population. Different waning rates are considered such that immunity wanes on average 1 year or 3 years between S and V classes, and consecutive S classes, which gives waning from full immunity to full susceptibility over 3 years or 9 years, respectively. For comparison, we also consider the case when immunity does not wane.
Finally, public health mitigation is incorporated into the model using modified contact matrices for home, school, work, and other types of contacts between age groups. Figure A2 plots the mitigation windows (gray bars) and the percent reduction in contacts gained from the contact matrix mitigation modifications in each window (x’s). The model also collectively accounts for compliance to social distancing and mask wearing, changes in testing and contact tracing rates, and changes in transmission due to VOCs and weather. This is done using a parameter κ that is determined from a model fit to daily incidence data.
Using our mathematical model, we track the distribution of immunity in the Canadian population over time, given different characteristics of protective capacity from the vaccine, and different assumptions with respect to the waning rate of immunity. There are nine scenarios that we consider altogether. We fit the model for each scenario from 25 January 2020 to 27 June 2021.
Projecting forward from the model fit, we modify the contact matrices corresponding to phase 1 in June to represent easing of lockdown restrictions across Canada with the schools remaining closed. In July, we increase contacts in public spaces, to reflect different reopening steps taken in Canadian jurisdictions. In September the contact restrictions are eased to phase 2, representing school reopening and easing of restrictions at work. For a detailed description of the mitigation phases see Table A3 and Table A4.
Behaviour is likely to relax over the summer and into the fall, increasing transmission. In July 2021 many jurisdictions relaxed social distancing rules. To reflect this, and the prevalence of the delta VOC, we increase κ by a factor of 2. Reductions in transmissibility due to weather should only affect the summer months. We suggest, therefore, that κ can be increased in September 2021. We compare model results with no change in κ in September 2021, a 20% increase to reflect increased transmissibility due to changes in the weather [20], and a 38% increase (to allow for consideration of introduction of a new VOC). Figure A2 provides an example of the contact rate change assuming an increase by a factor of 2 in the summer months following by a 20% increase to account for changes in weather.

3. Results

3.1. Model Fit

Figure 2 shows the daily incidence of severe, moderate + severe, and mild + moderate + severe infections for the model fitting to 27 June for all nine scenarios. Daily incidence data is shown in red and blue, with the last day of the fitting denoted by the red vertical line. The data indicated by the blue line shows appropriate model trend. We note that the model fit is similar between each of these subplots. We also note that one-year waning seems to match the trend more closely than the other waning rates.
The fitted values of κ are shown in Figure A2 (bottom panel) considering each vaccination and waning scenario. A description of the model parameters and the fitting algorithm is included in the Appendices. In Figure A2 (top panel), we plot the percent reduction in contacts from contact matrix modifications and κ (red +’s) for Vaccine 1 with a waning rate ω = 1 /year (see Table 1 and dashed blue line in Figure A2 bottom panel).
We note that Figures S1 and S4 are analogous, showing the same fitting results as presented in Figure 2. Large differences between these figures lie only from September 2021, when school opens, and variable increases in κ are implemented. From September 2021 we consider three different scenarios for κ : no predicted change from the last estimated value (Figure S1), a 20 % increase in κ (Figure 2), and a 38 % increase (Figure S4).

3.2. Seroprevalence

Given the model fitting above, we can now determine estimates of population seroprevalence. Figure 3, Figures S2 and S5 provide measurements of seroprevalence in the Canadian population given all vaccine types (Table 1) and changes in κ from September 2021 (assuming κ values increase by 20%, no change, and 38%, respectively). These figures show all classes that confer some immunity to SARS-CoV-2, over all age groups. The colours denote the sum of the one-dose vaccinated sub-populations ( i V i 1 ), the two-dose vaccinated sub-populations ( i V i 2 ), and the susceptible classes with partial and full immunity ( S 2 + S 3 + S 4 ), over 10-year age groups.The shading is related to the age classes, with lighter to darker coinciding with younger to older ages. In all nine scenarios, 60 to 80% of the Canadian population has some immunity against the pathogen (derived from infection and/or vaccination) by the time schools re-open in September 2021.
While the seroprevalence predicted by our model is higher than that projected by [5], there is little difference in the vaccine derived immunity between our results. The difference thus lies in the estimates of immunity derived from infections. While our model includes age structure, we must note that the system of ordinary differential equations assumes that the population is mixing at a higher level than in reality. It is therefore expected that the model will estimate higher levels of immunity. We note, however, that as vaccination coverage increases and becomes dominant in the population, the difference in seroprevalence estimates from our model and [5] should reduce. We also note that [5] does not include seroreversion (related to waning immunity). Inclusion of seroreversion will boost seroprevalence estimates in their work, so the difference between our results and theirs would reduce further.
For interest, in Figure 4, Figures S3 and S6 (assuming increases in κ values of 20%, no change, and 38%, respectively), we plot the corresponding distribution of protection against the virus of the different types of immunity in the population, over 10-year age groups, including no protection ( S 1 ), some protection ( S 2 + V 1 1 ), a higher level of protection ( S 3 + V 2 1 ), and full protection ( S 4 + i = 3 , 4 V i 1 + i = 1 , , 3 V i 2 ). The white area denotes the fraction of the population that resides in the exposed and infected classes E j k , I j , with j ( 2 j 4 ) and k ( 1 k 3 ).
Figure 3, Figure 4, Figures S2, S3, S5 and S6 show model predictions of the dominant level of immunity for each age class for each of scenario. We note that these outcomes are related to vaccine availability for each age group, model assumptions related to infection induced immunity, and the assumed vaccine acquired immunity characteristics (Table 1 which affect the circulation of the virus in the population. It is obvious in all of these figures that the younger ages, 0–10 year of age, have large levels of susceptibility (light red bars). This is due to the fact that vaccines are not yet available for these ages. Additionally, the mild infections that predominantly occur in these age groups confer low levels of immunity (light blue bars) that can wane quickly back to full susceptibility. The figures also show large levels of susceptibility in age groups 10–29 despite the availability of COVID-19 to 12–29 year olds. Finally, these figures show the effects of waning immunity. Immunity decays from higher immune classes to lower immune classes over time. Given that the early stages of the COVID-19 vaccine rollout in Canada centred on the older age groups, we observe increases in full susceptibility in the older age groups as time since vaccination increases (see dark red bars). Overall, the results in these figures suggest that increased vaccination coverage of age groups 12–29 should be pursued in government vaccination campaigns. Additionally, these results point to the need for a booster dose of vaccine in older age groups.

3.3. Herd Immunity

Herd immunity refers to a level of immunity, p c , that is needed in a population in order for that population to be resistant to further infection (a small number of infections may occur, but the disease will die out). Briefly, in a population with long-term effective immunity, the herd immunity threshold can be approximated by p c = 1 1 / R 0 , where p c is the critical fraction of the population that needs to have neutralizing immunity in order for the entire population to be protected, R 0 is the basic reproduction number of the pathogen in circulation, and homogeneous mixing in the population is assumed. Given a vaccine with vaccine efficacy 0 < v e < 1 , this approximation can be modified to be p c = ( 1 1 / R 0 ) / v e . We note that with the assumption of homogeneous mixing, this approximation will not provide a best estimate of the herd immunity threshold in age-based mixing models like ours (using contact matrices between age groups, and assuming preferential and proportional mixing), specifically considering different age-specific rates of vaccination [21,22,23,24,25,26,27,28,29]. Additionally, we note that the approximation assumes long-term or life-long immunity to be gained after infection or vaccination. If immunity wanes over time, and if the pathogen evolves, which decays effective immunity against infection over time, this will underestimate the immune population needed to provide herd immunity [29,30]. Nevertheless, we now consider the approximation with homogeneous mixing to estimate a herd immunity threshold for discussion. Currently, with the δ -variant in circulation, p c = 94 % , given a reproduction number of 6.5 (between 5 and 8 [31]) and a vaccine efficacy of 90% [13,14,15,16,17,18,19]. Our model predicts that 60 to 80% of the Canadian population will have some immunity to SARS-CoV-2 by the end of the vaccination campaign in late Summer 2021. We also observe that approximately 20 to 50% of the population will have neutralizing immunity, depending on the assumed waning rate and the characteristics of the vaccine (see Figure 4, Figures S3 and S6, green-shaded areas). This value of 20 to 50%, from our model (which incorporates age structure, waning immunity from infection and vaccination, and differential immunity gained after infection into our model structure) is far from the estimated herd immunity threshold of 94% (assuming homogeneous mixing and long-term immunity).
We note that benefits of vaccination programs are in their provision of protection from infection, but also in their protection from severe disease. Given 20–50% neutralizing immunity is not sufficient to protect the population from further infection, we now focus our attention to identifying populations that are at risk of severe infections in COVID-19 resurgence, given waning immunity from vaccination and infection, and protective capacities of immunity against infection and severe disease.

3.4. Resurgence

Figure 2, Figures S1 and S4 show no significant difference between model fits, and through the summer trend, but we do observe significant divergence past September, coinciding with the reopening of schools (Figure S1), and increases in κ by 20% and 38% (Figure 2 and Figure S4, respectively). When seroprevalence is high, corresponding to lower rates of waning immunity, we see that resurgence in infections, denoted by the white space, occurs in a reduced fashion. However, over all scenarios but one (when immunity does not wane, and vaccine protection is high, bottom right subplot), the model projects that the Fall 2021 resurgence will reach levels far greater than any wave of COVID-19 infection previously experienced.
We note that the resurgence in Fall 2021 is mainly driven by the virus circulation in Summer 2021. Figure 2 shows that our model captures the general trend in increased COVID-19 transmission in Summer 2021, and therefore, the Fall resurgence levels projected here are reasonable projections. Model results do however also show that if relaxation in July 2021 had not occurred, limiting the transmission of the delta VOC, resurgence would be greatly delayed and would have much lower infection levels in Summer and Fall 2021 (results not shown).
Individuals who received the vaccine early in the vaccination program will likely have experienced some effects of waning immunity by September 2021 (see Figure 3, Figure 4, Figures S2, S3, S5 and S6). Resurgence in infections will thus be stronger in these age groups. Figure 5 plots the proportion of daily I 4 infection incidence by age for each vaccine scenario, when school is reopened in September and κ is increased by 20%. Here, we see that resurgent cases are predominantly observed in age groups 50+ or 40+ when immunity wanes over 1 year or 3 years between consecutive immune classes, respectively (top and middle rows, respectively). Considering that COVID-19 infection fatality rates increase by age [5], we suggest that a vaccine booster campaign be considered in older age groups.
The degree of resurgence is affected by increases in κ and school re-opening. Given the higher levels of average daily contacts in younger age groups (see Figure A1, bottom right panel), and given that there is no vaccine yet approved for these ages, we recommend continued use of protective measures in schools, including mask wearing and social distancing. We also recommend that vaccination begin for these ages immediately after a vaccine is approved.

4. Discussion

Our model predicts that 60 to 80% of the Canadian population will have some immunity to SARS-CoV-2 by the end of the vaccination campaign in late Summer 2021. The population is vulnerable to virus resurgence, however, given the relaxation of non-pharmaceutical measures over Summer 2021 which allowed spread of the delta VOC. Model results pinpoint the need for increased vaccine coverage in ages 12–29, and booster doses in ages 50+.
The mathematical modelling study presented here is the first to address quantification of COVID-19 seroprevelance and distributions of immunity in a population. The model includes immunity gains from infection and vaccination, and also includes immunity decline due to waning. It is structured by 5-year age groups. As immunity is lost, model results can be used to identify age groups requiring vaccine booster campaigns, or increased vaccination coverage. It can also be used to determine the age groups most affected by COVID-19 resurgence.
The timing and severity of any resurgence is sensitive to the distribution of immunity. It is also sensitive to the introduction of VOCs, the cessation of personal protective measures and public health mitigation. This sensitivity is highlighted through the comparison of Figure 2 and Figure S4. We note that the possible effects of imported cases are not included in our model, and therefore cannot be gauged as to their effect on COVID-19 resurgence. This will require further evaluation, however, if the virus remains endemic in the population, importation is unlikely to substantially alter the model outputs according to recent modelling by members of our group [32].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/vaccines10010017/s1, Figure S1: Daily Incidence, with a no relaxation of κ . Showing daily incidence of I 4 (solid line), I 3 + I 4 (dashed line) and I 2 + I 3 + I 4 (dotted line). Smoothed incidence data is shown with a solid red line. The vertical blue, red and green lines indicate the start of vaccination, the beginning of the prediction phase, the start of relaxation, respectively. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. Figure S2: Seroprevalence, with no relaxation of κ . Showing seroprevalence as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The red region is the sum of a susceptible classes that have been exposed to the virus, either from natural infection or through waning from the vaccinated classes. The blue and green regions show the populations of the first and second dose vaccinated classes respectively. The total population with some immunity (the top of the red region) is equal to the vertical sum of the three Blue, Red, and Green regions. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase. Figure S3: Distribution of Immunity, with no relaxation of κ . Showing immune status of uninfected individuals as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The colours represent the immunity of the population from green, the fully immune population, through yellow and blue to red, the fully susceptible population. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase. Figure S4: Daily Incidence, with a 38 % relaxation of κ starting in September 2021. Showing daily incidence of I 4 (solid line), I 3 + I 4 (dashed line) and I 2 + I 3 + I 4 (dotted line). Smoothed incidence data is shown with a solid red line. The vertical blue, red and green lines indicate the start of vaccination, the beginning of the prediction phase, the start of relaxation, respectively. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. Figure S5: Seroprevalence, with a 38 % relaxation of κ starting in September 2021. Showing seroprevalence as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The red region is the sum of a susceptible classes that have been exposed to the virus, either from natural infection or through waning from the vaccinated classes. The blue and green regions show the populations of the first and second dose vaccinated classes respectively. The total population with some immunity (the top of the red region) is equal to the vertical sum of the three Blue, Red, and Green regions. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase. Figure S6: Distribution of Immunity, with a 38 % relaxation of κ starting in September 2021. Showing immune status of uninfected individuals as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The colours represent the immunity of the population from green, the fully immune population, through yellow and blue to red, the fully susceptible population. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase.

Author Contributions

Conceptualization, Z.F., J.L., G.R. and J.M.H.; Data curation, D.W.D. and J.M.H.; Formal analysis, D.W.D., L.C., Z.F. and J.M.H.; Funding acquisition, J.M.H.; Investigation, D.W.D., L.C., Z.F., J.L., G.R. and J.M.H.; Methodology, J.L., G.R. and J.M.H.; Project administration, J.M.H.; Supervision, J.M.H.; Validation, D.W.D., L.C., J.L., G.R. and J.M.H.; Visualization, D.W.D., L.C. and J.M.H.; Writing—original draft, D.W.D. and J.M.H.; Writing—review & editing, D.W.D., L.C., Z.F., J.L., G.R., D.L.B., N.H.O. and J.M.H. All authors have read and agreed to the published version of the manuscript.

Funding

Lauren Childs was supported by NSF Grant #2029262. J.M.H. is supported by NSERC and CIHR. D.W.D. is funded by the Canadian Immunization Research Network (CIRN). We gratefully acknowledge John Glasser, for fruitful discussions.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Canadian incidence data is taken from an online data portal. Please see [7] for details.

Conflicts of Interest

J.M.H. is a consultant to the health and vaccine industry. All contracts are not related to this study. All authors declare that they have no competing interest.

Appendix A

Appendix A.1. Model Equations

Our model tracks age, infection and immune status. Susceptible individuals of status i and age k are denoted by S i k ; similarly, infectious individuals by I i k . Infected but not-yet-infectious individuals of immune status i, age k and stage j are denoted by E i k j . Vaccinated individuals of initial immune status i, age k and dose k are denoted by V i k j . Parameter descriptions are found in Table A1. The system of ODEs for age group k is given by the following set of equations:
Susceptible compartments:
d d t S 1 k = j = 2 4 p 1 k j Λ 1 k S 1 k + ω 2 k S 2 k σ 1 k 1 ( t ) ρ S 1 k + ω 2 k V 1 k 1 , d d t S 2 k = j = 2 4 p 2 k j Λ 2 k S 2 k + ω 3 k S 3 k ω 2 k S 2 k σ 2 k 1 ( t ) ρ S 2 k + γ 2 k I 2 k + ω 3 k V 2 k 1 , d d t S 3 k = j = 3 4 p 3 k j Λ 3 k S 3 k + ω 4 k S 4 k ω 3 k S 3 k σ 3 k 1 ( t ) ρ S 3 k + γ 3 k I 3 k + ω 4 k j = 3 4 V j k 1 + j = 1 4 V j k 2 , d d t S 4 k = ω 4 k S 4 k σ 4 k 1 ( t ) ρ S 4 k + γ 4 k I 4 k ,
Vaccinated compartments:
d d t V 1 k 1 = σ 1 k 1 ( t ) ρ S 1 k σ 1 k 2 ( t ) ρ V 1 k 1 j = 2 4 p 2 k j Λ 2 k V 1 k 1 ω 2 k V 1 k 1 , d d t V 2 k 1 = σ 2 k 1 ( t ) ρ S 2 k σ 2 k 2 ( t ) ρ V 2 k 1 j = 3 4 p 3 k j Λ 3 k V 2 k 1 ω 3 k V 2 k 1 , d d t V 3 k 1 = σ 3 k 1 ( t ) ρ S 3 k σ 3 k 2 ( t ) ρ V 3 k 1 ω 4 k V 3 k 1 , d d t V 4 k 1 = σ 4 k 1 ( t ) ρ S 4 k σ 4 k 2 ( t ) ρ V 4 k 1 ω 4 k V 4 k 1 , d d t V 1 k 2 = σ 1 k 2 ( t ) ρ V 1 k 1 ω 4 k V 1 k 2 , d d t V 2 k 2 = σ 2 k 2 ( t ) ρ V 2 k 1 ω 4 k V 2 k 2 , d d t V 3 k 2 = σ 3 k 2 ( t ) ρ V 3 k 1 ω 4 k V 3 k 2 , d d t V 4 k 2 = σ 4 k 2 ( t ) ρ V 4 k 1 ω 4 k V 4 k 2 ,
Infected compartments:
d d t E 2 k 1 = p 1 k 2 Λ 1 k S 1 k + p 2 k 2 Λ 2 k S 2 k + p 2 k 2 Λ 2 k V 1 k 1 ζ 2 k 1 E 2 k 1 , d d t E 3 k 1 = p 1 k 3 Λ 1 k S 1 k + p 2 k 3 Λ 2 k S 2 k + p 2 k 3 Λ 2 k V 1 k 1 + p 3 k 3 Λ 3 k S 3 k + p 3 k 3 Λ 3 k V 2 k 1 ζ 3 k 1 E 3 k 1 , d d t E 4 k 1 = p 1 k 4 Λ 1 k S 1 k + p 2 k 4 Λ 2 k S 2 k + p 2 k 4 Λ 2 k V 1 k 1 + p 3 k 4 Λ 3 k S 3 k + p 3 k 4 Λ 3 k V 2 k 1 ζ 4 k 1 E 4 k 1 , d d t E 2 k 2 = ζ 2 k 1 E 2 k 1 ζ 2 k 2 E 2 k 2 , d d t E 3 k 2 = ζ 3 k 1 E 3 k 1 ζ 3 k 2 E 3 k 2 , d d t E 4 k 2 = ζ 4 k 1 E 4 k 1 ζ 4 k 2 E 4 k 2 , d d t E 2 k 3 = ζ 2 k 2 E 2 k 2 ζ 2 k 3 E 2 k 3 , d d t E 3 k 3 = ζ 3 k 2 E 3 k 2 ζ 3 k 3 E 3 k 3 , d d t E 4 k 3 = ζ 4 k 2 E 4 k 2 ζ 4 k 3 E 4 k 3 , d d t I 2 k = ζ 2 k 3 E 2 k 3 γ 2 k I 2 k δ 2 k I 2 k , d d t I 3 k = ζ 3 k 3 E 3 k 3 γ 3 k I 3 k δ 3 k I 3 k , d d t I 4 k = ζ 4 k 3 E 4 k 3 γ 4 k I 4 k δ 4 k I 41 ,
where for 1 k N ,
Λ i k ( t ) = α i k A k λ i k ( t ) , λ i k ( t ) = m = 1 N c k m j = 2 4 β j m I j m ( t ) j = 1 4 T j m ( t ) , T j m ( t ) = S j m ( t ) + V j m 1 ( t ) + V j m 2 ( t ) + E j m 1 ( t ) + E j m 2 ( t ) + E j m 3 ( t ) + I j m ( t ) .

Appendix A.2. Reproduction Number

The basic reproduction number is calculated using the Next Generation Method. It is calculated numerically, due to the dimensionality of the model system. Briefly, we determine the spectral radius of the matrix G = F V 1 , where matrix F is related to·new infections, and matrix V is related to transfers between classes infected classes. We refer the reader to [24,33,34,35] for important references on this method.
We assume a basic reproduction number of R 0 = 2.6 [36]. This is used to determine κ in the model fitting routine.

Appendix B. Parameters

Table A1. Parameter definitions for the ODE model.
Table A1. Parameter definitions for the ODE model.
ParameterDefinition
α i n susceptibility of individuals from S i n (i immunity status, n age group)
β j m infectivity of infected individuals from I j m (j immunity status, m age group)
γ j m recovery rate of infected individuals from I j m (j immunity status, m age group)
ζ j m k rates of progress through the pre-infectious period of infeciton (i immunity status, n age group, k stage)
δ j m disease-induced mortality rate of infected individuals from I j m (j immunity status, m age group)
ω i n waning rate of immunity of individuals from S i n (i immunity status, n age group)
p i n j proportion of S i going to I j upon infection with i = 1 , 2 , 3 and j = 2 , 3 , 4 ( i , j immunity status, n age group)
ρ i n vaccine efficacy (i immunity status, n age group)
σ k vaccination rate for first k = 1 and second k = 2 dose
A n per capita activity counts of individuals in age group n
c a n mixing matrix between individuals in age group a and age group n
Susceptibility: assume α 1 n = 1 , α 2 n = 2 3 , and α 3 n = 1 3 for any age group n; We assume that susceptibility decreases with increasing immunity status, but does not depend on age.
Infectivity: β 3 n = β = 0.08 , β 2 n = 0.5 β , β 4 n = 0.1 β for any age group n; Values are determined by fitting the model R 0 value to 2.6 [36]; We assume that infectivity varies by severity of infection and severity of disease—milder infections, with lower viral load, have lower infectivity, and more severe infections will have lower infectivity as such individuals will stay home, or will go to hospital.
Recovery: γ 1 n < γ 2 n < γ 3 n , with γ 1 n = 1 5 , γ 2 n = 1 10 , γ 3 n = 1 15 for any age group n; Recovery is assumed to depend on on disease severity, with milder disease associated with shorter infectious periods.
Exposed period: We assume an erlang distribution, erlang (3,2/3), with mean of 3 days.
Waning Immunity: ω i n = ω = 1 y r 365 for i = 1 , 2 , 3 and any age group n; We assume immunity gained from infection or vaccination lasts on average 1 to 3 years between successive immunity stages and is independent of age and immunity status. This assumption is varied in the current manuscript, and we also include model outcomes when it is assumed that immunity does not wane.
Probability of severe infection: see Table A2; We assume that severity of infection is associated with the prevalence of comorbidities in each age group. Precise values are determined from [9].
Table A2. Proportion of susceptible individuals from S i n going to I j infectious class.
Table A2. Proportion of susceptible individuals from S i n going to I j infectious class.
Age Group p i n ( 2 ) p i n ( 3 ) p i n ( 4 )
0–40.9791876250.0195323530.001280022
5–90.9703406740.027954180.001705146
10–140.9713547250.0269626030.001682673
15–190.9637904650.0338364240.002373111
20–240.947364080.0486183850.004017534
25–290.9237439930.0692897740.006966234
30–340.8972038020.091513710.011282488
35–390.8656053110.1167605190.01763417
40–440.8069848270.1625428010.030472373
45–490.7565879980.1978617180.045550284
50–540.6902451340.2415013710.068253495
55–590.6001904150.2963455920.103463993
60–640.5032450460.3468046340.14995032
65–690.4095054080.3839600710.206534521
70–740.3246640920.4040301630.271305745
75+0.2151502550.3737792070.411070537
Vaccination: the vaccination rate of the first dose is determined by the coverage detailed in the NACI document. The rate is determined each week, for particular age groups, following the NACI schedule. Vaccine efficacy for each dose is assumed to be 0.9. Individuals receive a second dose N d days following the first dose and acquire immunity after an additional 14 days. This is determined by vaccination guidelines in Canada [8].
Demographics: Given the short period of examination, we assume the absence of birth, natural mortality and aging. In the current study, we ignore disease induced death.
Contacts between age groups: Perturbations of the contact matrices are implemented to reflect public health mitigation strategies over time. A summary of the perturbations is given in Table A3 and Table A4. A graphical depiction of each public health mitigation phase listed in Table A3 and Table A4 is shown in Figure A1.
Table A3. Contact mitigation phase description (% reduction in contacts).
Table A3. Contact mitigation phase description (% reduction in contacts).
PhaseSchool ContactsOther ContactsWork Contacts
095%90% under 65, 95% over 6575% under 65, 95% over
195%75%70% under 65, 95% over
295%40% under 65, 65% over70% under 65, 95% over
315% under 20, 25% over40% under 65, 65% over35% under 65, 95% over
no mitigation0%0%0%
Table A4. Contact mitigation phase.
Table A4. Contact mitigation phase.
StartEndPhase
2020-02-052020-03-18no mitigation
2020-03-182020-03-313
2020-03-312020-04-251
2020-04-252020-06-190
2020-06-192020-09-011
2020-09-012021-01-083
2021-01-082021-01-221
2021-01-222021-02-190
2021-02-192021-04-091
2021-04-092021-05-140
2021-05-142021-07-151
2021-07-152021-09-012
2021-09-012021-12-313
2021-09-012021-12-313
Figure A1. Contact rates and activity under different phases of mitigation. Bar plots show contact rates between age groups under different mitigation strategies, with compliance (k) assumed to be 100%. Taller bars, by subplot, have lighter colors. All bar plots have identical y-axis scale for comparison (but are not identical in colour scale). The total activity is shown in the line plot (bottom right).
Figure A1. Contact rates and activity under different phases of mitigation. Bar plots show contact rates between age groups under different mitigation strategies, with compliance (k) assumed to be 100%. Taller bars, by subplot, have lighter colors. All bar plots have identical y-axis scale for comparison (but are not identical in colour scale). The total activity is shown in the line plot (bottom right).
Vaccines 10 00017 g0a1
Figure A2. Changes in contact levels over time. (top panel) Average overall contacts compared to no mitigation from modified contact matrices (x’s) plus parameter κ (red +’s). (bottom panel) Value of parameter κ from the model fit (before black vertical line), and into the future (after black vertical line). A 20% increase in κ in September 2021 is shown here. In both panels, mitigation phases are shown with background shading where darker shading represents more restrictive mitigation, see Table A3 and Table A4 for details.
Figure A2. Changes in contact levels over time. (top panel) Average overall contacts compared to no mitigation from modified contact matrices (x’s) plus parameter κ (red +’s). (bottom panel) Value of parameter κ from the model fit (before black vertical line), and into the future (after black vertical line). A 20% increase in κ in September 2021 is shown here. In both panels, mitigation phases are shown with background shading where darker shading represents more restrictive mitigation, see Table A3 and Table A4 for details.
Vaccines 10 00017 g0a2

Appendix C. Model Fitting

The model parameter κ was fit to daily incidence data over monthly windows using a least squares objective function. The objective function included a regularization hyper-parameter which penalized variation in κ to control over fitting. The objective function was minimized by comparison against the model solution for 3 5 I 3 + I 4 using the genetic algorithm as implemented in MATLAB. Convergence of the solution was confirmed by analyzing the the variance of the final population. The final fit for each scenario was fit for all time windows and independently of other scenarios. All additional parameters where informed by literature and chosen such that the basic reproduction number was 2.9 and are detailed above in Appendix B.

References

  1. Lewin, A.; Therrien, R.; De Serres, G.; Grégoire, Y.; Perreault, J.; Drouin, M.; Fournier, M.J.; Tremblay, T.; Beaudoin, J.; Beaudoin-Bussières, G.; et al. SARS-CoV-2 seroprevalence among blood donors in Québec, and analysis of symptoms associated with seropositivity: A nested case-control study. Can. J. Public Health 2021, 112, 576–586. [Google Scholar] [CrossRef]
  2. Saeed, S.; Drews, S.J.; Pambrun, C.; Yi, Q.L.; Osmond, L.; O’Brien, S.F. SARS-CoV-2 seroprevalence among blood donors after the first COVID-19 wave in Canada. Transfusion 2021, 61, 862–872. [Google Scholar] [CrossRef]
  3. Serotracker. Available online: https://serotracker.com/en/Explore (accessed on 12 July 2021).
  4. CITF. Task Force Funded Research; Technical Report. 2021. Available online: https://www.covid19immunitytaskforce.ca/task-force-research/ (accessed on 12 July 2021).
  5. CITF. Immunity Monitoring Report: Cumulative SARS-CoV-2 Seropositivity in Canada. Technical Report. 2021. Available online: https://www.covid19immunitytaskforce.ca (accessed on 12 July 2021).
  6. Childs, L.; Dick, D.W.; Feng, Z.; Heffernan, J.M.; Li, J.; Röst, G. Modeling waning and boosting of COVID-19 in Canada with vaccination. medRxiv 2021. [Google Scholar] [CrossRef]
  7. Berry, I.; Soucy, J.P.R.; Tuite, A.; Fisman, D. Open access epidemiologic data and an interactive dashboard to monitor the COVID-19 outbreak in Canada. Can. Med. Assoc. J. 2020, 192, E420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. National Advisory Council on Immunzation, Canada. High Consequence Infectious Diseases Working Group on COVID-19 Vaccines: Modelling Scenarios and Assumptions. Available online: https://www.medrxiv.org/content/10.1101/2021.05.18.21257426v1.full.pdf (accessed on 6 November 2020).
  9. Clark, A.; Jit, M.; Warren-Gash, C.; Guthrie, B.; Wang, H.H.X.; Mercer, S.W.; Sanderson, C.; McKee, M.; Troeger, C.; Ong, K.L.; et al. Global, regional, and national estimates of the population at increased risk of severe COVID-19 due to underlying health conditions in 2020: A modelling study. Lancet Glob. Health 2020, 8, E1003–E1017. [Google Scholar] [CrossRef]
  10. Robbiani, D.F.; Gaebler, C.; Muecksch, F.; Lorenzi, J.C.; Wang, Z.; Cho, A.; Agudelo, M.; Barnes, C.O.; Gazumyan, A.; Finkin, S.; et al. Convergent antibody responses to SARS-CoV-2 in convalescent individuals. Nature 2020, 584, 437–442. [Google Scholar] [CrossRef] [PubMed]
  11. Piccoli, L.; Park, Y.J.; Tortorici, M.A.; Czudnochowski, N.; Walls, A.C.; Beltramello, M.; Silacci-Fregni, C.; Pinto, D.; Rosen, L.E.; Bowen, J.E.; et al. Mapping neutralizing and immunodominant sites on the SARS-CoV-2 spike receptor-binding domain by structure-guided high-resolution serology. Cell 2020, 183, 1024–1042. [Google Scholar] [CrossRef]
  12. Dan, J.M.; Mateus, J.; Kato, Y.; Hastie, K.M.; Yu, E.D.; Faliti, C.E.; Grifoni, A.; Ramirez, S.I.; Haupt, S.; Frazier, A.; et al. Immunological memory to SARS-CoV-2 assessed for up to 8 months after infection. Science 2021, 371. [Google Scholar] [CrossRef]
  13. Abu-Raddad, L.J.; Chemaitelly, H.; Butt, A.A. Effectiveness of the BNT162b2 Covid-19 Vaccine against the B. 1.1. 7 and B. 1.351 Variants. N. Engl. J. Med. 2021, 385, 187–189. [Google Scholar] [CrossRef]
  14. Nasreen, S.; He, S.; Chung, H.; Brown, K.A.; Gubbay, J.B.; Buchan, S.A.; Wilson, S.E.; Sundaram, M.E.; Fell, D.B.; Chen, B.; et al. Effectiveness of COVID-19 vaccines against variants of concern, Canada. medRxiv 2021. [Google Scholar] [CrossRef]
  15. Emary, K.R.; Golubchik, T.; Aley, P.K.; Ariani, C.V.; Angus, B.; Bibi, S.; Blane, B.; Bonsall, D.; Cicconi, P.; Charlton, S.; et al. Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B. 1.1. 7): An exploratory analysis of a randomised controlled trial. Lancet 2021, 397, 1351–1362. [Google Scholar] [CrossRef]
  16. Liu, C.; Ginn, H.M.; Dejnirattisai, W.; Supasa, P.; Wang, B.; Tuekprakhon, A.; Nutalai, R.; Zhou, D.; Mentzer, A.J.; Zhao, Y.; et al. Reduced neutralization of SARS-CoV-2 B.1.617 by vaccine and convalescent serum. Cell 2021, 184, 4220–4236. [Google Scholar] [CrossRef]
  17. Public Health England. Vaccines Highly Effective against B.1.617.2 Variant after 2 Doses; Technical Report 2021. Available online: https://www.gov.uk/government/news/vaccines-highly-effective-against-b-1-617-2-variant-after-2-doses (accessed on 23 November 2021).
  18. Bernal, J.L.; Andrews, N.; Gower, C.; Gallagher, E.; Simmons, R.; Thelwall, S.; Stowe, J.; Tessier, E.; Groves, N.; Dabrera, G.; et al. Effectiveness of COVID-19 Vaccines against the B.1.617.2 Variant. N. Engl. J. Med. 2021, 385, 585–594. [Google Scholar]
  19. Sheikh, A.; McMenamin, J.; Taylor, B.; Robertson, C. SARS-CoV-2 Delta VOC in Scotland: Demographics, risk of hospital admission, and vaccine effectiveness. Lancet 2021, 397, 2461–2462. [Google Scholar] [CrossRef]
  20. Collin, A.; Hejblum, B.P.; Vignals, C.; Lehot, L.; Thiébaut, R.; Moireau, P.; Prague, M. Using Population Based Kalman Estimator to Model COVID-19 Epidemic in France: Estimating the Effects of Non-Pharmaceutical Interventions on the Dynamics of Epidemic. medRxiv 2021. [Google Scholar] [CrossRef]
  21. Knpil, D.H.; Rost, G. Modelling the strategies for age specific vaccination scheduling during influenza pandemic outbreaks. Math. Biosci. Eng. 2010, 8, 123–139. [Google Scholar]
  22. Glasser, J.; Feng, Z.; Moylan, A.; Del Valle, S.; Castillo-Chavez, C. Mixing in age-structured population models of infectious diseases. Math. Biosci. 2012, 235, 1–7. [Google Scholar] [CrossRef]
  23. Feng, Z.; Hill, A.N.; Smith, P.J.; Glasser, J.W. An elaboration of theory about preventing outbreaks in homogeneous populations to include heterogeneity or preferential mixing. J. Theor. Biol. 2015, 386, 177–187. [Google Scholar] [CrossRef] [Green Version]
  24. Feng, Z.; Han, Q.; Qiu, Z.; Hill, A.N.; Glasser, J.W. Computation of R in age-structured epidemiological models with maternal and temporary immunity. Discret. Contin. Dyn. Syst. Ser. B 2016, 21, 399. [Google Scholar] [CrossRef] [PubMed]
  25. Glasser, J.W.; Feng, Z.; Omer, S.B.; Smith, P.J.; Rodewald, L.E. The effect of heterogeneity in uptake of the measles, mumps, and rubella vaccine on the potential for outbreaks of measles: A modelling study. Lancet Infect. Dis. 2016, 16, 599–605. [Google Scholar] [CrossRef] [Green Version]
  26. Feng, Z.; Hill, A.N.; Curns, A.T.; Glasser, J.W. Evaluating targeted interventions via meta-population models with multi-level mixing. Math. Biosci. 2017, 287, 93–104. [Google Scholar] [CrossRef] [Green Version]
  27. Poghotanyan, G.; Feng, Z.; Glasser, J.W.; Hill, A.N. Constrained minimization problems for the reproduction number in meta-population models. J. Math. Biol. 2018, 77, 1795–1831. [Google Scholar] [CrossRef]
  28. Feng, Z.; Feng, Y.; Glasser, J.W. Influence of demographically-realistic mortality schedules on vaccination strategies in age-structured models. Theor. Popul. Biol. 2020, 132, 24–32. [Google Scholar] [CrossRef]
  29. Carlsson, R.M.; Childs, L.M.; Feng, Z.; Glasser, J.W.; Heffernan, J.M.; Li, J.; Röst, G. Modeling the waning and boosting of immunity from infection or vaccination. J. Theor. Biol. 2020, 497, 110265. [Google Scholar] [CrossRef]
  30. Heffernan, J.M.; Keeling, M.J. Implications of vaccination and waning immunity. Proc. R. Soc. B Biol. Sci. 2009, 276, 2071–2080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Chhatwal, J.; Mueller, P.C.P.; Ayer, T.; Adee, M.G.; Dalgic, O.; Ladd, M.A.; Linas, B. Changing Dynamics of COVID-19 in the US with the Emergence of the Delta Variant: Projections of the COVID-19 Simulator. medRxiv 2021. [Google Scholar] [CrossRef]
  32. Public Health Agency of Canada/National Microbiology Lab. COVID-19: PHAC Modelling Group Report. Public Health Agency of Canada/National Microbiology Lab. Available online: https://nccid.ca/wp-content/uploads/sites/2/2021/01/Modelling-Group-report-2020_01_14_Final.pdf (accessed on 25 February 2021).
  33. Heffernan, J.M.; Smith, R.J.; Wahl, L.M. Perspectives on the basic reproductive ratio. J. R. Soc. Interface 2005, 2, 281–293. [Google Scholar] [CrossRef]
  34. Van den Driessche, P.; Watmough, J. Further notes on the basic reproduction number. In Mathematical Epidemiology; Springer: Berlin/Heidelberg, Germany, 2008; pp. 159–178. [Google Scholar]
  35. Diekmann, O.; Heesterbeek, J.; Roberts, M.G. The construction of next-generation matrices for compartmental epidemic models. J. R. Soc. Interface 2010, 7, 873–885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Knight, J.; Mishra, S. Estimating effective reproduction number using generation time versus serial interval, with application to COVID-19 in the greater Toronto area, Canada. Infect. Dis. Model. 2020, 5, 889–896. [Google Scholar] [CrossRef]
Figure 1. Schematic of the model for one age group. Here, S 1 , S 2 , S 3 , and S 4 (purple shaded boxes) represent susceptible individuals who are immunologically naive, have some, moderate, and full immunity, respectively. I 2 , I 3 , and I 4 (red boxes) represent infected individuals with mild, moderate and severe symptoms, respectively, who will develop some, moderate, and full immunity once recovered (teal solid line), respectively. V i j ( i = 1 , 2 , 3 , 4 , j = 1 , 2 ) represent vaccinated individuals from the S i classes ( i = 1 , 2 , 3 , 4 ) after j = 1 , 2 doses of vaccine given a two-dose schedule. E i k ( i = 2 , 3 , 4 ; k = 1 , 2 , 3 ) represent exposed individuals (infected, asymptomatic, not infectious) with progressive stages k = 2 , 3 , 4 that will experience mild I 2 , moderate I 3 , and severe I 4 symptoms. Susceptible and vaccinated individuals can be infected and move to the exposed classes (red lines). Susceptible and vaccinated classes at the same location on the immunity continuum have similar characteristics. Immunity gained from infection and vaccination can wane (black lines). The bottom of the figure lists all the legends.
Figure 1. Schematic of the model for one age group. Here, S 1 , S 2 , S 3 , and S 4 (purple shaded boxes) represent susceptible individuals who are immunologically naive, have some, moderate, and full immunity, respectively. I 2 , I 3 , and I 4 (red boxes) represent infected individuals with mild, moderate and severe symptoms, respectively, who will develop some, moderate, and full immunity once recovered (teal solid line), respectively. V i j ( i = 1 , 2 , 3 , 4 , j = 1 , 2 ) represent vaccinated individuals from the S i classes ( i = 1 , 2 , 3 , 4 ) after j = 1 , 2 doses of vaccine given a two-dose schedule. E i k ( i = 2 , 3 , 4 ; k = 1 , 2 , 3 ) represent exposed individuals (infected, asymptomatic, not infectious) with progressive stages k = 2 , 3 , 4 that will experience mild I 2 , moderate I 3 , and severe I 4 symptoms. Susceptible and vaccinated individuals can be infected and move to the exposed classes (red lines). Susceptible and vaccinated classes at the same location on the immunity continuum have similar characteristics. Immunity gained from infection and vaccination can wane (black lines). The bottom of the figure lists all the legends.
Vaccines 10 00017 g001
Figure 2. Daily Incidence, with a 20 % relaxation of κ starting in September 2021. Showing daily incidence of I 4 (solid line), I 3 + I 4 (dashed line) and I 2 + I 3 + I 4 (dotted line). Smoothed incidence data is shown with a solid red line. The vertical blue, red and green lines indicate the start of vaccination, the beginning of the prediction phase, the start of relaxation, respectively. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively.
Figure 2. Daily Incidence, with a 20 % relaxation of κ starting in September 2021. Showing daily incidence of I 4 (solid line), I 3 + I 4 (dashed line) and I 2 + I 3 + I 4 (dotted line). Smoothed incidence data is shown with a solid red line. The vertical blue, red and green lines indicate the start of vaccination, the beginning of the prediction phase, the start of relaxation, respectively. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively.
Vaccines 10 00017 g002
Figure 3. Seroprevalence, with a 20 % relaxation of κ starting in September 2021. Showing seroprevalence as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The red region is the sum of a susceptible classes that have been exposed to the virus, either from natural infection or through waning from the vaccinated classes. The blue and green regions show the populations of the first and second dose vaccinated classes respectively. The total population with some immunity (the top of the red region) is equal to the vertical sum of the three Blue, Red, and Green regions. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase.
Figure 3. Seroprevalence, with a 20 % relaxation of κ starting in September 2021. Showing seroprevalence as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The red region is the sum of a susceptible classes that have been exposed to the virus, either from natural infection or through waning from the vaccinated classes. The blue and green regions show the populations of the first and second dose vaccinated classes respectively. The total population with some immunity (the top of the red region) is equal to the vertical sum of the three Blue, Red, and Green regions. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase.
Vaccines 10 00017 g003
Figure 4. Distribution of Immunity, with a 20 % relaxation of κ starting in September 2021. Showing immune status of uninfected individuals as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The colours represent the immunity of the population from green, the fully immune population, through yellow and blue to red, the fully susceptible population. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase.
Figure 4. Distribution of Immunity, with a 20 % relaxation of κ starting in September 2021. Showing immune status of uninfected individuals as a percent of the total population for 10 year age classes, with colour intensity corresponding to age class. The colours represent the immunity of the population from green, the fully immune population, through yellow and blue to red, the fully susceptible population. The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively. The vertical bar denotes the beginning of the prediction phase.
Vaccines 10 00017 g004
Figure 5. Age structure of daily incidence of severe disease, I 4 , over 7 different age groups (see legend and corresponding colours). The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively.
Figure 5. Age structure of daily incidence of severe disease, I 4 , over 7 different age groups (see legend and corresponding colours). The top row is waning of immunity by one year between consecutive classes; the middle row is waning of immunity of three years; the bottom row is no waning of immunity. Columns left to right represent vaccines 1 to 3, respectively.
Vaccines 10 00017 g005
Table 1. Vaccine efficacy. Three considered vaccine’s first and two-dose efficacy against infection and against severe disease.
Table 1. Vaccine efficacy. Three considered vaccine’s first and two-dose efficacy against infection and against severe disease.
Against InfectionAgainst Severe Disease
Two DosesFirst Dose
Vaccine 170%50%75%
Vaccine 280%70%80%
Vaccine 390%70%92%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dick, D.W.; Childs, L.; Feng, Z.; Li, J.; Röst, G.; Buckeridge, D.L.; Ogden, N.H.; Heffernan, J.M. COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network Study. Vaccines 2022, 10, 17. https://doi.org/10.3390/vaccines10010017

AMA Style

Dick DW, Childs L, Feng Z, Li J, Röst G, Buckeridge DL, Ogden NH, Heffernan JM. COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network Study. Vaccines. 2022; 10(1):17. https://doi.org/10.3390/vaccines10010017

Chicago/Turabian Style

Dick, David W., Lauren Childs, Zhilan Feng, Jing Li, Gergely Röst, David L. Buckeridge, Nick H. Ogden, and Jane M. Heffernan. 2022. "COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network Study" Vaccines 10, no. 1: 17. https://doi.org/10.3390/vaccines10010017

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop