Vaccine nationalism and the dynamics and control of SARS-CoV-2

Vaccines provide powerful tools to mitigate the enormous public health and economic costs that the ongoing SARS-CoV-2 pandemic continues to exert globally, yet vaccine distribution remains unequal among countries. To examine the potential epidemiological and evolutionary impacts of ‘vaccine nationalism’, we extend previous models to include simple scenarios of stockpiling between two regions. In general, when vaccines are widely available and the immunity they confer is robust, sharing doses minimizes total cases across regions. A number of subtleties arise when the populations and transmission rates in each region differ, depending on evolutionary assumptions and vaccine availability. When the waning of natural immunity contributes most to evolutionary potential, sustained transmission in low access regions results in an increased potential for antigenic evolution, which may result in the emergence of novel variants that affect epidemiological characteristics globally. Overall, our results stress the importance of rapid equitable vaccine distribution for global control of the pandemic.

The SARS-CoV-2 pandemic has led to more than 180 million infections and nearly 4 million fatalities to date (1). Effective vaccines [e.g., (2)(3)(4)] have now been approved and are actively being deployed, but numerous important questions remain. Eventually, community immunity may be attained through the deployment of vaccines; however if and when this occurs will be contingent on the characteristics of natural and vaccinal immunity (5)(6)(7). As illustrated by the rapid spread and high transmissibility of the 'delta' variant (8), SARS-CoV-2's evolutionary potential is a major potential obstacle for control (9).
Due to strong public and political pressures and fear of waning immunity, some countries with high vaccine availability are currently resorting to 'vaccine nationalism': stockpiling vaccines to prioritize rapid access to their citizenry (10). Indeed, at the time of writing, 99 and 117 doses per 100 individuals have been administered in the United States and United Kingdom, respectively, while an average of 25 and 3.8 doses per 100 individuals have been administered in India and across Africa, respectively (11). Recently, the World Health Organization recognized that delayed access to vaccines in countries with low vaccine availability may lead to more evolutionary potential (12), which could result in immune escape or other phenotypic changes of interest (e.g., increases in transmission). The emergence of future variants capable of evading natural or vaccinal immune responses could threaten containment efforts globally. These concepts underlie the development of a number of policy tools, including the existing COVAX initiative. Furthermore, to ensure that vaccine distribution is ethically-sound and equitable, the "Fair Priority Model" has been proposed (13-15) as a potential replacement to the currently planned proportional allocation (by population size) from COVAX.
Prior work has explored optimal prophylactic vaccine allocation for minimizing the final epidemic size of a fully immunizing infection (i.e., one that can be modeled using a susceptible-infected-recovered (SIR) framework) (16); when interaction between communities (or countries) is considered, equal vaccine distribution is increasingly advantageous in terms of minimizing case numbers (16). Modeling studies have also shown that coordinated influenza vaccine sharing would reduce the financial and infection burden of influenza outbreaks globally (17). Similar problems related to optimizing vaccine allocation have also been explored in networks Vaccines provide powerful tools to mitigate the enormous public health and economic costs that the ongoing SARS-CoV-2 pandemic continues to exert globally, yet vaccine distribution remains unequal among countries. To examine the potential epidemiological and evolutionary impacts of 'vaccine nationalism', we extend previous models to include simple scenarios of stockpiling between two regions. In general, when vaccines are widely available and the immunity they confer is robust, sharing doses minimizes total cases across regions. A number of subtleties arise when the populations and transmission rates in each region differ, depending on evolutionary assumptions and vaccine availability. When the waning of natural immunity contributes most to evolutionary potential, sustained transmission in low access regions results in an increased potential for antigenic evolution, which may result in the emergence of novel variants that affect epidemiological characteristics globally. Overall, our results stress the importance of rapid equitable vaccine distribution for global control of the pandemic. with community structure (18), as well as in the face of economic constraints (19) and vaccination coalition formation (20), and for SARS-CoV-2 with age-(and contact-) heterogeneity (21,22).
We have recently shown that the strength and duration of immunity elicited following infection or one or two doses of a vaccine will have a crucial impact on the medium-term epidemiological and potential evolutionary outcomes (5,6).
Here we extend these analyses to address potential epidemiological and evolutionary consequences of policies of vaccine nationalism or equitable access for a range of assumptions regarding the robustness of host immune responses. In reality, vaccine distribution is a public goods problem (23), and the optimal "global" allocation projected based on evolutionary and immunological uncertainties may differ from national optima due to the actual economic landscape of each country. In all cases, every nation, however, has a shared interest in reducing the potential for novel strains to arise, achievable by minimizing the global infection burden.
We consider a trans-national extension of our model, comprising two countries with possibly different population sizes and seasonal transmission patterns. One country, the high access region (HAR), chooses to allocate a fraction f of the total vaccine supply to the low access region (LAR). The underlying immuno-epidemiological models for both countries account for both the duration of natural and vaccinal immunity and the residual decrease in host susceptibility to infection (relative to immunologically naive individuals) after full natural or vaccinal immunity has waned; these models are described in detail in (5,6), where the more detailed structure (6) accounts for immunity after one or two vaccine doses. In the first "decoupled" framework (top panel of Fig. 1), we assume that the epidemiological dynamics of both countries are entirely independent, with the exception of their respective vaccination rates; we also compute a measure for the global potential for viral evolution of immune escape (6). In the second "coupled" framework (bottom panel of Fig. 1), we allow for immigration of infected individuals between the countries at rate η (19). Additionally, we approximate the stochastic occurrence of potential transmission increases (PTIs) in each country: briefly, if the 'potential net viral adaptation rate' [see supplementary materials, fig. S3, and (6)] exceeds a threshold, then there is a nonzero probability that the transmission rate in both countries increases. This follows evidence of enhanced binding of the SARS-CoV-2 spike protein receptor binding domain (RBD) with the ACE2 receptor in more contagious SARS-CoV-2 variants, as well as potentially higher viral loads (24). Mutational changes to the spike protein furin cleavage site (e.g., at site 681) may also contribute to increased viral transmissibility (25). In this way, this assumption represents a pessimistic scenario where the evolution of pathogen immune escape is inevitably accompanied by increases in transmission, although we also compare our results with the more optimistic scenario where transmission increases do not occur (note that we are not modeling the complexities of variant dynamics and evolution explicitly). In reality, transmission increases may plateau (26), and viral evolution may have more subtle effects on disease transmission including modulating the susceptibility of partially immune hosts. The full mathematical details for both frameworks are described in the supplementary materials.
We begin with the decoupled framework and the simpler underlying 'one-dose' vaccination model from (5) to compute the long-term equilibrium fraction of infections in both countries under a range of epidemiological and immunological scenarios. Then, with specific dosing regimes (6), we examine the short-and medium-term epidemiological dynamics and the global potential for evolution with the sharing of vaccines. Next, using the coupled framework, we compute national and combined case numbers in the medium term given different degrees of vaccine allocation from the HAR to the LAR for different immigration rates and average relative reproduction numbers. We do so for total as well as severe cases, with the expectation that the number of severe cases may be indicative of infections requiring hospitalization or the clinical burden of Covid-19, while the number of total cases reflects all infections regardless of severity. Finally, we compare the results of the coupled and decoupled frameworks for specific scenarios.

