Proliferation in malaria parasites: how resource limitation can prevent evolution of greater virulence

For parasites, robust proliferation within hosts is crucial for establishing the infection and creating opportu- nities for onward transmission. While faster proliferation enhances transmission rates, it is often assumed to curtail transmission duration by killing the host (virulence), a tradeoff constraining parasite evolution. Yet in many diseases, including malaria, the preponderance of infections with mild or absent symptoms suggests that host mortality is not a sufficient constraint, raising the question of what restrains evolution towards faster proliferation. In malaria infections, the maximum rate of proliferation is determined by the burst size, the number of daughter parasites produced per infected red blood cell. Larger burst sizes should expand the pool of infected red blood cells that can be used to produce the specialized transmission forms needed to infect mosquitoes. We use a within-host model parameterized for rodent malaria parasites (Plasmodium chabaudi ) to project the transmission consequences of burst size, focusing on initial acute infection where re- source limitation and risk of host mortality are greatest. We find that resource limitation restricts evolution towards higher burst sizes below the level predicted by host mortality alone. Our results suggest resource lim- itation could represent a more general constraint than virulence-transmission tradeoffs, preventing evolution towards faster proliferation.

Within-host dynamics of SARS-CoV-2 infection: A systematic review and meta-analysis

Within-host model specified by viral dynamic parameters is a mainstream tool to understand SARS-CoV-2 replication cycle in infected patients. The parameter uncertainty further affects the output of the model, such as the efficacy of potential antiviral drugs. However, gathering empirical data on these parameters is challenging. Here, we aim to conduct a systematic review of viral dynamic parameters used in within-host models by calibrating the model to the viral load data measured from upper respiratory specimens. We searched the PubMed, Embase and Web of Science databases (between 1 December 2019 and 10 February 2022) for within-host modelling studies. We identified seven independent within-host models from the above nine studies, including Type I interferon, innate response, humoral immune response or cell-mediated immune response. From these models, we extracted and analyse seven widely used viral dynamic parameters including the viral load at the point of infection or symptom onset, the rate of viral particles infecting susceptible cells, the rate of infected cells releasing virus, the rate of virus particles cleared, the rate of infected cells cleared and the rate of cells in the eclipse phase can become productively infected. We identified seven independent within-host models from nine eligible studies. The viral load at symptom onset is 4.78 (95% CI:2.93, 6.62) log(copies/ml), and the viral load at the point of infection is -1.00 (95% CI:-1.94, -0.05) log(copies/ml). The rate of viral particles infecting susceptible cells and the rate of infected cells cleared have the pooled estimates as -6.96 (95% CI:-7.66, -6.25) log([copies/ml]^-1 day^-1 ) and 0.92 (95% CI:-0.09, 1.93) day^-1 , respectively. We found that the rate of infected cells cleared was associated with the reported model in the meta-analysis by including the model type as a categorical variable (p < .01). Joint viral dynamic parameters estimates when parameterizing within-host models have been published for SARS-CoV-2. The reviewed viral dynamic parameters can be used in the same within-host model to understand SARS-CoV-2 replication cycle in infected patients and assess the impact of pharmaceutical interventions.

Linking functional and molecular mechanisms of host resilience to malaria infection

It remains challenging to understand why some hosts suffer severe illnesses, while others are unscathed by the same infection. We fitted a mathematical model to longitudinal measurements of parasite and red blood cell density in murine hosts from diverse genetic backgrounds to identify aspects of within-host interactions that explain variation in host resilience and survival during acute malaria infection. Among eight mouse strains that collectively span 90% of the common genetic diversity of laboratory mice, we found that high host mortality was associated with either weak parasite clearance, or a strong, yet imprecise response that inadvertently removes uninfected cells in excess. Subsequent cross-sectional cytokine assays revealed that the two distinct functional mechanisms of poor survival were underpinned by low expression of either pro- or anti-inflammatory cytokines, respectively. By combining mathematical modelling and molecular immunology assays, our study uncovered proximate mechanisms of diverse infection outcomes across multiple host strains and biological scales.

