Recombinant transmissible vaccines will be intrinsically contained despite the ability to superinfect.

INTRODUCTION: Transmissible vaccines offer a novel approach to suppressing viruses in wildlife populations, with possible applications against viruses that infect humans as zoonoses - Lassa, Ebola, rabies. To ensure safety, current designs propose a recombinant vector platform in which the vector is isolated from the target wildlife population. Because using an endemic vector creates the potential for preexisting immunity to block vaccine transmission, these designs focus on vector viruses capable of superinfection, spreading throughout the host population following vaccination of few individuals. AREAS COVERED: We present original theoretical arguments that, regardless of its R0 value, a recombinant vaccine using a superinfecting vector is not expected to expand its active infection coverage when released into a wildlife population that already carries the vector. However, if superinfection occurs at a high rate such that individuals are repeatedly infected throughout their lives, the immunity footprint in the population can be high despite a low incidence of active vaccine infections. Yet we provide reasons that the above expectation is optimistic. EXPERT OPINION: High vaccine coverage will typically require repeated releases or release into a population lacking the vector, but careful attention to vector choice and vaccine engineering should also help improve transmissible vaccine utility.

Predicting the fine-scale spatial distribution of zoonotic reservoirs using computer vision.

Zoonotic diseases threaten human health worldwide and are often associated with anthropogenic disturbance. Predicting how disturbance influences spillover risk is critical for effective disease intervention but difficult to achieve at fine spatial scales. Here, we develop a method that learns the spatial distribution of a reservoir species from aerial imagery. Our approach uses neural networks to extract features of known or hypothesized importance from images. The spatial distribution of these features is then summarized and linked to spatially explicit reservoir presence/absence data using boosted regression trees. We demonstrate the utility of our method by applying it to the reservoir of Lassa virus, Mastomys natalensis, within the West African nations of Sierra Leone and Guinea. We show that, when trained using reservoir trapping data and publicly available aerial imagery, our framework learns relationships between environmental features and reservoir occurrence and accurately ranks areas according to the likelihood of reservoir presence.

Optimizing the delivery of self-disseminating vaccines in fluctuating wildlife populations.

Zoonotic pathogens spread by wildlife continue to spill into human populations and threaten human lives. A potential way to reduce this threat is by vaccinating wildlife species that harbor pathogens that are infectious to humans. Unfortunately, even in cases where vaccines can be distributed en masse as edible baits, achieving levels of vaccine coverage sufficient for pathogen elimination is rare. Developing vaccines that self-disseminate may help solve this problem by magnifying the impact of limited direct vaccination. Although models exist that quantify how well these self-disseminating vaccines will work when introduced into temporally stable wildlife populations, how well they will perform when introduced into populations with pronounced seasonal population dynamics remains unknown. Here we develop and analyze mathematical models of fluctuating wildlife populations that allow us to study how reservoir ecology, vaccine design, and vaccine delivery interact to influence vaccine coverage and opportunities for pathogen elimination. Our results demonstrate that the timing of vaccine delivery can make or break the success of vaccination programs. As a general rule, the effectiveness of self-disseminating vaccines is optimized by introducing after the peak of seasonal reproduction when the number of susceptible animals is near its maximum.

How Interactions during Viral-Viral Coinfection Can Shape Infection Kinetics.

Respiratory viral infections are a leading global cause of disease with multiple viruses detected in 20-30% of cases, and several viruses simultaneously circulating. Some infections with unique viral copathogens result in reduced pathogenicity, while other viral pairings can worsen disease. The mechanisms driving these dichotomous outcomes are likely variable and have only begun to be examined in the laboratory and clinic. To better understand viral-viral coinfections and predict potential mechanisms that result in distinct disease outcomes, we first systematically fit mathematical models to viral load data from ferrets infected with respiratory syncytial virus (RSV), followed by influenza A virus (IAV) after 3 days. The results suggest that IAV reduced the rate of RSV production, while RSV reduced the rate of IAV infected cell clearance. We then explored the realm of possible dynamics for scenarios that had not been examined experimentally, including a different infection order, coinfection timing, interaction mechanisms, and viral pairings. IAV coinfection with rhinovirus (RV) or SARS-CoV-2 (CoV2) was examined by using human viral load data from single infections together with murine weight-loss data from IAV-RV, RV-IAV, and IAV-CoV2 coinfections to guide the interpretation of the model results. Similar to the results with RSV-IAV coinfection, this analysis shows that the increased disease severity observed during murine IAV-RV or IAV-CoV2 coinfection was likely due to the slower clearance of IAV-infected cells by the other viruses. The improved outcome when IAV followed RV, on the other hand, could be replicated when the rate of RV infected cell clearance was reduced by IAV. Simulating viral-viral coinfections in this way provides new insights about how viral-viral interactions can regulate disease severity during coinfection and yields testable hypotheses ripe for experimental evaluation.

