Vector-Borne Pathogen and Host Evolution in a Structured Immuno-Epidemiological System

The impact of blood on vector-borne diseases with emphasis on mosquitoes and sand flies

The impact of blood and its factors on vector-borne diseases is significant and multifaceted yet understudied. While blood is expected to play a central role in transmission, pathogen development, vector behavior, and vector competence, in experimental settings, most studies are developed in the frame of a single, infected blood meal. To effectively combat vector-borne diseases, we need to determine what is the influence of insect blood-feeding behavior on transmission and development of pathogens, toward translation to natural field settings. This review summarizes current findings, highlights key gaps, and outlines future research directions to enhance our understanding of the role of blood in vector-borne disease transmission.

Contact-number-driven virus evolution: A multi-level modeling framework for the evolution of acute or persistent RNA virus infection

Viruses evolve in infected host populations, and host population dynamics affect viral evolution. RNA viruses with a short duration of infection and a high peak viral load, such as SARS-CoV-2, are maintained in human populations. By contrast, RNA viruses characterized by a long infection duration and a low peak viral load (e.g., borna disease virus) can be maintained in nonhuman populations, and the process of the evolution of persistent viruses has rarely been explored. Here, using a multi-level modeling approach including both individual-level virus infection dynamics and population-scale transmission, we consider virus evolution based on the host environment, specifically, the effect of the contact history of infected hosts. We found that, with a highly dense contact history, viruses with a high virus production rate but low accuracy are likely to be optimal, resulting in a short infectious period with a high peak viral load. In contrast, with a low-density contact history, viral evolution is toward low virus production but high accuracy, resulting in long infection durations with low peak viral load. Our study sheds light on the origin of persistent viruses and why acute viral infections but not persistent virus infection tends to prevail in human society.

A nested model with boosting and waning of immunity from Tilapia Lake Virus infection with distributed resistance to pathogens carrier-state

This paper proposes and analyzes an immune-structured population model of tilapia subject to Tilapia Lake Virus (TiLV) disease. The model incorporates within-host dynamics, used to describe the interaction between the pathogen, the immune system and the waning of immunity. Individuals infected with a low dose acquire a low immunity level and those infected with a high dose acquire a high level of immunity. Since individuals' immune status plays an important role in the spread of infectious diseases at the population level, the within-host dynamics are connected to the between-host dynamics in the population. We define an explicit formula for the reproductive number [Formula: see text] and show that the disease-free equilibrium is locally asymptotically stable when [Formula: see text], while it is unstable when [Formula: see text]. Furthermore, we prove that an endemic equilibrium exists. We also study the influence of the initial distribution of host resistance on the spread of the disease, and find that hosts' initial resistance plays a crucial role in the disease dynamics. This suggests that the genetic selection aiming to improve hosts' initial resistance to TiLV could help fight the disease. The results also point out the crucial role played by the inoculum size. We find that the higher the initial inoculum size, the faster the dynamics of infection. Moreover, if the initial inoculum size is below a certain threshold, it may not result in an outbreak at the between-host level. Finally, the model shows that there is a strong negative correlation between heterogeneity and the probability of pathogen invasion.

A mathematical model for tilapia lake virus transmission with waning immunity

The goal of this paper is to investigate the influence of the waning immunity on the dynamics of Tilapia Lake Virus (TiLV) transmission in wild and farmed tilapia within freshwater. We formulate a model for which susceptible individuals can contract the disease in two ways: (i) direct mode caused by contact with infected individuals; (ii) indirect mode due to the presence of pathogenic agents in the water. We obtain an age-structured model which combines both age since infection and age since recovery. We derive an explicit formula for the reproductive number R0 and show that the disease-free equilibrium is locally asymptotically stable when, R0<1. We discuss on the form of the waning immunity parameter and show numerically that a Hopf bifurcation may occur for suitable immunity parameter values, which means that there is a periodic solution around the endemic equilibrium when, R0>1.

Mechanistic models of Rift Valley fever virus transmission: A systematic review

Rift Valley fever (RVF) is a zoonotic arbovirosis which has been reported across Africa including the northernmost edge, South West Indian Ocean islands, and the Arabian Peninsula. The virus is responsible for high abortion rates and mortality in young ruminants, with economic impacts in affected countries. To date, RVF epidemiological mechanisms are not fully understood, due to the multiplicity of implicated vertebrate hosts, vectors, and ecosystems. In this context, mathematical models are useful tools to develop our understanding of complex systems, and mechanistic models are particularly suited to data-scarce settings. Here, we performed a systematic review of mechanistic models studying RVF, to explore their diversity and their contribution to the understanding of this disease epidemiology. Researching Pubmed and Scopus databases (October 2021), we eventually selected 48 papers, presenting overall 49 different models with numerical application to RVF. We categorized models as theoretical, applied, or grey, depending on whether they represented a specific geographical context or not, and whether they relied on an extensive use of data. We discussed their contributions to the understanding of RVF epidemiology, and highlighted that theoretical and applied models are used differently yet meet common objectives. Through the examination of model features, we identified research questions left unexplored across scales, such as the role of animal mobility, as well as the relative contributions of host and vector species to transmission. Importantly, we noted a substantial lack of justification when choosing a functional form for the force of infection. Overall, we showed a great diversity in RVF models, leading to important progress in our comprehension of epidemiological mechanisms. To go further, data gaps must be filled, and modelers need to improve their code accessibility.

