key: cord-0004009-zs8urh6c
authors: Nuismer, Scott L.; Basinski, Andrew; Bull, James J.
title: Evolution and containment of transmissible recombinant vector vaccines
date: 2019-06-12
journal: Evol Appl
DOI: 10.1111/eva.12806
sha: bafab6b3dd88dcdefe111698d02f81998c9accdb
doc_id: 4009
cord_uid: zs8urh6c

Transmissible vaccines offer a revolutionary approach for controlling infectious disease and may provide one of the few feasible methods for eliminating pathogens from inaccessible wildlife populations. Current efforts to develop transmissible vaccines use recombinant vector technology whereby pathogen antigens are engineered to be expressed from innocuous infectious viral vectors. The resulting vaccines can transmit from host to host, amplifying the number of vaccine‐protected individuals beyond those initially vaccinated directly through parenteral inoculation. One main engineering challenge is the potential for natural selection to favor vaccine mutants that eliminate or reduce expression of antigenic inserts, resulting in immunogenic decay of the vaccine over time. Here, we study a mathematical model of vector mutation whereby continuous elimination of the antigenic insert results in reversion of the vaccine back into the insert‐free vector. We use this model to quantify the maximum allowable rate of reversion that can be tolerated for a transmissible vaccine to maintain a critical threshold level of immunogenicity against a target pathogen. Our results demonstrate that even for transmissible vaccines where reversion is frequent, performance will often substantially exceed that of conventional, directly administered vaccines. Further, our results demonstrate the feasibility of designing transmissible vaccines that yield desired levels of immunogenicity, yet degrade at a rate sufficient for persistence of the recombinant vaccine within the environment to be minimized.

large public health impacts. For instance, the 2014 West African Ebola epidemic killed more than 11,000 people over the course of the 30-month epidemic (Kaner & Schaack, 2016) , and outbreaks of SARS and MERS have killed hundreds (de Wit, Doremalen, Falzarano, & Munster, 2016) . Although large-scale vaccination campaigns of wild animals have shown the ability to reduce spillover events in temperate regions (e.g., Europe and North America for rabies in red foxes and raccoons), similar efforts cannot provide the coverage needed for control of most emerging zoonoses because of difficulties in reaching wildlife populations with conventional vaccines in the more remote and inaccessible environments where these pathogens frequently emerge.

One promising innovation for the vaccination of wildlife is the use of transmissible, rather than individually administered, vaccines.

Conventional, nontransmissible vaccines protect only those individuals directly vaccinated (e.g., administered parenterally or through bait). In contrast, transmissible vaccines have the potential to reach more of the targeted wildlife population through animal-to-animal transmission. A transmissible vaccine works on the same principles as other live vaccines-following administration to an individual, the vaccine leads to stimulation of a protective immune response-but with its capacity for transmission, the vaccine is able to spread and immunologically protect other susceptible individuals beyond those directly vaccinated.

Results from mathematical modeling indicate that transmissible vaccines substantially increase vaccine coverage and thereby the scope for pathogen elimination or eradication Nuismer et al., 2016; . To date, transmissible vaccines remain largely untested. However, a transmissible vaccine comprised of an attenuated myxoma virus encoding an antigen from rabbit hemorrhagic fever virus (RHFV) was able to transmit between rabbits and proved protective against both myxoma and rabbit hemorrhagic fever (Angulo & Barcena, 2007; Bárcena et al., 2000) .

Efforts are now underway to develop transmissible vaccines targeting other pathogens in wildlife populations including Ebola and Lassa fever viruses in their natural animal reservoirs (Marzi et al., 2016; Murphy, Redwood, & Jarvis, 2016; Tsuda et al., 2011 Tsuda et al., , 2015 .

Of the two most common types of live vaccines, attenuated and recombinant vectored, the latter appears to be the most promising for development as a transmissible vaccine platform (Bull, Smithson, & Nuismer, 2018; Murphy et al., 2016) . Recombinant vector vaccines (RVVs) consist of an antigen-expressing pathogen gene(s) engineered into a competent but avirulent viral vector (Bull et al., 2018) . The possible combinations afforded by this approach appear almost unlimited, and the ability to choose vectors from a wide range of possibilities should permit rates and pathways of vaccine transmission to be tuned to at least some degree. However, recombinant vaccines face regulatory issues due to their genetically modified status, and the ability to transmit-a big advantage in pathogen control-may provide a hurdle to regulatory approval. There are multiple means by which transmissible RVVs may overcome potential regulatory issues for their use in wildlife populations. If the vaccine is only weakly transmissible (e.g., R 0 < 1), the vaccine will inevitably approach extinction once it is no longer being directly administered to the target population (Nuismer et al., 2016) .

Similarly, if the vector occurs naturally in the host population, vectorspecific immunity may lead to vaccine extinction .

Finally, if the RVV is genetically unstable and progressively evolves to delete or downregulate the inserted transgene, the vaccine may revert over time, leaving behind only the original viral vector.

Although previous work has studied the efficacy of weakly transmissible vaccines and explored how cross-immunity between vector and vaccine can influence vaccine performance Nuismer et al., 2016) , the consequences of genetic instability for transmissible vaccines have not yet been rigorously evaluated. As a result, it is currently unclear what level of instability can be tolerated for a transmissible vaccine to remain effective, or how unstable a transmissible vaccine needs to be to guarantee its ultimate extinction. This balance between immunogenicity and target antigen loss will dictate whether vaccine decay is slow enough to enable broad protection against a target pathogen but fast enough to ensure vaccine extinction in a reasonable time frame.