How Interactions During Viral-Viral Coinfection Can Shape Infection Kinetics.

Respiratory virus infections are a leading cause of disease worldwide with multiple viruses detected in 20-30% of cases and several viruses simultaneously circulating. Some infections with viral copathogens have been shown to result in reduced pathogenicity while other virus pairings can worsen disease. The mechanisms driving these dichotomous outcomes are likely variable and have only begun to be examined in the laboratory and clinic. To better understand viral-viral coinfections and predict potential mechanisms that result in distinct disease outcomes, we first systematically fit mathematical models to viral load data from ferrets infected with respiratory syncytial virus (RSV) followed by influenza A virus (IAV) after 3 days. The results suggested that IAV reduced the rate of RSV production while RSV reduced the rate of IAV infected cell clearance. We then explored the realm of possible dynamics for scenarios not examined experimentally, including different infection order, coinfection timing, interaction mechanisms, and viral pairings. IAV coinfection with rhinovirus (RV) or SARS-CoV-2 (CoV2) was examined by using human viral load data from single infections together with murine weight loss data from IAV-RV, RV-IAV, and IAV-CoV2 coinfections to guide the interpretation of the model results. Similar to the results with RSV-IAV coinfection, this analysis showed that the increased disease severity observed during murine IAV-RV or IAV-CoV2 coinfection was likely due to slower clearance of IAV infected cells by the other viruses. On the contrary, the improved outcome when IAV followed RV could be replicated when the rate of RV infected cell clearance was reduced by IAV. Simulating viral-viral coinfections in this way provides new insights about how viral-viral interactions can regulate disease severity during coinfection and yields testable hypotheses ripe for experimental evaluation.

Correction: Microbial phenotypic heterogeneity in response to a metabolic toxin: Continuous, dynamically shifting distribution of formaldehyde tolerance in Methylobacterium extorquens populations.

[This corrects the article DOI: 10.1371/journal.pgen.1008458.].

Identifying the genetic basis of viral spillover using Lassa virus as a test case.

The rate at which zoonotic viruses spill over into the human population varies significantly over space and time. Remarkably, we do not yet know how much of this variation is attributable to genetic variation within viral populations. This gap in understanding arises because we lack methods of genetic analysis that can be easily applied to zoonotic viruses, where the number of available viral sequences is often limited, and opportunistic sampling introduces significant population stratification. Here, we explore the feasibility of using patterns of shared ancestry to correct for population stratification, enabling genome-wide association methods to identify genetic substitutions associated with spillover into the human population. Using a combination of phylogenetically structured simulations and Lassa virus sequences collected from humans and rodents in Sierra Leone, we demonstrate that existing methods do not fully correct for stratification, leading to elevated error rates. We also demonstrate, however, that the Type I error rate can be substantially reduced by confining the analysis to a less-stratified region of the phylogeny, even in an already-small dataset. Using this method, we detect two candidate single-nucleotide polymorphisms associated with spillover in the Lassa virus polymerase gene and provide generalized recommendations for the collection and analysis of zoonotic viruses.

Reservoir population ecology, viral evolution and the risk of emerging infectious disease.

The ecology and life history of wild animals influences their potential to harbour infectious disease. This observation has motivated studies identifying empirical relationships between traits of wild animals and historical patterns of spillover and emergence into humans. Although these studies have identified compelling broad-scale patterns, they are generally agnostic with respect to underlying mechanisms. Here, we develop mathematical models that couple reservoir population ecology with viral epidemiology and evolution to clarify existing verbal arguments and pinpoint the conditions that favour spillover and emergence. Our results support the idea that average lifespan influences the likelihood of an animal serving as a reservoir for human infectious disease. At the same time, however, our results show that the magnitude of this effect is sensitive to the rate of viral mutation. Our results also demonstrate that viral pathogens causing persistent infections or a transient immune response within the reservoir are more likely to fuel emergence. Genetically explicit stochastic simulations enrich these mathematical results by identifying relationships between the genetic basis of transmission and the risk of spillover and emergence. Together, our results clarify the scope of applicability for existing hypotheses and refine our understanding of emergence risk.