Equilibrium analyses
To obtain analytical intuition for the effect of vaccine nationalism, we first examine the simplest model of vaccine sharing, where both regions are only coupled through their vaccination rates. The dynamics of prophylactic vaccine distribution strategies are well understood when infections lead to recovery and lifelong immune protection (SIR) (16). However, natural and vaccinal immunity to SARS-CoV-2 is likely not lifelong, yet complete re-susceptibility after the waning of immunity, as is assumed in susceptible-infected-recovered-susceptible (SIRS) frameworks, is also unlikely. Instead, we examine the role of both the strength and duration of immune responses with more general SIR(S) models (5,6). We first ignore the complexities of dosing regimes and extend the model in (5) to consider vaccine sharing in the decoupled framework (top panel of Fig. 1) assuming that a single immune category exists for vaccinated individuals, and that vaccinal immunity may wane at a rate distinct from natural immunity. Since we assume that the infection dynamics in both countries are only coupled through their respective vaccination rates, a unique equilibrium of total infections exists [either disease-free or endemic; see supplementary materials and (5) for details]. To examine the long-term ν f , respectively. We examine four immunity scenarios that range from poor to robust natural and vaccinal immune responses (Fig. 2).
When the characteristics of both countries are the same, sharing vaccines always decreases or maintains the total fraction of infections at the long-term equilibrium [see Keeling (16) for the SIR extreme with a focus on two differently-sized populations]. The intuition for this result is apparent from examining the underlying values for each country (Fig. 2C). The total fraction of infections are minimized whenever one of the countries does not vaccinate beyond the rate needed for herd immunity. Additionally, sharing does not have an appreciable impact on the total fraction of infections at equilibrium when vaccination rates are too low ( Fig. 2A, top panel), or overall transmission rates are more elevated and host immune responses are poorer (Fig. 2B, top panel). Because of nonpharmaceutical interventions (NPIs) or intrinsic factors [e.g., demographics (27), population density (28), or vulnerabilities (29)], transmission rates in the two countries may be asymmetric. If there is less disease transmission in the HAR (modeled as a reduction in the transmission rate), then the 'optimal' fraction of vaccines shared to minimize the combined equilibrium fraction of infections crucially depends on the magnitude of the vaccination rate. If vaccine supplies are low and immune responses poor, then sharing only a very small fraction of the vaccine supply is epidemiologically beneficial in terms of decreasing the overall burden. For stronger immune responses, augmenting vaccine sharing rates becomes increasingly beneficial from an epidemiological perspective, as the protective effects of the vaccine are maintained for longer within the population (compare the columns of the middle row of Fig. 2A). Similarly, as vaccine supplies increase (compare the colored curves in the middle row of Fig. 2A), the minimum value of infections occurs for increasingly large values of f, or fractions of vaccines shared. Eventually, when global vaccination rates are high, even for poor host immune responses, this minimum is attained when more than half of the vaccine supply is allocated to the LAR (leftmost panel of middle row of Fig. 2A). By symmetry, the opposite occurs if there is less transmission in the LAR. These trends are further magnified if overall transmission rates are increased (Fig. 2B). To further emphasize these effects, we present the long-term equilibrium of each country under representative scenarios in Fig. 2, C to E. In particular, the comparison between Fig. 2, D and E, illustrates the importance, and indirect benefit, of increasing vaccine supply. The relative sizes of the HAR and LAR populations can have important consequences for the fraction of vaccine allocation that minimizes the weighted fraction of infections (figs. S1 and S2). Overall, these results highlight the importance of continued NPIs that decrease transmission, such as rapidtesting and physical distancing, in conjunction with ramping up vaccination and sharing vaccine supplies equitably to decrease overall burden.

