A hybrid Lagrangian-Eulerian model for vector-borne diseases.

In this paper, a multi-patch and multi-group vector-borne disease model is proposed to study the effects of host commuting (Lagrangian approach) and/or vector migration (Eulerian approach) on disease spread. We first define the basic reproduction number of the model, R 0 , which completely determines the global dynamics of the model system. Namely, if R 0 ≤ 1 , then the disease-free equilibrium is globally asymptotically stable, and if R 0 > 1 , then there exists a unique endemic equilibrium which is globally asymptotically stable. Then, we show that the basic reproduction number has lower and upper bounds which are independent of the host residence times matrix and the vector migration matrix. In particular, nonhomogeneous mixing of hosts and vectors in a homogeneous environment generally increases disease persistence and the basic reproduction number of the model attains its minimum when the distributions of hosts and vectors are proportional. Moreover, R 0 can also be estimated by the basic reproduction numbers of disconnected patches if the environment is homogeneous. The optimal vector control strategy is obtained for a special scenario. In the two-patch and two-group case, we numerically analyze the dependence of the basic reproduction number and the total number of infected people on the host residence times matrix and illustrate the optimal vector control strategy in homogeneous and heterogeneous environments.

Vector-borne disease models with Lagrangian approach.

We develop a multi-group and multi-patch model to study the effects of population dispersal on the spatial spread of vector-borne diseases across a heterogeneous environment. The movement of host and/or vector is described by Lagrangian approach in which the origin or identity of each individual stays unchanged regardless of movement. The basic reproduction number [Formula: see text] of the model is defined and the strong connectivity of the host-vector network is succinctly characterized by the residence times matrices of hosts and vectors. Furthermore, the definition and criterion of the strong connectivity of general infectious disease networks are given and applied to establish the global stability of the disease-free equilibrium. The global dynamics of the model system are shown to be entirely determined by its basic reproduction number. We then obtain several biologically meaningful upper and lower bounds on the basic reproduction number which are independent or dependent of the residence times matrices. In particular, the heterogeneous mixing of hosts and vectors in a homogeneous environment always increases the basic reproduction number. There is a substantial difference on the upper bound of [Formula: see text] between Lagrangian and Eulerian modeling approaches. When only host movement between two patches is concerned, the subdivision of hosts (more host groups) can lead to a larger basic reproduction number. In addition, we numerically investigate the dependence of the basic reproduction number and the total number of infected hosts on the residence times matrix of hosts, and compare the impact of different vector control strategies on disease transmission.

Effects of behaviour change on HFMD transmission.

We propose a hand, foot and mouth disease (HFMD) transmission model for children with behaviour change and imperfect quarantine. The symptomatic and quarantined states obey constant behaviour change while others follow variable behaviour change depending on the numbers of new and recent infections. The basic reproduction number R0 of the model is defined and shown to be a threshold for disease persistence and eradication. Namely, the disease-free equilibrium is globally asymptotically stable if R0≤1 whereas the disease persists and there is a unique endemic equilibrium otherwise. By fitting the model to weekly HFMD data of Shanghai in 2019, the reproduction number is estimated at 2.41. Sensitivity analysis for R0 shows that avoiding contagious contacts and implementing strict quarantine are essential to lower HFMD persistence. Numerical simulations suggest that strong behaviour change not only reduces the peak size and endemic level dramatically but also impairs the role of asymptomatic transmission.

Relative prevalence-based dispersal in an epidemic patch model.

In this paper, we propose a two-patch SIRS model with a nonlinear incidence rate: [Formula: see text] and nonconstant dispersal rates, where the dispersal rates of susceptible and recovered individuals depend on the relative disease prevalence in two patches. In an isolated environment, the model admits Bogdanov-Takens bifurcation of codimension 3 (cusp case) and Hopf bifurcation of codimension up to 2 as the parameters vary, and exhibits rich dynamics such as multiple coexistent steady states and periodic orbits, homoclinic orbits and multitype bistability. The long-term dynamics can be classified in terms of the infection rates [Formula: see text] (due to single contact) and [Formula: see text] (due to double exposures). In a connected environment, we establish a threshold [Formula: see text] between disease extinction and uniform persistence under certain conditions. We numerically explore the effect of population dispersal on disease spread when [Formula: see text] and patch 1 has a lower infection rate, our results indicate: (i) [Formula: see text] can be nonmonotonic in dispersal rates and [Formula: see text] ([Formula: see text] is the basic reproduction number of patch i) may fail; (ii) the constant dispersal of susceptible individuals (or infective individuals) between two patches (or from patch 2 to patch 1) will increase (or reduce) the overall disease prevalence; (iii) the relative prevalence-based dispersal may reduce the overall disease prevalence. When [Formula: see text] and the disease outbreaks periodically in each isolated patch, we find that: (a) small unidirectional and constant dispersal can lead to complex periodic patterns like relaxation oscillations or mixed-mode oscillations, whereas large ones can make the disease go extinct in one patch and persist in the form of a positive steady state or a periodic solution in the other patch; (b) relative prevalence-based and unidirectional dispersal can make periodic outbreak earlier.

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.

