Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Antibody escape, the risk of serotype formation, and rapid immune waning: Modeling the implications of SARS-CoV-2 immune evasion

  • Catherine Albright,

    Roles Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Boston University, Boston, MA, United States of America

  • Debra Van Egeren,

    Roles Conceptualization, Writing – review & editing

    Affiliation Stanford University School of Medicine, Stanford, CA, United States of America

  • Aditya Thakur,

    Roles Methodology, Writing – review & editing

    Affiliation Boston University, Boston, MA, United States of America

  • Arijit Chakravarty,

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation Fractal Therapeutics, Cambridge, MA, United States of America

  • Laura F. White,

    Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

    Affiliation Boston University School of Public Health, Boston, MA, United States of America

  • Madison Stoddard

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    madison.stoddard@fractaltx.com

    Affiliation Fractal Therapeutics, Cambridge, MA, United States of America

Abstract

As the COVID-19 pandemic progresses, widespread community transmission of SARS-CoV-2 has ushered in a volatile era of viral immune evasion rather than the much-heralded stability of “endemicity” or “herd immunity.” At this point, an array of viral strains has rendered essentially all monoclonal antibody therapeutics obsolete and strongly undermined the impact of vaccinal immunity on SARS-CoV-2 transmission. In this work, we demonstrate that antibody escape resulting in evasion of pre-existing immunity is highly evolutionarily favored and likely to cause waves of short-term transmission. In the long-term, invading strains that induce weak cross-immunity against pre-existing strains may co-circulate with those pre-existing strains. This would result in the formation of serotypes that increase disease burden, complicate SARS-CoV-2 control, and raise the potential for increases in viral virulence. Less durable immunity does not drive positive selection as a trait, but such strains may transmit at high levels if they establish. Overall, our results draw attention to the importance of inter-strain cross-immunity as a driver of transmission trends and the importance of early immune evasion data to predict the trajectory of the pandemic.

Introduction

Some early commentators bullishly predicted the end of the COVID-19 pandemic [14], with the build-up of vaccine and natural immunity eventually curtailing SARS-CoV-2 transmission. However, the pandemic is now entering its fourth year despite a vast burden of prior infection, over 13 billion vaccination doses globally [5] and high prevalence of anti-SARS-CoV-2 antibodies [6,7]. Consistent with early-pandemic warnings [813], the pace of immune evasion has proven rapid [14,15], and transmission has continued robustly in the post-vaccine era [16]. Warnings about insufficient vaccine acceptance [17], rapid waning of vaccine and post-infection immune protection [18,19], and antibody evasion [810,20] have all materialized at this point [2123], leading to the high levels of SARS-CoV-2 transmission.

The ability of SARS-CoV-2 to evade immunity through mutations that degrade neutralizing antibody binding has been a major driver of ongoing viral transmission, as neutralizing antibodies have been demonstrated to be the correlate of immune protection for SARS-CoV-2 [2426]. Indeed, the post-omicron era of the pandemic has been marked by successive waves driven by immune-evading strains including BA.1, BA.5, XBB, and XBB.1.5 [14,22,2729]. These immune-evading strains acquire an evolutionary advantage in the context of widespread immunity through mutations that degrade the binding of neutralizing antibodies induced by infection with prior strains or by vaccines (antibody escape through “antigenic drift”) [30]. As neutralizing antibodies mediate sterilizing immunity to SARS-CoV-2–that is, they block infection upon exposure—evasion of these antibodies promotes reinfection [31]. This advantage has allowed these immune-evading strains to achieve dominance, drive spikes in transmission, and replace (succeed) pre-existing strains [20]. Between December 2021 (initial dissemination of the original BA.1 omicron) and December 2022, several strain succession events were documented resulting in an approximately 35-fold loss in neutralizing titer [32].

Understanding the potential of emerging strains to drive waves of infection, persist in circulation, and co-circulate with pre-existing strains is vital for understanding and reacting to this volatile phase of the pandemic. Immune evasion has implications for short-term and long-term transmission levels [33], possible changes in disease severity [34], and the efficacy of vaccines and therapeutics, especially monoclonal antibodies [22,35]. Anticipating the behavior of viral strains has tremendous practical significance for designing nonpharmaceutical and biomedical interventions.

In this study, we use an epidemiological modeling framework to build a quantitative understanding of the role of immune evasion in inter-strain competition and selection dynamics under endemic conditions. To this end, we developed a two-strain Susceptible- Infectious- Recovered- Susceptible (SIRS) model accounting for variable cross-immunity between a pre-existing and an invading strain. This paper explores viral evolutionary strategies by simulating a few relevant immunological scenarios: antigenic drift, which we surmise may result in symmetric or unilateral antibody escape, and induction of less durable immunity.

Antigenic drift results in reduced cross-immunity (immunity induced by one strain against another) compared to homologous immunity (immunity induced by a strain against itself) [36]. If the impact of antigenic drift is symmetric, the invading strain’s cross immunity against the original strain will equal the original strain’s cross immunity against the invading strain. The plausibility of this scenario is supported by the tolerance of SARS-CoV-2’s spike for a wide variety of mutations [8,37,38]. However, omicron BA.1 appeared to benefit from essentially unilateral antibody escape: while BA.1 strongly evaded pre-existing immunity to delta, delta was impeded by immunity induced by BA.1 [39]. The final scenario regards the possibility of viral strains with reduced durability of immunological protection from reinfection. Possibly exemplifying this scenario is omicron, which appears to exert weaker protection against homologous reinfection than delta (prior omicron reduces risk of omicron reinfection by 59.3%; prior delta infection reduces risk of delta reinfection by 92.3% [40].)