Medium-term dynamics Epidemiological considerations
To consider the near-and medium-term dynamic epidemiological effects of vaccine sharing, in Fig. 3 we explore the landscapes of immunity and infections across multiple scenarios for otherwise-symmetric countries (i.e., population size and seasonal transmission rates). In all scenarios, vaccine supply is assumed to be limited initially in the HAR [modeled as a one dose policy, a lower maximal rate of administration of the first dose 0,tot ν , and no sharing (f = 0)] and then is assumed to increase. In conjunction with an increase in 0,tot ν , we allow for a transition to the recommended two-dose strategy ( 1 4 weeks ω = , right two columns of Fig. 3), and/or the initiation of equal sharing (f = 0.5) with the LAR (which is assumed to distribute vaccines using the same strategy as the HAR) may be initiated (second and fourth columns of Fig. 3). Intuitively, if one-dose immunity is robust (bottom panel of Fig. 3), then transitioning to a two-dose strategy leads to fewer individuals with robust vaccinal immunity, in turn giving rise to increases in infections in both in the short-and medium-term [compare the corresponding scenarios of the bottom panel of Fig. 3 and see also (6)]. In such a situation, 'one-dose' strategies (i.e. either the first dose of a 2-dose vaccine or the unique dose of a 1-dose vaccine) with equal sharing between countries suppress overall burden. On the other hand, if one-dose immunity is poor, switching to the recommended two-dose regimen prevents the accumulation of individuals with waned one-dose immunity and thus potentially larger infections peaks in the longer term (top panel of Fig. 3). If poor one-dose immunity nevertheless reduces severity of infection after waning (unlike our pessimistic assumption), then the predicted clinical burden of severe cases would likely be lower. Finally, if one-dose immunity is poor and a one-dose policy is pursued, the first infection peak after ramping up vaccination in the HAR may be higher without sharing (top left panel of Fig. 3). This counter-intuitive finding arises due to the large accumulation of individuals with waned one-dose immunity who experience infection. This highlights the important role for population-level susceptibility (modulated by natural and vaccinal immune responses) and its dynamical interplay with transmission in determining the timing and burden of infections. The effects of different NPI scenarios, transmission patterns, and vaccination rates in either the HAR or LAR can be further explored with the online application (30).

