Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Canadian Covid-19 pandemic public health mitigation measures at the province level.

The Covid-19 pandemic has prompted governments across the world to enforce a range of public health interventions. We introduce the Covid-19 Policy Response Canadian tracker (CPRCT) database that tracks and records implemented public health measures in every province and territory in Canada. The implementations are recorded on a four-level ordinal scale (0-3) for three domains, (Schools, Work, and Other), capturing differences in degree of response. The data-set allows the exploration of the effects of public health mitigation on the spread of Covid-19, as well as provides a near-real-time record in an accessible format that is useful for a diverse range of modeling and research questions.

Seasonality as a driver of pH1N12009 influenza vaccination campaign impact.

Although the most recent respiratory virus pandemic was triggered by a Coronavirus, sustained and elevated prevalence of highly pathogenic avian influenza viruses able to infect mammalian hosts highlight the continued threat of pandemics of influenza A virus (IAV) to global health. Retrospective analysis of pandemic outcomes, including comparative investigation of intervention efficacy in different regions, provide important contributions to the evidence base for future pandemic planning. The swine-origin IAV pandemic of 2009 exhibited regional variation in onset, infection dynamics and annual infection attack rates (IARs). For example, the UK experienced three severe peaks of infection over two influenza seasons, whilst Australia experienced a single severe wave. We adopt a seasonally forced 2-subtype model for the transmission of pH1N12009 and seasonal H3N2 to examine the role vaccination campaigns may play in explaining differences in pandemic trajectories in temperate regions. Our model differentiates between the nature of vaccine- and infection-acquired immunity. In particular, we assume that immunity triggered by infection elicits heterologous cross-protection against viral shedding in addition to long-lasting neutralising antibody, whereas vaccination induces imperfect reduction in susceptibility. We employ an Approximate Bayesian Computation (ABC) framework to calibrate the model using data for pH1N12009 seroprevalence, relative subtype dominance, and annual IARs for Australia and the UK. Heterologous cross-protection substantially suppressed the pandemic IAR over the posterior, with the strength of protection against onward transmission inversely correlated with the initial reproduction number. We show that IAV pandemic timing relative to the usual seasonal influenza cycle influenced the size of the initial waves of pH1N12009 in temperate regions and the impact of vaccination campaigns.

Pharmaceutical and non-pharmaceutical interventions for controlling the COVID-19 pandemic.

Disease spread can be affected by pharmaceutical interventions (such as vaccination) and non-pharmaceutical interventions (such as physical distancing, mask-wearing and contact tracing). Understanding the relationship between disease dynamics and human behaviour is a significant factor to controlling infections. In this work, we propose a compartmental epidemiological model for studying how the infection dynamics of COVID-19 evolves for people with different levels of social distancing, natural immunity and vaccine-induced immunity. Our model recreates the transmission dynamics of COVID-19 in Ontario up to December 2021. Our results indicate that people change their behaviour based on the disease dynamics and mitigation measures. Specifically, they adopt more protective behaviour when mandated social distancing measures are in effect, typically concurrent with a high number of infections. They reduce protective behaviour when vaccination coverage is high or when mandated contact reduction measures are relaxed, typically concurrent with a reduction of infections. We demonstrate that waning of infection and vaccine-induced immunity are important for reproducing disease transmission in autumn 2021.

Immunogenicity of COVID-19 vaccines and their effect on HIV reservoir in older people with HIV.

Older individuals and people with HIV (PWH) were prioritized for COVID-19 vaccination, yet comprehensive studies of the immunogenicity of these vaccines and their effects on HIV reservoirs are not available. Our study on 68 PWH and 23 HIV-negative participants aged 55 and older post-three vaccine doses showed equally strong anti-spike IgG responses in serum and saliva through week 48 from baseline, while PWH salivary IgA responses were low. PWH had diminished live-virus neutralization responses after two vaccine doses, which were 'rescued' post-booster. Spike-specific T cell immunity was enhanced in PWH with normal CD4+ T cell count, suggesting Th1 imprinting. The frequency of detectable HIV viremia increased post-vaccination, but vaccines did not affect the size of the HIV reservoir in most PWH, except those with low-level viremia. Thus, older PWH require three doses of COVID-19 vaccine for maximum protection, while individuals with unsuppressed viremia should be monitored for adverse reactions from HIV reservoirs.

Immunogenicity of COVID-19 vaccines and their effect on the HIV reservoir in older people with HIV.