A generalized within-host model of dengue infection with a non-constant monocyte production rate

In this paper, we generalize a previous model of within-host dengue infection with a nonconstant monocyte production rate. We establish the existence of three equilibria and give some local stability results. We then estimate three parameters in the model from clinical data for dengue virus serotype 2. It is then shown that the model can exhibit behaviours that are not possible under the assumption of constant monocyte production. Lastly, we perform a sensitivity analysis of the model in two contexts, antiviral treatment and immunostimulatory treatment. The results predict that antiviral treatments that reduce the viral replication rate in infected monocytes are the most effective, while immunostimulatory treatments that increase the rate at which infected monocytes are removed are best.

The impact of within-vector parasite development on the extrinsic incubation period

Mosquito-borne diseases, in particular malaria, have a significant burden worldwide leading to nearly half a million deaths each year. The malaria parasite requires a vertebrate host, such as a human, and a vector host, the Anopheles mosquito, to complete its full life cycle. Here, we focus on the parasite dynamics within the vector to examine the first appearance of sporozoites in the salivary glands, which indicates a first time of infectiousness of mosquitoes. The timing of this period of pathogen development in the mosquito until transmissibility, known as the extrinsic incubation period, remains poorly understood. We develop compartmental models of within-mosquito parasite dynamics fitted with experimental data on oocyst and sporozoite counts. We find that only a fraction of oocysts burst to release sporozoites and bursting must be delayed either via a time-dependent function or a gamma-distributed set of compartments. We use Bayesian inference to estimate distributions of parameters and determine that bursting rate is a key epidemiological parameter. A better understanding of the factors impacting the extrinsic incubation period will aid in the development of interventions to slow or stop the spread of malaria.

Accelerated viral dynamics in bat cell lines, with implications for zoonotic emergence

Bats host virulent zoonotic viruses without experiencing disease. A mechanistic understanding of the impact of bats’ virus hosting capacities, including uniquely constitutive immune pathways, on cellular-scale viral dynamics is needed to elucidate zoonotic emergence. We carried out virus infectivity assays on bat cell lines expressing induced and constitutive immune phenotypes, then developed a theoretical model of our in vitro system, which we fit to empirical data. Best fit models recapitulated expected immune phenotypes for representative cell lines, supporting robust antiviral defenses in bat cells that correlated with higher estimates for within-host viral propagation rates. In general, heightened immune responses limit pathogen-induced cellular morbidity, which can facilitate the establishment of rapidly-propagating persistent infections within-host. Rapidly-transmitting viruses that have evolved with bat immune systems will likely cause enhanced virulence following emergence into secondary hosts with immune systems that diverge from those unique to bats.

Bats can carry viruses that are deadly to other mammals without themselves showing serious symptoms. In fact, bats are natural reservoirs for viruses that have some of the highest fatality rates of any viruses that people acquire from wild animals – including rabies, Ebola and the SARS coronavirus.

Bats have a suite of antiviral defenses that keep the amount of virus in check. For example, some bats have an antiviral immune response called the interferon pathway perpetually switched on. In most other mammals, having such a hyper-vigilant immune response would cause harmful inflammation. Bats, however, have adapted anti-inflammatory traits that protect them from such harm, include the loss of certain genes that normally promote inflammation. However, no one has previously explored how these unique antiviral defenses of bats impact the viruses themselves.

Now, Brook et al. have studied this exact question using bat cells grown in the laboratory. The experiments made use of cells from one bat species – the black flying fox – in which the interferon pathway is always on, and another – the Egyptian fruit bat – in which this pathway is only activated during an infection. The bat cells were infected with three different viruses, and then Brook et al. observed how the interferon pathway helped keep the infections in check, before creating a computer model of this response.

The experiments and model helped reveal that the bats’ defenses may have a potential downside for other animals, including humans. In both bat species, the strongest antiviral responses were countered by the virus spreading more quickly from cell to cell. This suggests that bat immune defenses may drive the evolution of faster transmitting viruses, and while bats are well protected from the harmful effects of their own prolific viruses, other creatures like humans are not.