Determining the immunological properties likely to be selected for is crucial for predicting the trajectory of SARS-CoV-2 under widespread transmission. Although the rapid pace of SARS-CoV-2 evolution and the simultaneous emergence of immune-evading multiple strains paints a complex picture [41,42], this simplified analysis provides a basis for understanding the inter-strain competition and selection that underpin these dynamics. Identifying the characteristics of strains likely to be successful and drive significant waves of transmission is crucial to support early-warning systems. Co-circulation of viral serotypes–that is, viral strains sufficiently antigenically distinct to coexist [43]–is an emergent threat that requires greater understanding and may lead to dramatically worse outcomes in the long-term trajectory of the pandemic.

Methods

To evaluate short-term and long-term transmission of novel strains of SARS-CoV-2, we built a two-strain SIRS (susceptible-infectious-recovered-susceptible) model represented by a series of ordinary differential equations (ODEs). We implemented the model in Python using Scipy’s solve_ivp ODE solver [44]. Relevant code has been uploaded to a Github repository.

The model has two sets of parallel infected and recovered compartments, representing those infected with and recovered from the original strain and the invading strain. While an individual may possess immunity to both strains simultaneously, we excluded the possibility of co-infection. The model considers disease transmission, waning immunity, and induction of immunity due to exposure to either strain or both strains. We assumed a constant population size, such that birth rate is equal to death rate. We did not account for disease-related deaths. Fig 1 is a visual representation of the model compartments and transitions between compartments through infection, recovery, and immune waning.

thumbnail
Fig 1. Schematic representing the two-strain SIRS model.

Light blue represents fully susceptible individuals, dark blue represents infected individuals, and gray represents recovered individuals.

https://doi.org/10.1371/journal.pone.0292099.g001

The model parameters that are fixed among all model iterations are listed in Table 1A; parameters varied during the analysis and their ranges are described in Table 1B. The US population size and birth rate are used to inform the model. However, birth rate is expected to have minimal impact on immunological dynamics for frequent reinfection pathogen such as SARS-CoV-2, and the SIRS model scales linearly with population size [45]. As a result, the per-capita infection rates predicted by our model are expected to extend well to other populations.

thumbnail
Table 1. A: Model parameters fixed between scenarios.

B: Variable parameters and their values by scenario.

https://doi.org/10.1371/journal.pone.0292099.t001

In Table 1B, variable parameter ranges were chosen to reflect the range of plausible values. For example, commonly cited estimates for the intrinsic reproductive number (R0) of measles–which is highly contagious–reach up to 18 [52]. As a result, we explored a range of R0 values from 0 to 20 to address a wide range of possible scenarios. We varied cross-immunity values from none (0) to complete (1). The waning factor represents a fold-increase in the waning rate of immunity, which we varied from 1-fold to 10-fold to represent a wide range of possibilities.

The effective infectiousness for the original and invading strains, βO and βI, respectively, are derived from the relationship between R0,O and RO,I and the recovery rate, γ. βO and βI are calculated as follows:

The effective infectiousness for reinfection, βRO, and βRI can be derived from the relationship between the effective infectiousness and cross-immunity parameter, C:

To represent symmetric antibody escape, the cross-immunity of the original strain against the invading strain (COvI) is set equal to the cross-immunity of the invading strain against the original strain (CIvO). In the shorter duration of immunity scenario, there is assumed to be complete cross-immunity between the two strains; as a result, no individuals are simultaneously immune to both strains (no individuals in the ROI compartment).

The ODEs for the model are summarized below: Where S is the susceptible population, I is the infectious population, and R is the recovered population. Subscripts O and I represent the original strain and the invading strain, respectively. Individuals in the ROI compartment have immunity against both strains. IrO represents individuals who are infected with the original strain while immune to the invading strain; IrI represents individuals who are infected with the invading strain while immune to the original strain.

Initial conditions for the model are set such that the original strain is at an endemic steady-state, meaning the number of active infections is constant. The steady-state distribution of individuals in each compartment was determined by running the simulation for 10,000 days in the absence of the invading strain. Using this method, we found that the steady-state number of active original strain infections was 5.27 x 106, with 2.84 x 108 individuals recovered. To simulate invasion, we set the initial value of II−infections with the invading strain–to one.

We used the model to assess the invading strain’s ability to invade and drive a wave of transmission in the short-term as well as the impact of invading strains on long-term transmission trends. We define an outbreak as successful if the invading strain does not go extinct. We evaluated short-term transmission over a six-month period to reflect the period of enhanced transmission before immunity induced by the novel strain constrains its spread. To determine the long-term impact of strain invasion on steady-state SARS-CoV-2 transmission levels, we ran the two-strain model for 10 years, after which point immunological equilibrium is reached. Yearly infections are read out over the last year of the simulation.