A Model for Brucellosis Disease Incorporating Age of Infection and Waning Immunity

This paper proposes a model for brucellosis transmission. The model takes into account the age of infection and waning immunity, that is, the progressive loss of immunity after recovery. Three routes of transmissions are considered: vertical transmission, and both direct and indirect routes of horizontal transmission. According to the well-posedness results, we provide explicit formulas for the equilibria. Next, we derive the basic reproduction number R0 and prove some stability results depending on the basic reproductive number. Finally, we perform numerical simulations using model parameters estimated from biological data to confirm our theoretical results. The results of these simulations suggest that for certain values of parameters, there will be periodic outbreaks of epidemics, and the disease will not be eradicated from the population. Our results also highlight the fact that the birth rate of cattle significantly influences the dynamics of the disease. The proposed model can be of a good use in studying the effects of vaccination on the cattle population.

Sensitivity Analysis in an Immuno-Epidemiological Vector-Host Model

Sensitivity Analysis (SA) is a useful tool to measure the impact of changes in model parameters on the infection dynamics, particularly to quantify the expected efficacy of disease control strategies. SA has only been applied to epidemic models at the population level, ignoring the effect of within-host virus-with-immune-system interactions on the disease spread. Connecting the scales from individual to population can help inform drug and vaccine development. Thus the value of understanding the impact of immunological parameters on epidemiological quantities. Here we consider an age-since-infection structured vector-host model, in which epidemiological parameters are formulated as functions of within-host virus and antibody densities, governed by an ODE system. We then use SA for these immuno-epidemiological models to investigate the impact of immunological parameters on population-level disease dynamics such as basic reproduction number, final size of the epidemic or the infectiousness at different phases of an outbreak. As a case study, we consider Rift Valley Fever Disease utilizing parameter estimations from prior studies. SA indicates that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\%$$\end{document}1% increase in within-host pathogen growth rate can lead up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal R_0,$$\end{document}R0, up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1 \%$$\end{document}1% increase in steady-state infected host abundance, and up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4\%$$\end{document}4% increase in infectiousness of hosts when the reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal R_0$$\end{document}R0 is larger than one. These significant increases in population-scale disease quantities suggest that control strategies that reduce the within-host pathogen growth can be important in reducing disease prevalence.

A Network Immuno-Epidemiological HIV Model

In this paper we formulate a multi-scale nested immuno-epidemiological model of HIV on complex networks. The system is described by ordinary differential equations coupled with a partial differential equation. First, we prove the existence and uniqueness of the immunological model and then establish the well-posedness of the multi-scale model. We derive an explicit expression of the basic reproduction number [Formula: see text] of the immuno-epidemiological model. The system has a disease-free equilibrium and an endemic equilibrium. The disease-free equilibrium is globally stable when [Formula: see text] and unstable when [Formula: see text]. Numerical simulations suggest that [Formula: see text] increases as the number of nodes in the network increases. Further, we find that for a scale-free network the number of infected individuals at equilibrium is a hump-like function of the within-host reproduction number; however, the dependence becomes monotone if the network has predominantly low connectivity nodes or high connectivity nodes.

AN IMMUNO-EPIDEMIOLOGICAL VECTOR–HOST MODEL WITH WITHIN-VECTOR VIRAL KINETICS

A current challenge for disease modeling and public health is understanding pathogen dynamics across scales since their ecology and evolution ultimately operate on several coupled scales. This is particularly true for vector-borne diseases, where within-vector, within-host, and between vector–host populations all play crucial roles in diversity and distribution of the pathogen. Despite recent modeling efforts to determine the effect of within-host virus-immune response dynamics on between-host transmission, the role of within-vector viral dynamics on disease spread is overlooked. Here, we formulate an age-since-infection-structured epidemic model coupled to nonlinear ordinary differential equations describing within-host immune-virus dynamics and within-vector viral kinetics, with feedbacks across these scales. We first define the within-host viral-immune response and within-vector viral kinetics-dependent basic reproduction number [Formula: see text] Then we prove that whenever [Formula: see text] the disease-free equilibrium is locally asymptotically stable, and under certain biologically interpretable conditions, globally asymptotically stable. Otherwise, if [Formula: see text] it is unstable and the system has a unique positive endemic equilibrium. In the special case of constant vector to host inoculum size, we show the positive equilibrium is locally asymptotically stable and the disease is weakly uniformly persistent. Furthermore, numerical results suggest that within-vector-viral kinetics and dynamic inoculum size may play a substantial role in epidemics. Finally, we address how the model can be utilized to better predict the success of control strategies such as vaccination and drug treatment.