In the present study, we develop and analyze mathematical models and computer simulations of transmissible recombinant vector vaccines. Our models explicitly integrate genetic instability by allowing the vaccine to revert to viral vector. We consider scenarios where cross-immunity between vector and vaccine is absent as well as scenarios where cross-immunity is strong. We use these models to identify the amount of genetic instability that can be tolerated for a transmissible vaccine to remain effective, and to evaluate the scope for tuning genetic instability for maximal efficacy while ensuring vaccine extinction.

Recombinant vector vaccines are assumed to exist in large well-mixed populations, where mutation combined with within-host evolution continuously causes reversion of the vaccine to the insert-free vector at a rate, μ, either through insert deletion or through downregulation of antigenic insert expression. Although we represent reversion as a simple single-step process, in reality it depends on a complex progression of processes (e.g., mutation, selection, sampling) playing out within-host individuals. We assume the vaccine is continuously administered directly to a fraction, σ, of newborn animals. Although this may be difficult to accomplish in many wild animal populations, it may be possible using the distribution of baits or by directly vaccinating and releasing animals reared in captivity. Our model assumes directly vaccinated animals become infected with the vaccine and recover from vaccine infection at rate, V . Upon recovery, previously infected animals develop lifelong immunity against the target pathogen (hence against the vaccine). The degree of cross-immunity between viral vector and intact vaccine has been shown to influence the efficacy of recombinant vector transmissible vaccines . Here, we study two models that capture the extreme ends of the cross-immunity spectrum. In the first model, we assume a complete absence of cross-immunity; the dynamics of vaccine and vector are independent of each other (Figure 1a ). In the second model, we investigate the opposite extreme, where complete cross-immunity between vector and vaccine exists, leading to strong competition between them for susceptible hosts (Figure 1b ).

An absence of cross-immunity might accrue to a RVV if prior vaccine exposure and the resultant host adaptive immune response does not biologically impact vaccine re-infection as reported for cytomegalovirus- (Gandhi & Khanna, 2004) or if immunity to the foreign antigen is highly immunodominant, minimizing other immunological responses to the vector. In the absence of cross-immunity, the dynamics of the transmissible vaccine and pathogen can be described without explicitly tracking the density of vector generated by reversion: When a reversion transforms vaccine to vector, the host that was vaccine-infected merely moves back into the pool of susceptible hosts. Consequently, we need follow only the densities of hosts susceptible to vaccine and pathogen (S), infected with vaccine (V), infected with pathogen (P), or recovered from vaccine infection (R V ) or pathogen infection (R P ) and immune to infection by both vaccine and pathogen. The following system of differential equations applies:

where all parameters and variables are defined in Table 1 . For simplicity, we assume immunity to the pathogen and vaccine is perfect and lifelong.

Although some viral vectors are expected to elicit little if any crossimmunity (i.e., CMV), in other cases at least some degree of crossimmunity is likely. To bracket the range of possible scenarios, our second model assumes the opposite extreme of the first: complete cross-immunity between vector and vaccine. Integrating cross-immunity necessitates a more complex model where the dynamics of the vector are tracked explicitly. In this model, when a reversion transforms the vaccine into vector, the revertants are removed from the vaccine-infected class (V) and moved to the vector-infected class (M). Because the vector no longer carries or expresses pathogen genes, vector and pathogen share no cross-immunity and can co-infect. In contrast, because we assume the vaccine carries both pathogen and vector genes that elicit an immune response, crossimmunity exists between vaccine and pathogen and vaccine and vector. We assume this cross-immunity is complete and precludes co-infection . These assumptions give rise to the following system of differential equations:

Flowcharts of the model without (a) and with (b) cross-immunity.

All classes die at rate d (not shown). In the model with cross-immunity, shading indicates classes that are infected with pathogen (pink) and vector (blue)

where P = P S + P M + C and M = M P + M S + C are the total number of animals infected with pathogen and revertant vaccine, respectively, and all other parameters and variables are as defined in Table 2 .

Transmission rates and recovery rates of vector and vaccine are identical, and immunity is perfect and lifelong (as before).

We begin our analysis by identifying conditions where the level of transmissible vaccine administration precludes invasion by the pathogen and the transmissible vaccine can thus be used to protect a naïve population from a pathogen. The pathogen will be characterized by its basic reproductive number (denoted here as 

In the absence of cross-immunity between vaccine and viral vector, identifying conditions where a transmissible vaccine precludes pathogen invasion is straightforward (Appendix 1). A population can be protected from pathogen invasion if the transmissible vaccine is administered to newborn individuals at a rate that exceeds the following threshold:

where R 0,V is the basic reproductive number of the vaccine. This result generalizes our previous work for a transmissible vaccine that does not decay (Nuismer et al., 2016) .

This threshold can be used to define the conditions under which a transmissible vaccine will outperform a standard vaccine. Specifically, dividing the right-hand side of (3) by the vaccination rate required to achieve prophylaxis using a conventional vaccine 1 − 1 R 0,P , yields the following expression for the relative performance of a transmissible vaccine that experiences reversion to viral vector:

If ρ is less than one, a transmissible vaccine will outperform a conventional, directly administered vaccine. If, in contrast, ρ is greater than one, a transmissible vaccine will perform worse than a traditional vaccine. This result seems to defy intuition-that transmission could be worse than no transmission. However, it stems simply from the assumption that individuals infected with the transmissible vaccine are allowed to decay into susceptibles (at rate μ) without first gaining immunity to the target pathogen. Thus, individuals that are directly Inspection of (4) demonstrates that the relative performance of a transmissible vaccine depends on its rate of reversion to insert-free vector ( Figure 2 ). Specifically, if the reversion rate is small relative to host death rate, d, and recovery rate from vaccine/vector infection, V , reversion has only a very slight impact on vaccine performance ( Figure 2 ). In contrast, if the reversion rate is large relative to death and recovery rates, (4) demonstrates that reversion can have a substantial adverse effect on the performance of the transmissible vaccine ( Figure 2 ). This implies that, all else being equal, genetically unstable transmissible vaccines that generate short-lived infections (i.e., large recovery rate, V ) will be more effective than those that generate long-term infections from which the host recovers only very slowly (i.e., small recovery rate, V ). The reason the rate of recovery matters, of course, is that long-lived infections offer substantially greater scope for reversion to vector and thus may undermine the effectiveness of genetically unstable transmissible vaccines.

When cross-immunity between vector and vaccine exists, identifying conditions under which a transmissible vaccine protects a population from pathogen invasion becomes considerably more complex (Appendix 2). Consequently, although we successfully solved for the critical level of direct vaccination required to prevent pathogen invasion, the expression becomes unwieldy. Instead, we report the key results gleaned from numerical investigations. Plotting the relative performance of a transmissible vaccine as a function of its basic reproductive number, R 0,V , and reversion rate to vector, μ, demonstrates that any cross-immunity between vector and vaccine substantially reduces the scope for a transmissible vaccine to outperform a traditional vaccine ( Figure 3 , compare to Figure 2 ). This result generalizes previous work by showing that the rate of reversion to vector plays an important role. Specifically, in the presence of cross-immunity, even rates of reversion as low as 1.5 × 10 −2 day −1 (one reversion every 66.7 days, on average), can cause a highly transmissible vaccine to be outperformed by a conventional and directly administered vaccine ( Figure 3 ). Here too, we find that, all else being equal, transmissible vaccines with a rapid host recovery rate are more effective than those with a slow recovery rate. The benefit of rapid host recovery is, once again, the reduction in opportunities for within-host evolution to cause reversion of genetically unstable transmissible vaccines and the production of insert-free vector with which the vaccine must effectively compete.

These results demonstrate that when cross-immunity between vector and vaccine exists, the development of effective transmissible vaccines will depend on engineering highly stable recombinant vector vaccines with low rates of reversion to vector.

The foregoing results identify conditions where a transmissible vaccine outperforms a conventional, nontransmissible vaccine with respect to protecting a naïve population from pathogen invasion.

Although these results identify a baseline set of conditions for determining whether a transmissible vaccine will be more effective than 

The vaccination effort (relative to a nontransmissible vaccine) required for pathogen eradication using a recombinant vector transmissible vaccine characterized by particular combinations of reversion rate to insert-free vector (y axis) and vector R 0,V (x axis). Cross-immunity between vector and vaccine was assumed to be absent. The gray shaded area indicates parameter combinations for which more vaccination effort is required for a transmissible vaccine than a nontransmissible vaccine. In the top panel, the vaccine is constructed from a vector characterized by a recovery rate that is large relative to the reversion rate ( V = 0.05), whereas in the bottom panel, the vector has a recovery rate that is low relative to the reversion rate ( V = 0.002809). In both panels, the pathogen had an R 0,P = 2, with remaining background parameters given by: b = 10, d = 0.002809, and P = 0.05 

The vaccination effort (relative to a nontransmissible vaccine) required for pathogen eradication using a recombinant vector transmissible vaccine characterized by particular combinations of reversion rate to insert-free vector (y axis) and vector R 0,V (x axis). Cross-immunity between vector and vaccine was assumed to be complete. The gray shaded area indicates parameter combinations for which more vaccination effort is required for a transmissible vaccine than a nontransmissible vaccine. In the top panel, the vaccine is constructed from a vector characterized by a recovery rate that is large relative to the reversion rate ( V = 0.05), whereas in the bottom panel, the vector has a recovery rate that is low relative to the reversion rate ( V = 0.002809). In both panels, the pathogen had an R 0,P = 2, with remaining background parameters given by: b = 10, d = 0.002809, and P = 0.05 if a recombinant vector transmissible vaccine is constructed by inserting pathogen genes with antigenic activity into an endemic viral vector, mutations that delete the inserted genes transform vaccine back into vector and may ultimately return all vaccine descendants to their pre-engineered state. Thus, mutation/reversion may enable self-extinguishing transmissible vaccines that eradicate a pathogen, but then decay. In this section, our goal is to identify the conditions under which a recombinant vector transmissible vaccine will decay, making the idea of a "clean" or "self-extinguishing" transmissible vaccine a real possibility.

In the absence of the pathogen and without cross-immunity between vector and vaccine, the only limit to the spread of a recombinant vector transmissible vaccine is reversion. In Appendix 3, we

show that a recombinant vector transmissible vaccine will ultimately be displaced by vector anytime the following condition holds:

Thus, the more transmissible the vaccine (increasing values of R 0,V ), the greater the rate of reversion must be for the vaccine to be self-extinguishing (Figure 4 ). Because these conditions are in direct opposition to those that increase the performance of a transmissible vaccine, an effective yet self-extinguishing transmissible vaccine must walk a tightrope between having a low enough reversion rate to be effective against the pathogen and yet a high enough reversion rate to self-extinguish once direct vaccination ceases. Later, we will explore the scope for a transmissible vaccine to navigate this tightrope using stochastic simulations.

In the presence of complete cross-immunity, both reversion and competition with viral vector act in opposition to the spread of a transmissible vaccine. Consequently, results derived in Appendix 4 demonstrate that the only condition for such a vaccine to be self-extinguishing is that reversion occurs. Thus, even though we have ignored potential costs of carrying genetic cargo, the presence of cross-immunity guarantees that even a slight rate of reversion from vaccine to vector tips the balance in favor of vector and guarantees the ultimate extinction of the transmissible vaccine.

This process may, however, be slow if rates of reversion are not appreciable.

Given the preceding results, the scope for engineering a transmissible vaccine that is highly effective, yet guaranteed to self-extinguish once direct vaccination ceases, is unclear. Here, we use stochastic simulations to explore the scope for a transmissible vaccine to eradicate a target pathogen, but ultimately extinguish itself from the population. Specifically, we use the Gillespie algorithm (Gillespie, 1977) to simulate the systems of differential Equations (1) and (2) in finite populations where stochastic extinction of the pathogen and vaccine is possible. The general scenario we focus on is one where

The range of reversion rates to insert-free vector (y axis) and vector R 0,V (x axis) for which a recombinant vector transmissible vaccine is expected to be self-extinguishing in the absence of cross-immunity between vector and vaccine. The gray shaded area indicates parameter combinations for which the vaccine is self-extinguishing. In the top panel, the vaccine is constructed from a vector characterized by a recovery rate that is large relative to the reversion rate ( V = 0.05), whereas in the bottom panel, the vector has a recovery rate that is low relative to the reversion rate ( V = 0.002809). 

In the absence of cross-immunity, our results show that a transmissible vaccine can be an effective tool for eradicating/eliminating an endemic pathogen even when significant reversion occurs ( Figure 5 ). At the same time, however, our results show that for such a vaccine to self-extinguish, reversion must increase in proportion to the vaccine R 0,V . This sets up a conundrum where the better able the vaccine is to immunologically protect against the pathogen, the harder it is for the vaccine to self-extinguish. We investigated the scope for resolving this conflict by identifying parameter combinations that allowed the vaccine to first eradicate the pathogen and then subsequently self-extinguish ( Figure 5 ). This analysis confirmed the challenge of engineering an effective, yet self-extinguishing, transmissible vaccine and demonstrates that reversion rate, μ, and vaccine/vector basic reproductive number, R 0,V , must fall within a narrow band of combinations ( Figure 5 ).

In addition, all else being equal, as the recovery rate ( V ) increases, the range of suitable parameter combinations also increases. Thus, stochastic simulations support the deterministic results in demonstrating that developing an effective, yet self-extinguishing, transmissible vaccine in the absence of cross-immunity is theoretically possible but requires challenging engineering and reliable estimates of key vaccine parameters.

Stochastic simulations demonstrate that cross-immunity between vector and vaccine reduces the scope for pathogen eradication but has little overall impact on opportunities for designing an effective and self-extinguishing transmissible vaccine ( Figure 6 ). The reason is that, although cross-immunity makes pathogen eradication more challenging, it practically guarantees that a transmissible vaccine self-extinguishes from competition with vector, although this process may be slow when rates of reversion are not appreciable. On balance, then, cross-immunity does little to reduce opportunities for designing an effective, self-extinguishing, vaccine. Cross-immunity does, however, affect the parameter combinations that enable success, requiring increased transmission rates and reduced reversion rates (to see this, compare the bottom rows of Figures 5 and 6 ). Just as in the absence of cross-immunity, the scope for engineering an effective and self-extinguishing transmissible vaccine is substantially greater when recovery from vector/vaccine infection is relatively rapid (i.e., V is relatively large). Also, as we saw before, successful engineering will require that key parameters-such as rates of reversion and transmission-can be accurately measured and finely tuned.

The rate at which a recombinant vector transmissible vaccine reverts to insert-free vector (μ) is critical for predicting both vaccine efficacy and its loss from the population when vaccination ceases. There are as yet no direct estimates of this parameter, and because it involves the interplay of multiple within-host processes (i.e., mutation, selection, sampling), it also cannot be easily inferred from existing data on those constituent processes. What we do know, however, is that the stability of inserted transgenes, in general, varies substantially (Duch, Carrasco, Jespersen, Hansen, & Pedersen, 2004; Knowles et al., 2009; Logg, Logg, Tai, Cannon, & Kasahara, 2001; Malczyk et al., 2015; Schmerer et al., 2014; Yu et al., 2014) , emphasizing the critical importance of estimating this parameter for any candidate vaccine.

Here, we develop a brief description of how this parameter might be estimated. Our example is deliberately simplified and should not be considered representative of most or even many ways that such an estimation protocol might be conducted. The example is offered to provide a sense of what is involved in the estimation.

A straightforward approach to estimating μ for a recombinant vector vaccine is through replicated serial passage experiments (Figure 7) . would be scored for their infection status (also observing individuals F I G U R E 5 The impact of the reversion rate to insert-free vector (y axis) and vector R 0,V (x axis) on the probability the pathogen is eradicated (top panels), the probability that the vaccine is self-extinguishing (middle panels), and the probability that the pathogen is eradicated, and the vaccine has been entirely degraded/extinguished (bottom panels) for the case where vector and vaccine have no cross-immunity. In the left-hand column, the vaccine is constructed from a vector characterized by recovery rate that is large relative to the reversion rate ( V = 0.05), whereas in the right-hand column, the vector has a recovery rate that is low relative to the reversion rate ( V = 0.002809). In all panels, the pathogen had an R 0,P = 2, with remaining background parameters given by: b = 10, d = 0.002809, and P = 0.05 . In this figure, the vaccine was administered to 30% of newborn animals over a period of 720 days that could affect the rate of reversion (e.g., when to manually transmit, how much virus to transmit, which tissues to inoculate).

The experimental design presented above allows the rate of reversion to be estimated using maximum likelihood. We note that any known design could be subjected to its own likelihood analysis, and ours was chosen merely to keep the illustration tractable. The data required are the fraction of the initial lines still carrying vaccine at each point in time, adjusted to exclude lines in which no transmission occurred. This approach neglects issues such as viral polymorphism (or dichotomizes it). The probability that a lineage carries intact vaccine at time t in the experiment is given by:

where t is the number of days since the experiment began and all individuals were infected with intact vaccine at the start of the experiment. The probability of observing X lineages with vaccine given a total of N replicate lineages is then given by: and the maximum likelihood solution for the reversion rate is given by: If PCR and/or serological tests have been conducted at multiple times over the course of the experiment, more refined estimates of the reversion rate can be achieved by integrating this information.

Specifically, if N 1 , N 2 , …, N τ lineages have been assayed on days, t 1 , t 2 , … t τ , and X 1 , X 2 , …, X τ of these lineages found to be PCR and/or serologically positive for the insert, the likelihood function is:

Numerically maximizing this function, or the equivalent log-likelihood expression, with respect to yields an estimate for the key reversion rate parameter of our model.

Our goal has been to investigate the consequences of transmissible vaccine instability on vaccine efficacy and fate. We have focused on recombinant vector vaccines, for which instability is the loss or downregulation of an antigenic insert and loss of vaccine function.

Our results demonstrate that recombinant vector vaccines can remain effective, even when genetically unstable and prone to reversion to insert-free vector. Not surprisingly, this is particularly true for transmissible vaccines that experience little cross-immunity with the viral vector. In addition to confirming the importance of cross-immunity, our results show that, all else being equal, those vaccines that transmit rapidly and generate only short-lived infection will be the most tolerant to genetic instability. Short-lived infections are favorable because they reduce opportunities for mutation and strain selection within hosts to increase the abundance of revertant vaccine.

In addition to demonstrating that transmissible vaccines can be robust to genetic instability, our results suggest that it may be possible to capitalize on this instability to develop self-extinguishing transmissible vaccines. The most immediate reason for ensuring that transmissible vaccines self-extinguish is for regulatory reasons. In this respect, deletion of the insert would restore the vaccine to the status of viral vector, an entity that may be naturally circulating in the wild population. Fortunately, deletion is also the process most easily controlled by engineering. For example, deletions are well known to occur between direct repeats in DNA (Cunningham et al., 1997; Springman, Molineux, Duong, Bull, & Bull, 2012) , and modification or elimination of those repeats reduces the spontaneous rate of deletions (Sleight, Bartley, Lieviant, & Sauro, 2010) . Genome engineering may also afford the possibility of designing accelerating rates of genomic instability. Our model assumed a constant rate of decay to revertant; a rate that accelerated over time would have many advantages. Such methods are not obviously ready for implementation, but they are at least suggested by "counting mechanisms" that have been employed in bacteria (reviewed in Bull & Barrick, 2017) .

Our study contributes to an increasing awareness of the many issues facing the choice of an appropriate vector for a transmissible recombinant vector vaccine. In most cases, the choice of a vector may be based on considerations that lie outside model results, and (6)

The impact of the reversion rate to insert-free vector (y axis) and vector R 0,P (x axis) on the probability the pathogen is eradicated (top panels), the probability that the vaccine is self-extinguishing (middle panels), and the probability that the pathogen is eradicated, and the vaccine has been entirely degraded/extinguished (bottom panels) for the case where vector and vaccine have complete cross-immunity. In the left-hand column, the vaccine is constructed from a vector characterized by a recovery rate that is large relative to the reversion rate ( V = 0.05), whereas in the right-hand column, the vector has a recovery rate that is low relative to the reversion rate ( V = 0.002809). In all panels, the pathogen had an R 0,P = 2, with remaining background parameters given by: b = 10, d = 0.002809, and P = 0.05 . In this figure, the vaccine was administered to 30% of newborn animals over a period of 720 days. Although the vaccine would ultimately be eliminated from the population, as predicted by our analytical results, for very low rate of reversion, this process is very slow and requires more time than our simulations were run the model results will then inform the possibilities and avenues for success. Consider cytomegalovirus (CMV), a vector currently being used in the development of transmissible vaccines targeting Ebola and Lassa fever (Marzi et al., 2016; Murphy et al., 2016; Tsuda et al., 2011 Tsuda et al., , 2015 . CMV appears promising because of its high level of species specificity, and ability to readily reinfect previously infected hosts , a property shown to be critically important here, and in a previous study . In contrast, our results also illuminate a potential weakness of CMV, which tends to generate long-lived infections that our results suggest increase opportunities for vaccine reversion. Ultimately, a rigorous assessment of vector quality will require robust estimates for key parameters such as transmission rates and rates of reversion to vector, something that should be a central focus of ongoing efforts to engineer successful transmissible vaccines.

As with any model, ours make simplifying assumptions that may influence our conclusions. One important assumption is that reversion from vaccine to vector has no impact on rates of transmission, despite the clear reasons to expect greater rates of vector transmission after the genetic cargo has been unloaded. Although we did not explicitly model this scenario, we expect doing so would support our current results. For instance, if the vector transmits more rapidly than the vaccine, we expect cross-immunity to become even more important and the efficacy of a transmissible vaccine to be further reduced. As a consequence, the difference between vectors that generate cross-immunity and those that do not would become even more significant. Another important assumption of our model is that it is possible to directly vaccinate newborn animals and that these newborn animals have the same rate of transmission as adult animals. If recovery from vaccine infection is rapid and newborn animals have limited encounters with the population at large, this could substantially limit the efficacy of a transmissible vaccine. Evaluating how this important assumption influences our results will require more complex models that explicitly integrate age structure. Finally, our modeling framework focuses on a single well-mixed population and thus ignores potential impacts of spatial heterogeneity and population structure (Varrelman, Basinski, Remien, & Nuismer, 2019) .

As a whole, our results suggest transmissible vaccines may be relatively robust to modest levels of genetic instability, particularly in cases where cross-immunity between vector and vaccine is weak. It may also prove possible to capitalize on genetic instability to develop transmissible vaccines that are both effective and self-extinguishing.

Identifying the quantitative impact transmissible vaccines are likely to have on real pathogen populations, however, will require considerable progress in estimating key parameters such as vector transmission rates, recovery rates, and rates of reversion to insert-free vector. Only when reliable estimates become available for these parameters will it become possible to critically evaluate the promise of various vectors and the scope for designing effective and self-extinguishing vaccines.

Funding for this work was provided by NIH R01GM122079 to SLN.

We thank Michael Jarvis for helpful comments on this work.

None declared.

Data sharing not applicable-no new data generated. All simulation code is available upon request. 

We identify the vaccination rate that, in the absence of any crossimmunity between vector and vaccine, preempts the invasion of a pathogen. If the host population is at steady state and a small number of pathogen-infected hosts are introduced into the population, the pathogen is unable to invade the population when

Here, S * is the steady state density of susceptible hosts, and P is the initial (small) density of pathogen-infected hosts. Because p > 0, inequality (A1-1) can be rearranged as

where R 0,p is the pathogen's basic reproduction number. Inequality A1-2 shows that the host population is protected from pathogen invasion so long as the density of susceptible hosts at steady state is below the threshold b d 1 R 0,p . Next, we find the proportion σ of newborns that must be vaccinated to maintain the susceptible population at this prophylactic threshold.

In the absence of the pathogen, the densities of susceptible, vaccine-infected, and vaccine-recovered hosts are described by a subset of system (1) with the P state variable set to 0, and the equations for pathogen-infected and pathogen-recovered hosts (P,R P ) omitted.

The submodel consists of Equations 1a, 1b, and 1d:

In the supplementary Mathematica file (available on request), we

show that system A1-3 possesses a single biologically relevant (i.e., positive state variables) equilibrium that is locally stable.

As a result, if continual vaccination is carried out in perpetuity, the densities of susceptible, vaccine-infected, and vaccine-recovered hosts each approach steady state values described by this equilibrium. In particular, the density of susceptible individuals approaches which was found by setting the derivatives in system A1-3 to 0 and solving for S, V, and R V . The steady state values V * and R V * are contained in the Mathematica file, but are omitted here.

We derive the threshold vaccination proportion by plugging (A1-4) into the threshold condition for prophylaxis, inequality A1-2, and solving for σ. The critical proportion is APPENDIX 2

Cross-immunity between the vector and vaccine will hinder the spread of the vaccine. Consequently, a greater proportion of the newborns will have to be vaccinated to prevent pathogen invasion. In Appendix The total rate at which an invading pathogen spreads is described by the summed rates of spread in the pathogen-infected classes. In the Mathematica supplementary, we show that, if the population is at steady state, the total rate at which an invading pathogen spreads is where P = P S + P M + C and S * denotes the steady state level of S . The pathogen will not be able to invade the population so long as equation A2-2 is <0 for an initial (small) P , which is the case so long as Inequality A2-3 implies that prophylaxis is achieved when the sum of the steady state densities of S, M S , and R M is below b d 1 R 0,P . Next, we find the values of these state variables at steady state and determine the proportion of newborns that must be vaccinated to preemptively bring the density of S * below threshold A2-3.

In the absence of the pathogen, the levels of vector-infected, vaccine-infected, and susceptible hosts are described by a simplified version of System 2 with all instances of pathogen-infected and pathogen-recovered state variables set to 0 (P = 0, C = 0, R P = 0, M P = 0), and their associated equations omitted. Because the recovered classes R V and R M do not influence the steady state level of S * , they can be omitted from the analysis that follows. The reduced system is comprised of Equations 2a, 2b, 2c, and 2i of the main text:

(A1-1) P = p S * P − p + d P < 0.

The steady states of System A2 are obtained by setting the lefthand sides to 0, and solving the resulting system of equations for S, M S , and R M . The expressions that result are omitted here due to their complexity but presented in the supplementary Mathematica file.

We use the resulting steady states to calculate the effective density of susceptible hosts, S * , and then solve for the critical value of σ that brings S * below the threshold that prevents pathogen invasion, inequality A2-3. The formula that results is complex and omitted here, but displayed in the Mathematica file.

After using a transmissible vaccine to block an invading pathogen, it may be desirable for the vaccine infection process to self-extinguish once direct vaccination stops. We show that this behavior is possible if the rate at which the vaccine reverts into the vector is sufficiently large. Upon cessation of direct vaccination, the levels of susceptible, vaccine-infected, and vaccine-recovered hosts are described by

Equations 1a, 1b, and 1d of the main text, with the P state variable and vaccination parameter σ both set to 0:

To be self-extinguishing, we require that any solution with nonnegative initial conditions approaches the equilibrium for which only susceptible individuals are present. In the supplementary Mathematica file, we show that this occurs when

If inequality A3-2 is satisfied, the susceptible-only equilibrium is the only biologically relevant equilibrium (i.e., all state variables are non-negative) and, in addition, is locally stable if condition A3-2 is satisfied, but not otherwise. As a result, condition A3-2 is necessary for solutions of system A3-1 to approach an equilibrium with zero vaccine-infected hosts.

We use the relationship R 0,V = V b d(d+ V ) to rewrite condition A3-2 in terms of the vector R 0,V . Condition A3-2 can be rearranged to describe a threshold rate of reversion for a self-extinguishing vaccine:

We check for conditions for which, following pathogen eradication, the levels of vaccine-infected hosts fall to zero when vaccination ceases. We use a subset of System 2, which describes the scenario in which the vaccine and vector share cross-immunity. Because the pathogen is not present in this scenario, we only work with Equation 

Towards a unique and transmissible vaccine against myxomatosis and rabbit haemorrhagic disease for rabbit populations

Horizontal transmissible protection against myxomatosis and rabbit hemorrhagic disease by using a recombinant myxoma virus

A little goes a long way: Weak vaccine transmission facilitates oral vaccination campaigns against zoonotic pathogens

Evaluating the promise of recombinant transmissible vaccines

Arresting evolution

Transmissible viral vaccines

Parallel molecular evolution of deletions and nonsense mutations in bacteriophage T7

SARS and MERS: Recent insights into emerging coronaviruses

Transgene stability for three replication-competent murine leukemia virus vectors

Human cytomegalovirus: Clinical aspects, immune regulation, and emerging treatments

Exact stochastic simulation of coupled chemicalreactions

Understanding Ebola: The 2014 epidemic

In vitro and in vivo genetic stability studies of a human adenovirus type 5 recombinant rabies glycoprotein vaccine (ONRAB)

Genomic stability of murine leukemia viruses containing insertions at the Env-3 ' untranslated region boundary

A highly immunogenic and protective middle east respiratory syndrome coronavirus vaccine based on a recombinant measles virus vaccine platform

Cytomegalovirus-based vaccine expressing Ebola virus glycoprotein protects nonhuman primates from Ebola virus infection

Self-disseminating vaccines for emerging infectious diseases

Eradicating infectious disease using weakly transmissible vaccines

Controlling epidemics with transmissible vaccines

Pathways to zoonotic spillover

Lassa fever: Epidemiology, clinical features, and social consequences

Challenges in predicting the evolutionary maintenance of a phage transgene

Designing and engineering evolutionary robust genetic circuits

Evolutionary stability of a refactored phage genome

Global epidemiology of canine rabies: Past, present, and future prospects

A replicating cytomegalovirus-based vaccine encoding a single Ebola virus nucleoprotein CTL epitope confers protection against Ebola virus

A cytomegalovirus-based vaccine provides long-lasting protection against lethal Ebola virus challenge after a single dose

Transmissible vaccines in heterogeneous populations: Implications for vaccine design

Construction and in vitro evaluation of a recombinant live attenuated PRRSV expressing GM-CSF

Evolution and containment of transmissible recombinant vector vaccines