Re-assessing thermal response of schistosomiasis transmission risk: evidence for a higher thermal optimum than previously predicted.

The geographical range of schistosomiasis is affected by the ecology of schistosome parasites and their obligate host snails, including their response to temperature. Previous models predicted schistosomiasis' thermal optimum at 21.7 °C, which is not compatible with the temperature in sub-Saharan Africa (SSA) regions where schistosomiasis is hyperendemic. We performed an extensive literature search for empirical data on the effect of temperature on physiological and epidemiological parameters regulating the free-living stages of S. mansoni and S. haematobium and their obligate host snails, i.e., Biomphalaria spp. and Bulinus spp., respectively. We derived nonlinear thermal responses fitted on these data to parameterize a mechanistic, process-based model of schistosomiasis. We then re-cast the basic reproduction number and the prevalence of schistosome infection as functions of temperature. We found that the thermal optima for transmission of S. mansoni and S. haematobium range between 23.1-27.3 °C and 23.6-27.9 °C (95 % CI) respectively. We also found that the thermal optimum shifts toward higher temperatures as the human water contact rate increases with temperature. Our findings align with an extensive dataset of schistosomiasis prevalence in SSA. The refined nonlinear thermal-response model developed here suggests a more suitable current climate and a greater risk of increased transmission with future warming for more than half of the schistosomiasis suitable regions with mean annual temperature below the thermal optimum.

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.

Parameter identifiability of a within-host SARS-CoV-2 epidemic model.

Parameter identification involves the estimation of undisclosed parameters within a system based on observed data and mathematical models. In this investigation, we employ DAISY to meticulously examine the structural identifiability of parameters of a within-host SARS-CoV-2 epidemic model, taking into account an array of observable datasets. Furthermore, Monte Carlo simulations are performed to offer a comprehensive practical analysis of model parameters. Lastly, sensitivity analysis is employed to ascertain that decreasing the replication rate of the SARS-CoV-2 virus and curbing the infectious period are the most efficacious measures in alleviating the dissemination of COVID-19 amongst hosts.

The impact of temperature, humidity and closing school on the mumps epidemic: a case study in the mainland of China.

BACKGROUND: To control resurging infectious diseases like mumps, it is necessary to resort to effective control and preventive measures. These measures include increasing vaccine coverage, providing the community with advice on how to reduce exposure, and closing schools. To justify such intervention, it is important to understand how well each of these measures helps to limit transmission. METHODS: In this paper, we propose a simple SEILR (susceptible-exposed-symptomatically infectious-asymptomatically infectious-recovered) model by using a novel transmission rate function to incorporate temperature, humidity, and closing school factors. This new transmission rate function allows us to verify the impact of each factor either separately or combined. Using reported mumps cases from 2004 to 2018 in the mainland of China, we perform data fitting and parameter estimation to evaluate the basic reproduction number  R 0 . As a wide range of one-dose measles, mumps, and rubella (MMR) vaccine programs in China started only in 2008, we use different vaccination proportions for the first Stage I period (from 2004 to 2008) and the second Stage II period (from 2009 to 2018). This allows us to verify the importance of higher vaccine coverage with a possible second dose of MMR vaccine. RESULTS: We find that the basic reproduction number  R 0  is generally between 1 and 3. We then use the Akaike Information Criteria to assess the extent to which each of the three factors contributed to the spread of mumps. The findings suggest that the impact of all three factors is substantial, with temperature having the most significant impact, followed by school opening and closing, and finally humidity. CONCLUSION: We conclude that the strategy of increasing vaccine coverage, changing micro-climate (temperature and humidity), and closing schools can greatly reduce mumps transmission.

Mathematical analysis and prediction of future outbreak of dengue on time-varying contact rate using machine learning approach.

This article introduces a novel mathematical model analyzing the dynamics of Dengue in the recent past, specifically focusing on the 2023 outbreak of this disease. The model explores the patterns and behaviors of dengue fever in Bangladesh. Incorporating a sinusoidal function reveals significant mid-May to Late October outbreak predictions, aligning with the government's exposed data in our simulation. For different amplitudes (A) within a sequence of values (A = 0.1 to 0.5), the highest number of infected mosquitoes occurs in July. However, simulations project that when ßM = 0.5 and A = 0.1, the peak of human infections occurs in late September. Not only the next-generation matrix approach along with the stability of disease-free and endemic equilibrium points are observed, but also a cutting-edge Machine learning (ML) approach such as the Prophet model is explored for forecasting future Dengue outbreaks in Bangladesh. Remarkably, we have fitted our solution curve of infection with the reported data by the government of Bangladesh. We can predict the outcome of 2024 based on the ML Prophet model situation of Dengue will be detrimental and proliferate 25 % compared to 2023. Finally, the study marks a significant milestone in understanding and managing Dengue outbreaks in Bangladesh.