Results

Mutual immune evasion favors successful invasion and co-circulation

In Figs 2 and 3, we explore epidemiological outcomes after the introduction of an invader strain with mutually reduced cross-immunity with respect to the original strain. We assumed immune evasion is perfectly symmetric: that is, the invader strain’s sensitivity to immunity induced by the original strain is equal to the original strain’s sensitivity to new immunity induced by the invader strain. Fig 2 demonstrates that if immune evasion is symmetric, invasion is successful under a wide range of transmissibility and immune evasion conditions. Even a poorly transmissible strain relative to currently circulating strains may successfully invade if cross-immunity is low enough. This suggests that symmetric immune evasion is highly evolutionarily favorable as a trait, and selection for symmetric immune evasion may overcome deficits in intrinsic reproductive number.

thumbnail
Fig 2. Invasion outcomes for invader strains with varying transmissibility (R0) and symmetric cross-immunity with the original strain.

Outcomes are classified based on whether either strain becomes extinct or the strains co-circulate. Invasion is successful unless the invader strain becomes extinct (yellow).

https://doi.org/10.1371/journal.pone.0292099.g002

thumbnail
Fig 3. Short-term and long-term transmission trends after introduction of invader strains with varying transmissibility (R0) and symmetric cross-immunity properties.

The total number of infections is expressed as a fraction of the population size (represented by the colorbar). Fractions greater than 1 indicate reinfections. Total infections over the first six months after invader strain introduction due to A) the original strain, B) the invader strain, or C) both summed together. Total infections over one year at steady-state (SS) for D) the original strain, E) the invader strain, or F) both summed together. Colormaps are scale-matched across subpanels A-C and D-F.

https://doi.org/10.1371/journal.pone.0292099.g003

Additionally, if immune evasion is strong and symmetric, the success of the invader strain may not come at the expense of the original strain. Succession (extinction of the original strain) is expected only in scenarios where cross-immunity between the strains is high, and the invading strain is more transmissible than the original strain.

Invasion increases transmission in the short-term and co-circulation increases transmission in the long-term

Fig 3 explores short-term and long-term strain transmission of the original and invader strains after various invasion scenarios. As shown in Fig 3A, invasion of a novel strain may only modestly reduce transmission of the original strain in the short-term, while extensive transmission of the invading strain may occur (Fig 3B). This results in a large increase in overall transmission in the short-term (a “wave” or a “spike”) for successful invasion scenarios compared to failed invasions, which reflect baseline transmission of the original strain (Fig 3C). Interestingly, short-term reductions in transmission of the original strain upon invasion (Fig 3A) do not always result in extinction of the original strain (Fig 2). In general, higher transmissibility of the invader strain and greater mutual immune evasion (lower cross-immunity) promote larger waves of transmission.

In the long-term, when steady-state is reached, low symmetric cross-immunity that results in co-circulation of strains causes the highest overall transmission (Fig 3D–3F.) This is because as cross-immunity levels decrease, competition between the strains also decreases. We also note that increases in transmissibility in the absence of immune evasion–represented by the uppermost section of each plot where cross-immunity approaches one–are unlikely to drive significant increases in transmission in the short- or long-term. However, overall transmission remains high in these scenarios, nearing 60% of the population infected annually (Fig 3F).

Unilateral immune evasion favors succession

Fig 4 explores invasion scenarios in which cross-immunity is unilateral: the invader strain evades the original strain’s immunity but exerts full cross-immunity against the original strain. While Fig 2 demonstrates that succession (extinction of the original strain) is rare if immune evasion is perfectly symmetric, Fig 4 indicates that unilateral immune evasion favors succession, especially if the R0 of the invader strain is equal to or greater than the original strain’s R0. However, the conditions under which invasion succeeds–represented by both co-circulation and original strain extinction outcomes–are identical to the symmetric immune evasion scenario.

thumbnail
Fig 4. Invasion outcomes for novel strains with varying transmissibility (R0) and degrees of unilateral immune evasion.

https://doi.org/10.1371/journal.pone.0292099.g004

In the short-term, transmission dynamics after emergence of an invading strain in the unilateral immune evasion scenario are similar to the symmetric scenario (Fig 5A–5C). This is intuitive: the difference between the two scenarios is in the nature of the invading strain’s immunity, which is not yet present in the short-term. In the long-term, immunity created by the invading strain constrains spread of the original strain, in many cases leading to its extinction (Fig 5D–5F). This reduces transmission of the original strain with limited benefit to the invading strain, except in cases where the invading strain’s dominance is weaker. As a result, long-term, overall transmission is not substantially impacted by successful invasion when cross-immunity induced by the invading strain against the original strain is strong.

thumbnail
Fig 5. Short-term and long-term transmission trends after introduction of invader strains with varying transmissibility (R0) and unilateral cross-immunity properties.

In these scenarios, the invader strain induces strong immunity against the original strain, while the original strain induces weaker immunity against the invader strain. The total number of infections is expressed as a fraction of the population size. Fractions greater than 1 indicate reinfections. Total infections over six months due to A) the original strain, B) the invader strain, or C) both summed together. Total infections over one year at steady-state for D) the original strain, E) the invader strain, or F) both summed together. Colormaps are scale-matched across subpanels A-C and D-F.