Older individuals and people with HIV (PWH) were prioritized for COVID-19 vaccination, yet comprehensive studies of the immunogenicity of these vaccines and their effects on HIV reservoirs are not available. We followed 68 PWH aged 55 and older and 23 age-matched HIV-negative individuals for 48 weeks from the first vaccine dose, after the total of three doses. All PWH were on antiretroviral therapy (cART) and had different immune status, including immune responders (IR), immune non-responders (INR), and PWH with low-level viremia (LLV). We measured total and neutralizing Ab responses to SARS-CoV-2 spike and RBD in sera, total anti-spike Abs in saliva, frequency of anti-RBD/NTD B cells, changes in frequency of anti-spike, HIV gag/nef-specific T cells, and HIV reservoirs in peripheral CD4 + T cells. The resulting datasets were used to create a mathematical model for within-host immunization. Various regimens of BNT162b2, mRNA-1273, and ChAdOx1 vaccines elicited equally strong anti-spike IgG responses in PWH and HIV - participants in serum and saliva at all timepoints. These responses had similar kinetics in both cohorts and peaked at 4 weeks post-booster (third dose), while half-lives of plasma IgG also dramatically increased post-booster in both groups. Salivary spike IgA responses were low, especially in INRs. PWH had diminished live virus neutralizing titers after two vaccine doses which were 'rescued' after a booster. Anti-spike T cell immunity was enhanced in IRs even in comparison to HIV - participants, suggesting Th1 imprinting from HIV, while in INRs it was the lowest. Increased frequency of viral 'blips' in PWH were seen post-vaccination, but vaccines did not affect the size of the intact HIV reservoir in CD4 + T cells in most PWH, except in LLVs. Thus, older PWH require three doses of COVID-19 vaccine to maximize neutralizing responses against SARS-CoV-2, although vaccines may increase HIV reservoirs in PWH with persistent viremia.

Monkeypox: a review of epidemiological modelling studies and how modelling has led to mechanistic insight.

Human monkeypox (mpox) virus is a viral zoonosis that belongs to the Orthopoxvirus genus of the Poxviridae family, which presents with similar symptoms as those seen in human smallpox patients. Mpox is an increasing concern globally, with over 80,000 cases in non-endemic countries as of December 2022. In this review, we provide a brief history and ecology of mpox, its basic virology, and the key differences in mpox viral fitness traits before and after 2022. We summarize and critique current knowledge from epidemiological mathematical models, within-host models, and between-host transmission models using the One Health approach, where we distinguish between models that focus on immunity from vaccination, geography, climatic variables, as well as animal models. We report various epidemiological parameters, such as the reproduction number, R0, in a condensed format to facilitate comparison between studies. We focus on how mathematical modelling studies have led to novel mechanistic insight into mpox transmission and pathogenesis. As mpox is predicted to lead to further infection peaks in many historically non-endemic countries, mathematical modelling studies of mpox can provide rapid actionable insights into viral dynamics to guide public health measures and mitigation strategies.

Determination of significant immunological timescales from mRNA-LNP-based vaccines in humans.

A compartment model for an in-host liquid nanoparticle delivered mRNA vaccine is presented. Through non-dimensionalisation, five timescales are identified that dictate the lifetime of the vaccine in-host: decay of interferon gamma, antibody priming, autocatalytic growth, antibody peak and decay, and interleukin cessation. Through asymptotic analysis we are able to obtain semi-analytical solutions in each of the time regimes which allows us to predict maximal concentrations and better understand parameter dependence in the model. We compare our model to 22 data sets for the BNT162b2 and mRNA-1273 mRNA vaccines demonstrating good agreement. Using our analysis, we estimate the values for each of the five timescales in each data set and predict maximal concentrations of plasma B-cells, antibody, and interleukin. Through our comparison, we do not observe any discernible differences between vaccine candidates and sex. However, we do identify an age dependence, specifically that vaccine activation takes longer and that peak antibody occurs sooner in patients aged 55 and greater.

COVID-19 Vaccination and Healthcare Demand.

