Lyme Disease Models of Tick-Mouse Dynamics with Seasonal Variation in Births, Deaths, and Tick Feeding.

Lyme disease is the most common vector-borne disease in the United States impacting the Northeast and Midwest at the highest rates. Recently, it has become established in southeastern and south-central regions of Canada. In these regions, Lyme disease is caused by Borrelia burgdorferi, which is transmitted to humans by an infected Ixodes scapularis tick. Understanding the parasite-host interaction is critical as the white-footed mouse is one of the most competent reservoir for B. burgdorferi. The cycle of infection is driven by tick larvae feeding on infected mice that molt into infected nymphs and then transmit the disease to another susceptible host such as mice or humans. Lyme disease in humans is generally caused by the bite of an infected nymph. The main aim of this investigation is to study how diapause delays and demographic and seasonal variability in tick births, deaths, and feedings impact the infection dynamics of the tick-mouse cycle. We model tick-mouse dynamics with fixed diapause delays and more realistic Erlang distributed delays through delay and ordinary differential equations (ODEs). To account for demographic and seasonal variability, the ODEs are generalized to a continuous-time Markov chain (CTMC). The basic reproduction number and parameter sensitivity analysis are computed for the ODEs. The CTMC is used to investigate the probability of Lyme disease emergence when ticks and mice are introduced, a few of which are infected. The probability of disease emergence is highly dependent on the time and the infected species introduced. Infected mice introduced during the summer season result in the highest probability of disease emergence.

Modeling the Immune Response for Pathogenic and Nonpathogenic Orthohantavirus Infections in Human Lung Microvasculature Endothelial Cells.

Hantaviruses, genus Orthohantavirus, family Hantaviridae, order Bunyavirales, are negative-sense, single-stranded, tri-segmented RNA viruses that persistently infect rodents, shrews, and moles. Of these, only certain virus species harbored by rodents are pathogenic to humans. Infection begins with inhalation of virus particles into the lung and trafficking to the lung microvascular endothelial cells (LMVEC). The reason why certain rodent-borne hantavirus species are pathogenic has long been hypothesized to be related to their ability to downregulate and dysregulate the immune response as well as increase vascular permeability of infected endothelial cells. We set out to study the temporal dynamics of host immune response modulation in primary human LMVECs following infection by Prospect Hill (nonpathogenic), Andes (pathogenic), and Hantaan (pathogenic) viruses. We measured the level of RNA transcripts for genes representing antiviral, proinflammatory, anti-inflammatory, and metabolic pathways from 12 to 72 h with time points every 12 h. Gene expression analysis in conjunction with mathematical modeling revealed a similar profile for all three viruses in terms of upregulated genes that partake in interferon signaling (TLR3, IRF7, IFNB1), host immune cell recruitment (CXCL10, CXCL11, and CCL5), and host immune response modulation (IDO1). We examined secreted protein levels of IFN-ß, CXCL10, CXCL11, CCL5, and IDO in two male and two female primary HLMVEC donors at 48 and 60 h post infection. All three viruses induced similar levels of CCL5, CXCL10, and CXCL11 within a particular donor, and the levels were similar in three of the four donors. All three viruses induced different protein secretion levels for both IFN-ß and IDO and secretion levels differed between donors. In conclusion, we show that there was no difference in the transcriptional profiles of key genes in primary HLMVECs following infection by pathogenic and nonpathogenic hantaviruses, with protein secretion levels being more donor-specific than virus-specific.

A Hybrid Epidemic Model to Explore Stochasticity in COVID-19 Dynamics.

The dynamic nature of the COVID-19 pandemic has demanded a public health response that is constantly evolving due to the novelty of the virus. Many jurisdictions in the USA, Canada, and across the world have adopted social distancing and recommended the use of face masks. Considering these measures, it is prudent to understand the contributions of subpopulations-such as "silent spreaders"-to disease transmission dynamics in order to inform public health strategies in a jurisdiction-dependent manner. Additionally, we and others have shown that demographic and environmental stochasticity in transmission rates can play an important role in shaping disease dynamics. Here, we create a model for the COVID-19 pandemic by including two classes of individuals: silent spreaders, who either never experience a symptomatic phase or remain undetected throughout their disease course; and symptomatic spreaders, who experience symptoms and are detected. We fit the model to real-time COVID-19 confirmed cases and deaths to derive the transmission rates, death rates, and other relevant parameters for multiple phases of outbreaks in British Columbia (BC), Canada. We determine the extent to which SilS contributed to BC's early wave of disease transmission as well as the impact of public health interventions on reducing transmission from both SilS and SymS. To do this, we validate our model against an existing COVID-19 parameterized framework and then fit our model to clinical data to estimate key parameter values for different stages of BC's disease dynamics. We then use these parameters to construct a hybrid stochastic model that leverages the strengths of both a time-nonhomogeneous discrete process and a stochastic differential equation model. By combining these previously established approaches, we explore the impact of demographic and environmental variability on disease dynamics by simulating various scenarios in which a COVID-19 outbreak is initiated. Our results demonstrate that variability in disease transmission rate impacts the probability and severity of COVID-19 outbreaks differently in high- versus low-transmission scenarios.