https://doi.org/10.1371/journal.pone.0292099.g005

Duration of invader strain’s immunity does not impact its potential to invade

We also simulated invasion attempts by strains that exert shorter-term immunity. We found that in the absence of neutralizing antibody evasion, the duration of immunity induced by the invading strain does not impact its ability to establish (Fig 6). This is consistent with the intuition that immunity induced by the invader strain does not yet exist at the time of invasion and likewise does not impact invasion success. Additionally, a strain that induces less durable immunity achieves no selective advantage in the absence of immune evasion because the shorter immunity increases the reinfection potential of both strains equally. Additionally, we note that variation in this property in the absence of immune evasion does not support coexistence of strains. However, strains with shorter duration of immunity may transmit at much higher levels in the long term (S1 Fig).

thumbnail
Fig 6. Outcomes for invasion scenarios in which the invader strain induces less durable immunity than the original strain.

The duration of immunity against wild-type SARS-CoV-2 infection is estimated to be 550 days [48].

https://doi.org/10.1371/journal.pone.0292099.g006

Discussion

The epidemiological modeling analyses in this work demonstrate that cross-immunity between circulating and emergent strains is expected to be a major driver of the ability of invading strains to successfully establish as well as short-term and long-term transmission trends. This includes the possibility of stable co-circulation of antigenically distinct strains (serotypes) [Figs 2 & 4]. Under conditions of widespread anti-SARS-CoV-2 immunity, neutralizing antibody evasion likely provides an evolutionary advantage. Increases in immune evasiveness can be expected to drive selection and transmission to a greater extent than further increases in intrinsic transmissibility (R0) alone [Figs 3 & 5].

In this work, we showed that the intrinsic transmissibility of the emerging strain and the immunological landscape a new strain encounters (by definition, imposed by pre-existing strains) shape its invasion potential and extent of transmission in the short-term. This is the cross-immunity of the original strain against the invader strain. Strains with antibody escape that weakens pre-existing immunity are likely to be strongly selected for and drive significant increases in short-term transmission [Figs 3A–3C and 5A–5C]. In this way, ongoing viral evolution is likely to contribute to long-term high and volatile levels of SARS-CoV-2 transmission. Notably, our work shows that even poorly transmissible (low R0) strains may successfully invade if they sufficiently evade pre-existing neutralizing antibodies.

Our analyses further showed that while evasion of pre-existing immunity dictates which strains will successfully emerge, the properties of immunity induced by successful invading strains contribute to their long-term impact on overall transmission levels. For example, we have shown that invader strains inducing less durable immunity can sustain greater long-term transmission if invasion is successful [S1 Fig].

Additionally, we have shown here that weak cross-immunity of the invader strain against the original strain can permit the invader and original strains to coexist with each other and continue to co-circulate [Figs 2 & 4]. This situation, referred to as serotype formation [53], may significantly increase long-term transmission levels. The existence of serotypes for other human pathogens (e.g. dengue) has made their control more complicated, as tests and vaccines have to be created to match each serotype [54]. Serotypes have already been observed for human coronaviruses OC43 and NL63 [55,56]. Co-circulation of SARS-CoV-2 serotypes would increase demand on already-strained global sequencing systems for tracking viral evolution and may exacerbate the challenge of matching vaccines and therapeutics to circulating strains. Serotypes are also thought to increase the risk of antibody-dependent enhancement (ADE), as has been observed for dengue [57]. Although ADE is not thought to be a risk for SARS-CoV-2 at present, the co-circulation of two or more SARS-CoV-2 serotypes would create conditions known to contribute to ADE in other viral diseases [58].

These findings highlight the importance of measuring cross-immunities between pre-existing and emergent strains to forecast the likelihood of invasion success and short-term and long-term implications of a successful invasion. Although real-world cross-immunity is challenging to measure, neutralizing potency of post-immunization sera has proven to be a practical correlate [59]. Many studies assess the neutralization potency induced by vaccines against emerging strains [60,61], while some studies have also assessed the neutralization potency raised by circulating strains against emerging strains [6264]. Using models linking neutralizing antibody potency to extent of protection, this data can be leveraged to estimate cross-immunity against invading strains [25]. This is important for predicting which strains are likely to establish successfully and significantly increase short-term transmission.

However, data is often lacking regarding the neutralization potency induced by emergent strains against pre-existing strains. Measuring the neutralizing potency of invading strain infectee sera against pre-existing strains could reveal the risk of long-term co-circulation. In this work, we have demonstrated that co-circulation may follow apparent declines in the transmission of the original strain during high post-emergence transmission of the invader strain. This underlines the importance of in vitro or clinical data regarding cross-immunity of invading strains against pre-existing strains. The likelihood of co-circulation may not be obvious in early epidemiological data, at which point interventions are most likely to succeed [65].

