Minimal within-host dengue models highlight the specific roles of the immune response in primary and secondary dengue infections

In recent years, the within-host viral dynamics of dengue infections have been increasingly characterized, and the relationship between aspects of these dynamics and the manifestation of severe disease has been increasingly probed. Despite this progress, there are few mathematical models of within-host dengue dynamics, and the ones that exist focus primarily on the general role of immune cells in the clearance of infected cells, while neglecting other components of the immune response in limiting viraemia. Here, by considering a suite of mathematical within-host dengue models of increasing complexity, we aim to isolate the critical components of the innate and the adaptive immune response that suffice in the reproduction of several well-characterized features of primary and secondary dengue infections. By building up from a simple target cell limited model, we show that only the innate immune response is needed to recover the characteristic features of a primary symptomatic dengue infection, while a higher rate of viral infectivity (indicative of antibody-dependent enhancement) and infected cell clearance by T cells are further needed to recover the characteristic features of a secondary dengue infection. We show that these minimal models can reproduce the increased risk of disease associated with secondary heterologous infections that arises as a result of a cytokine storm, and, further, that they are consistent with virological indicators that predict the onset of severe disease, such as the magnitude of peak viraemia, time to peak viral load, and viral clearance rate. Finally, we show that the effectiveness of these virological indicators to predict the onset of severe disease depends on the contribution of T cells in fuelling the cytokine storm.     

Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity

Dengue, a mosquito-borne virus of humans, infects over 50 million people annually. Infection with any of the four dengue serotypes induces protective immunity to that serotype, but does not confer long-term protection against infection by other serotypes. The immunological interactions between serotypes are of central importance in understanding epidemiological dynamics and anticipating the impact of dengue vaccines. We analysed a 38-year time series with 12 197 serotyped dengue infections from a hospital in Bangkok, Thailand. Using novel mechanistic models to represent different hypothesized immune interactions between serotypes, we found strong evidence that infection with dengue provides substantial short-term cross-protection against other serotypes (approx. 1–3 years). This is the first quantitative evidence that short-term cross-protection exists since human experimental infection studies performed in the 1950s. These findings will impact strategies for designing dengue vaccine studies, future multi-strain modelling efforts, and our understanding of evolutionary pressures in multi-strain disease systems.     

Dynamics of killer T cell inflation in viral infections

Upon acute viral infection, a typical cytotoxic T lymphocyte (CTL) response is characterized by a phase of expansion and contraction after which it settles at a relatively stable memory level. Recently, experimental data from mice infected with murine cytomegalovirus (MCMV) showed different and unusual dynamics. After acute infection had resolved, some antigen specific CTL started to expand over time despite the fact that no replicative virus was detectable. This phenomenon has been termed as "CTL memory inflation". In order to examine the dynamics of this system further, we developed a mathematical model analysing the impact of innate and adaptive immune responses. According to this model, a potentially important contributor to CTL inflation is competition between the specific CTL response and an innate natural killer (NK) cell response. Inflation occurs most readily if the NK cell response is more efficient than the CTL at reducing virus load during acute infection, but thereafter maintains a chronic virus load which is sufficient to induce CTL proliferation. The model further suggests that weaker NK cell mediated protection can correlate with more pronounced CTL inflation dynamics over time. We present experimental data from mice infected with MCMV which are consistent with the theoretical predictions. This model provides valuable information and may help to explain the inflation of CMV specific CD8+T cells seen in humans as they age.     

一类带有特异性免疫细胞钟形增殖率的慢性病毒感染模型的全局动力学性态

特异性免疫应答对控制宿主体内的病毒感染起着非常重要的作用.本文提出并研究了一类具有特异性免疫细胞钟形增殖率的慢性病毒感染模型.这里免疫细胞的钟形增殖意指当病毒载量足够大时其繁殖率会降低.病毒对免疫应答的损害也在本文的模型中被考虑.在找到该模型免疫应答基本再生数的同时,完整分析了其局部动力学行为.为了确定其全局动力学性态,应用中心流型理论对一些临界情形进行了分析,并通过构造适当的Dulac函数排除了该模型周期解的存在性.本文得到的结果显示在一定条件下模型会出现后向分支,这意味着模型的动力学性质会因初始状态的不同而改变.最后的数值模拟说明最终的单调和持续震荡对病毒种群和免疫应答都是有可能发生的.     

Effect of the initial dose of foot-and-mouth disease virus on the early viral dynamics within pigs