Data-driven models for replication kinetics of Orthohantavirus infections.

The Hantaviridae constitute a family of viruses harbored by mice, rats, shrews, voles, moles and bats. Intriguingly, only viruses harbored by mice and rats may cause disease in humans with up to 40% case fatality rate in the Americas. Transmission of virus from rodents to humans occurs via the respiratory route and results in replication of the virus in the microvascular endothelial cells of the lung or kidney. Understanding the replication kinetics of these viruses in various cell types and how replication is abrogated by the host is critical to the development of effective therapeutics for treatment for which there are none. We formulate several new ordinary differential equation (ODE) models to examine the replication kinetics of Prospect Hill orthohantavirus (PHV). The models are distinguished by the distribution of the viral replication delay. A new threshold, RGE, the genome equivalent replication number, is defined in terms of the model parameters. New final density relations are derived that associate RGE to the asymptotic number of virions in each model. All models are fit to real time (qRT)-PCR data of genomic RNA from PHV released from Vero E6 cells over 192 h. A sensitivity analysis of the parameters is performed and models are tested for best fit. Our findings provide a basis for future research into formulating more complex mathematical models for evaluation of the replication of hantaviruses in various cell types and sources.

Stochastic models of infectious diseases in a periodic environment with application to cholera epidemics.

Seasonal variation affects the dynamics of many infectious diseases including influenza, cholera and malaria. The time when infectious individuals are first introduced into a population is crucial in predicting whether a major disease outbreak occurs. In this investigation, we apply a time-nonhomogeneous stochastic process for a cholera epidemic with seasonal periodicity and a multitype branching process approximation to obtain an analytical estimate for the probability of an outbreak. In particular, an analytic estimate of the probability of disease extinction is shown to satisfy a system of ordinary differential equations which follows from the backward Kolmogorov differential equation. An explicit expression for the mean (resp. variance) of the first extinction time given an extinction occurs is derived based on the analytic estimate for the extinction probability. Our results indicate that the probability of a disease outbreak, and mean and standard derivation of the first time to disease extinction are periodic in time and depend on the time when the infectious individuals or free-living pathogens are introduced. Numerical simulations are then carried out to validate the analytical predictions using two examples of the general cholera model. At the end, the developed theoretical results are extended to more general models of infectious diseases.

Effects of environmental variability on superspreading transmission events in stochastic epidemic models.

Superspreaders (individuals with a high propensity for disease spread) have played a pivotal role in recent emerging and re-emerging diseases. In disease outbreak studies, host heterogeneity based on demographic (e.g. age, sex, vaccination status) and environmental (e.g. climate, urban/rural residence, clinics) factors are critical for the spread of infectious diseases, such as Ebola and Middle East Respiratory Syndrome (MERS). Transmission rates can vary as demographic and environmental factors are altered naturally or due to modified behaviors in response to the implementation of public health strategies. In this work, we develop stochastic models to explore the effects of demographic and environmental variability on human-to-human disease transmission rates among superspreaders in the case of Ebola and MERS. We show that the addition of environmental variability results in reduced probability of outbreak occurrence, however the severity of outbreaks that do occur increases. These observations have implications for public health strategies that aim to control environmental variables.

Probability of a zoonotic spillover with seasonal variation.

