Global sensitivity analysis in epidemiological modeling

Machine Learning Models to Predict Early Breakthrough of Recalcitrant Organic Micropollutants in Granular Activated Carbon Adsorbers

Granular activated carbon (GAC) adsorption is frequently used to remove recalcitrant organic micropollutants (MPs) from water. The overarching aim of this research was to develop machine learning (ML) models to predict GAC performance from adsorbent, adsorbate, and background water matrix properties. For model calibration, MP breakthrough curves were compiled and analyzed to determine the bed volumes of water that can be treated until MP breakthrough reaches ten percent of the influent MP concentration (BV10). Over 400 data points were split into training, validation, and testing sets. Seventeen variables describing MP, background water matrix, and GAC properties were explored in ML models to predict log10-transformed BV10 values. Using the ML models on the testing set, predicted BV10 values exhibited mean absolute errors of ∼0.12 log units and were highly correlated with experimentally determined values (R2 ≥ 0.88). The top three drivers influencing BV10 predictions were the air-hexadecane partition coefficient and hydrogen bond acidity (Abraham parameters L and A) of the MPs and the dissolved organic carbon concentration of the GAC influent water. The model can be used to rapidly estimate the GAC bed life, select effective GAC products for a given treatment scenario, and explore the suitability of GAC treatment for remediating emerging MPs.

Quantifying trade-offs between therapeutic efficacy and resistance dissemination for enrofloxacin dose regimens in cattle

The use of antimicrobial drugs in food-producing animals contributes to the selection pressure on pathogenic and commensal bacteria to become resistant. This study aims to evaluate the existence of trade-offs between treatment effectiveness, cost, and the dynamics of resistance in gut commensal bacteria. We developed a within-host ordinary differential equation model to track the dynamics of antimicrobial drug concentrations and bacterial populations in the site of infection (lung) and the gut. The model was parameterized to represent enrofloxacin treatment for bovine respiratory disease (BRD) caused by Pastereulla multocida in cattle. Three approved enrofloxacin dosing regimens were compared for their effects on resistance on P. multocida and commensal E. coli: 12.5 mg/kg and 7.5 mg/kg as a single dose, and 5 mg/kg as three doses. Additionally, we explored non-FDA-approved regimes. Our results indicated that both 12.5 mg/kg and 7.5 mg/kg as a single dose scenario increased the most the treatment costs and prevalence of P. multocida resistance in the lungs, while 5 mg/kg as three doses increased resistance in commensal E. coli bacteria in the gut the most out of the approved scenarios. A proposed non-FDA-approved scenario (7.5 mg/kg, two doses 24 h apart) showed low economic costs, minimal P. multocida, and moderate effects on resistant E. coli. Overall, the scenarios that decrease P. multocida, including resistant P. multocida did not coincide with those that decrease resistant E. coli the most, suggesting a trade-off between both outcomes. The sensitivity analysis suggests that bacterial populations were the most sensitive to drug conversion factors into plasma ( β ), elimination of the drug from the colon ( ϑ ), fifty percent sensitive bacteria (P. multocida) killing effect ( L s50 ), fifty percent of bacteria (E. coli) above ECOFF killing effect ( C r50 ), and net drug transfer rate in the lung ( γ ) parameters.

Predicting Tuberculosis Incidence and Its Trend in Tigray, Ethiopia: A Reality-Counterfactual Modeling Approach

Background

The Tigray region of Ethiopia, which has been affected by civil war from 2020 to 2022, is facing an increase in tuberculosis in the damaged health system. Our study employed mathematical modeling to predict the incidence of tuberculosis and its trends during the war and in the post-conflict setting of Tigray, Northern Ethiopia.

Methods

We predicted the incidence of tuberculosis from 2020 to 2025 in Tigray using the SEIRD model in the context of the recent war and compared it with its counterfactual trend in the absence of war. The counterfactual trend was forecasted using an autoregressive integrated moving average (ARIMA) model for stationary time-series data. We performed rolling origin cross-validation for ARIMA and sensitivity analysis for the SEIRD model. The initial tuberculosis data and model parameters were obtained from the Institute for Health Metrics and Evaluation and the literature, respectively.

Results