Infection severity across scales in multi-strain immuno-epidemiological Dengue model structured by host antibody level

Infection by distinct Dengue virus serotypes and host immunity are intricately linked. In particular, certain levels of cross-reactive antibodies in the host may actually enhance infection severity leading to Dengue hemorrhagic fever (DHF). The coupled immunological and epidemiological dynamics of Dengue calls for a multi-scale modeling approach. In this work, we formulate a within-host model which mechanistically recapitulates characteristics of antibody dependent enhancement in Dengue infection. The within-host scale is then linked to epidemiological spread by a vector-host partial differential equation model structured by host antibody level. The coupling allows for dynamic population-wide antibody levels to be tracked through primary and secondary infections by distinct Dengue strains, along with waning of cross-protective immunity after primary infection. Analysis of both the within-host and between-host systems are conducted. Stability results in the epidemic model are formulated via basic and invasion reproduction numbers as a function of immunological variables. Additionally, we develop numerical methods in order to simulate the multi-scale model and assess the influence of parameters on disease spread and DHF prevalence in the population.

A general method for multiscale modelling of vector-borne disease systems

The inability to develop multiscale models which can describe vector-borne disease systems in terms of the complete pathogen life cycle which represents multiple targets for control has hindered progress in our efforts to control, eliminate and even eradicate these multi-host infections. This is because it is currently not easy to determine precisely where and how in the life cycles of vector-borne disease systems the key constrains which are regarded as crucial in regulating pathogen population dynamics in both the vertebrate host and vector host operate. In this article, we present a general method for development of multiscale models of vector-borne disease systems which integrate the within-host and between-host scales for the two hosts (a vertebrate host and a vector host) that are implicated in vector-borne disease dynamics. The general multiscale modelling method is an extension of our previous work on multiscale models of infectious disease systems which established a basic science and accompanying theory of how pathogen population dynamics at within-host scale scales up to between-host scale and in turn how it scales down from between-host scale to within-host scale. Further, the general method is applied to multiscale modelling of human onchocerciasis-a vector-borne disease system which is sometimes called river blindness as a case study.

A two-stage experimental design for dilution assays

Dilution assays to determine solute concentration have found wide use in biomedical research. Many dilution assays return imprecise concentration estimates because they are only done to orders of magnitude. Previous statistical work has focused on how to design efficient experiments that can return more precise estimates, however this work has not considered the practical difficulties of implementing these designs in the laboratory. We developed a two-stage experiment with a first stage that obtains an order of magnitude estimate and a second stage that concentrates effort on the most informative dilution to increase estimator precision. We show using simulations and an empirical example that the best two-stage experimental designs yield estimates that are remarkably more accurate than standard methods with equivalent effort. This work demonstrates how to utilize previous advances in experimental design in a manner consistent with current laboratory practice. We expect that multi-stage designs will prove to be useful for obtaining precise estimates with minimal experimental effort.

Structural and Practical Identifiability Issues of Immuno-Epidemiological Vector-Host Models with Application to Rift Valley Fever

In this article, we discuss the structural and practical identifiability of a nested immuno-epidemiological model of arbovirus diseases, where host-vector transmission rate, host recovery, and disease-induced death rates are governed by the within-host immune system. We incorporate the newest ideas and the most up-to-date features of numerical methods to fit multi-scale models to multi-scale data. For an immunological model, we use Rift Valley Fever Virus (RVFV) time-series data obtained from livestock under laboratory experiments, and for an epidemiological model we incorporate a human compartment to the nested model and use the number of human RVFV cases reported by the CDC during the 2006-2007 Kenya outbreak. We show that the immunological model is not structurally identifiable for the measurements of time-series viremia concentrations in the host. Thus, we study the non-dimensionalized and scaled versions of the immunological model and prove that both are structurally globally identifiable. After fixing estimated parameter values for the immunological model derived from the scaled model, we develop a numerical method to fit observable RVFV epidemiological data to the nested model for the remaining parameter values of the multi-scale system. For the given (CDC) data set, Monte Carlo simulations indicate that only three parameters of the epidemiological model are practically identifiable when the immune model parameters are fixed. Alternatively, we fit the multi-scale data to the multi-scale model simultaneously. Monte Carlo simulations for the simultaneous fitting suggest that the parameters of the immunological model and the parameters of the immuno-epidemiological model are practically identifiable. We suggest that analytic approaches for studying the structural identifiability of nested models are a necessity, so that identifiable parameter combinations can be derived to reparameterize the nested model to obtain an identifiable one. This is a crucial step in developing multi-scale models which explain multi-scale data.