Evolutionary considerations
The accumulation of individuals with various states of immunity (i.e., waned one-dose immunity or immunity following natural infection) may also lead to different evolutionary outcomes depending on the vaccine sharing scheme pursued. Current evidence suggests that adaptive immune responses following natural infection with SARS-CoV-2 are fairly robust and long-lasting (31,32). This protection may be less certain in the context of subsequent infection with variant strains (33,34). Encouragingly, studies indicate that previously-infected hosts are also largely protected (clinically and against breakthrough infections) against emerging variants after a single vaccine dose (35). However, recent work indicates poorer protection against the rapidly spreading 'delta' variant of concern particularly after a single dose (36), although protection against severe disease still appears fairly robust (37)(38)(39). Interestingly, the protection conferred by a single vaccine dose in previously infected individuals may be even greater than that provided by two vaccine doses in naÃ¯ve hosts (40). However, there are still many immunological uncertainties, e.g., the duration and longer-term strength of this protection against existing strains and potential emerging variants remain unknown. In fig. S3, we use the evolutionary framework from (6) to project the potential net viral adaptation rate (see supplementary materials for details and online application for additional scenarios). Overall, we find that uncertainties in evolutionary outcomes dominate our projections, echoing previous findings (6) ( fig. S3 and online application). However, if the evolutionary potential for immune escape is highest among infections in hosts with natural immunity, the framework predicts that sharing vaccines always decreases global evolutionary potential (green curves of fig. S3) in the decoupled framework. Overall, when immunity after a single vaccine dose is robust, natural and vaccine-derived immunity will limit damaging pathogen evolution relative to the scenario with poor single dose immunity (compare the top and bottom panels of fig. S3). However, if immunity is partial or waning, ongoing transmission might accelerate adaption, supporting the need for continued monitoring of variants and their interaction with natural and vaccine-derived immunity.
Another intuitive result is that, in the LAR, sharing vaccines leads to increases in population immunity (for any dosing regime) and thus a decrease in infections and burden in the short term, even with poor one-dose immunity. In general, in the HAR, sharing decreases population immunity and increases infections in the short term. However, these changes are minimal and likely acceptable given the combined decrease in infections, illustrating the long term benefits of vaccine sharing. In particular, while local cases in the HAR may increase in the short-term, the longer-term disease risk would be lowered by sharing due to a decrease in the potential for the evolution of more transmissible or immuneescape viral variants.

Coupled framework
So far, we have assumed that the countries have decoupled disease dynamics. This simplification for tractability ignores infection importation as well as the possible emergence of variants from regions with more persistent infections. The issues that arise from the global circulation of SARS-CoV-2, particularly the variants of concern, are of considerable public health importance. Thus, we next explore these effects using the coupled framework presented in the bottom panel of Fig. 1.

Cumulative cases and pathogen evolution in the medium-term
In Fig. 4, we plot the cumulative number of total and severe cases (see supplementary materials for details) assuming equal population sizes in both countries from the time of vaccine introduction until 5 years after the pandemic onset in the HAR, LAR, and combined, as well as the projected number of PTIs to have occurred in both regions by the end of the 5 year period. We do so for various vaccine allocation fractions between the HAR and LAR, as well as a range of immigration rates assuming symmetric transmission rates (Fig. 4, A and C) and relative mean reproduction numbers assuming constant immigration rates (Fig. 4, B and D; see Materials and methods). In Fig. 4, A and B, we assume that infection following waned natural immunity contributes the most to viral adaptation, while in Fig. 4, C and D, we assume that infection following waned vaccinal immunity, and one-dose immunity in particular, contributes the most.

