Modeling insect growth regulators for pest management.

Insect growth regulators (IGRs) have been developed as effective control measures against harmful insect pests to disrupt their normal development. This study is to propose a mathematical model to evaluate the cost-effectiveness of IGRs for pest management. The key features of the model include the temperature-dependent growth of insects and realistic impulsive IGRs releasing strategies. The impulsive releases are carefully modeled by counting the number of implements during an insect's temperature-dependent development duration, which introduces a surviving probability determined by a product of terms corresponding to each release. Dynamical behavior of the model is illustrated through dynamical system analysis and a threshold-type result is established in terms of the net reproduction number. Further numerical simulations are performed to quantitatively evaluate the effectiveness of IGRs to control populations of harmful insect pests. It is interesting to observe that the time-changing environment plays an important role in determining an optimal pest control scheme with appropriate release frequencies and time instants.

Getting Jab or Regular Test: Observations from an Impulsive Epidemic COVID-19 Model.

Several safe and effective vaccines are available to prevent individuals from experiencing severe illness or death as a result of COVID-19. Widespread vaccination is widely regarded as a critical tool in the fight against the disease. However, some individuals may choose not to vaccinate due to vaccine hesitancy or other medical conditions. In some sectors, regular compulsory testing is required for such unvaccinated individuals. Interestingly, different sectors require testing at various frequencies, such as weekly or biweekly. As a result, it is essential to determine the optimal testing frequency and identify underlying factors. This study proposes a population-based model that can accommodate different personal decision choices, such as getting vaccinated or undergoing regular tests, as well as vaccine efficacies and uncertainties in epidemic transmission. The model, formulated as impulsive differential equations, uses time instants to represent the reporting date for the test result of an unvaccinated individual. By employing well-accepted indices to measure transmission risk, including the basic reproduction number, the peak time, the final size, and the number of severe infections, the study shows that an optimal testing frequency is highly sensitive to parameters involved in the transmission process, such as vaccine efficacy, disease transmission rate, test accuracy, and existing vaccination coverage. The testing frequency should be appropriately designed with the consideration of all these factors, as well as the control objectives measured by epidemiological quantities of great concern.

Evaluation of Effectiveness of Global COVID-19 Vaccination Campaign.

To model estimated deaths averted by COVID-19 vaccines, we used state-of-the-art mathematical modeling, likelihood-based inference, and reported COVID-19 death and vaccination data. We estimated that >1.5 million deaths were averted in 12 countries. Our model can help assess effectiveness of the vaccination program, which is crucial for curbing the COVID-19 pandemic.

A Zika Endemic Model for the Contribution of Multiple Transmission Routes.

Zika virus disease is a viral disease primarily transmitted to humans through the bite of infected female mosquitoes. Recent evidence indicates that the virus can also be sexually transmitted in hosts and vertically transmitted in vectors. In this paper, we propose a Zika model with three transmission routes, that is, vector-borne transmission between humans and mosquitoes, sexual transmission within humans and vertical transmission within mosquitoes. The basic reproduction number [Formula: see text] is computed and shown to be a sharp threshold quantity. Namely, the disease-free equilibrium is globally asymptotically stable as [Formula: see text], whereas there exists a unique endemic equilibrium which is globally asymptotically stable as [Formula: see text]. The relative contributions of each transmission route on the reproduction number, and the short- and long-term host infections are analyzed. Numerical simulations confirm that vectorial transmission contributes the most to the initial and subsequent transmission. The role of sexual transmission in the early phase of a Zika outbreak is greater than the long term, while vertical transmission is the opposite. Reducing mosquito bites is the most effective measure in lowering the risk of Zika virus infection.

Dynamics of a periodic tick-borne disease model with co-feeding and multiple patches.