There are a few limitations to this analysis. Firstly, we did not simulate the impact of vaccines in this analysis. In our prior work [66], we have explored the impact of vaccines on SARS-CoV-2 evolutionary dynamics; we determined that the impact of vaccines on evolutionary dynamics is limited due to minimal vaccinal protection against infection. Secondly, our SIRS model implementation categorizes individuals into four immunological groups: fully susceptible, immune to original strain and susceptible to invader strain, immune to invader strain and susceptible to original strain, and immune to both strains. As we have explored in other analyses using agent-based models, the immunity landscape is more complex than this SIRS model can capture. This complexity arises due to heterogeneity in individual exposure to infection and vaccination, interindividual variability in neutralizing antibody durability, and neutralizing antibody build-up over successive infections and vaccinations [48,67]. In particular, we surmise that the build-up of neutralizing potency over successive infections may cause the apparent cross-immunities between strains to not be fixed over time. This may result in additional competitive dynamics in the long term. Lastly, our analysis addresses competition between two strains for hosts in short and long timeframes. In reality, multiple strains may be in competition simultaneously, and the emergence of further invader strains may occur before the “long-term” two-strain dynamics explored here are fully realized [68,69]. We also acknowledge the possibility of short-term “pseudo” co-circulation, in which strains may appear to co-circulate in the short-term before sufficient infections have occurred to impose immunological constraints. Nevertheless, our work provides a basis for understanding basic competition dynamics between any two strains.

With these limitations in mind, the purpose of this analysis is to broadly identify strain properties conducive to invasion and to highlight the risk and implications of co-circulation of SARS-CoV-2 strains. Our work demonstrates that the emergence of immune-evading strains is a key driver of transmission rates at this stage of the pandemic, as invasion events are likely to drive extensive transmission in the short-term. We emphasize the value of early data characterizing the antibody evasion of and immune responses to emergent strains to inform SARS-CoV-2 vaccination and other mitigation strategies.

Our work also suggests the emergence of SARS-CoV-2 serotypes as a further potential risk of the current public health strategy. Stable co-existence of strains has not yet been demonstrated at the time of writing, perhaps due to strong cross-immunity against original strains as demonstrated with Delta and Omicron BA.1 [39]. In sum, the results presented here point to the likelihood of high baseline SARS-CoV-2 transmission, with the possibility of transient increases in transmission (“waves”) driven by successive (or co-circulating) waves of immune-evasive strains. These predictions are reflected in the real-world emergence and widespread circulation of many competing immune-evasive SARS-CoV-2 strains [41,42]. We identify inter-strain cross-immunity as an important variable for predicting outcomes and gauging risk in this seemingly unpredictable scenario.

Supporting information

S1 Fig. Short-term and long-term transmission trends after introduction of invader strains with varying transmissibility (R0) and durability of immunity.

The total number of infections is expressed as a fraction of the population size. Fractions greater than 1 indicate reinfection. Total infections over six months due to A) the original strain, B) the invader strain, or C) both summed together. Total infections over one year at steady-state for D) the original strain, E) the invader strain, or F) both summed together. Colormaps are matched across subpanels A-C and D-F.

https://doi.org/10.1371/journal.pone.0292099.s001

(ZIP)