Dependence on immigration rate
For equal population sizes and symmetric transmission rates, a weak dependence of total and severe cases as well as PTIs on the immigration rate η is observed, regardless of whether infection after waned natural or vaccinal immunity is assumed to contribute more to viral evolution (Fig. 4, A and C). Additionally, more equitable vaccine distribution minimizes total PTIs and combined cases in both scenarios. When natural infections contribute more to evolution (Fig. 4A) and η is low, the HAR must retain an increasing share of the vaccines to minimize local total cases as the immigration rate increases, but this is done at the expense of more cases in the LAR and PTIs. Notably, for f ≈ 0, severe cases are minimized regardless of the immigration rate in the HAR and maximized in the LAR, which may have important clinical consequences. When infection following waned vaccinal immunity contributes more to viral evolution, this approach is no longer advantageous for the HAR, and the retained vaccine fraction sets the observed case numbers nearly independently of the assumed immigration rate. Further, large asymmetries in vaccine sharing (i.e., f = 0 or f = 1) result in much more marked relative numbers of PTIs in this scenario. For the same total vaccine availability, the most realistic population asymmetry is that the LAR has a larger population. Under this condition ( fig. S4) and when infection after waned natural immunity contributes more to evolution, total cases in the LAR are relatively insensitive to the amount of vaccine allocated, except for very large f; however severe cases can be substantially reduced with vaccine sharing. Here, combined total cases and PTIs are minimized by minimizing cases in the HAR. When a greater number of total vaccines are available (i.e., a larger 0,tot ν ; fig. S5), total and severe cases in the more populous LAR decrease approximately monotonically with increasing vaccine allocation f and more equitable vaccine allocation once again minimizes combined total and severe cases. When infection after waned vaccinal immunity contributes more to evolution, the trends are more similar to those for symmetric population sizes, and more equitable vaccine sharing is favored for both vaccination rates 0,tot ν given a larger population in the LAR (figs. S4 and S5).
When the LAR has a smaller population ( fig. S6), a relatively weak dependence of total and severe case numbers and PTIs on η is still observed for both evolutionary scenarios, particularly for higher immigration rates. However, and particularly when infection after waned natural immunity contributes more to evolution, the minima in combined cases and PTIs are now observed for f < 0.5, i.e., when the HAR retains more than half of the available vaccines. Importantly, we note that additional booster doses may further change these landscapes of immunity, and consequently the projected burdens of total and severe cases.

Dependence on relative transmission rate
The number of total and severe cases and PTIs show a greater sensitivity to the average reproduction number ratio 0,LAR 0,HAR / R R between the two countries for a fixed immigration rate. Intuitively, for equal population sizes (Fig. 4) and when 0 R in the LAR country is lower, having the HAR retain more than half of the vaccines (f < 0.5) is a good strategy for minimizing total PTIs and cases, regardless of the evolutionary scenario. Indeed, the optimal vaccine allocation shifts closer and closer to equal sharing as the 0 R values of both countries approach each other, along with an increase in cases. These trends are similar when the LAR has a larger or smaller population than the HAR (figs. S4 and S6, respectively).
When 0 R in the LAR is higher, trends are more complex. In general, regardless of the relative population size or evolutionary scenario, for a given 0,LAR 0,HAR / R R , cases in the HAR decrease with increasing vaccine retention (smaller f), while cases in the LAR increase ( Fig. 4 and figs. S4 to S6). Severe cases in each region in particular are strongly reduced by increased vaccine availability. The increase in cases in the LAR is increasingly large at higher 0,LAR 0,HAR / R R . When infection after waned natural immunity contributes more to evolution and vaccine supply is sufficiently high to reduce case numbers in the LAR, PTIs are numerous when the HAR retains a large fraction of the vaccines ( Fig. 4B and figs. S5B and S6B). This is due to sustained elevated case numbers in unvaccinated individuals in the LAR. On the other hand, when infection following waned vaccinal immunity contributes more to evolution, then having the HAR (with lower transmission rates) retain a larger fraction of the vaccines minimizes PTIs for any relative population size, since the high 0 R in the LAR would result in large subsequent peaks containing individuals whose vaccinal immunity has waned with sharing ( Fig.  4D and figs. S4D, S5D, and S6D). However, this strategy also leads to highly elevated case numbers, including severe cases, in the LAR. Since LARs may also have more fragile healthcare systems, the elevated clinical burden under this scenario may be particularly problematic.