Total biomass of a single population in two-patch environments.

For the two-patch logistic model, we study the effect of dispersal intensity and dispersal asymmetry on the total population abundance and its distribution. Two complete classifications of the model parameter space are given: one concerning when dispersal causes smaller or larger total biomass than no dispersal, and the other addressing how the total biomass changes with dispersal intensity and dispersal asymmetry. The dependencies of the population abundance of each individual patch on dispersal intensity and dispersal asymmetry are also fully characterized. In addition, the maximal and minimal total population sizes induced by dispersal are determined for the logistic model with an arbitrary number of patches, and a weak order-preserving result correlated the local population abundances with and without dispersal is established.

Modelling COVID-19 outbreak on the Diamond Princess ship using the public surveillance data.

The novel coronavirus disease 2019 (COVID-19) outbreak on the Diamond Princess (DP) ship has caused over 634 cases as of February 20, 2020. We model the transmission process on DP ship as a stochastic branching process, and estimate the reproduction number at the innitial phase of 2.9 (95%CrI: 1.7-7.7). The epidemic doubling time is 3.4 days, and thus timely actions on COVID-19 control were crucial. We estimate the COVID-19 transmissibility reduced 34% after the quarantine program on the DP ship which was implemented on February 5. According to the model simulation, relocating the population at risk may sustainably decrease the epidemic size, postpone the timing of epidemic peak, and thus relieve the tensive demands in the healthcare. The lesson learnt on the ship should be considered in other similar settings.

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.

Shrinkage in serial intervals across transmission generations of COVID-19.

One of the key epidemiological characteristics that shape the transmission of coronavirus disease 2019 (COVID-19) is the serial interval (SI). Although SI is commonly considered following a probability distribution at a population scale, recent studies reported a slight shrinkage (or contraction) of the mean of effective SI across transmission generations or over time. Here, we develop a likelihood-based statistical inference framework with truncation to explore the change in SI across transmission generations after adjusting the impacts of case isolation. The COVID-19 contact tracing surveillance data in Hong Kong are used for exemplification. We find that for COVID-19, the mean of individual SI is likely to shrink with a factor at 0.72 per generation (95%CI: 0.54, 0.96) as the transmission generation increases, where a threshold may exist as the lower boundary of this shrinking process. We speculate that one of the probable explanations for the shrinkage in SI might be an outcome due to the competition among multiple candidate infectors within the same case cluster. Thus, the nonpharmaceutical interventive strategies are crucially important to block the transmission chains, and mitigate the COVID-19 epidemic.

Impact of State-Dependent Dispersal on Disease Prevalence.

Based on a susceptible-infected-susceptible patch model, we study the influence of dispersal on the disease prevalence of an individual patch and all patches at the endemic equilibrium. Specifically, we estimate the disease prevalence of each patch and obtain a weak order-preserving result that correlated the patch reproduction number with the patch disease prevalence. Then we assume that dispersal rates of the susceptible and infected populations are proportional and derive the overall disease prevalence, or equivalently, the total infection size at no dispersal or infinite dispersal as well as the right derivative of the total infection size at no dispersal. Furthermore, for the two-patch submodel, two complete classifications of the model parameter space are given: one addressing when dispersal leads to higher or lower overall disease prevalence than no dispersal, and the other concerning how the overall disease prevalence varies with dispersal rate. Numerical simulations are performed to further investigate the effect of movement on disease prevalence.

Comparing COVID-19 and the 1918-19 influenza pandemics in the United Kingdom.

We compared the COVID-19 and 1918-19 influenza pandemics in the United Kingdom. We found that the ongoing COVID-19 wave of infection matched the major wave of the 1918-19 influenza pandemic surprisingly well, with both reaching similar magnitudes (in terms of estimated weekly new infections) and spending the same duration with over five cases per 1000 inhabitants over the previous two months. We also discussed the similarities in epidemiological characteristics between these two pandemics.

Imitation dynamics in the mitigation of the novel coronavirus disease (COVID-19) outbreak in Wuhan, China from 2019 to 2020.

BACKGROUND: The coronavirus disease 2019 (COVID-19) was first identified in Wuhan, China on December 2019 in patients presenting with atypical pneumonia. Although 'city-lockdown' policy reduced the spatial spreading of the COVID-19, the city-level outbreaks within each city remain a major concern to be addressed. The local or regional level disease control mainly depends on individuals self-administered infection prevention actions. The contradiction between choice of taking infection prevention actions or not makes the elimination difficult under a voluntary acting scheme, and represents a clash between the optimal choice of action for the individual interest and group interests. METHODS: We develop a compartmental epidemic model based on the classic susceptible-exposed-infectious-recovered model and use this to fit the data. Behavioral imitation through a game theoretical decision-making process is incorporated to study and project the dynamics of the COVID-19 outbreak in Wuhan, China. By varying the key model parameters, we explore the probable course of the outbreak in terms of size and timing under several public interventions in improving public awareness and sensitivity to the infection risk as well as their potential impact. RESULTS: We estimate the basic reproduction number, R 0, to be 2.5 (95% CI: 2.4-2.7). Under the current most realistic setting, we estimate the peak size at 0.28 (95% CI: 0.24-0.32) infections per 1,000 population. In Wuhan, the final size of the outbreak is likely to infect 1.35% (95% CI: 1.00-2.12%) of the population. The outbreak will be most likely to peak in the first half of February and drop to daily incidences lower than 10 in June 2020. Increasing sensitivity to take infection prevention actions and the effectiveness of infection prevention measures are likely to mitigate the COVID-19 outbreak in Wuhan. CONCLUSIONS: Through an imitating social learning process, individual-level behavioral change on taking infection prevention actions have the potentials to significantly reduce the COVID-19 outbreak in terms of size and timing at city-level. Timely and substantially resources and supports for improving the willingness-to-act and conducts of self-administered infection prevention actions are recommended to reduce to the COVID-19 associated risks.