By extending a mechanistic model for the tick-borne pathogen systemic transmission with the consideration of seasonal climate impacts, host movement as well as the co-feeding transmission route, this paper proposes a novel modeling framework for describing the spatial dynamics of tick-borne diseases. The net reproduction number for tick growth and basic reproduction number for disease transmission are derived, which predict the global dynamics of tick population growth and disease transmission. Numerical simulations not only verify the analytical results, but also characterize the contribution of co-feeding transmission route on disease prevalence in a habitat and the effect of host movement on the spatial spreading of the pathogen.

Modelling epidemics with fractional-dose vaccination in response to limited vaccine supply.

The control strategies of emergency infectious diseases are constrained by limited medical resources. The fractional dose vaccination strategy as one of feasible strategies was proposed in response to global shortages of vaccine stockpiles. Although a variety of epidemic models have been developed under the circumstances of limited resources in treatment, few models particularly investigated vaccination strategies in resource-limited settings. In this paper, we develop a two-group SIR model with incorporation of proportionate mixing patterns and n-fold fractional dose vaccination related parameters to evaluate the efficiency of fractional dose vaccination on disease control at the population level. The existence and uniqueness of the final size of the two-group SIR epidemic model, the formulation of the basic reproduction number and the relationship between them are established. Moreover, numerical simulations are performed based on this two-group vector-free model to investigate the effectiveness of n-fold fractional dose vaccination by using the emergency outbreaks of yellow fever in Angola in 2016. By employing linear and nonlinear dose-response relationships, we compare the resulting fluctuations of four characteristics of the epidemics, which are the outbreak size, the peak time of the outbreak, the basic reproduction number and the infection attack rate (IAR). For both types of dose-response relationships, dose-fractionation takes positive effects in lowering the outbreak size, delay the peak time of the outbreak, reducing the basic reproduction number and the IAR of yellow fever only when the vaccine efficacy is high enough. Moreover, five-fold fractional dose vaccination strategy may not be the optimal vaccination strategy as proposed by the World Health Organization if the dose-response relationship is nonlinear.

Low dispersion in the infectiousness of COVID-19 cases implies difficulty in control.

The individual infectiousness of coronavirus disease 2019 (COVID-19), quantified by the number of secondary cases of a typical index case, is conventionally modelled by a negative-binomial (NB) distribution. Based on patient data of 9120 confirmed cases in China, we calculated the variation of the individual infectiousness, i.e., the dispersion parameter k of the NB distribution, at 0.70 (95% confidence interval: 0.59, 0.98). This suggests that the dispersion in the individual infectiousness is probably low, thus COVID-19 infection is relatively easy to sustain in the population and more challenging to control. Instead of focusing on the much fewer super spreading events, we also need to focus on almost every case to effectively reduce transmission.

Modelling diapause in mosquito population growth.

Diapause, a period of arrested development caused by adverse environmental conditions, serves as a key survival mechanism for insects and other invertebrate organisms in temperate and subtropical areas. In this paper, a novel modelling framework, motivated by mosquito species, is proposed to investigate the effects of diapause on seasonal population growth, where the diapause period is taken as an independent growth process, during which the population dynamics are completely different from that in the normal developmental and post-diapause periods. More specifically, the annual growth period is divided into three intervals, and the population dynamics during each interval are described by different sets of equations. We formulate two models of delay differential equations (DDE) to explicitly describe mosquito population growth with a single diapausing stage, either immature or adult. These two models can be further unified into one DDE model, on which the well-posedness of the solutions and the global stability of the trivial and positive periodic solutions in terms of an index [Formula: see text] are analysed. The seasonal population abundances of two temperate mosquito species with different diapausing stages are simulated to identify the essential role on population persistence that diapause plays and predict that killing mosquitoes during the diapause period can lower but fail to prevent the occurrence of peak abundance in the following season. Instead, culling mosquitoes during the normal growth period is much more efficient to decrease the outbreak size. Our modelling framework may shed light on the diapause-induced variations in spatiotemporal distributions of different mosquito species.