Zoonotic infectious diseases are spread from animals to humans. It is estimated that over 60% of human infectious diseases are zoonotic and 75% of them are emerging zoonoses. The majority of emerging zoonotic infectious diseases are caused by viruses including avian influenza, rabies, Ebola, coronaviruses and hantaviruses. Spillover of infection from animals to humans depends on a complex transmission pathway, which is influenced by epidemiological and environmental processes. In this investigation, the focus is on direct transmission between animals and humans and the effects of seasonal variations on the transmission and recovery rates. Fluctuations in transmission and recovery, besides being influenced by physiological processes and behaviors of pathogen and host, are driven by seasonal variations in temperature, humidity or rainfall. A new time-nonhomogeneous stochastic process is formulated for infectious disease spread from animals to humans when transmission and recovery rates are time-periodic. A branching process approximation is applied near the disease-free state to predict the probability of the first spillover event from animals to humans. This probability is a periodic function of the time when infection is introduced into the animal population. It is shown that the highest risk of a spillover depends on a combination of animal to human transmission, animal to animal transmission and animal recovery. The results are applied to a stochastic model for avian influenza with spillover from domestic poultry to humans.

The effect of demographic and environmental variability on disease outbreak for a dengue model with a seasonally varying vector population.

Seasonal changes in temperature, humidity, and rainfall affect vector survival and emergence of mosquitoes and thus impact the dynamics of vector-borne disease outbreaks. Recent studies of deterministic and stochastic epidemic models with periodic environments have shown that the average basic reproduction number is not sufficient to predict an outbreak. We extend these studies to time-nonhomogeneous stochastic dengue models with demographic variability wherein the adult vectors emerge from the larval stage vary periodically. The combined effects of variability and periodicity provide a better understanding of the risk of dengue outbreaks. A multitype branching process approximation of the stochastic dengue model near the disease-free periodic solution is used to calculate the probability of a disease outbreak. The approximation follows from the solution of a system of differential equations derived from the backward Kolmogorov differential equation. This approximation shows that the risk of a disease outbreak is also periodic and depends on the particular time and the number of the initial infected individuals. Numerical examples are explored to demonstrate that the estimates of the probability of an outbreak from that of branching process approximations agree well with that of the continuous-time Markov chain. In addition, we propose a simple stochastic model to account for the effects of environmental variability on the emergence of adult vectors from the larval stage.

Disease Emergence in Multi-Patch Stochastic Epidemic Models with Demographic and Seasonal Variability.

Factors such as seasonality and spatial connectivity affect the spread of an infectious disease. Accounting for these factors in infectious disease models provides useful information on the times and locations of greatest risk for disease outbreaks. In this investigation, stochastic multi-patch epidemic models are formulated with seasonal and demographic variability. The stochastic models are used to investigate the probability of a disease outbreak when infected individuals are introduced into one or more of the patches. Seasonal variation is included through periodic transmission and dispersal rates. Multi-type branching process approximation and application of the backward Kolmogorov differential equation lead to an estimate for the probability of a disease outbreak. This estimate is also periodic and depends on the time, the location, and the number of initial infected individuals introduced into the patch system as well as the magnitude of the transmission and dispersal rates and the connectivity between patches. Examples are given for seasonal transmission and dispersal in two and three patches.

Coinfections by noninteracting pathogens are not independent and require new tests of interaction.

If pathogen species, strains, or clones do not interact, intuition suggests the proportion of coinfected hosts should be the product of the individual prevalences. Independence consequently underpins the wide range of methods for detecting pathogen interactions from cross-sectional survey data. However, the very simplest of epidemiological models challenge the underlying assumption of statistical independence. Even if pathogens do not interact, death of coinfected hosts causes net prevalences of individual pathogens to decrease simultaneously. The induced positive correlation between prevalences means the proportion of coinfected hosts is expected to be higher than multiplication would suggest. By modelling the dynamics of multiple noninteracting pathogens causing chronic infections, we develop a pair of novel tests of interaction that properly account for nonindependence between pathogens causing lifelong infection. Our tests allow us to reinterpret data from previous studies including pathogens of humans, plants, and animals. Our work demonstrates how methods to identify interactions between pathogens can be updated using simple epidemic models.

Modelling Vector Transmission and Epidemiology of Co-Infecting Plant Viruses.

Co-infection of plant hosts by two or more viruses is common in agricultural crops and natural plant communities. A variety of models have been used to investigate the dynamics of co-infection which track only the disease status of infected and co-infected plants, and which do not explicitly track the density of inoculative vectors. Much less attention has been paid to the role of vector transmission in co-infection, that is, acquisition and inoculation and their synergistic and antagonistic interactions. In this investigation, a general epidemiological model is formulated for one vector species and one plant species with potential co-infection in the host plant by two viruses. The basic reproduction number provides conditions for successful invasion of a single virus. We derive a new invasion threshold which provides conditions for successful invasion of a second virus. These two thresholds highlight some key epidemiological parameters important in vector transmission. To illustrate the flexibility of our model, we examine numerically two special cases of viral invasion. In the first case, one virus species depends on an autonomous virus for its successful transmission and in the second case, both viruses are unable to invade alone but can co-infect the host plant when prevalence is high.