Quantifying the improvement in confirmation efficiency of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the early phase of the outbreak in Hong Kong in 2020.

BACKGROUNDS: The emerging virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), caused a large outbreak of coronavirus disease, COVID-19, in Wuhan, China, since December 2019. COVID-19 soon spread to other regions of China and overseas. In Hong Kong, local mitigation measures were implemented since the first imported case was confirmed on January 23, 2020. Here we evaluated the temporal variation of detection delay from symptoms onset to laboratory confirmation of SARS-CoV-2 in Hong Kong. METHODS: A regression model is adopted to quantify the association between the SARS-CoV-2 detection delay and calendar time. The association is tested and further validated by a Cox proportional hazard model. FINDINGS: The estimated median detection delay was 9.5 days (95%CI: 6.5-11.5) in the second half of January, reduced to 6.0 days (95%CI: 5.5-9.5) in the first half of February 2020. We estimate that SARS-CoV-2 detection efficiency improved at a daily rate of 5.40% (95%CI: 2.54-8.33) in Hong Kong. CONCLUSIONS: The detection efficiency of SARS-CoV-2 was likely being improved substantially in Hong Kong since the first imported case was detected. Sustaining enforcement in timely detection and other effective control measures are recommended to prevent the SARS-CoV-2 infection.

Mathematical Analysis of the Ross-Macdonald Model with Quarantine.

People infected with malaria may receive less mosquito bites when they are treated in well-equipped hospitals or follow doctors' advice for reducing exposure to mosquitoes at home. This quarantine-like intervention measure is especially feasible in countries and areas approaching malaria elimination. Motivated by mathematical models with quarantine for directly transmitted diseases, we develop a mosquito-borne disease model where imperfect quarantine is considered to mitigate the disease transmission from infected humans to susceptible mosquitoes. The basic reproduction number [Formula: see text] is computed and the model equilibria and their stabilities are analyzed when the incidence rate is standard or bilinear. In particular, the model system may undergo a subcritical (backward) bifurcation at [Formula: see text] when standard incidence is adopted, whereas the disease-free equilibrium is globally asymptotically stable as [Formula: see text] and the unique endemic equilibrium is locally asymptotically stable as [Formula: see text] when the infection incidence is bilinear. Numerical simulations suggest that the quarantine strategy can play an important role in decreasing malaria transmission. The success of quarantine mainly relies on the reduction of bites on quarantined individuals.

A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action.

The ongoing coronavirus disease 2019 (COVID-19) outbreak, emerged in Wuhan, China in the end of 2019, has claimed more than 2600 lives as of 24 February 2020 and posed a huge threat to global public health. The Chinese government has implemented control measures including setting up special hospitals and travel restriction to mitigate the spread. We propose conceptual models for the COVID-19 outbreak in Wuhan with the consideration of individual behavioural reaction and governmental actions, e.g., holiday extension, travel restriction, hospitalisation and quarantine. We employe the estimates of these two key components from the 1918 influenza pandemic in London, United Kingdom, incorporated zoonotic introductions and the emigration, and then compute future trends and the reporting ratio. The model is concise in structure, and it successfully captures the course of the COVID-19 outbreak, and thus sheds light on understanding the trends of the outbreak.

Modelling triatomine bug population and Trypanosoma rangeli transmission dynamics: Co-feeding, pathogenic effect and linkage with chagas disease.

Trypanosoma rangeli (T. rangeli), a parasite, is not pathogenic to human but pathogenic to some vector species to induce the behavior changes of infected vectors and subsequently impact the transmission dynamics of other diseases such as Chagas disease which shares the same vector species. Here we develop a mathematical model and conduct qualitative analysis for the transmission dynamics of T. rangeli. We incorporate both systemic and co-feeding transmission routes, and account for the pathogenic effect using infection-induced fecundity and fertility change of the triatomine bugs. We derive two thresholds Rv (the triatomine bug basic reproduction number) and R0 (the T. rangeli basic reproduction number) to delineate the dynamical behaviors of the ecological and epidemiological systems. We show that when Rv>1 and R0>1, a unique parasite positive equilibrium E* appears. We find that E* can be unstable and periodic oscillations can be observed where the pathogenic effect plays a significant role. Implications of the qualitative analysis and numerical simulations suggest the need of an integrative vector-borne disease prevention and control strategy when multiple vector-borne diseases are transmitted by the same set of vector species.