Between 2000 and 2017, the incidence of tuberculosis in Tigray decreased at an annual rate of 3.0%. Shortly before the war, the incidence of tuberculosis in the region was 178 per 100,000 people. In a counterfactual scenario where there was no war, the incidence was projected to decrease to 144.3 in 2022 and 126.3 in 2025. However, owing to the war and siege, the SEIRD-projected incidence of tuberculosis would have increased to 965.5 (95% CI: 958.5-972.7) in 2022 and 372.4 (95% CI: 367.7-376.6) in 2025. Over 800 cases of tuberculosis per 100,000 people were attributed to the war in 2022. In the postwar period, the incidence is projected to decrease by 30% by 2023.

Conclusion

The Tigray War reversed a two-decade decline in tuberculosis cases, causing a five-fold increase compared to the no-war scenario. Urgent interventions are needed to support tuberculosis prevention, testing, and treatment, particularly in key and vulnerable populations.

A modeling analysis of whole body potassium regulation on a high-potassium diet: proximal tubule and tubuloglomerular feedback effects

Potassium (K+) is an essential electrolyte that plays a key role in many physiological processes, including mineralcorticoid action, systemic blood-pressure regulation, and hormone secretion and action. Indeed, maintaining K+ balance is critical for normal cell function, as too high or too low K+ levels can have serious and potentially deadly health consequences. K+ homeostasis is achieved by an intricate balance between the intracellular and extracellular fluid as well as balance between K+ intake and excretion. This is achieved via the coordinated actions of regulatory mechanisms such as the gastrointestinal feedforward effect, insulin and aldosterone upregulation of Na+-K+-ATPase uptake, and hormone and electrolyte impacts on renal K+ handling. We recently developed a mathematical model of whole body K+ regulation to unravel the individual impacts of these regulatory mechanisms. In this study, we extend our mathematical model to incorporate recent experimental findings that showed decreased fractional proximal tubule reabsorption under a high-K+ diet. We conducted model simulations and sensitivity analyses to investigate how these renal alterations impact whole body K+ regulation. Model predictions quantify the sensitivity of K+ regulation to various levels of proximal tubule K+ reabsorption adaptation and tubuloglomerular feedback. Our results suggest that the reduced proximal tubule K+ reabsorption under a high-K+ diet could achieve K+ balance in isolation, but the resulting tubuloglomerular feedback reduces filtration rate and thus K+ excretion.NEW & NOTEWORTHY Potassium homeostasis is maintained in the body by a complex system of regulatory mechanisms. This system, when healthy, maintains a small extracellular potassium concentration, despite large fluctuations of dietary potassium. The complexities of the system make this problem well suited for investigation with mathematical modeling. In this study, we extend our mathematical model to consider recent experimental results on renal potassium handling on a high potassium diet and investigate the impacts from a whole body perspective.

Quantifying trade-offs between therapeutic efficacy and resistance dissemination for enrofloxacin dose regimens in cattle

The use of antimicrobial drugs in food-producing animals increases the selection pressure on pathogenic and commensal bacteria to become resistant. This study aims to evaluate the existence of trade-offs between treatment effectiveness, cost, and the dissemination of resistance in gut commensal bacteria. We developed a within-host ordinary differential equation model to track the dynamics of antimicrobial drug concentrations and bacterial populations in the site of infection (lung) and the gut. The model was parameterized to represent enrofloxacin treatment for bovine respiratory disease (BRD) caused by Pastereulla multocida in cattle. Three approved enrofloxacin dosing regimens were compared for their effects on resistance on P. multocida and commensal E. coli: 12.5 mg/kg and 7.5 mg/kg as a single dose, and 5 mg/kg as three doses. Additionally, we explored non-approved regimes. Our results indicated that both 12.5 mg/kg and 7.5 mg/kg as a single dose scenario increased the most the treatment costs and prevalence of P. multocida resistance in the lungs, while 5 mg/kg as three doses increased resistance in commensal E. coli bacteria in the gut the most out of the approved scenarios. A proposed scenario (7.5 mg/kg, two doses 24 hours apart) showed low economic costs, minimal P. multocida, and moderate effects on resistant E. coli. Overall, the scenarios that decrease P. multocida, including resistant P. multocida did not coincide with the scenarios that decrease resistant E. coli the most, suggesting a trade-off between both outcomes. The sensitivity analysis indicates that bacterial populations were the most sensitive to drug conversion factors into plasma (β), elimination of the drug from the colon (υ), fifty percent sensitive bacteria (P. multocida) killing effect (Ls50), fifty percent of bacteria (E. coli) above ECOFF killing effect (Cr50), and net drug transfer rate in the lung (γ) parameters.