Approximating reproduction numbers: a general numerical method for age-structured models.

In this paper, we introduce a general numerical method to approximate the reproduction numbers of a large class of multi-group, age-structured, population models with a finite age span. To provide complete flexibility in the definition of the birth and transition processes, we propose an equivalent formulation for the age-integrated state within the extended space framework. Then, we discretize the birth and transition operators via pseudospectral collocation. We discuss applications to epidemic models with continuous and piecewise continuous rates, with different interpretations of the age variable (e.g., demographic age, infection age and disease age) and the transmission terms (e.g., horizontal and vertical transmission). The tests illustrate that the method can compute different reproduction numbers, including the basic and type reproduction numbers as special cases.

Prediction and control of cholera outbreak: Study case of Cameroon.

This paper deals with the problem of the prediction and control of cholera outbreak using real data of Cameroon. We first develop and analyze a deterministic model with seasonality for the cholera, the novelty of which lies in the incorporation of undetected cases. We present the basic properties of the model and compute two explicit threshold parameters R¯0 and R_0 that bound the effective reproduction number R0, from below and above, that is R_0≤R0≤R¯0. We prove that cholera tends to disappear when R¯0≤1, while when R_0>1, cholera persists uniformly within the population. After, assuming that the cholera transmission rates and the proportions of newly symptomatic are unknown, we develop the EnKf approach to estimate unmeasurable state variables and these unknown parameters using real data of cholera from 2014 to 2022 in Cameroon. We use this result to estimate the upper and lower bound of the effective reproduction number and reconstructed active asymptomatic and symptomatic cholera cases in Cameroon, and give a short-term forecasts of cholera in Cameroon until 2024. Numerical simulations show that (i) the transmission rate from free Vibrio cholerae in the environment is more important than the human transmission and begin to be high few week after May and in October, (ii) 90% of newly cholera infected cases that present the symptoms of cholera are not diagnosed and (iii) 60.36% of asymptomatic are detected at 14% and 86% of them recover naturally. The future trends reveals that an outbreak appeared from July to November 2023 with the number of cases reported monthly peaked in October 2023. An impulsive control strategy is incorporated in the model with the aim to avoid or prevent the cholera outbreak. In the first year of monitoring, we observed a reduction of more than 75% of incidences and the disappearance of the peaks when no control are available in Cameroon. A second monitoring of control led to a further reduction of around 60% of incidences the following year, showing how impulse control could be an effective means of eradicating cholera.

Threshold dynamics of a reaction-advection-diffusion schistosomiasis epidemic model with seasonality and spatial heterogeneity.

Most water-borne disease models ignore the advection of water flows in order to simplify the mathematical analysis and numerical computation. However, advection can play an important role in determining the disease transmission dynamics. In this paper, we investigate the long-term dynamics of a periodic reaction-advection-diffusion schistosomiasis model and explore the joint impact of advection, seasonality and spatial heterogeneity on the transmission of the disease. We derive the basic reproduction number R 0 and show that the disease-free periodic solution is globally attractive when R 0 < 1 whereas there is a positive endemic periodic solution and the system is uniformly persistent in a special case when R 0 > 1 . Moreover, we find that R 0 is a decreasing function of the advection coefficients which offers insights into why schistosomiasis is more serious in regions with slow water flows.

Modelling phytoplankton-virus interactions: phytoplankton blooms and lytic virus transmission.

A dynamic reaction-diffusion model of four variables is proposed to describe the spread of lytic viruses among phytoplankton in a poorly mixed aquatic environment. The basic ecological reproductive index for phytoplankton invasion and the basic reproduction number for virus transmission are derived to characterize the phytoplankton growth and virus transmission dynamics. The theoretical and numerical results from the model show that the spread of lytic viruses effectively controls phytoplankton blooms. This validates the observations and experimental results of Emiliana huxleyi-lytic virus interactions. The studies also indicate that the lytic virus transmission cannot occur in a low-light or oligotrophic aquatic environment.

Genotypic-phenotypic landscape computation based on first principle and deep learning.

The relationship between genotype and fitness is fundamental to evolution, but quantitatively mapping genotypes to fitness has remained challenging. We propose the Phenotypic-Embedding theorem (P-E theorem) that bridges genotype-phenotype through an encoder-decoder deep learning framework. Inspired by this, we proposed a more general first principle for correlating genotype-phenotype, and the P-E theorem provides a computable basis for the application of first principle. As an application example of the P-E theorem, we developed the Co-attention based Transformer model to bridge Genotype and Fitness model, a Transformer-based pre-train foundation model with downstream supervised fine-tuning that can accurately simulate the neutral evolution of viruses and predict immune escape mutations. Accordingly, following the calculation path of the P-E theorem, we accurately obtained the basic reproduction number (${R}_0$) of SARS-CoV-2 from first principles, quantitatively linked immune escape to viral fitness and plotted the genotype-fitness landscape. The theoretical system we established provides a general and interpretable method to construct genotype-phenotype landscapes, providing a new paradigm for studying theoretical and computational biology.