One of the driving concerns during any epidemic is the strain on the healthcare system. As we have seen many times over the globe with the COVID-19 pandemic, hospitals and ICUs can quickly become overwhelmed by cases. While strict periods of public health mitigation have certainly helped decrease incidence and thus healthcare demand, vaccination is the only clear long-term solution. In this paper, we develop a two-module model to forecast the effects of relaxation of non-pharmaceutical intervention and vaccine uptake on daily incidence, and the cascade effects on healthcare demand. The first module is a simple epidemiological model which incorporates non-pharmaceutical intervention, the relaxation of such measures and vaccination campaigns to predict caseloads into the Fall of 2021. This module is then fed into a healthcare module which can forecast the number of doctor visits, the number of occupied hospital beds, number of occupied ICU beds and any excess demand of these. From this module, we can also estimate the length of stay of individuals in ICU. For model verification and forecasting, we use the four most populous Canadian provinces as a case study.

Multiple cohort study of hospitalized SARS-CoV-2 in-host infection dynamics: Parameter estimates, identifiability, sensitivity and the eclipse phase profile.

Within-host SARS-CoV-2 modelling studies have been published throughout the COVID-19 pandemic. These studies contain highly variable numbers of individuals and capture varying timescales of pathogen dynamics; some studies capture the time of disease onset, the peak viral load and subsequent heterogeneity in clearance dynamics across individuals, while others capture late-time post-peak dynamics. In this study, we curate multiple previously published SARS-CoV-2 viral load data sets, fit these data with a consistent modelling approach, and estimate the variability of in-host parameters including the basic reproduction number, R0, as well as the best-fit eclipse phase profile. We find that fitted dynamics can be highly variable across data sets, and highly variable within data sets, particularly when key components of the dynamic trajectories (e.g. peak viral load) are not represented in the data. Further, we investigated the role of the eclipse phase time distribution in fitting SARS-CoV-2 viral load data. By varying the shape parameter of an Erlang distribution, we demonstrate that models with either no eclipse phase, or with an exponentially-distributed eclipse phase, offer significantly worse fits to these data, whereas models with less dispersion around the mean eclipse time (shape parameter two or more) offered the best fits to the available data across all data sets used in this work. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

A mathematical model of protein subunits COVID-19 vaccines.

We consider a general mathematical model for protein subunit vaccine with a focus on the MF59-adjuvanted spike glycoprotein-clamp vaccine for SARS-CoV-2, and use the model to study immunological outcomes in the humoral and cell-mediated arms of the immune response from vaccination. The mathematical model is fit to vaccine clinical trial data. We elucidate the role of Interferon-γ and Interleukin-4 in stimulating the immune response of the host. Model results, and results from a sensitivity analysis, show that a balance between the TH1 and TH2 arms of the immune response is struck, with the TH1 response being dominant. The model predicts that two-doses of the vaccine at 28 days apart will result in approximately 85% humoral immunity loss relative to peak immunity approximately 6 months post dose 1.

An HIV stochastic model with cell-to-cell infection, B-cell immune response and distributed delay.

In this study, a delayed HIV stochastic model with virus-to-cell infection, cell-to-cell transmission and B-cell immune response is proposed. We first transform the stochastic differential equation with distributed delay into a high-dimensional degenerate stochastic differential equation, and then theoretically analyze the dynamic behaviour of the degenerate model. The unique global solution of the model is given by rigorous analysis. By formulating suitable Lyapunov functions, the existence of the stationary Markov process is obtained if the stochastic B-cell-activated reproduction number is greater than one. We also use the law of large numbers theorem and the spectral radius analysis method to deduce that the virus can be cleared if the stochastic B-cell-inactivated reproduction number is less than one. Through uncertainty and sensitivity analysis, we obtain key parameters that determine the value of the stochastic B-cell-activated reproduction number. Numerically, we examine that low level noise can maintain the number of the virus and B-cell populations at a certain range, while high level noise is helpful for the elimination of the virus. Furthermore, the effect of the cell-to-cell infection on model behaviour, and the influence of the key parameters on the size of the stochastic B-cell-activated reproduction number are also investigated.

Long-term durability of immune responses to the BNT162b2 and mRNA-1273 vaccines based on dosage, age and sex.