This paper investigates the early viral dynamics of foot-and-mouth disease (FMD) within infected pigs. Using an existing within-host model, we investigate whether individual variation can be explained by the effect of the initial dose of FMD virus. To do this, we consider the experimental data on the concentration of FMD virus genomes in the blood (viral load). In this experiment, 12 pigs were inoculated with one of three different doses of FMD virus: low; medium; or high. Measurements of the viral load were recorded over a time course of approximately 11 days for every 8 hours. The model is a set of deterministic differential equations with the following variables: viral load; virus in the interstitial space; and the proportion of epithelial cells available for infection, infected and uninfected. The model was fitted to the data for each animal individually and also simultaneously over all animals varying only the initial dose. We show that the general trend in the data can be explained by varying only the initial dose. The higher the initial dose the earlier the development of a detectable viral load.     

Dissecting the determinants of malaria chronicity: why within-host models struggle to reproduce infection dynamics

The duration of infection is fundamental to the epidemiological behaviour of any infectious disease, but remains one of the most poorly understood aspects of malaria. In endemic areas, the malaria parasite Plasmodium falciparum can cause both acute, severe infections and asymptomatic, chronic infections through its interaction with the host immune system. Frequent superinfection and massive parasite genetic diversity make it extremely difficult to accurately measure the distribution of infection lengths, complicating the estimation of basic epidemiological parameters and the prediction of the impact of interventions. Mathematical models have qualitatively reproduced parasite dynamics early during infection, but reproducing long-lived chronic infections remains much more challenging. Here, we construct a model of infection dynamics to examine the consequences of common biological assumptions for the generation of chronicity and the impact of co-infection. We find that although a combination of host and parasite heterogeneities are capable of generating chronic infections, they do so only under restricted parameter choices. Furthermore, under biologically plausible assumptions, co-infection of parasite genotypes can alter the course of infection of both the resident and co-infecting strain in complex non-intuitive ways. We outline the most important puzzles for within-host models of malaria arising from our analysis, and their implications for malaria epidemiology and control.     

Perspectives on the basic reproductive ratio.

The basic reproductive ratio, R0, is defined as the expected number of secondary infections arising from a single individual during his or her entire infectious period, in a population of susceptibles. This concept is fundamental to the study of epidemiology and within-host pathogen dynamics. Most importantly, R0 often serves as a threshold parameter that predicts whether an infection will spread. Related parameters which share this threshold behaviour, however, may or may not give the true value of R0. In this paper we give a brief overview of common methods of formulating R0 and surrogate threshold parameters from deterministic, non-structured models. We also review common means of estimating R0 from epidemiological data. Finally, we survey the recent use of R0 in assessing emerging diseases, such as severe acute respiratory syndrome and avian influenza, a number of recent livestock diseases, and vector-borne diseases malaria, dengue and West Nile virus.     

Modelling dynamics of the type I interferon response to in vitro viral infection.

Innate immunity is crucial in the early stages of resistance to novel viral infection. The family of cytokines known as the interferons (IFNs) forms an essential component of this system: they are responsible for signalling that an infection is underway and for promoting an antiviral response in susceptible cells. We construct a spatial stochastic model, parameterized by experimental data and informed by analytic approximation, to capture the dynamics of virus-IFN interaction during in vitro infection of Madin-Darby bovine kidney cell monolayers by Herpes simplex virus 1. The dose dependence of infection progression, subsequent monolayer destruction and IFN-beta production are investigated. Implications for in vivo infections, in particular the priming of susceptible cells by IFN-beta during infection, are considered.     

Modelling HIV-1 2-LTR dynamics following raltegravir intensification

A model of reservoir activation and viral replication is introduced accounting for the production of 2-LTR HIV-1 DNA circles following antiviral intensification with the HIV integrase inhibitor raltegravir, considering contributions of de novo infection events and exogenous sources of infected cells, including quiescent infected cell activation. The model shows that a monotonic increase in measured 2-LTR concentration post intensification is consistent with limited de novo infection primarily maintained by sources of infected cells unaffected by raltegravir, such as quiescent cell activation, while a transient increase in measured 2-LTR concentration is consistent with significant levels of efficient (R₀ > 1) de novo infection. The model is validated against patient data from the INTEGRAL study and is shown to have a statistically significant fit relative to the null hypothesis of random measurement variation about a mean. We obtain estimates and confidence intervals for the model parameters, including 2-LTR half-life. Seven of the 13 patients with detectable 2-LTR concentrations from the INTEGRAL study have measured 2-LTR dynamics consistent with significant levels of efficient replication of the virus prior to treatment intensification.     