The effect of delay in viral production in within-host models during early infection.

Delay in viral production may have a significant impact on the early stages of infection. During the eclipse phase, the time from viral entry until active production of viral particles, no viruses are produced. This delay affects the probability that a viral infection becomes established and timing of the peak viral load. Deterministic and stochastic models are formulated with either multiple latent stages or a fixed delay for the eclipse phase. The deterministic model with multiple latent stages approaches in the limit the model with a fixed delay as the number of stages approaches infinity. The deterministic model framework is used to formulate continuous-time Markov chain and stochastic differential equation models. The probability of a minor infection with rapid viral clearance as opposed to a major full-blown infection with a high viral load is estimated from a branching process approximation of the Markov chain model and the results are confirmed through numerical simulations. In addition, parameter values for influenza A are used to numerically estimate the time to peak viral infection and peak viral load for the deterministic and stochastic models. Although the average length of the eclipse phase is the same in each of the models, as the number of latent stages increases, the numerical results show that the time to viral peak and the peak viral load increase.

Models of cytokine dynamics in the inflammatory response of viral zoonotic infectious diseases.

Inflammatory responses to an infection from a zoonotic pathogen, such as avian influenza viruses, hantaviruses and some coronaviruses, are distinctly different in their natural reservoir versus human host. While not as well studied in the natural reservoirs, the pro-inflammatory response and viral replication appear controlled and show no obvious pathology. In contrast, infection in humans results in an initial high viral load marked by an aggressive pro-inflammatory response known as a cytokine storm. The key difference in the course of the infection between the reservoir and human host is the inflammatory response. In this investigation, we apply a simple two-component differential equation model for pro-inflammatory and anti-inflammatory responses and a detailed mathematical analysis to identify specific regions in parameter space for single stable endemic equilibrium, bistability or periodic solutions. The extensions of the deterministic model to two stochastic models account for variability in responses seen at the cell (local) or tissue (global) levels. Numerical solutions of the stochastic models exhibit outcomes that are typical of a chronic infection in the natural reservoir or a cytokine storm in human infection. In the chronic infection, occasional flare-ups between high and low responses occur when model parameters are in a region of bistability or periodic solutions. The cytokine storm with a vigorous pro-inflammatory response and less vigorous anti-inflammatory response occurs in the parameter region for a single stable endemic equilibrium with a strong pro-inflammatory response. The results of the model analyses and the simulations are interpreted in terms of the functional role of the cytokines and the inflammatory responses seen in infection of the natural reservoir or of the human host.

Stochastic two-group models with transmission dependent on host infectivity or susceptibility.

Stochastic epidemic models with two groups are formulated and applied to emerging and re-emerging infectious diseases. In recent emerging diseases, disease spread has been attributed to superspreaders, highly infectious individuals that infect a large number of susceptible individuals. In some re-emerging infectious diseases, disease spread is attributed to waning immunity in susceptible hosts. We apply a continuous-time Markov chain (CTMC) model to study disease emergence or re-emergence from different groups, where the transmission rates depend on either the infectious host or the susceptible host. Multitype branching processes approximate the dynamics of the CTMC model near the disease-free equilibrium and are used to estimate the probability of a minor or a major epidemic. It is shown that the probability of a major epidemic is greater if initiated by an individual from the superspreader group or by an individual from the highly susceptible group. The models are applied to Severe Acute Respiratory Syndrome and measles.

Duration of a minor epidemic.

Disease outbreaks in stochastic SIR epidemic models are characterized as either minor or major. When ℛ 0 < 1 , all epidemics are minor, whereas if ℛ 0 > 1 , they can be minor or major. In 1955, Whittle derived formulas for the probability of a minor or a major epidemic. A minor epidemic is distinguished from a major one in that a minor epidemic is generally of shorter duration and has substantially fewer cases than a major epidemic. In this investigation, analytical formulas are derived that approximate the probability density, the mean, and the higher-order moments for the duration of a minor epidemic. These analytical results are applicable to minor epidemics in stochastic SIR, SIS, and SIRS models with a single infected class. The probability density for minor epidemics in more complex epidemic models can be computed numerically applying multitype branching processes and the backward Kolmogorov differential equations. When ℛ 0 is close to one, minor epidemics are more common than major epidemics and their duration is significantly longer than when ℛ 0 ≪ 1 or ℛ 0 ≫ 1 .