The lipid nanoparticle (LNP)-formulated mRNA vaccines BNT162b2 and mRNA-1273 are a widely adopted multi vaccination public health strategy to manage the COVID-19 pandemic. Clinical trial data has described the immunogenicity of the vaccine, albeit within a limited study time frame. Here, we use a within-host mathematical model for LNP-formulated mRNA vaccines, informed by available clinical trial data from 2020 to September 2021, to project a longer term understanding of immunity as a function of vaccine type, dosage amount, age, and sex. We estimate that two standard doses of either mRNA-1273 or BNT162b2, with dosage times separated by the company-mandated intervals, results in individuals losing more than 99% humoral immunity relative to peak immunity by 8 months following the second dose. We predict that within an 8 month period following dose two (corresponding to the original CDC time-frame for administration of a third dose), there exists a period of time longer than 1 month where an individual has lost more than 99% humoral immunity relative to peak immunity, regardless of which vaccine was administered. We further find that age has a strong influence in maintaining humoral immunity; by 8 months following dose two we predict that individuals aged 18-55 have a four-fold humoral advantage compared to aged 56-70 and 70+ individuals. We find that sex has little effect on the immune response and long-term IgG counts. Finally, we find that humoral immunity generated from two low doses of mRNA-1273 decays at a substantially slower rate relative to peak immunity gained compared to two standard doses of either mRNA-1273 or BNT162b2. Our predictions highlight the importance of the recommended third booster dose in order to maintain elevated levels of antibodies.

Modeling waning and boosting of COVID-19 in Canada with vaccination.

SARS-CoV-2, the causative agent of COVID-19, has caused devastating health and economic impacts around the globe since its appearance in late 2019. The advent of effective vaccines leads to open questions on how best to vaccinate the population. To address such questions, we developed a model of COVID-19 infection by age that includes the waning and boosting of immunity against SARS-CoV-2 in the context of infection and vaccination. The model also accounts for changes to infectivity of the virus, such as public health mitigation protocols over time, increases in the transmissibility of variants of concern, changes in compliance to mask wearing and social distancing, and changes in testing rates. The model is employed to study public health mitigation and vaccination of the COVID-19 epidemic in Canada, including different vaccination programs (rollout by age), and delays between doses in a two-dose vaccine. We find that the decision to delay the second dose of vaccine is appropriate in the Canadian context. We also find that the benefits of a COVID-19 vaccination program in terms of reductions in infections is increased if vaccination of 15-19 year olds are included in the vaccine rollout.

Estimating the basic reproduction number at the beginning of an outbreak.

We compare several popular methods of estimating the basic reproduction number, R0, focusing on the early stages of an epidemic, and assuming weekly reports of new infecteds. We study the situation when data is generated by one of three standard epidemiological compartmental models: SIR, SEIR, and SEAIR; and examine the sensitivity of the estimators to the model structure. As some methods are developed assuming specific epidemiological models, our work adds a study of their performance in both a well-specified (data generating model and method model are the same) and miss-specified (data generating model and method model differ) settings. We also study R0 estimation using Canadian COVID-19 case report data. In this study we focus on examples of influenza and COVID-19, though the general approach is easily extendable to other scenarios. Our simulation study reveals that some estimation methods tend to work better than others, however, no singular best method was clearly detected. In the discussion, we provide recommendations for practitioners based on our results.

Longitudinal Assessment of SARS-CoV-2-Specific T Cell Cytokine-Producing Responses for 1 Year Reveals Persistence of Multicytokine Proliferative Responses, with Greater Immunity Associated with Disease Severity.

Cell-mediated immunity is critical for long-term protection against most viral infections, including coronaviruses. We studied 23 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected survivors over a 1-year post-symptom onset (PSO) interval by ex vivo cytokine enzyme-linked immunosorbent spot assay (ELISpot) assay. All subjects demonstrated SARS-CoV-2-specific gamma interferon (IFN-γ), interleukin 2 (IL-2), and granzyme B (GzmB) T cell responses at presentation, with greater frequencies in severe disease. Cytokines, mainly produced by CD4+ T cells, targeted all structural proteins (nucleocapsid, membrane, and spike) except envelope, with GzmB and IL-2 greater than IFN-γ. Mathematical modeling predicted that (i) cytokine responses peaked at 6 days for IFN-γ, 36 days for IL-2, and 7 days for GzmB, (ii) severe illness was associated with reduced IFN-γ and GzmB but increased IL-2 production rates, and (iii) males displayed greater production of IFN-γ, whereas females produced more GzmB. Ex vivo responses declined over time, with persistence of IL-2 in 86% and of IFN-γ and GzmB in 70% of subjects at a median of 336 days PSO. The average half-life of SARS-CoV-2-specific cytokine-producing cells was modeled to be 139 days (~4.6 months). Potent T cell proliferative responses persisted throughout observation, were CD4 dominant, and were capable of producing all 3 cytokines. Several immunodominant CD4 and CD8 epitopes identified in this study were shared by seasonal coronaviruses or SARS-CoV-1 in the nucleocapsid and membrane regions. Both SARS-CoV-2-specific CD4+ and CD8+ T cell clones were able to kill target cells, though CD8 tended to be more potent. IMPORTANCE Our findings highlight the relative importance of SARS-CoV-2-specific GzmB-producing T cell responses in SARS-CoV-2 control and shared CD4 and CD8 immunodominant epitopes in seasonal coronaviruses or SARS-CoV-1, and they indicate robust persistence of T cell memory at least 1 year after infection. Our findings should inform future strategies to induce T cell vaccines against SARS-CoV-2 and other coronaviruses.