Effects of fish-human transmission and different life stages of fish on Clonorchiasis: A novel mathematical model.

Clonorchiasis is a zoonotic disease mainly caused by eating raw fish and shrimp, and there is no vaccine to prevent it. More than 30 million people are infected worldwide, of which China alone accounts for about half, and is one of the countries most seriously affected by Clonorchiasis. In this work, we formulate a novel Ordinary Differential Equation (ODE) model to discuss the biological attributes of fish within authentic ecosystems and the complex lifecycle of Clonorchis sinensis. This model includes larval fish, adult fish, infected fish, humans, and cercariae. We derive the basic reproduction number and perform a rigorous stability analysis of the proposed model. Numerically, we use data from 2016 to 2021 in Guangxi, China, to discuss outbreaks of Clonorchiasis and obtain the basic reproduction number R0=1.4764. The fitted curve appropriately reflects the overall trend and replicates a low peak in the case number of Clonorchiasis. By reducing the release rate of cercariae in 2018, the fitted values of Clonorchiasis cases dropped rapidly and almost disappeared. If we decrease the transmission rate from infected fish to humans, Clonorchiasis can be controlled. Our studies also suggest that strengthening publicity education and cleaning water quality can effectively control the transmission of Clonorchiasis in Guangxi, China.

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.

Evaluation of age-structured vaccination strategies for curbing the disease spread.

Age structure is one of the crucial factors in characterizing the heterogeneous epidemic transmission. Vaccination is regarded as an effective control measure for prevention and control epidemics. Due to the shortage of vaccine capacity during the outbreak of epidemics, how to design vaccination policy has become an urgent issue in suppressing the disease transmission. In this paper, we make an effort to propose an age-structured SVEIHR model with the disease-caused death to take account of dynamics of age-related vaccination policy for better understanding disease spread and control. We present an explicit expression of the basic reproduction number R 0 , which determines whether or not the disease persists, and then establish the existence and stability of endemic equilibria under certain conditions. Numerical simulations are illustrated to show that the age-related vaccination policy has a tremendous influence on curbing the disease transmission. Especially, vaccination of people over 65 is better than for people aged 21-65 in terms of rapid eradication of the disease in Italy.

Transmission dynamics of a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods.

In this paper, we propose a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods. Firstly, we establish the well-posedness of the model. Secondly, we define the basic reproduction number ℜ 0 for this model and show that ℜ 0 is a threshold parameter: if ℜ 0 < 1 , then the disease-free periodic solution is globally attractive; if ℜ 0 > 1 , the system is uniformly persistent. Thirdly, we study the global attractivity of the positive steady state when the spatial environment is homogeneous and the advection of mosquitoes is ignored. As an example, we use the model to investigate the dengue fever transmission case in Guangdong Province, China, and explore the impact of model parameters on ℜ 0 . Our findings indicate that ignoring seasonality may underestimate ℜ 0 . Additionally, the spatial heterogeneity of transmission may increase the risk of disease transmission, while the increase of seasonal developmental durations, intrinsic incubation periods and advection rates can all reduce the risk of disease transmission.

The role of long-lived plasma cells in viral clearance.

The adaptive immune system has two types of plasma cells (PC), long-lived plasma cells (LLPC) and short-lived plasma cells (SLPC), that differ in their lifespan. In this paper, we propose that LLPC is crucial to the clearance of viral particles in addition to reducing the viral basic reproduction number in secondary infections. We use a sequence of within-host mathematical models to show that, CD8 T cells, SLPC and memory B cells cannot achieve full viral clearance, and the viral load will reach a low positive equilibrium level because of a continuous replenishment of target cells. However, the presence of LLPC is crucial for viral clearance.

Dynamical Analysis of an Improved Bidirectional Immunization SIR Model in Complex Network.