Modelling the skip-and-resurgence of Japanese encephalitis epidemics in Hong Kong.

Japanese encephalitis virus (JEV) is a zoonotic mosquito-borne virus, persisting in pigs, Ardeid birds and Culex mosquitoes. It is endemic to China and Southeastern Asia. The case-fatality ratio (CFR) or the rate of permanent psychiatric sequelae is 30% among symptomatic patients. There were no reported local JEV human cases between 2006 and 2010 in Hong Kong, but it was followed by a resurgence of cases from 2011 to 2017. The mechanism behind this "skip-and-resurgence" patterns is unclear. This work aims to reveal the mechanism behind the "skip-and-resurgence" patterns using mathematical modelling and likelihood-based inference techniques. We found that pig-to-pig transmission increases the size of JEV epidemics but is unlikely to maintain the same level of transmission among pigs. The disappearance of JEV human cases in 2006-2010 could be explained by a sudden reduction of the population of farm pigs as a result of the implementation of the voluntary "pig-rearing licence surrendering" policy. The resurgence could be explained by of a new strain in 2011, which increased the transmissibility of the virus or the spill-over ratio from reservoir to host or both.

Behavioral synchronization induced by epidemic spread in complex networks.

During the spread of an epidemic, individuals in realistic networks may exhibit collective behaviors. In order to characterize this kind of phenomenon and explore the correlation between collective behaviors and epidemic spread, in this paper, we construct several mathematical models (including without delay, with a coupling delay, and with double delays) of epidemic synchronization by applying the adaptive feedback motivated by real observations. By using Lyapunov function methods, we obtain the conditions for local and global stability of these epidemic synchronization models. Then, we illustrate that quenched mean-field theory is more accurate than heterogeneous mean-field theory in the prediction of epidemic synchronization. Finally, some numerical simulations are performed to complement our theoretical results, which also reveal some unexpected phenomena, for example, the coupling delay and epidemic delay influence the speed of epidemic synchronization. This work makes further exploration on the relationship between epidemic dynamics and synchronization dynamics, in the hope of being helpful to the study of other dynamical phenomena in the process of epidemic spread.

Local immunization program for susceptible-infected-recovered network epidemic model.

The immunization strategies through contact tracing on the susceptible-infected-recovered framework in social networks are modelled to evaluate the cost-effectiveness of information-based vaccination programs with particular focus on the scenario where individuals belonging to a specific set can get vaccinated due to the vaccine shortages and other economic or humanity constraints. By using the block heterogeneous mean-field approach, a series of discrete-time dynamical models is formulated and the condition for epidemic outbreaks can be established which is shown to be not only dependent on the network structure but also closely related to the immunization control parameters. Results show that increasing the immunization strength can effectively raise the epidemic threshold, which is different from the predictions obtained through the susceptible-infected-susceptible network framework, where epidemic threshold is independent of the vaccination strength. Furthermore, a significant decrease of vaccine use to control the infectious disease is observed for the local vaccination strategy, which shows the promising applications of the local immunization programs to disease control while calls for accurate local information during the process of disease outbreak.

Impact of biodiversity and seasonality on Lyme-pathogen transmission.

Lyme disease imposes increasing global public health challenges. To better understand the joint effects of seasonal temperature variation and host community composition on the pathogen transmission, a stage-structured periodic model is proposed by integrating seasonal tick development and activity, multiple host species and complex pathogen transmission routes between ticks and reservoirs. Two thresholds, one for tick population dynamics and the other for Lyme-pathogen transmission dynamics, are identified and shown to fully classify the long-term outcomes of the tick invasion and disease persistence. Seeding with the realistic parameters, the tick reproduction threshold and Lyme disease spread threshold are estimated to illustrate the joint effects of the climate change and host community diversity on the pattern of Lyme disease risk. It is shown that climate warming can amplify the disease risk and slightly change the seasonality of disease risk. Both the "dilution effect" and "amplification effect" are observed by feeding the model with different possible alternative hosts. Therefore, the relationship between the host community biodiversity and disease risk varies, calling for more accurate measurements on the local environment, both biotic and abiotic such as the temperature and the host community composition.