References

  1. 1. Powell A. Anthony Fauci offers a timeline for ending COVID-19 pandemic. In: Harvard Gazette [Internet]. 10 Dec 2020 [cited 7 Oct 2022]. Available: https://news.harvard.edu/gazette/story/2020/12/anthony-fauci-offers-a-timeline-for-ending-covid-19-pandemic/.
  2. 2. The covid-19 pandemic will be over by the end of 2021, says Bill Gates. The Economist. Available: https://www.economist.com/international/2020/08/18/the-covid-19-pandemic-will-be-over-by-the-end-of-2021-says-bill-gates. Accessed 22 Jan 2023.
  3. 3. Pence M. There Isn’t a Coronavirus ‘Second Wave.’ Wall Street Journal. 16 Jun 2020. Available: https://www.wsj.com/articles/there-isnt-a-coronavirus-second-wave-11592327890. Accessed 22 Jan 2023.
  4. 4. Marais BJ, Sorrell TC. Pathways to COVID-19 ‘community protection.’ Int J Infect Dis. 2020;96: 496–499. pmid:32425642
  5. 5. Mathieu E, Ritchie H, Rodés-Guirao L, Appel C, Giattino C, Hasell J, et al. Coronavirus Pandemic (COVID-19). Our World in Data. 2020 [cited 22 Jan 2023]. Available: https://ourworldindata.org/covid-vaccinations.
  6. 6. Coronavirus (COVID-19) latest insights: Antibodies. In: Office for National Statistics [Internet]. [cited 22 Jan 2023]. Available: https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/articles/coronaviruscovid19latestinsights/antibodies.
  7. 7. COVID-19 SeroHub. In: NIH [Internet]. 21 Nov 2022 [cited 22 Jan 2023]. Available: https://covid19serohub.nih.gov/.
  8. 8. Egeren DV, Novokhodko A, Stoddard M, Tran U, Zetter B, Rogers M, et al. Risk of rapid evolutionary escape from biomedical interventions targeting SARS-CoV-2 spike protein. PLOS ONE. 2021;16: e0250780. pmid:33909660
  9. 9. Greaney AJ, Starr TN, Gilchuk P, Zost SJ, Binshtein E, Loes AN, et al. Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain that Escape Antibody Recognition. Cell Host Microbe. 2021;29: 44–57.e9. pmid:33259788
  10. 10. Thomson EC, Rosen LE, Shepherd JG, Spreafico R, Filipe A da S, Wojcechowskyj JA, et al. Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity. Cell. 2021;184: 1171–1187.e20. pmid:33621484
  11. 11. Nabel KG, Clark SA, Shankar S, Pan J, Clark LE, Yang P, et al. Structural basis for continued antibody evasion by the SARS-CoV-2 receptor binding domain. Science. 2021;375: eabl6251. pmid:34855508
  12. 12. Stoddard M, Sarkar S, Yuan L, Nolan RP, White DE, White LF, et al. Beyond the new normal: Assessing the feasibility of vaccine-based suppression of SARS-CoV-2. PLoS One. 2021;16: e0254734. pmid:34270597
  13. 13. Egeren DV, Stoddard M, Novokhodko A, Rogers M, Joseph-McCarthy D, Zetter B, et al. The specter of Manaus: the risks of a rapid return to pre-pandemic conditions after COVID-19 vaccine rollout. medRxiv. 2021; 2021.05.14.21257250.
  14. 14. Yue C, Song W, Wang L, Jian F, Chen X, Gao F, et al. Enhanced transmissibility of XBB.1.5 is contributed by both strong ACE2 binding and antibody evasion. Immunology; 2023 Jan.
  15. 15. Aggarwal A, Akerman A, Milogiannakis V, Silva MR, Walker G, Stella AO, et al. SARS-CoV-2 Omicron BA.5: Evolving tropism and evasion of potent humoral responses and resistance to clinical immunotherapeutics relative to viral variants of concern. eBioMedicine. 2022;84. pmid:36130476
  16. 16. Coronavirus (COVID-19) latest insights: Infections. In: Office for National Statistics [Internet]. [cited 22 Jan 2023]. Available: https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/articles/coronaviruscovid19latestinsights/infections.
  17. 17. Stoddard M, Van Egeren D, Johnson KE, Rao S, Furgeson J, White DE, et al. Individually optimal choices can be collectively disastrous in COVID-19 disease control. BMC Public Health. 2021;21: 832. pmid:33931055
  18. 18. Bottino D, Hather G, Yuan L, Stoddard M, White L, Chakravarty A. Using mixed-effects modeling to estimate decay kinetics of response to SARS-CoV-2 infection. Antib Ther. 2021;4: 144–148. pmid:34286216
  19. 19. Ibarrondo FJ, Hofmann C, Fulcher JA, Goodman-Meza D, Mu W, Hausner MA, et al. Primary, Recall, and Decay Kinetics of SARS-CoV-2 Vaccine Antibody Responses. ACS Nano. 2021. pmid:34159781
  20. 20. Fischer W, Giorgi EE, Chakraborty S, Nguyen K, Bhattacharya T, Theiler J, et al. HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens. Cell Host Microbe. 2021;29: 1093–1110. pmid:34242582
  21. 21. Cuadros DF, Gutierrez JD, Moreno CM, Escobar S, Miller FD, Musuka G, et al. Impact of healthcare capacity disparities on the COVID-19 vaccination coverage in the United States: A cross-sectional study. The Lancet Regional Health–Americas. 2023;18. pmid:36536782
  22. 22. Wang Q, Iketani S, Li Z, Liu L, Guo Y, Huang Y, et al. Alarming antibody evasion properties of rising SARS-CoV-2 BQ and XBB subvariants. Cell. 2023;186: 279–286.e8. pmid:36580913
  23. 23. Goldberg Y, Mandel M, Bar-On YM, Bodenheimer O, Freedman LS, Ash N, et al. Protection and Waning of Natural and Hybrid Immunity to SARS-CoV-2. New England Journal of Medicine. 2022;386: 2201–2212. pmid:35613036
  24. 24. Morales-Núñez JJ, Muñoz-Valle JF, Torres-Hernández PC, Hernández-Bello J. Overview of Neutralizing Antibodies and Their Potential in COVID-19. Vaccines (Basel). 2021;9: 1376. pmid:34960121
  25. 25. Cromer D, Steain M, Reynaldi A, Schlub TE, Wheatley AK, Juno JA, et al. Neutralising antibody titres as predictors of protection against SARS-CoV-2 variants and the impact of boosting: a meta-analysis. Lancet Microbe. 2022;3: e52–e61. pmid:34806056
  26. 26. Khoury DS, Cromer D, Reynaldi A, Schlub TE, Wheatley AK, Juno JA, et al. Neutralizing antibody levels are highly predictive of immune protection from symptomatic SARS-CoV-2 infection. Nat Med. 2021;27: 1205–1211. pmid:34002089
  27. 27. Willett BJ, Grove J, MacLean OA, Wilkie C, De Lorenzo G, Furnon W, et al. SARS-CoV-2 Omicron is an immune escape variant with an altered cell entry pathway. Nat Microbiol. 2022;7: 1161–1179. pmid:35798890
  28. 28. Cao Y, Wang J, Jian F, Xiao T, Song W, Yisimayi A, et al. Omicron escapes the majority of existing SARS-CoV-2 neutralizing antibodies. Nature. 2022;602: 657–663. pmid:35016194
  29. 29. Ke H, Chang MR, Marasco WA. Immune Evasion of SARS-CoV-2 Omicron Subvariants. Vaccines. 2022;10: 1545. pmid:36146623
  30. 30. How COVID-19 variants evade immune response. In: National Institutes of Health (NIH) [Internet]. 7 Jun 2021 [cited 22 Jan 2023]. Available: https://www.nih.gov/news-events/nih-research-matters/how-covid-19-variants-evade-immune-response.
  31. 31. Focosi D, Maggi F, Casadevall A. Mucosal Vaccines, Sterilizing Immunity, and the Future of SARS-CoV-2 Virulence. Viruses. 2022;14: 187. pmid:35215783
  32. 32. Cao Y, Jian F, Wang J, Yu Y, Song W, Yisimayi A, et al. Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution. bioRxiv; 2022. p. 2022.09.15.507787. pmid:36535326
  33. 33. Sun K, Tempia S, Kleynhans J, von Gottberg A, McMorrow ML, Wolter N, et al. SARS-CoV-2 transmission, persistence of immunity, and estimates of Omicron’s impact in South African population cohorts. Science Translational Medicine. 2022;14: eabo7081. pmid:35638937
  34. 34. Stoddard M, Yuan L, Sarkar S, Egeren DV, White LF, Chakravarty A. The gray swan: model-based assessment of the risk of sudden failure of hybrid immunity to SARS-CoV-2. medRxiv; 2023. p. 2023.02.26.
  35. 35. Cox M, Peacock TP, Harvey WT, Hughes J, Wright DW, Willett BJ, et al. SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies. Nat Rev Microbiol. 2023;21: 112–124. pmid:36307535
  36. 36. CDC. How Flu Viruses Can Change. In: Centers for Disease Control and Prevention [Internet]. 12 Dec 2022 [cited 22 Jan 2023]. Available: https://www.cdc.gov/flu/about/viruses/change.htm.
  37. 37. Starr TN, Greaney AJ, Hilton SK, Ellis D, Crawford KHD, Dingens AS, et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell. 2020;182: 1295–1310.e20. pmid:32841599
  38. 38. Starr TN, Greaney AJ, Hannon WW, Loes AN, Hauser K, Dillen JR, et al. Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolution. Science. 2022;377: 420–424. pmid:35762884
  39. 39. Khan K, Karim F, Cele S, Reedoy K, San JE, Lustig G, et al. Omicron infection enhances Delta antibody immunity in vaccinated persons. Nature. 2022;607: 356–359. pmid:35523247
  40. 40. Powell AA, Kirsebom F, Stowe J, Ramsay ME, Lopez-Bernal J, Andrews N, et al. Protection against symptomatic infection with delta (B.1.617.2) and omicron (B.1.1.529) BA.1 and BA.2 SARS-CoV-2 variants after previous infection and vaccination in adolescents in England, August, 2021–March, 2022: a national, observational, test-negative, case-control study. The Lancet Infectious Diseases. 2022;0. pmid:36436536
  41. 41. Callaway E. COVID ‘variant soup’ is making winter surges hard to predict. Nature. 2022;611: 213–214. pmid:36307585
  42. 42. Focosi D, Quiroga R, McConnell SA, Johnson MC, Casadevall A. Convergent evolution in SARS-CoV-2 Spike creates a variant soup that causes new COVID-19 waves. bioRxiv; 2022. p. 2022.12.05.518843.
  43. 43. Simon-Loriere E, Schwartz O. Towards SARS-CoV-2 serotypes? Nat Rev Microbiol. 2022;20: 187–188. pmid:35181769
  44. 44. scipy.integrate.solve_ivp—SciPy v1.11.2 Manual. [cited 1 Sep 2023]. Available: https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.solve_ivp.html.
  45. 45. Morris SE, Pitzer VE, Viboud C, Metcalf CJE, Bjørnstad ON, Grenfell BT. Demographic buffering: titrating the effects of birth rate and imperfect immunity on epidemic dynamics. J R Soc Interface. 2015;12: 20141245. pmid:25589567
  46. 46. Martin J, Hamilton B, Osterman M, Driscoll A. Births: Final Data for 2018. National Vital Statistics Reports. 2018;68: 47.
  47. 47. U.S. Census Bureau QuickFacts: United States. [cited 22 Jan 2023]. Available: https://www.census.gov/quickfacts/fact/table/US/PST045221.
  48. 48. Stoddard M, Yuan L, Sarkar S, Mazewski M, Egeren DV, Mangalaganesh S, et al. Shielding under endemic SARS-CoV-2 conditions is easier said than done: a model-based analysis. medRxiv; 2023. p. 2023.01.22.23284884.
  49. 49. Rocklöv J, Sjödin H, Wilder-Smith A. COVID-19 outbreak on the Diamond Princess cruise ship: estimating the epidemic potential and effectiveness of public health countermeasures. J Travel Med. 2020 [cited 10 Aug 2020]. pmid:32109273
  50. 50. Liu Y, Rocklöv J. The effective reproductive number of the Omicron variant of SARS-CoV-2 is several times relative to Delta. J Travel Med. 2022;29: taac037. pmid:35262737
  51. 51. Du Z, Hong H, Wang S, Ma L, Liu C, Bai Y, et al. Reproduction Number of the Omicron Variant Triples That of the Delta Variant. Viruses. 2022;14: 821. pmid:35458551
  52. 52. Delamater PL, Street EJ, Leslie TF, Yang YT, Jacobsen KH. Complexity of the Basic Reproduction Number (R0). Emerg Infect Dis. 2019;25: 1–4. pmid:30560777
  53. 53. Bell SM, Katzelnick L, Bedford T. Dengue genetic divergence generates within-serotype antigenic variation, but serotypes dominate evolutionary dynamics. eLife. 8: e42496. pmid:31385805
  54. 54. Hernández-García E, Muñoz M de L, David RE, Pérez-Ramírez G, Navarrete-Espinosa J, Díaz-Badillo Á, et al. Epidemiological implications of the genetic diversification of dengue virus (DENV) serotypes and genotypes in Mexico. Infect Genet Evol. 2020;84: 104391. pmid:32502732
  55. 55. Vijgen L, Keyaerts E, Lemey P, Moës E, Li S, Vandamme A-M, et al. Circulation of genetically distinct contemporary human coronavirus OC43 strains. Virology. 2005;337: 85–92. pmid:15914223
  56. 56. Pyrc K, Dijkman R, Deng L, Jebbink MF, Ross HA, Berkhout B, et al. Mosaic structure of human coronavirus NL63, one thousand years of evolution. J Mol Biol. 2006;364: 964–973. pmid:17054987
  57. 57. Shukla R, Ramasamy V, Shanmugam RK, Ahuja R, Khanna N. Antibody-Dependent Enhancement: A Challenge for Developing a Safe Dengue Vaccine. Frontiers in Cellular and Infection Microbiology. 2020;10. Available: https://www.frontiersin.org/articles/10.3389/fcimb.2020.572681.
  58. 58. Ajmeriya S, Kumar A, Karmakar S, Rana S, Singh H. Neutralizing Antibodies and Antibody-Dependent Enhancement in COVID-19: A Perspective. J Indian Inst Sci. 2022;102: 671–687. pmid:35136306
  59. 59. Khoury DS, Cromer D, Reynaldi A, Schlub TE, Wheatley AK, Juno JA, et al. Neutralizing antibody levels are highly predictive of immune protection from symptomatic SARS-CoV-2 infection. Nat Med. 2021;27: 1205–1211. pmid:34002089
  60. 60. Miller J, Hachmann NP, Collier AY, Lasrado N, Mazurek CR, Patio RC, et al. Substantial Neutralization Escape by the SARS-CoV-2 Omicron Variant BQ.1.1. bioRxiv; 2022. p. 2022.11.01.514722.
  61. 61. Cheng SM, Mok CKP, Chan KC, Ng SS, Lam BH, Luk LL, et al. SARS-CoV-2 Omicron variant BA.2 neutralisation in sera of people with Comirnaty or CoronaVac vaccination, infection or breakthrough infection, Hong Kong, 2020 to 2022. Euro Surveill. 2022;27: 2200178. pmid:35514306
  62. 62. Cao Y, Song W, Wang L, Liu P, Yue C, Jian F, et al. Characterization of the enhanced infectivity and antibody evasion of Omicron BA.2.75. Cell Host & Microbe. 2022;30: 1527–1539.e5. pmid:36270286
  63. 63. Cao Y, Jian F, Wang J, Yu Y, Song W, Yisimayi A, et al. Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution. Nature. 2022; 1–3. pmid:36535326
  64. 64. Qu P, Evans JP, Faraone J, Zheng Y-M, Carlin C, Anghelina M, et al. Distinct Neutralizing Antibody Escape of SARS-CoV-2 Omicron Subvariants BQ.1, BQ.1.1, BA.4.6, BF.7 and BA.2.75.2. bioRxiv. 2022; 2022.10.19.512891. pmid:36299423
  65. 65. Zou Z, Fairley CK, Shen M, Scott N, Xu X, Li Z, et al. Critical timing and extent of public health interventions to control outbreaks dominated by SARS-CoV-2 variants in Australia: a mathematical modelling study. International Journal of Infectious Diseases. 2022;115: 154–165. pmid:34800688
  66. 66. Egeren DV, Stoddard M, White LF, Hochberg NS, Rogers MS, Zetter B, et al. Vaccines alone cannot slow the evolution of SARS-CoV-2. medRxiv; 2022. p. 2022.10.04.22280696.
  67. 67. Stoddard M, Yuan L, Sarkar S, Mangalaganesh S, Nolan RP, Bottino D, et al. Heterogeneity in vaccinal immunity to SARS-CoV-2 can be addressed by a personalized booster strategy. medRxiv; 2022. p. 2022.08.30.22279397.
  68. 68. Jaspe RC, Loureiro CL, Sulbaran Y, Moros ZC, D’Angelo P, Hidalgo M, et al. Description of a One-Year Succession of Variants of Interest and Concern of SARS-CoV-2 in Venezuela. Viruses. 2022;14: 1378. pmid:35891359
  69. 69. CoVariants: Per Country. [cited 22 Jan 2023]. Available: https://covariants.org/per-country.