Quantifying the effectiveness of betaherpesvirus-vectored transmissible vaccines.

Transmissible vaccines have the potential to revolutionize how zoonotic pathogens are controlled within wildlife reservoirs. A key challenge that must be overcome is identifying viral vectors that can rapidly spread immunity through a reservoir population. Because they are broadly distributed taxonomically, species specific, and stable to genetic manipulation, betaherpesviruses are leading candidates for use as transmissible vaccine vectors. Here we evaluate the likely effectiveness of betaherpesvirus-vectored transmissible vaccines by developing and parameterizing a mathematical model using data from captive and free-living mouse populations infected with murine cytomegalovirus (MCMV). Simulations of our parameterized model demonstrate rapid and effective control for a range of pathogens, with pathogen elimination frequently occurring within a year of vaccine introduction. Our results also suggest, however, that the effectiveness of transmissible vaccines may vary across reservoir populations and with respect to the specific vector strain used to construct the vaccine.

Structural identifiability of the generalized Lotka-Volterra model for microbiome studies.

Population dynamic models can be used in conjunction with time series of species abundances to infer interactions. Understanding microbial interactions is a prerequisite for numerous goals in microbiome research, including predicting how populations change over time, determining how manipulations of microbiomes affect dynamics and designing synthetic microbiomes to perform tasks. As such, there is great interest in adapting population dynamic theory for microbial systems. Despite the appeal, numerous hurdles exist. One hurdle is that the data commonly obtained from DNA sequencing yield estimates of relative abundances, while population dynamic models such as the generalized Lotka-Volterra model track absolute abundances or densities. It is not clear whether relative abundance data alone can be used to infer parameters of population dynamic models such as the Lotka-Volterra model. We used structural identifiability analyses to determine the extent to which a time series of relative abundances can be used to parametrize the generalized Lotka-Volterra model. We found that only with absolute abundance data to accompany relative abundance estimates from sequencing can all parameters be uniquely identified. However, relative abundance data alone do contain information on relative interaction strengths, which is sufficient for many studies where the goal is to estimate key interactions and their effects on dynamics. Using synthetic data of a simple community for which we know the underlying structure, local practical identifiability analysis showed that modest amounts of both process and measurement error do not fundamentally affect these identifiability properties.

Suppressing evolution in genetically engineered systems through repeated supplementation.

Genetically engineered organisms are prone to evolve in response to the engineering. This evolution is often undesirable and can negatively affect the purpose of the engineering. Methods that maintain the stability of engineered genomes are therefore critical to the successful design and use of genetically engineered organisms. One potential method to limit unwanted evolution is by taking advantage of the ability of gene flow to counter local adaption, a process of supplementation. Here, we investigate the feasibility of supplementation as a mechanism to offset the evolutionary degradation of a transgene in three model systems: a bioreactor, a gene drive, and a transmissible vaccine. In each model, continual introduction from a stock is used to balance mutation and selection against the transgene. Each system has its unique features. The bioreactor system is especially tractable and has a simple answer: The level of supplementation required to maintain the transgene at a frequency p ^ is approximately p ^ s , where s is the selective disadvantage of the transgene. Supplementation is also feasible in the transmissible vaccine case but is probably not practical to prevent the evolution of resistance against a gene drive. We note, however, that the continual replacement of even a small fraction of a large population can be challenging, limiting the usefulness of supplementation as a means of controlling unwanted evolution.

Vector dynamics influence spatially imperfect genetic interventions against disease.

BACKGROUND AND OBJECTIVES: Genetic engineering and similar technologies offer promising new approaches to controlling human diseases by blocking transmission from vectors. However, in spatially structured populations, imperfect coverage of the vector will leave pockets in which the parasite may persist. Movement by humans may disrupt this local persistence and facilitate eradication when these pockets are small, spreading parasite reproduction outside unprotected areas and into areas that block its reproduction. Here, we consider the sensitivity of this process to biological details: do simple generalities emerge that may facilitate interventions? METHODOLOGY: We develop formal mathematical models of this process similar to standard Ross-Macdonald models, but (i) specifying spatial structure of two patches, with vector transmission blocked in one patch but not in the other, (ii) allowing temporary human movement (travel instead of migration) and (iii) considering two different modes of mosquito biting. RESULTS: We find that there is no invariant effect of disrupting spatial structure with travel. For both biting models, travel out of the unprotected patch has different consequences than travel by visitors into the patch, but the effects are reversed between the two biting models. CONCLUSIONS AND IMPLICATIONS: Overall, the effect of human travel on the maintenance of vector-borne diseases in structured habitats must be considered in light of the actual biology of mosquito abundances, biting dynamics and human movement patterns. Lay summary: Genetic interventions against pathogens transmitted by insect vectors are promising methods of controlling infectious diseases. These interventions may be imperfect, leaving pockets where the parasite persists. How will human movement between protected and unprotected areas affect persistence? Mathematical models developed here show that the answer is ecology-dependent, depending on vector biting behavior.