Glycyrrhizic‐Acid‐Based Carbon Dots with High Antiviral Activity by Multisite Inhibition Mechanisms

With the gradual usage of carbon dots (CDs) in the area of antiviral research, attempts have been stepped up to develop new antiviral CDs with high biocompatibility and antiviral effects. In this study, a kind of highly biocompatible CDs (Gly‐CDs) is synthesized from active ingredient (glycyrrhizic acid) of Chinese herbal medicine by a hydrothermal method. Using the porcine reproductive and respiratory syndrome virus (PRRSV) as a model, it is found that the Gly‐CDs inhibit PRRSV proliferation by up to 5 orders of viral titers. Detailed investigations reveal that Gly‐CDs can inhibit PRRSV invasion and replication, stimulate antiviral innate immune responses, and inhibit the accumulation of intracellular reactive oxygen species (ROS) caused by PRRSV infection. Proteomics analysis demonstrates that Gly‐CDs can stimulate cells to regulate the expression of some host restriction factors, including DDX53 and NOS3, which are directly related to PRRSV proliferation. Moreover, it is found that Gly‐CDs also remarkably suppress the propagation of other viruses, such as pseudorabies virus (PRV) and porcine epidemic diarrhea virus (PEDV), suggesting the broad antiviral activity of Gly‐CDs. The integrated results demonstrate that Gly‐CDs possess extraordinary antiviral activity with multisite inhibition mechanisms, providing a promising candidate for alternative therapy for PRRSV infection.     

An explicit immunogenetic model of gastrointestinal nematode infection in sheep

Gastrointestinal nematodes are a global cause of disease and death in humans, wildlife and livestock. Livestock infection has historically been controlled with anthelmintic drugs, but the development of resistance means that alternative controls are needed. The most promising alternatives are vaccination, nutritional supplementation and selective breeding, all of which act by enhancing the immune response. Currently, control planning is hampered by reliance on the faecal egg count (FEC), which suffers from low accuracy and a nonlinear and indirect relationship with infection intensity and host immune responses. We address this gap by using extensive parasitological, immunological and genetic data on the sheep–Teladorsagia circumcincta interaction to create an immunologically explicit model of infection dynamics in a sheep flock that links host genetic variation with variation in the two key immune responses to predict the observed parasitological measures. Using our model, we show that the immune responses are highly heritable and by comparing selective breeding based on low FECs versus high plasma IgA responses, we show that the immune markers are a much improved measure of host resistance. In summary, we have created a model of host–parasite infections that explicitly captures the development of the adaptive immune response and show that by integrating genetic, immunological and parasitological understanding we can identify new immune-based markers for diagnosis and control.     

Modelling the impact of household size distribution on the transmission dynamics of COVID-19

Under the implementation of non-pharmaceutical interventions such as social distancing and lockdowns, household transmission has been shown to be significant for COVID-19, posing challenges for reducing incidence in settings where people are asked to self-isolate at home and to spend increasing amounts of time at home due to distancing measures. Accordingly, characteristics of households in a region have been shown to relate to transmission heterogeneity of the virus. We introduce a discrete-time stochastic epidemiological model to examine the impact of the household size distribution in a region on the transmission dynamics. We choose parameters to reflect incidence in two health regions of the Greater Vancouver area in British Columbia and simulate the impact of distancing measures on transmission, with household size distribution the only different parameter between simulations for the two regions. Our result suggests that the dissimilarity in household size distribution alone can cause significant differences in incidence of the two regions, and the distributions drive distinct dynamics that match reported cases. Furthermore, our model suggests that offering individuals a place to isolate outside their household can speed the decline in cases, and does so more effectively where there are more larger households.     

Multi-scale modelling reveals that early super-spreader events are a likely contributor to novel variant predominance

The emergence of new SARS-CoV-2 variants of concern (VOC) has hampered international efforts to contain the COVID-19 pandemic. VOCs have been characterized to varying degrees by higher transmissibility, worse infection outcomes and evasion of vaccine and infection-induced immunologic memory. VOCs are hypothesized to have originated from animal reservoirs, communities in regions with low surveillance and/or single individuals with poor immunologic control of the virus. Yet, the factors dictating which variants ultimately predominate remain incompletely characterized. Here we present a multi-scale model of SARS-CoV-2 dynamics that describes population spread through individuals whose viral loads and numbers of contacts (drawn from an over-dispersed distribution) are both time-varying. This framework allows us to explore how super-spreader events (SSE) (defined as greater than five secondary infections per day) contribute to variant emergence. We find stochasticity remains a powerful determinant of predominance. Variants that predominate are more likely to be associated with higher infectiousness, an SSE early after variant emergence and ongoing decline of the current dominant variant. Additionally, our simulations reveal that most new highly infectious variants that infect one or a few individuals do not achieve permanence in the population. Consequently, interventions that reduce super-spreading may delay or mitigate emergence of VOCs.     