A PERIODIC ROSS-MACDONALD MODEL IN A PATCHY ENVIRONMENT.

Based on the classical Ross-Macdonald model, in this paper we propose a periodic malaria model to incorporate the effects of temporal and spatial heterogeneity on disease transmission. The temporal heterogeneity is described by assuming that some model coefficients are time-periodic, while the spatial heterogeneity is modeled by using a multi-patch structure and assuming that individuals travel among patches. We calculate the basic reproduction number [Formula: see text] and show that either the disease-free periodic solution is globally asymptotically stable if [Formula: see text] or the positive periodic solution is globally asymptotically stable if [Formula: see text]. Numerical simulations are conducted to confirm the analytical results and explore the effect of travel control on the disease prevalence.

Effectiveness evaluation of mosquito suppression strategies on dengue transmission under changing temperature and precipitation.

Widespread resurgence of dengue outbreaks has seriously threatened the global health. Due to lack of treatments and vaccines, one key strategy in dengue control is to reduce the vector population size. As an environment-friendly mosquito control approach, releasing male mosquitoes transinfected with specific Wolbachia strain into the field to suppress the wild mosquito population size has become wildly accepted. The current study evaluates the effectiveness of this suppression strategy on dengue control under changing temperature and precipitation profiles. We formulate a mathematical model which includes larval intra-specific competition, the maturation period for mosquitoes, the extrinsic incubation period (EIP) and intrinsic incubation period (IIP). The persistence of mosquitoes and disease is discussed in terms of two basic reproduction numbers (RM and R0) and the release ratio pw. Further numerical simulations are carried out to not only validate theoretical results, but also provide interesting quantitative observations. Sensitivity analysis on the reproduction numbers, peak size, peak time and the final epidemic size is performed with respect to model parameters, which highlights effective control measures against dengue transmission. Moreover, by assuming temperature and precipitation dependent mosquito-related parameters, the model can be used to project the effectiveness of releasing Wolbachia-carrying males under climatic variations. It is shown that the effectiveness of various control strategies is highly dependent on the changing temperature and precipitation profiles. In particular, the model projects that it is most challenging to control the disease at the favorable temperature (around 27∼30∘C) and precipitation (5∼8mm/day) range, during which the basic reproduction number R0 is very high and more Wolbachia-infected males should be released.

A story of viral co-infection, co-transmission and co-feeding in ticks: how to compute an invasion reproduction number.

With a single circulating vector-borne virus, the basic reproduction number incorporates contributions from tick-to-tick (co-feeding), tick-to-host and host-to-tick transmission routes. With two different circulating vector-borne viral strains, resident and invasive, and under the assumption that co-feeding is the only transmission route in a tick population, the invasion reproduction number depends on whether the model system of ordinary differential equations possesses the property of neutrality. We show that a simple model, with two populations of ticks infected with one strain, resident or invasive, and one population of co-infected ticks, does not have Alizon's neutrality property. We present model alternatives that are capable of representing the invasion potential of a novel strain by including populations of ticks dually infected with the same strain. The invasion reproduction number is analysed with the next-generation method and via numerical simulations.

Developing a temperature-driven map of the basic reproductive number of the emerging tick vector of Lyme disease Ixodes scapularis in Canada.