In order to investigate the impact of two immunization strategies-vaccination targeting susceptible individuals to reduce their infection rate and clinical medical interventions targeting infected individuals to enhance their recovery rate-on the spread of infectious diseases in complex networks, this study proposes a bilinear SIR infectious disease model that considers bidirectional immunization. By analyzing the conditions for the existence of endemic equilibrium points, we derive the basic reproduction numbers and outbreak thresholds for both homogeneous and heterogeneous networks. The epidemic model is then reconstructed and extensively analyzed using continuous-time Markov chain (CTMC) methods. This analysis includes the investigation of transition probabilities, transition rate matrices, steady-state distributions, and the transition probability matrix based on the embedded chain. In numerical simulations, a notable concordance exists between the outcomes of CTMC and mean-field (MF) simulations, thereby substantiating the efficacy of the CTMC model. Moreover, the CTMC-based model adeptly captures the inherent stochastic fluctuation in the disease transmission, which is consistent with the mathematical properties of Markov chains. We further analyze the relationship between the system's steady-state infection density and the immunization rate through MCS. The results suggest that the infection density decreases with an increase in the immunization rate among susceptible individuals. The current research results will enhance our understanding of infectious disease transmission patterns in real-world scenarios, providing valuable theoretical insights for the development of epidemic prevention and control strategies.

Stability analysis of a SAIR epidemic model on scale-free community networks.

The presence of asymptomatic carriers, often unrecognized as infectious disease vectors, complicates epidemic management, particularly when inter-community migrations are involved. We introduced a SAIR (susceptible-asymptomatic-infected-recovered) infectious disease model within a network framework to explore the dynamics of disease transmission amid asymptomatic carriers. This model facilitated an in-depth analysis of outbreak control strategies in scenarios with active community migrations. Key contributions included determining the basic reproduction number, $ R_0 $, and analyzing two equilibrium states. Local asymptotic stability of the disease-free equilibrium is confirmed through characteristic equation analysis, while its global asymptotic stability is investigated using the decomposition theorem. Additionally, the global stability of the endemic equilibrium is established using the Lyapunov functional theory.

Final epidemic size of a two-community SIR model with asymmetric coupling.

Communities are commonly not isolated but interact asymmetrically with each other, allowing the propagation of infectious diseases within the same community and between different communities. To reveal the impact of asymmetrical interactions and contact heterogeneity on disease transmission, we formulate a two-community SIR epidemic model, in which each community has its contact structure while communication between communities occurs through temporary commuters. We derive an explicit formula for the basic reproduction number R 0 , give an implicit equation for the final epidemic size z, and analyze the relationship between them. Unlike the typical positive correlation between R 0 and z in the classic SIR model, we find a negatively correlated relationship between counterparts of our model deviating from homogeneous populations. Moreover, we investigate the impact of asymmetric coupling mechanisms on R 0 . The results suggest that, in scenarios with restricted movement of susceptible individuals within a community, R 0 does not follow a simple monotonous relationship, indicating that an unbending decrease in the movement of susceptible individuals may increase R 0 . We further demonstrate that network contacts within communities have a greater effect on R 0 than casual contacts between communities. Finally, we develop an epidemic model without restriction on the movement of susceptible individuals, and the numerical simulations suggest that the increase in human flow between communities leads to a larger R 0 .

Transmission dynamics of Zika virus with multiple infection routes and a case study in Brazil.

The Zika virus (ZIKV) is a serious global public health crisis. A major control challenge is its multiple transmission modes. This paper aims to simulate the transmission patterns of ZIKV using a dynamic process-based epidemiological model written in ordinary differential equations, which incorporates the human-to-mosquito infection by bites and sewage, mosquito-to-human infection by bites, and human-to-human infection by sex. Mathematical analyses are carried out to calculate the basic reproduction number and backward bifurcation, and prove the existence and stability of the equilibria. The model is validated with infection data by applying it to the 2015-2016 ZIKV epidemic in Brazil. The results indicate that the reproduction number is estimated to be 2.13, in which the contributions by mosquito bite, sex and sewage account for 85.7%, 3.5% and 10.8%, respectively. This number and the morbidity rate are most sensitive to parameters related to mosquito ecology, rather than asymptomatic or human-to-human transmission. Multiple transmission routes and suitable temperature exacerbate ZIKV infection in Brazil, and the vast majority of human infection cases were prevented by the intervention implemented. These findings may provide new insights to improve the risk assessment of ZIKV infection.

Optimal control of a multi-scale HIV-opioid model.

In this study, we apply optimal control theory to an immuno-epidemiological model of HIV and opioid epidemics. For the multi-scale model, we used four controls: treating the opioid use, reducing HIV risk behaviour among opioid users, entry inhibiting antiviral therapy, and antiviral therapy which blocks the viral production. Two population-level controls are combined with two within-host-level controls. We prove the existence and uniqueness of an optimal control quadruple. Comparing the two population-level controls, we find that reducing the HIV risk of opioid users has a stronger impact on the population who is both HIV-infected and opioid-dependent than treating the opioid disorder. The within-host-level antiviral treatment has an effect not only on the co-affected population but also on the HIV-only infected population. Our findings suggest that the most effective strategy for managing the HIV and opioid epidemics is combining all controls at both within-host and between-host scales.