Role of high-dose exposure in transmission hot zones as a driver of SARS-CoV-2 dynamics

Epidemiological data about SARS-CoV-2 spread indicate that the virus is not transmitted uniformly in the population. The transmission tends to be more effective in select settings that involve exposure to relatively high viral dose, such as in crowded indoor settings, assisted living facilities, prisons or food processing plants. To explore the effect on infection dynamics, we describe a new mathematical model where transmission can occur (i) in the community at large, characterized by low-dose exposure and mostly mild disease, and (ii) in so-called transmission hot zones, characterized by high-dose exposure that can be associated with more severe disease. The model yields different types of epidemiological dynamics, depending on the relative importance of hot zone and community transmission. Interesting dynamics occur if the rate of virus release/deposition from severely infected people is larger than that of mildly infected individuals. Under this assumption, we find that successful infection spread can hinge upon high-dose hot zone transmission, yet the majority of infections are predicted to occur in the community at large with mild disease. In this regime, residual hot zone transmission can account for continued virus spread during community lockdowns, and the suppression of hot zones after community interventions are relaxed can cause a prolonged lack of infection resurgence following the reopening of society. This gives rise to the notion that targeted interventions specifically reducing virus transmission in the hot zones have the potential to suppress overall infection spread, including in the community at large. Epidemiological trends in the USA and Europe are interpreted in light of this model.     

Towards a quantitative understanding of the within-host dynamics of influenza A infections

Although the influenza A virus has been extensively studied, a quantitative understanding of the infection dynamics is still lacking. To make progress in this direction, we designed several mathematical models and compared them with data from influenza A infections of mice. We find that the immune response (IR) plays an important part in the infection dynamics. Both an innate and an adaptive IR are required to provide adequate explanation of the data. In contrast, regrowth of epithelial cells did not seem to be an important mechanism on the time scale of the infection. We also find that different model variants for both innate and adaptive responses fit the data well, indicating the need for additional data to allow further model discrimination.     

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.     

Hyperendemicity associated with increased dengue burden

Over 105 million dengue infections are estimated to occur annually. Understanding the disease dynamics of dengue is often difficult due to multiple strains circulating within a population. Interactions between dengue serotype dynamics may result in complex cross-immunity dynamics at the population level and create difficulties in terms of formulating intervention strategies for the disease. In this study, a nationally representative 16-year time series with over 43 000 serotyped dengue infections was used to infer the long-run effects of between and within strain interactions and their impacts on past outbreaks. We used a novel identification strategy incorporating sign-identified Bayesian vector autoregressions, using structural impulse responses, historical decompositions and counterfactual analysis to conduct inference on dengue dynamics post-estimation. We found that on the population level: (i) across-serotype interactions on the population level were highly persistent, with a one time increase in any other serotype associated with long run decreases in the serotype of interest (range: 0.5–2.5 years) and (ii) over 38.7% of dengue cases of any serotype were associated with across-serotype interactions. The findings in this paper will substantially impact public health policy interventions with respect to dengue.     

Demographic buffering: titrating the effects of birth rate and imperfect immunity on epidemic dynamics

Host demography can alter the dynamics of infectious disease. In the case of perfectly immunizing infections, observations of strong sensitivity to demographic variation have been mechanistically explained through analysis of the susceptible–infected–recovered (SIR) model that assumes lifelong immunity following recovery from infection. When imperfect immunity is incorporated into this framework via the susceptible–infected–recovered–susceptible (SIRS) model, with individuals regaining full susceptibility following recovery, we show that rapid loss of immunity is predicted to buffer populations against the effects of demographic change. However, this buffering is contrary to the dependence on demography recently observed for partially immunizing infections such as rotavirus and respiratory syncytial virus. We show that this discrepancy arises from a key simplification embedded in the SIR(S) framework, namely that the potential for differential immune responses to repeat exposures is ignored. We explore the minimum additional immunological information that must be included to reflect the range of observed dependencies on demography. We show that including partial protection and lower transmission following primary infection is sufficient to capture more realistic reduced levels of buffering, in addition to changes in epidemic timing, across a range of partially and fully immunizing infections. Furthermore, our results identify key variables in this relationship, including R₀.