Modeling Virus Coinfection to Inform Management of Maize Lethal Necrosis in Kenya.

Maize lethal necrosis (MLN) has emerged as a serious threat to food security in sub-Saharan Africa. MLN is caused by coinfection with two viruses, Maize chlorotic mottle virus and a potyvirus, often Sugarcane mosaic virus. To better understand the dynamics of MLN and to provide insight into disease management, we modeled the spread of the viruses causing MLN within and between growing seasons. The model allows for transmission via vectors, soil, and seed, as well as exogenous sources of infection. Following model parameterization, we predict how management affects disease prevalence and crop performance over multiple seasons. Resource-rich farmers with large holdings can achieve good control by combining clean seed and insect control. However, crop rotation is often required to effect full control. Resource-poor farmers with smaller holdings must rely on rotation and roguing, and achieve more limited control. For both types of farmer, unless management is synchronized over large areas, exogenous sources of infection can thwart control. As well as providing practical guidance, our modeling framework is potentially informative for other cropping systems in which coinfection has devastating effects. Our work also emphasizes how mathematical modeling can inform management of an emerging disease even when epidemiological information remains scanty. [Formula: see text] Copyright © 2017 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

The evolution of parasitic and mutualistic plant-virus symbioses through transmission-virulence trade-offs.

Virus-plant interactions range from parasitism to mutualism. Viruses have been shown to increase fecundity of infected plants in comparison with uninfected plants under certain environmental conditions. Increased fecundity of infected plants may benefit both the plant and the virus as seed transmission is one of the main virus transmission pathways, in addition to vector transmission. Trade-offs between vertical (seed) and horizontal (vector) transmission pathways may involve virulence, defined here as decreased fecundity in infected plants. To better understand plant-virus symbiosis evolution, we explore the ecological and evolutionary interplay of virus transmission modes when infection can lead to an increase in plant fecundity. We consider two possible trade-offs: vertical seed transmission vs infected plant fecundity, and horizontal vector transmission vs infected plant fecundity (virulence). Through mathematical models and numerical simulations, we show (1) that a trade-off between virulence and vertical transmission can lead to virus extinction during the course of evolution, (2) that evolutionary branching can occur with subsequent coexistence of mutualistic and parasitic virus strains, and (3) that mutualism can out-compete parasitism in the long-run. In passing, we show that ecological bi-stability is possible in a very simple discrete-time epidemic model. Possible extensions of this study include the evolution of conditional (environment-dependent) mutualism in plant viruses.

A primer on stochastic epidemic models: Formulation, numerical simulation, and analysis.

Some mathematical methods for formulation and numerical simulation of stochastic epidemic models are presented. Specifically, models are formulated for continuous-time Markov chains and stochastic differential equations. Some well-known examples are used for illustration such as an SIR epidemic model and a host-vector malaria model. Analytical methods for approximating the probability of a disease outbreak are also discussed.

The evolution of plant virus transmission pathways.

The evolution of plant virus transmission pathways is studied through transmission via seed, pollen, or a vector. We address the questions: under what circumstances does vector transmission make pollen transmission redundant? Can evolution lead to the coexistence of multiple virus transmission pathways? We restrict the analysis to an annual plant population in which reproduction through seed is obligatory. A semi-discrete model with pollen, seed, and vector transmission is formulated to investigate these questions. We assume vector and pollen transmission rates are frequency-dependent and density-dependent, respectively. An ecological stability analysis is performed for the semi-discrete model and used to inform an evolutionary study of trade-offs between pollen and seed versus vector transmission. Evolutionary dynamics critically depend on the shape of the trade-off functions. Assuming a trade-off between pollen and vector transmission, evolution either leads to an evolutionarily stable mix of pollen and vector transmission (concave trade-off) or there is evolutionary bi-stability (convex trade-off); the presence of pollen transmission may prevent evolution of vector transmission. Considering a trade-off between seed and vector transmission, evolutionary branching and the subsequent coexistence of pollen-borne and vector-borne strains is possible. This study contributes to the theory behind the diversity of plant-virus transmission patterns observed in nature.