Bridging the gap: Using reservoir ecology and human serosurveys to estimate Lassa virus spillover in West Africa.

Forecasting the risk of pathogen spillover from reservoir populations of wild or domestic animals is essential for the effective deployment of interventions such as wildlife vaccination or culling. Due to the sporadic nature of spillover events and limited availability of data, developing and validating robust, spatially explicit, predictions is challenging. Recent efforts have begun to make progress in this direction by capitalizing on machine learning methodologies. An important weakness of existing approaches, however, is that they generally rely on combining human and reservoir infection data during the training process and thus conflate risk attributable to the prevalence of the pathogen in the reservoir population with the risk attributed to the realized rate of spillover into the human population. Because effective planning of interventions requires that these components of risk be disentangled, we developed a multi-layer machine learning framework that separates these processes. Our approach begins by training models to predict the geographic range of the primary reservoir and the subset of this range in which the pathogen occurs. The spillover risk predicted by the product of these reservoir specific models is then fit to data on realized patterns of historical spillover into the human population. The result is a geographically specific spillover risk forecast that can be easily decomposed and used to guide effective intervention. Applying our method to Lassa virus, a zoonotic pathogen that regularly spills over into the human population across West Africa, results in a model that explains a modest but statistically significant portion of geographic variation in historical patterns of spillover. When combined with a mechanistic mathematical model of infection dynamics, our spillover risk model predicts that 897,700 humans are infected by Lassa virus each year across West Africa, with Nigeria accounting for more than half of these human infections.

The role of alcohol consumption on acetaminophen induced liver injury: Implications from a mathematical model.

Acetaminophen (APAP) overdose is one of the predominant causes of drug induced acute liver injury in the U.S and U.K. Clinical studies show that ingestion of alcohol may increase the risk of APAP induced liver injury. Chronic alcoholism may potentiate APAP hepatotoxicity and this increased risk of APAP toxicity is observed when APAP is ingested even shortly after alcohol is cleared from the body. However, clinical reports also suggest that acute alcohol consumption may have a protective effect against hepatotoxicity by inhibiting microsomal acetaminophen oxidation and thereby reducing N-acetyl-p-benzoquinone imine (NAPQI) production. The aim of this study is to model this dual role of alcohol to determine how the timing of alcohol ingestion affects APAP metabolism and resulting liver injury and identify mechanisms of APAP induced liver injury. The mathematical model is developed to capture condition of a patient of single time APAP overdose who may be an acute or chronic alcohol user. The analysis suggests that the risk of APAP-induced hepatotoxicity is increased if APAP is ingested shortly after alcohol is cleared from the body in chronic alcohol users. A protective effect of acute consumption of alcohol is also observed in patients with APAP overdose. For example, simultaneous ingestion of alcohol and APAP overdose or alcohol intake after or before few hours of APAP overdose may result in less APAP-induced hepatotoxicity when compared to a single time APAP overdose. The rate of hepatocyte damage in APAP overdose patients depends on trade-off between induction and inhibition of CYP enzyme.

Single-Inclusion Kinetics of Chlamydia trachomatis Development.