A mechanistic model of the tick vector of Lyme disease, Ixodes scapularis, was adapted to a deterministic structure. Using temperature normals smoothed by Fourier analysis to generate seasonal temperature-driven development rates and host biting rates, and a next generation matrix approach, the model was used to obtain values for the basic reproduction number (R(0)) for I. scapularis at locations in southern Canada where the tick is established and emerging. The R(0) at Long Point, Point Pelee and Chatham sites where I. scapularis are established, was estimated at 1.5, 3.19 and 3.65, respectively. The threshold temperature conditions for tick population survival (R(0)=1) were shown to be the same as those identified using the mechanistic model (2800-3100 cumulative annual degree days >0°C), and a map of R(0) for I. scapularis, the first such map for an arthropod vector, was drawn for Canada east of the Rocky Mountains. This map supports current risk assessments for Lyme disease risk emergence in Canada. Sensitivity analysis identified host abundance, tick development rates and summer temperatures as highly influential variables in the model, which is consistent with our current knowledge of the biology of this tick. The development of a deterministic model for I. scapularis that is capable of providing values for R(0) is a key step in our evolving ability to develop tools for assessment of Lyme disease risk emergence and for development of public health policies on surveillance, prevention and control.

A Reaction-Diffusion Model with Spatially Inhomogeneous Delays.

Motivated by population growth in a heterogeneous environment, this manuscript builds a reaction-diffusion model with spatially dependent parameters. In particular, a term for spatially uneven maturation durations is included in the model, which puts the current investigation among the very few studies on reaction-diffusion systems with spatially dependent delays. Rigorous analysis is performed, including the well-posedness of the model, the basic reproduction ratio formulation and long-term behavior of solutions. Under mild assumptions on model parameters, extinction of the species is predicted when the basic reproduction ratio is less than one. When the birth rate is an increasing function and the basic reproduction ratio is greater than one, uniqueness and global attractivity of a positive equilibrium can be established with the help of a novel functional phase space. Permanence of the species is shown when the birth function is in a unimodal form and the basic reproduction ratio is greater than one. The synthesized approach proposed here is applicable to broader contexts of studies on the impact of spatial heterogeneity on population dynamics, in particular, when the delayed feedbacks are involved and the response time is spatially varying.

Threshold virus dynamics with impulsive antiretroviral drug effects.

The purposes of this paper are twofold: to develop a rigorous approach to analyze the threshold behaviors of nonlinear virus dynamics models with impulsive drug effects and to examine the feasibility of virus clearance following the Manuals of National AIDS Free Antiviral Treatment in China. An impulsive system of differential equations is developed to describe the within-host virus dynamics of both wild-type and drug-resistant strains when a combination of antiretroviral drugs is used to induce instantaneous drug effects at a sequence of dosing times equally spaced while drug concentrations decay exponentially after the dosing time. Threshold parameters are derived using the basic reproduction number of periodic epidemic models, and are used to depict virus clearance/persistence scenarios using the theory of asymptotic periodic systems and the persistence theory of discrete dynamical systems. Numerical simulations using model systems parametrized in terms of the antiretroviral therapy recommended in the aforementioned Manuals illustrate the theoretical threshold virus dynamics, and examine conditions under which the impulsive antiretroviral therapy leads to treatment success. In particular, our results show that only the drug-resistant strain can dominate (the first-line treatment program guided by the Manuals) or both strains may be rapidly eliminated (the second-line treatment program), thus the work indicates the importance of implementing the second-line treatment program as soon as possible.

Stage duration distributions and intraspecific competition: a review of continuous stage-structured models.

Stage structured models, by grouping individuals with similar demographic characteristics together, have proven useful in describing population dynamics. This manuscript starts from reviewing two widely used modeling frameworks that are in the form of integral equations and age-structured partial differential equations. Both modeling frameworks can be reduced to the same differential equation structures with/without time delays by applying Dirac and gamma distributions for the stage durations. Each framework has its advantages and inherent limitations. The net reproduction number and initial growth rate can be easily defined from the integral equation. However, it becomes challenging to integrate the density-dependent regulations on the stage distribution and survival probabilities in an integral equation, which may be suitably incorporated into partial differential equations. Further recent modeling studies, in particular those by Stephen A. Gourley and collaborators, are reviewed under the conditions of the stage duration distribution and survival probability being regulated by population density.