Comparison of coupled and decoupled frameworks
We note that assumptions of large 0 R in the LAR also result in very large initial infection peaks, which increase community immunity in the medium-term. These initial waves are not reflected in the total case counts, however, since these values are summations from the time of vaccine initiation through the end of the 5 year period after the onset of the pandemic. Further, in Fig. 5 we illustrate the temporal effect of the coupled framework on the infection dynamics in the LAR and HAR relative to the first model with no immigration or explicit effect of PTIs on the transmission rate. When moderate asymmetry in 0 R is assumed ( 0,LAR 0,HAR / 1.2 R R = ), simulations using the decoupled framework suggest that a strategy where the HAR retains all vaccines would be highly beneficial for that country (top panel of Fig. 5; no PTIs are projected to occur, and low case numbers are observed throughout). However, this occurs at the expense of PTIs and infection burden in the LAR, which are both substantially higher. With the more realistic coupled framework, immigration and increases in transmission illustrate that this strategy is far less beneficial to the HAR than the decoupled framework would suggest, as substantially higher case numbers and pathogen evolution are predicted in this region. Although total cases in the HAR increase slightly when vaccines are equally distributed under the coupled framework (lower panel of Fig. 5), substantial reductions in case numbers in the LAR result in fewer PTIs in that country, and total combined case numbers are also slightly lower. To untangle the effects of immigration and PTIs on dynamics in the coupled framework, we reproduce Figs. 4 and 5 allowing for immigration only in figs. S7 and S8, respectively. In other words, these figures represent a more optimistic evolutionary scenario in which the occurrence of a PTI does not increase transmission rates in either the HAR or LAR. Overall, we show that vaccines play an important role in minimizing cases (particularly severe cases) as well as potential viral adaptation in both regions. We also emphasize that imperfect vaccinal and natural immunity and asymmetries in population size and transmission rates add many nuances to this picture.