The findings may help to explain why bats are often the source for viruses that are deadly in humans. Learning more about bats' antiviral defenses and how they drive virus evolution may help scientists develop better ways to predict, prevent or limit the spread of viruses from bats to humans. More studies are needed in bats to help these efforts. In the meantime, the experiments highlight the importance of warning people to avoid direct contact with wild bats.

Could the human papillomavirus vaccines drive virulence evolution?

The human papillomavirus (HPV) vaccines hold great promise for preventing several cancers caused by HPV infections. Yet little attention has been given to whether HPV could respond evolutionarily to the new selection pressures imposed on it by the novel immunity response created by the vaccine. Here, we present and theoretically validate a mechanism by which the vaccine alters the transmission-recovery trade-off that constrains HPV's virulence such that higher oncogene expression is favoured. With a high oncogene expression strategy, the virus is able to increase its viral load and infected cell population before clearance by the vaccine, thus improving its chances of transmission. This new rapid cell-proliferation strategy is able to circulate between hosts with medium to high turnover rates of sexual partners. We also discuss the importance of better quantifying the duration of challenge infections and the degree to which a vaccinated host can shed virus. The generality of the models presented here suggests a wider applicability of this mechanism, and thus highlights the need to investigate viral oncogenicity from an evolutionary perspective.

A population model capturing dynamics of tuberculosis granulomas predicts host infection outcomes

Granulomas play a centric role in tuberculosis (TB) infection progression. Multiple granulomas usually develop within a single host. These granulomas are not synchronized in size or bacteria load, and will follow different trajectories over time. How the fate of individual granulomas influence overall infection outcome at host scale is not understood, although computational models have been developed to predict single granuloma behavior. Here we present a within-host population model that tracks granulomas in two key organs during Mycobacteria tuberculosis (Mtb) infection: lung and lymph nodes (LN). We capture various time courses of TB progression, including latency and reactivation. The model predicts that there is no steady state; rather it is a continuous process of progressing to active disease over differing time periods. This is consistent with recently posed ideas suggesting that latent TB exists as a spectrum of states and not a single state. The model also predicts a dual role for granuloma development in LNs during Mtb infection: in early phases of infection granulomas suppress infection by providing additional antigens to the site of immune priming; however, this induces a more rapid reactivation at later stages by disrupting immune responses. We identify mechanisms that strongly correlate with better host-level outcomes, including elimination of uncontained lung granulomas by inducing low levels of lung tissue damage and inhibition of bacteria dissemination within the lung.

Pathogen growth in insect hosts: inferring the importance of different mechanisms using stochastic models and response-time data

Pathogen population dynamics within individual hosts can alter disease epidemics and pathogen evolution, but our understanding of the mechanisms driving within-host dynamics is weak. Mathematical models have provided useful insights, but existing models have only rarely been subjected to rigorous tests, and their reliability is therefore open to question. Most models assume that initial pathogen population sizes are so large that stochastic effects due to small population sizes, so-called demographic stochasticity, are negligible, but whether this assumption is reasonable is unknown. Most models also assume that the dynamic effects of a host's immune system strongly affect pathogen incubation times or "response times," but whether such effects are important in real host-pathogen interactions is likewise unknown. Here we use data for a baculovirus of the gypsy moth to test models of within-host pathogen growth. By using Bayesian statistical techniques and formal model-selection procedures, we are able to show that the response time of the gypsy moth virus is strongly affected by both demographic stochasticity and a dynamic response of the host immune system. Our results imply that not all response-time variability can be explained by host and pathogen variability, and that immune system responses to infection may have important effects on population-level disease dynamics.

How sticky should a virus be? The impact of virus binding and release on transmission fitness using influenza as an example