The obligate intracellular bacterial pathogen Chlamydia trachomatis is reliant on a developmental cycle consisting of two cell forms, termed the elementary body (EB) and the reticulate body (RB). The EB is infectious and utilizes a type III secretion system and preformed effector proteins during invasion, but it does not replicate. The RB replicates in the host cell but is noninfectious. This developmental cycle is central to chlamydial pathogenesis. In this study, we developed mathematical models of the developmental cycle that account for potential factors influencing RB-to-EB cell type switching during infection. Our models predicted that two categories of regulatory signals for RB-to-EB development could be differentiated experimentally, an "intrinsic" cell-autonomous program inherent to each RB and an "extrinsic" environmental signal to which RBs respond. To experimentally differentiate between mechanisms, we tracked the expression of C. trachomatis development-specific promoters in individual inclusions using fluorescent reporters and live-cell imaging. These experiments indicated that EB production was not influenced by increased multiplicity of infection or by superinfection, suggesting the cycle follows an intrinsic program that is not directly controlled by environmental factors. Additionally, live-cell imaging revealed that EB development is a multistep process linked to RB growth rate and cell division. The formation of EBs followed a progression with expression from the euo and ihtA promoters evident in RBs, while expression from the promoter for hctA was apparent in early EBs/IBs. Finally, expression from the promoters for the true late genes, hctB, scc2, and tarp, was evident in the maturing EB.IMPORTANCE Chlamydia trachomatis is an obligate intracellular bacterium that can cause trachoma, cervicitis, urethritis, salpingitis, and pelvic inflammatory disease. To establish infection in host cells, Chlamydia must complete a multiple-cell-type developmental cycle. The developmental cycle consists of specialized cells, the EB cell, which mediates infection of new host cells, and the RB cell, which replicates and eventually produces more EB cells to mediate the next round of infection. By developing and testing mathematical models to discriminate between two competing hypotheses for the nature of the signal controlling RB-to-EB cell type switching, we demonstrate that RB-to-EB development follows a cell-autonomous program that does not respond to environmental cues. Additionally, we show that RB-to-EB development is a function of chlamydial growth and division. This study serves to further our understanding of the chlamydial developmental cycle that is central to the bacterium's pathogenesis.

Bayesian estimation of Lassa virus epidemiological parameters: Implications for spillover prevention using wildlife vaccination.

Lassa virus is a significant burden on human health throughout its endemic region in West Africa, with most human infections the result of spillover from the primary rodent reservoir of the virus, the natal multimammate mouse, M. natalensis. Here we develop a Bayesian methodology for estimating epidemiological parameters of Lassa virus within its rodent reservoir and for generating probabilistic predictions for the efficacy of rodent vaccination programs. Our approach uses Approximate Bayesian Computation (ABC) to integrate mechanistic mathematical models, remotely-sensed precipitation data, and Lassa virus surveillance data from rodent populations. Using simulated data, we show that our method accurately estimates key model parameters, even when surveillance data are available from only a relatively small number of points in space and time. Applying our method to previously published data from two villages in Guinea estimates the time-averaged R0 of Lassa virus to be 1.74 and 1.54 for rodent populations in the villages of Bantou and Tanganya, respectively. Using the posterior distribution for model parameters derived from these Guinean populations, we evaluate the likely efficacy of vaccination programs relying on distribution of vaccine-laced baits. Our results demonstrate that effective and durable reductions in the risk of Lassa virus spillover into the human population will require repeated distribution of large quantities of vaccine.

Microbial phenotypic heterogeneity in response to a metabolic toxin: Continuous, dynamically shifting distribution of formaldehyde tolerance in Methylobacterium extorquens populations.

While microbiologists often make the simplifying assumption that genotype determines phenotype in a given environment, it is becoming increasingly apparent that phenotypic heterogeneity (in which one genotype generates multiple phenotypes simultaneously even in a uniform environment) is common in many microbial populations. The importance of phenotypic heterogeneity has been demonstrated in a number of model systems involving binary phenotypic states (e.g., growth/non-growth); however, less is known about systems involving phenotype distributions that are continuous across an environmental gradient, and how those distributions change when the environment changes. Here, we describe a novel instance of phenotypic diversity in tolerance to a metabolic toxin within wild-type populations of Methylobacterium extorquens, a ubiquitous phyllosphere methylotroph capable of growing on the methanol periodically released from plant leaves. The first intermediate in methanol metabolism is formaldehyde, a potent cellular toxin that is lethal in high concentrations. We have found that at moderate concentrations, formaldehyde tolerance in M. extorquens is heterogeneous, with a cell's minimum tolerance level ranging between 0 mM and 8 mM. Tolerant cells have a distinct gene expression profile from non-tolerant cells. This form of heterogeneity is continuous in terms of threshold (the formaldehyde concentration where growth ceases), yet binary in outcome (at a given formaldehyde concentration, cells either grow normally or die, with no intermediate phenotype), and it is not associated with any detectable genetic mutations. Moreover, tolerance distributions within the population are dynamic, changing over time in response to growth conditions. We characterized this phenomenon using bulk liquid culture experiments, colony growth tracking, flow cytometry, single-cell time-lapse microscopy, transcriptomics, and genome resequencing. Finally, we used mathematical modeling to better understand the processes by which cells change phenotype, and found evidence for both stochastic, bidirectional phenotypic diversification and responsive, directed phenotypic shifts, depending on the growth substrate and the presence of toxin.