A machine learning approach to differentiate between COVID-19 and influenza infection using synthetic infection and immune response data.

Data analysis is widely used to generate new insights into human disease mechanisms and provide better treatment methods. In this work, we used the mechanistic models of viral infection to generate synthetic data of influenza and COVID-19 patients. We then developed and validated a supervised machine learning model that can distinguish between the two infections. Influenza and COVID-19 are contagious respiratory illnesses that are caused by different pathogenic viruses but appeared with similar initial presentations. While having the same primary signs COVID-19 can produce more severe symptoms, illnesses, and higher mortality. The predictive model performance was externally evaluated by the ROC AUC metric (area under the receiver operating characteristic curve) on 100 virtual patients from each cohort and was able to achieve at least AUC = 91% using our multiclass classifier. The current investigation highlighted the ability of machine learning models to accurately identify two different diseases based on major components of viral infection and immune response. The model predicted a dominant role for viral load and productively infected cells through the feature selection process.

In Elimination Settings, Measles Antibodies Wane After Vaccination but Not After Infection: A Systematic Review and Meta-Analysis.

BACKGROUND: We conducted a systematic review to assess whether measles humoral immunity wanes in previously infected or vaccinated populations in measles elimination settings. METHODS: After screening 16 822 citations, we identified 9 articles from populations exposed to wild-type measles and 16 articles from vaccinated populations that met our inclusion criteria. RESULTS: Using linear regression, we found that geometric mean titers (GMTs) decreased significantly in individuals who received 2 doses of measles-containing vaccine (MCV) by 121.8 mIU/mL (95% confidence interval [CI], -212.4 to -31.1) per year since vaccination over 1 to 5 years, 53.7 mIU/mL (95% CI, -95.3 to -12.2) 5 to 10 years, 33.2 mIU/mL (95% CI, -62.6 to -3.9), 10 to 15 years, and 24.1 mIU/mL (95% CI, -51.5 to 3.3) 15 to 20 years since vaccination. Decreases in GMT over time were not significant after 1 dose of MCV or after infection. Decreases in the proportion of seropositive individuals over time were not significant after 1 or 2 doses of MCV or after infection. CONCLUSIONS: Measles antibody waning in vaccinated populations should be considered in planning for measles elimination.

Commentary on the use of the reproduction number R during the COVID-19 pandemic.

Since the beginning of the COVID-19 pandemic, the reproduction number [Formula: see text] has become a popular epidemiological metric used to communicate the state of the epidemic. At its most basic, [Formula: see text] is defined as the average number of secondary infections caused by one primary infected individual. [Formula: see text] seems convenient, because the epidemic is expanding if [Formula: see text] and contracting if [Formula: see text]. The magnitude of [Formula: see text] indicates by how much transmission needs to be reduced to control the epidemic. Using [Formula: see text] in a naïve way can cause new problems. The reasons for this are threefold: (1) There is not just one definition of [Formula: see text] but many, and the precise definition of [Formula: see text] affects both its estimated value and how it should be interpreted. (2) Even with a particular clearly defined [Formula: see text], there may be different statistical methods used to estimate its value, and the choice of method will affect the estimate. (3) The availability and type of data used to estimate [Formula: see text] vary, and it is not always clear what data should be included in the estimation. In this review, we discuss when [Formula: see text] is useful, when it may be of use but needs to be interpreted with care, and when it may be an inappropriate indicator of the progress of the epidemic. We also argue that careful definition of [Formula: see text], and the data and methods used to estimate it, can make [Formula: see text] a more useful metric for future management of the epidemic.