Budding viruses face a trade-off: virions need to efficiently attach to and enter uninfected cells while newly generated virions need to efficiently detach from infected cells. The right balance between attachment and detachment-the right amount of stickiness-is needed for maximum fitness. Here, we design and analyse a mathematical model to study in detail the impact of attachment and detachment rates on virus fitness. We apply our model to influenza, where stickiness is determined by a balance of the haemagglutinin (HA) and neuraminidase (NA) proteins. We investigate how drugs, the adaptive immune response and vaccines impact influenza stickiness and fitness. Our model suggests that the location in the 'stickiness landscape' of the virus determines how well interventions such as drugs or vaccines are expected to work. We discuss why hypothetical NA enhancer drugs might occasionally perform better than the currently available NA inhibitors in reducing virus fitness. We show that an increased antibody or T-cell-mediated immune response leads to maximum fitness at higher stickiness. We further show that antibody-based vaccines targeting mainly HA or NA, which leads to a shift in stickiness, might reduce virus fitness above what can be achieved by the direct immunological action of the vaccine. Overall, our findings provide potentially useful conceptual insights for future vaccine and drug development and can be applied to other budding viruses beyond influenza.

The role of backward mutations on the within-host dynamics of HIV-1

The quality of life for patients infected with human immunodeficiency virus (HIV-1) has been positively impacted by the use of antiretroviral therapy (ART). However, the benefits of ART are usually halted by the emergence of drug resistance. Drug-resistant strains arise from virus mutations, as HIV-1 reverse transcription is prone to errors, with mutations normally carrying fitness costs to the virus. When ART is interrupted, the wild-type drug-sensitive strain rapidly out-competes the resistant strain, as the former strain is fitter than the latter in the absence of ART. One mechanism for sustaining the sensitive strain during ART is given by the virus mutating from resistant to sensitive strains, which is referred to as backward mutation. This is important during periods of treatment interruptions as prior existence of the sensitive strain would lead to replacement of the resistant strain. In order to assess the role of backward mutations in the dynamics of HIV-1 within an infected host, we analyze a mathematical model of two interacting virus strains in either absence or presence of ART. We study the effect of backward mutations on the definition of the basic reproductive number, and the value and stability of equilibrium points. The analysis of the model shows that, thanks to both forward and backward mutations, sensitive and resistant strains co-exist. In addition, conditions for the dominance of a viral strain with or without ART are provided. For this model, backward mutations are shown to be necessary for the persistence of the sensitive strain during ART.

Time and dose-dependent risk of pneumococcal pneumonia following influenza: a model for within-host interaction between influenza and Streptococcus pneumoniae

A significant fraction of seasonal and in particular pandemic influenza deaths are attributed to secondary bacterial infections. In animal models, influenza virus predisposes hosts to severe infection with both Streptococcus pneumoniae and Staphylococcus aureus. Despite its importance, the mechanistic nature of the interaction between influenza and pneumococci, its dependence on the timing and sequence of infections as well as the clinical and epidemiological consequences remain unclear. We explore an immune-mediated model of the viral-bacterial interaction that quantifies the timing and the intensity of the interaction. Taking advantage of the wealth of knowledge gained from animal models, and the quantitative understanding of the kinetics of pathogen-specific immunological dynamics, we formulate a mathematical model for immune-mediated interaction between influenza virus and S. pneumoniae in the lungs. We use the model to examine the pathogenic effect of inoculum size and timing of pneumococcal invasion relative to influenza infection, as well as the efficacy of antivirals in preventing severe pneumococcal disease. We find that our model is able to capture the key features of the interaction observed in animal experiments. The model predicts that introduction of pneumococcal bacteria during a 4-6 day window following influenza infection results in invasive pneumonia at significantly lower inoculum size than in hosts not infected with influenza. Furthermore, we find that antiviral treatment administered later than 4 days after influenza infection was not able to prevent invasive pneumococcal disease. This work provides a quantitative framework to study interactions between influenza and pneumococci and has the potential to accurately quantify the interactions. Such quantitative understanding can form a basis for effective clinical care, public health policies and pandemic preparedness.