Analysis of the dynamics of a vector-borne infection with the effect of imperfect vaccination from a fractional perspective

The burden of vector-borne infections is significant, particularly in low- and middle-income countries where vector populations are high and healthcare infrastructure may be inadequate. Further, studies are required to investigate the key factors of vector-borne infections to provide effective control measure. This study focuses on formulating a mathematical framework to characterize the spread of chikungunya infection in the presence of vaccines and treatments. The research is primarily dedicated to descriptive study and comprehension of dynamic behaviour of chikungunya dynamics. We use Banach's and Schaefer's fixed point theorems to investigate the existence and uniqueness of the suggested chikungunya framework resolution. Additionally, we confirm the Ulam-Hyers stability of the chikungunya system. To assess the impact of various parameters on the dynamics of chikungunya, we examine solution pathways using the Laplace-Adomian method of disintegration. Specifically, to visualise the impacts of fractional order, vaccination, bite rate and treatment computer algorithms are employed on the infection level of chikungunya. Our research identified the framework's essential input settings for managing chikungunya infection. Notably, the intensity of chikungunya infection can be reduced by lowering mosquito bite rates in the affected area. On the other hand, vaccination, memory index or fractional order, and treatment could be used as efficient controlling variables.

There is a need for more precise models to assess the determinants of health crises like COVID-19

The COVID-19 pandemic has had a significant impact on global mortality. While the causal relationship between SARS-CoV-2 and the anomalous increase in deaths is established, more precise and complex models are needed to determine the exact weight of epidemiological factors involved. Indeed, COVID-19 behavior is influenced by a wide range of variables, including demographic characteristics, population habits and behavior, healthcare performance, and environmental and seasonal risk factors. The bidirectional causality between impacted and impacting aspects, as well as confounding variables, complicates efforts to draw clear, generalizable conclusions regarding the effectiveness and cost-benefit ratio of non-pharmaceutical health countermeasures. Thus, it is imperative that the scientific community and health authorities worldwide develop comprehensive models not only for the current pandemic but also for future health crises. These models should be implemented locally to account for micro-differences in epidemiological characteristics that may have relevant effects. It is important to note that the lack of a universal model does not imply that local decisions have been unjustified, and the request to decrease scientific uncertainty does not mean denying the evidence of the effectiveness of the countermeasures adopted. Therefore, this paper must not be exploited to denigrate either the scientific community or the health authorities.

How containment policy and medical service impact COVID-19 transmission: A cross-national comparison among China, the USA, and Sweden

As COVID-19 shows a heterogeneous spreading process globally, investigating factors associated with COVID-19 spreading among different countries will provide information for containment strategy and medical service decisions. A significant challenge for analyzing how these factors impact COVID-19 transmission is assessing key epidemiological parameters and how they change under different containment strategies across different nations. This paper builds a COVID-19 spread simulation model to estimate the core COVID-19 epidemiological parameters. Then, the correlation between these core COVID-19 epidemiological parameters and the times of publicly announced interventions is analyzed, including three typical countries, China (strictly containment), the USA (moderately control), and Sweden (loose control). Results show that the recovery rate leads to a distinct COVID-19 transmission process in the three countries, as all three countries finally have similar and close to zero spreading rates in the third period of COVID-19 transmission. Then, an epidemic fundamental diagram between COVID-19 "active infections" and "current patients" is discovered, which could plan a country's COVID-19 medical capacity and containment strategies when combined with the COVID-19 spreading simulation model. Based on that, the hypothetical policies are proved effectively, which will give support for future infectious diseases.

Introduction to the special issue on the role of operational research in future epidemics/ pandemics

In this special issue, 23 research papers are published focusing on COVID-19 and operational research solution techniques. First, we detail the process from advertising the call for papers to the point where the best papers are accepted. Then, we provide a summary of each paper focusing on applications, solution techniques and insights for practitioners and policy makers. To provide a holistic view for readers, we have clustered the papers into different groups: transmission, propagation and forecasting, non-pharmaceutical intervention, healthcare network configuration, healthcare resource allocation, hospital operations, vaccine and testing kits, and production and manufacturing. Then, we introduce other possible subjects that can be considered for future research.