Spatial structure undermines parasite suppression by gene drive cargo.

Gene drives may be used in two ways to curtail vectored diseases. Both involve engineering the drive to spread in the vector population. One approach uses the drive to directly depress vector numbers, possibly to extinction. The other approach leaves intact the vector population but suppresses the disease agent during its interaction with the vector. This second application may use a drive engineered to carry a genetic cargo that blocks the disease agent. An advantage of the second application is that it is far less likely to select vector resistance to block the drive, but the disease agent may instead evolve resistance to the inhibitory cargo. However, some gene drives are expected to spread so fast and attain such high coverage in the vector population that, if the disease agent can evolve resistance only gradually, disease eradication may be feasible. Here we use simple models to show that spatial structure in the vector population can greatly facilitate persistence and evolution of resistance by the disease agent. We suggest simple approaches to avoid some types of spatial structure, but others may be intrinsic to the populations being challenged and difficult to overcome.

Gene-drive-mediated extinction is thwarted by population structure and evolution of sib mating.

BACKGROUND AND OBJECTIVES: Genetic engineering combined with CRISPR technology has developed to the point that gene drives can, in theory, be engineered to cause extinction in countless species. Success of extinction programs now rests on the possibility of resistance evolution, which is largely unknown. Depending on the gene-drive technology, resistance may take many forms, from mutations in the nuclease target sequence (e.g. for CRISPR) to specific types of non-random population structures that limit the drive (that may block potentially any gene-drive technology). METHODOLOGY: We develop mathematical models of various deviations from random mating to consider escapes from extinction-causing gene drives. A main emphasis here is sib mating in the face of recessive-lethal and Y-chromosome drives. RESULTS: Sib mating easily evolves in response to both kinds of gene drives and maintains mean fitness above 0, with equilibrium fitness depending on the level of inbreeding depression. Environmental determination of sib mating (as might stem from population density crashes) can also maintain mean fitness above 0. A version of Maynard Smith's haystack model shows that pre-existing population structure can enable drive-free subpopulations to be maintained against gene drives. CONCLUSIONS AND IMPLICATIONS: Translation of mean fitness into population size depends on ecological details, so understanding mean fitness evolution and dynamics is merely the first step in predicting extinction. Nonetheless, these results point to possible escapes from gene-drive-mediated extinctions that lie beyond the control of genome engineering. LAY SUMMARY: Recent gene drive technologies promise to suppress and even eradicate pests and disease vectors. Simple models of gene-drive evolution in structured populations show that extinction-causing gene drives can be thwarted both through the evolution of sib mating as well as from purely demographic processes that cluster drive-free individuals.

Transmissible vaccines in heterogeneous populations: Implications for vaccine design.

Transmissible vaccines may provide a promising solution for improving the control of infectious disease, particularly zoonotic pathogens with wildlife reservoirs. Although it is well known that heterogeneity in pathogen transmission impacts the spread of infectious disease, the effects of heterogeneity on vaccine transmission are largely unknown. Here we develop and analyze a mathematical model that quantifies the potential benefits of a transmissible vaccine in a population where transmission is heterogeneous between two subgroups. Our results demonstrate that the effect of heterogeneity on the benefit of vaccine transmission largely depends on the vaccine design and the pattern of vaccine administration across subgroups. Specifically, our results show that in most cases a transmissible vaccine designed to mirror the transmission of the pathogen is optimal. If the vaccination effort can be preferentially biased towards a given subgroup, a vaccine with a pattern of transmission opposite to that of the pathogen can become optimal in some cases. To better understand the consequences of heterogeneity on the effectiveness of a transmissible vaccine in the real world, we parameterized our model using data from Sin Nombre virus in deer mice (Peromyscus maniculatus). The results of this analysis reveal that when a vaccination campaign is limited in vaccine availability, a traditional vaccine must be administered primarily to males for the spread of Sin Nombre virus to be prevented. In contrast, a transmissible vaccine remains effective even when it cannot be preferentially administered to males.