Caveats
A full list of caveats and future directions is presented in the supplementary materials; we briefly summarize them below. First, building on prior work and in order to focus on qualitative features, we ignore heterogeneities within countries, such as due to age (41) or superspreading (42). We have also assumed simple scenarios for nonpharmaceutical interventions as in (6), and assumed that the seasonal transmission rates are similar in both countries. More granular, well-parameterized epidemiological models with these complexities would lead to more accurate quantitative predictions. Additionally, underlying differences in susceptibility between the HAR and LAR, through factors such as variations in age structure or degrees of pre-existing population-level immunity, may result in differential rates of severe cases and hospitalizations independent of vaccination. An explicit treatment of age-structure and serostatus [as in (21)] in modeling vaccine roll-out to estimate clinical burdens in HARs and LARs under different sharing schemes should be investigated further. We also note that our models focused on the situation where dose-administration is contingent on supply, so that 'sharing' vaccines from the HAR to the LAR decreases the vaccination rate in the HAR. However, certain HARs may have access to more doses than they can administer. In such a scenario, sharing would be even more beneficial: it would lead to no change in the HAR vaccination rate while increasing that of the LAR. Moving forward, additional booster doses may be administered that could alter population-level immune landscapes, and including these in future models will be important for qualitative and quantitative predictions. Additionally, we omit vaccine hesitancy (23), though simple extensions of our previous models with hesitancy (5,6) could examine the resulting interplay with vaccine nationalism. Furthermore, our equilibrium analyses highlight the long-term importance of vaccine sharing regardless of when sharing is initiated. However, our scenarios that examine short-and medium-term implications strongly emphasize the benefits of close to "parallel" or concurrent vaccine sharing. The effects of delaying the initiation of sharing, say until a certain proportion of individuals in the HAR are vaccinated, should be carefully examined.
Lastly, we have assumed the simplest evolutionary models, both for determining potential viral adaptation rates as well as for simulating potential increases in transmission rates. As more data become available, these should be refined (43)(44)(45)(46), with possible directions including extending the model to explicitly track the transmission of different strains, and accounting for potential reductions in the strength of vaccinal immunity (i.e., the parameters due to the emergence of novel immune escape variants (47).
In particular, prolonged infections in immunosuppressed hosts may be disproportionately important for evolution (48). The online interactive application (30) allows for an in-depth exploration of the effect of different climate-driven seasonal transmission rates as well as a broad range of assumptions related to NPIs and immuno-epidemiological parameters.

Conclusion
Even as vaccine production increases, a number of countries are choosing to share little or no vaccines with countries that have very low vaccine availability. Vaccine nationalism, dosing regimes, and host immune responses have important interactive effects, and these will substantially shape epidemiological dynamics and evolutionary potential in the medium term. Additionally, unstable vaccine supply will also increase variability in the timing or availability of first and second doses. Using extensions of our prior work (5, 6), we incorporated vaccine sharing scenarios in two countries whose infection dynamics are either otherwise independent or coupled through immigration of infectious individuals and evolutiondriven increases in transmission rates. When country profiles are symmetric, we find that sharing vaccines with countries that have low availability decreases overall infections and may also mitigate potential antigenic evolution. Asymmetries in population size or transmission rates introduce additional complexities, which are particularly marked when natural and vaccinal immunity is weak. Nevertheless, our models indicate that the prompt redistribution of vaccine surpluses is likely advantageous in terms of epidemiological and evolutionary outcomes in both countries and, by extension, globally. Ethical arguments also support this policy (13,14). Persistent elevated disease transmission in countries with low vaccine availability also substantially undermines attempts at infection control via stockpiling in the country with high vaccine availability, which is not accounted for when disease transmission in both countries is assumed to be decoupled. Overall, our work highlights the importance of continued efforts in quantifying the robustness of immunity following vaccination. Furthermore, reevaluation of stockpiling policies as vaccine supplies increase is imperative, and ramping up global vaccination efforts is crucial.

Determination of seasonal reproduction numbers
In order to reflect observed seasonal variation in transmission rates for respiratory infections arising from related coronaviruses (49), influenza (49) and respiratory syncytial virus (50), we base seasonal reproduction numbers in this work on those in (5), which were calculated in (49) based on the climate of New York City. Other seasonal patterns can be explored using the interactive online application. In all simulations, we modify these values to force a mean value for the basic reproduction number of ( )

Modeling of nonpharmaceutical interventions (NPIs)
In all simulations, we enforce periods of NPI adoption (arising from behaviors and policies such as lock downs, maskwearing, and social distancing) in which the transmission rate is reduced from its seasonal value described in the previous section.
For the simulations of medium-term dynamics using the decoupled framework (i.e., Fig. 3 and fig. S3), we used the same NPI scenarios as in (6). For all simulations using the coupled framework (i.e., Figs. 4 and 5 and figs. S4 to S8), we assume that NPIs are adopted between weeks 8 and 44 following the pandemic onset resulting in the transmission rate being reduced to 40% of its seasonal value. Between weeks 45 and 79, we assume that the transmission rate is 60% of its seasonal value; higher than during the previous period due to either behavioral changes following the introduction of the vaccine or the emergence of more transmissible strains. Finally, we assume that NPIs are completely relaxed beyond week 80.  ) through the end of the 5 year period for both countries (leftmost two columns), the HAR (third and fourth columns from the left), the LAR (fifth and sixth columns from the left), as well as the combined number of PTIs to have occurred in both countries at the end of 5 years (rightmost column). Each grid-point denotes the mean value of 100 simulations. The population of both countries is taken to be the same. Each area plot is internally normalized, such that the largest value in each plot is 1. The x-axis indicates the fraction of vaccines retained by the HAR (i.e., 1 − f); thus the far right of a plot is the scenario where the HAR retains all vaccines (f = 0). In (A) and (C), both countries have the same average transmission rate ( 0 R , see Materials and methods), and the immigration rate η is varied. In (B) and (D), the immigration rate is fixed at η = 0.01, and the relative mean transmission rate in the LAR, i.e.,