Search

Mathematical modelling of (re)emerging infectious respiratory diseases among humans poses multiple challenges for modellers, which can arise as a result of limited data and surveillance, uncertainty in the natural history of the disease, as well as public health and individual responses to outbreaks. Here, we propose a COVID-19-inspired health state diagram (HSD) to serve as a foundational framework for conceptualising the modelling process for (re)emerging respiratory diseases, and public health responses, in the early stages of their emergence. The HSD aims to serve as a starting point for reflection on the structure and parameterisation of a transmission model to assess the impact of the (re)emerging disease and the capacity of public health interventions to control transmission. We also explore the adaptability of the HSD to different (re)emerging diseases using the characteristics of three respiratory diseases of historical public health importance. We outline key questions to contemplate when applying and adapting this HSD to (re)emerging infectious diseases and provide reflections on adapting the framework for public health-related interventions. [ABSTRACT FROM AUTHOR]

Background: Recently, several novel RSV immunisation products that protect infants and older adults against RSV disease have been licensed in Europe. We estimated the effectiveness and efficiency of introducing these RSV immunisation strategies in Germany. Methods: We used a Bayesian framework to fit a deterministic age-structured dynamic transmission model of RSV to sentinel surveillance and RSV-specific hospitalisation data in Germany from 2015 to 2019. The calibrated model was used to evaluate different RSV intervention strategies over 5 years: long-acting, single-dose monoclonal antibodies (mAbs) in high-risk infants aged 1–5 months; long-acting mAbs in all infants aged 1–5 months; seasonal vaccination of pregnant women and one-time seasonal vaccination of older adults (75 + /65 + /55 + years). We performed sensitivity analysis on vaccine uptake, seasonal vs. year-round maternal vaccination, and the effect of under-ascertainment for older adults. Results: The model was able to match the various RSV datasets. Replacing the current short-acting mAB for high-risk infants with long-acting mAbs prevented 1.1% of RSV-specific hospitalisations in infants per year at the same uptake. Expanding the long-acting mAB programme to all infants prevented 39.3% of infant hospitalisations per year. Maternal vaccination required a larger number to be immunised to prevent one additional hospitalisation than a long-acting mAB for the same uptake. Vaccination of adults older than 75 years at an uptake of 40% in addition to Nirsevimab in all infants prevented an additional 4.5% of all RSV hospitalisations over 5 years, with substantial uncertainty in the correction for under-ascertainment of the RSV burden. Conclusions: Immunisation has the potential to reduce the RSV disease burden in Germany. [ABSTRACT FROM AUTHOR]

Climatic conditions are a key determinant of malaria transmission intensity, through their impacts on both the parasite and its mosquito vectors. Mathematical models relating climatic conditions to malaria transmission can be used to develop spatial maps of climatic suitability for malaria. These maps underpin efforts to quantify the distribution and burden of malaria in humans, enabling improved monitoring and control. Previous work has developed mathematical models and global maps for the suitability of temperature for malaria transmission. In this paper, existing temperature-based models are extended to include two other important bioclimatic factors: humidity and rainfall. This model is combined with fine spatial resolution climatic data to produce a more biologically-realistic global map of climatic suitability for Plasmodium falciparum malaria. The climatic suitability index developed corresponds more closely than previous temperature suitability indices with the global distribution of P. falciparum malaria. There is weak agreement between the Malaria Atlas Project estimates of P. falciparum prevalence in Africa and the estimates of suitability solely based on temperature (Spearman Correlation coefficient of ρ = 0.24 ). The addition of humidity and then rainfall improves the comparison ( ρ = 0.62 when humidity added; ρ = 0.70 when both humidity and rainfall added). By incorporating the impacts of humidity and rainfall, this model identifies arid regions that are not climatically suitable for transmission of P. falciparum malaria. Incorporation of this improved index of climatic suitability into geospatial models can improve global estimates of malaria prevalence and transmission intensity. [ABSTRACT FROM AUTHOR]

Background: Including structural determinants (e.g. criminalisation, stigma, inequitable gender norms) in dynamic HIV transmission models is important to help quantify their population-level impacts and guide implementation of effective interventions that reduce the burden of HIV and inequalities thereof. However, evidence-based modelling of structural determinants is challenging partly due to a limited understanding of their causal pathways and few empirical estimates of their effects on HIV acquisition and transmission. Methods: We conducted a scoping review of dynamic HIV transmission modelling studies that evaluated the impacts of structural determinants, published up to August 28, 2023, using Ovid Embase and Medline online databases. We appraised studies on how models represented exposure to structural determinants and causal pathways. Building on this, we developed a new methodological framework and recommendations to support the incorporation of structural determinants in transmission dynamics models and their analyses. We discuss the data and analyses that could strengthen the evidence used to inform these models. Results: We identified 17 HIV modelling studies that represented structural determinants and/or interventions, including incarceration of people who inject drugs (number of studies [n] = 5), violence against women (n = 3), HIV stigma (n = 1), and housing instability (n = 1), among others (n = 7). Most studies (n = 10) modelled exposures dynamically. Almost half (8/17 studies) represented multiple exposure histories (e.g. current, recent, non-recent exposure). Structural determinants were often assumed to influence HIV indirectly by influencing mediators such as contact patterns, condom use, and antiretroviral therapy use. However, causal pathways' assumptions were sometimes simple, with few mediators explicitly represented in the model, and largely based on cross-sectional associations. Although most studies calibrated models using HIV epidemiological data, less than half (7/17) also fitted or cross-validated to data on the prevalence, frequency, or effects of exposure to structural determinants. Conclusions: Mathematical models can play a crucial role in elucidating the population-level impacts of structural determinants and interventions on HIV. We recommend the next generation of models reflect exposure to structural determinants dynamically and mechanistically, and reproduce the key causal pathways, based on longitudinal evidence of links between structural determinants, mediators, and HIV. This would improve the validity and usefulness of predictions of the impacts of structural determinants and interventions. [ABSTRACT FROM AUTHOR]

Zoonotic diseases represent a significant societal challenge in terms of their health and economic impacts. One Health approaches to managing zoonotic diseases are becoming more prevalent, but require novel thinking, tools and cross-disciplinary collaboration. Bovine tuberculosis (bTB) is one example of a costly One Health challenge with a complex epidemiology involving humans, domestic animals, wildlife and environmental factors, which require sophisticated collaborative approaches. We undertook a scoping review of multi-host bTB epidemiology to identify trends in species publication focus, methodologies, and One Health approaches. We aimed to identify knowledge gaps where novel research could provide insights to inform control policy, for bTB and other zoonoses. The review included 532 articles. We found different levels of research attention across episystems, with a significant proportion of the literature focusing on the badger-cattle-TB episystem, with far less attention given to tropical multi-host episystems. We found a limited number of studies focusing on management solutions and their efficacy, with very few studies looking at modelling exit strategies. Only a small number of studies looked at the effect of human disturbances on the spread of bTB involving wildlife hosts. Most of the studies we reviewed focused on the effect of badger vaccination and culling on bTB dynamics with few looking at how roads, human perturbations and habitat change may affect wildlife movement and disease spread. Finally, we observed a lack of studies considering the effect of weather variables on bTB spread, which is particularly relevant when studying zoonoses under climate change scenarios. Significant technological and methodological advances have been applied to bTB episystems, providing explicit insights into its spread and maintenance across populations. We identified a prominent bias towards certain species and locations. Generating more high-quality empirical data on wildlife host distribution and abundance, high-resolution individual behaviours and greater use of mathematical models and simulations are key areas for future research. Integrating data sources across disciplines, and a "virtuous cycle" of well-designed empirical data collection linked with mathematical and simulation modelling could provide additional gains for policy-makers and managers, enabling optimised bTB management with broader insights for other zoonoses. [ABSTRACT FROM AUTHOR]

Background: Defining healthcare facility catchment areas is a key step in predicting future healthcare demand in epidemic settings. Forecasts of hospitalisations can be informed by leading indicators measured at the community level. However, this relies on the definition of so-called catchment areas or the geographies whose populations make up the patients admitted to a given hospital, which are often not well-defined. Little work has been done to quantify the impact of hospital catchment area definitions on healthcare demand forecasting. Methods: We made forecasts of local-level hospital admissions using a scaled convolution of local cases (as defined by the hospital catchment area) and delay distribution. Hospital catchment area definitions were derived from either simple heuristics (in which people are admitted to their nearest hospital or any nearby hospital) or historical admissions data (all emergency or elective admissions in 2019, or COVID-19 admissions), plus a marginal baseline definition based on the distribution of all hospital admissions. We evaluated predictive performance using each hospital catchment area definition using the weighted interval score and considered how this changed by the length of the predictive horizon, the date on which the forecast was made, and by location. We also considered the change, if any, in the relative performance of each definition in retrospective vs. real-time settings, or at different spatial scales. Results: The choice of hospital catchment area definition affected the accuracy of hospital admission forecasts. The definition based on COVID-19 admissions data resulted in the most accurate forecasts at both a 7- and 14-day horizon and was one of the top two best-performing definitions across forecast dates and locations. The "nearby" heuristic also performed well, but less consistently than the COVID-19 data definition. The marginal distribution baseline, which did not include any spatial information, was the lowest-ranked definition. The relative performance of the definitions was larger when using case forecasts compared to future observed cases. All results were consistent across spatial scales of the catchment area definitions. Conclusions: Using catchment area definitions derived from context-specific data can improve local-level hospital admission forecasts. Where context-specific data is not available, using catchment areas defined by carefully chosen heuristics is a sufficiently good substitute. There is clear value in understanding what drives local admissions patterns, and further research is needed to understand the impact of different catchment area definitions on forecast performance where case trends are more heterogeneous. [ABSTRACT FROM AUTHOR]

Background: The World Health Organization (WHO) has outlined a set of targets to achieve eliminating hepatitis C by 2030. In May 2022, Lithuanian health authorities initiated a hepatitis C virus (HCV) screening program to start working towards elimination. In the program, bonus was given to general practitioners (GPs) to promote and conduct anti-HCV tests for two situations: (1) one time testing for individuals born in 1945–1994 and (2) annual HCV testing for persons who inject drugs or are living with human immunodeficiency virus (HIV) regardless of age. This study aimed to model the current viral hepatitis C epidemiological status in Lithuania and to outline the requirements for WHO elimination targets using the first-year HCV screening results. Methods: Individuals were invited to participate in the anti-HCV screening by GPs during routine visits. Patients who tested positive were then referred to a gastroenterologist or infectious disease doctor for further confirmatory testing. If a patient received a positive RNA test and a fibrosis staging result of ≥ F2, the doctor prescribed direct-acting antivirals. Information on the patients screened, diagnosed, and treated was obtained from the National Health Insurance Fund. The Markov disease progression model, developed by the CDA Foundation, was used to evaluate the screening program results and HCV elimination progress in Lithuania. Results: Between May 2022 and April 2023, 790,070 individuals underwent anti-HCV testing, with 11,943 individuals (1.5%) receiving positive results. Anti-HCV seroprevalence was found to be higher among males than females, 1.9% and 1.2%, respectively. Within the risk population tested, 2087 (31.1%) seropositive individuals were identified. When comparing the screening program results to WHO elimination targets through modelling, 2180 patients still need to be treated annually until 2030, along with expanding fibrosis restrictions. If an elimination approach was implemented, 1000 new infections would be prevented, while saving 150 lives and averting 90 decompensated cirrhosis cases and 110 hepatocellular carcinoma cases. Conclusions: During the first year of the Lithuanian screening program, GPs were able to screen 44% of the target population. However, the country will not meet elimination targets as it currently stands without increasing treatment levels and lifting fibrosis restrictions. [ABSTRACT FROM AUTHOR]

Background: Since the emergence of SARS-CoV-2 (COVID-19), there have been multiple waves of infection and multiple rounds of vaccination rollouts. Both prior infection and vaccination can prevent future infection and reduce severity of outcomes, combining to form hybrid immunity against COVID-19 at the individual and population level. Here, we explore how different combinations of hybrid immunity affect the size and severity of near-future Omicron waves. Methods: To investigate the role of hybrid immunity, we use an agent-based model of COVID-19 transmission with waning immunity to simulate outbreaks in populations with varied past attack rates and past vaccine coverages, basing the demographics and past histories on the World Health Organization Western Pacific Region. Results: We find that if the past infection immunity is high but vaccination levels are low, then the secondary outbreak with the same variant can occur within a few months after the first outbreak; meanwhile, high vaccination levels can suppress near-term outbreaks and delay the second wave. Additionally, hybrid immunity has limited impact on future COVID-19 waves with immune-escape variants. Conclusions: Enhanced understanding of the interplay between infection and vaccine exposure can aid anticipation of future epidemic activity due to current and emergent variants, including the likely impact of responsive vaccine interventions. [ABSTRACT FROM AUTHOR]

Background: Highlighted by the rise of COVID-19, climate change, and conflict, socially vulnerable populations are least resilient to disaster. In infectious disease management, mathematical models are a commonly used tool. Researchers should include social vulnerability in models to strengthen their utility in reflecting real-world dynamics. We conducted a scoping review to evaluate how researchers have incorporated social vulnerability into infectious disease mathematical models. Methods: The methodology followed the Joanna Briggs Institute and updated Arksey and O'Malley frameworks, verified by the PRISMA-ScR checklist. PubMed, Clarivate Web of Science, Scopus, EBSCO Africa Wide Information, and Cochrane Library were systematically searched for peer-reviewed published articles. Screening and extracting data were done by two independent researchers. Results: Of 4075 results, 89 articles were identified. Two-thirds of articles used a compartmental model (n = 58, 65.2%), with a quarter using agent-based models (n = 24, 27.0%). Overall, routine indicators, namely age and sex, were among the most frequently used measures (n = 42, 12.3%; n = 22, 6.4%, respectively). Only one measure related to culture and social behaviour (0.3%). For compartmental models, researchers commonly constructed distinct models for each level of a social vulnerability measure and included new parameters or influenced standard parameters in model equations (n = 30, 51.7%). For all agent-based models, characteristics were assigned to hosts (n = 24, 100.0%), with most models including age, contact behaviour, and/or sex (n = 18, 75.0%; n = 14, 53.3%; n = 10, 41.7%, respectively). Conclusions: Given the importance of equitable and effective infectious disease management, there is potential to further the field. Our findings demonstrate that social vulnerability is not considered holistically. There is a focus on incorporating routine demographic indicators but important cultural and social behaviours that impact health outcomes are excluded. It is crucial to develop models that foreground social vulnerability to not only design more equitable interventions, but also to develop more effective infectious disease control and elimination strategies. Furthermore, this study revealed the lack of transparency around data sources, inconsistent reporting, lack of collaboration with local experts, and limited studies focused on modelling cultural indicators. These challenges are priorities for future research. [ABSTRACT FROM AUTHOR]

Background: The increased prevalence of insulin resistance is one of the major health risks in society today. Insulin resistance involves both short-term dynamics, such as altered meal responses, and long-term dynamics, such as the development of type 2 diabetes. Insulin resistance also occurs on different physiological levels, ranging from disease phenotypes to organ-organ communication and intracellular signaling. To better understand the progression of insulin resistance, an analysis method is needed that can combine different timescales and physiological levels. One such method is digital twins, consisting of combined mechanistic mathematical models. We have previously developed a model for short-term glucose homeostasis and intracellular insulin signaling, and there exist long-term weight regulation models. Herein, we combine these models into a first interconnected digital twin for the progression of insulin resistance in humans. Methods: The model is based on ordinary differential equations representing biochemical and physiological processes, in which unknown parameters were fitted to data using a MATLAB toolbox. Results: The interconnected twin correctly predicts independent data from a weight increase study, both for weight-changes, fasting plasma insulin and glucose levels, and intracellular insulin signaling. Similarly, the model can predict independent weight-change data in a weight loss study with the weight loss drug topiramate. The model can also predict non-measured variables. Conclusions: The model presented herein constitutes the basis for a new digital twin technology, which in the future could be used to aid medical pedagogy and increase motivation and compliance and thus aid in the prevention and treatment of insulin resistance. [ABSTRACT FROM AUTHOR]

Background: In a heterogeneous population of cells, individual cells can behave differently and respond variably to the environment. This cellular diversity can be assessed by measuring DNA methylation patterns. The loci with variable methylation patterns are informative of cellular heterogeneity and may serve as biomarkers of diseases and developmental progression. Cell-to-cell methylation heterogeneity can be evaluated through single-cell methylomes or computational techniques for pooled cells. However, the feasibility and performance of these approaches to precisely estimate methylation heterogeneity require further assessment. Results: Here, we proposed model-based methods adopted from a mathematical framework originally from biodiversity, to estimate genome-wide DNA methylation heterogeneity. We evaluated the performance of our models and the existing methods with feature comparison, and tested on both synthetic datasets and real data. Overall, our methods have demonstrated advantages over others because of their better correlation with the actual heterogeneity. We also demonstrated that methylation heterogeneity offers an additional layer of biological information distinct from the conventional methylation level. In the case studies, we showed that distinct profiles of methylation heterogeneity in CG and non-CG methylation can predict the regulatory roles between genomic elements in Arabidopsis. This opens up a new direction for plant epigenomics. Finally, we demonstrated that our score might be able to identify loci in human cancer samples as putative biomarkers for early cancer detection. Conclusions: We adopted the mathematical framework from biodiversity into three model-based methods for analyzing genome-wide DNA methylation heterogeneity to monitor cellular heterogeneity. Our methods, namely MeH, have been implemented, evaluated with existing methods, and are open to the research community. [ABSTRACT FROM AUTHOR]

Background: Mathematical models are increasingly used to inform HIV policy and planning. Comparing estimates obtained using different mathematical models can test the robustness of estimates and highlight research gaps. As part of a larger project aiming to determine the optimal allocation of funding for HIV services, in this study we compare projections from five mathematical models of the HIV epidemic in South Africa: EMOD-HIV, Goals, HIV-Synthesis, Optima, and Thembisa. Methods: The five modelling groups produced estimates of the total population, HIV incidence, HIV prevalence, proportion of people living with HIV who are diagnosed, ART coverage, proportion of those on ART who are virally suppressed, AIDS-related deaths, total deaths, and the proportion of adult males who are circumcised. Estimates were made under a "status quo" scenario for the period 1990 to 2040. For each output variable we assessed the consistency of model estimates by calculating the coefficient of variation and examining the trend over time. Results: For most outputs there was significant inter-model variability between 1990 and 2005, when limited data was available for calibration, good consistency from 2005 to 2025, and increasing variability towards the end of the projection period. Estimates of HIV incidence, deaths in people living with HIV, and total deaths displayed the largest long-term variability, with standard deviations between 35 and 65% of the cross-model means. Despite this variability, all models predicted a gradual decline in HIV incidence in the long-term. Projections related to the UNAIDS 95–95-95 targets were more consistent, with the coefficients of variation below 0.1 for all groups except children. Conclusions: While models produced consistent estimates for several outputs, there are areas of variability that should be investigated. This is important if projections are to be used in subsequent cost-effectiveness studies. [ABSTRACT FROM AUTHOR]

Background: Influenza is a major year-round cause of respiratory illness in Kenya, particularly in children under 5. Current influenza vaccines result in short-term, strain-specific immunity and were found in a previous study not to be cost-effective in Kenya. However, next-generation vaccines are in development that may have a greater impact and cost-effectiveness profile. Methods: We expanded a model previously used to evaluate the cost-effectiveness of seasonal influenza vaccines in Kenya to include next-generation vaccines by allowing for enhanced vaccine characteristics and multi-annual immunity. We specifically examined vaccinating children under 5 years of age with improved vaccines, evaluating vaccines with combinations of increased vaccine effectiveness, cross-protection between strains (breadth) and duration of immunity. We evaluated cost-effectiveness using incremental cost-effectiveness ratios (ICERs) and incremental net monetary benefits (INMBs) for a range of values for the willingness-to-pay (WTP) per DALY averted. Finally, we estimated threshold per-dose vaccine prices at which vaccination becomes cost-effective. Results: Next-generation vaccines can be cost-effective, dependent on the vaccine characteristics and assumed WTP thresholds. Universal vaccines (assumed to provide long-term and broad immunity) are most cost-effective in Kenya across three of four WTP thresholds evaluated, with the lowest median value of ICER per DALY averted ($263, 95% Credible Interval (CrI): $ − 1698, $1061) and the highest median INMBs. At a WTP of $623, universal vaccines are cost-effective at or below a median price of $5.16 per dose (95% CrI: $0.94, $18.57). We also show that the assumed mechanism underlying infection-derived immunity strongly impacts vaccine outcomes. Conclusions: This evaluation provides evidence for country-level decision makers about future next-generation vaccine introduction, as well as global research funders about the potential market for these vaccines. Next-generation vaccines may offer a cost-effective intervention to reduce influenza burden in low-income countries with year-round seasonality like Kenya. [ABSTRACT FROM AUTHOR]

Background: India had an estimated 2.9 million tuberculosis cases and 506 thousand deaths in 2021. Novel vaccines effective in adolescents and adults could reduce this burden. M72/AS01E and BCG-revaccination have recently completed phase IIb trials and estimates of their population-level impact are needed. We estimated the potential health and economic impact of M72/AS01E and BCG-revaccination in India and investigated the impact of variation in vaccine characteristics and delivery strategies. Methods: We developed an age-stratified compartmental tuberculosis transmission model for India calibrated to country-specific epidemiology. We projected baseline epidemiology to 2050 assuming no-new-vaccine introduction, and M72/AS01E and BCG-revaccination scenarios over 2025–2050 exploring uncertainty in product characteristics (vaccine efficacy, mechanism of effect, infection status required for vaccine efficacy, duration of protection) and implementation (achieved vaccine coverage and ages targeted). We estimated reductions in tuberculosis cases and deaths by each scenario compared to the no-new-vaccine baseline, as well as costs and cost-effectiveness from health-system and societal perspectives. Results: M72/AS01E scenarios were predicted to avert 40% more tuberculosis cases and deaths by 2050 compared to BCG-revaccination scenarios. Cost-effectiveness ratios for M72/AS01E vaccines were around seven times higher than BCG-revaccination, but nearly all scenarios were cost-effective. The estimated average incremental cost was US$190 million for M72/AS01E and US$23 million for BCG-revaccination per year. Sources of uncertainty included whether M72/AS01E was efficacious in uninfected individuals at vaccination, and if BCG-revaccination could prevent disease. Conclusions: M72/AS01E and BCG-revaccination could be impactful and cost-effective in India. However, there is great uncertainty in impact, especially given the unknowns surrounding the mechanism of effect and infection status required for vaccine efficacy. Greater investment in vaccine development and delivery is needed to resolve these unknowns in vaccine product characteristics. [ABSTRACT FROM AUTHOR]

Background : Binge drinking, inequitable gender norms and sexual risk behaviour are closely interlinked. This study aims to model the potential effect of alcohol counselling interventions (in men and women) and gender-transformative interventions (in men) as strategies to reduce HIV transmission. Methods : We developed an agent-based model of HIV and other sexually transmitted infections, allowing for effects of binge drinking on sexual risk behaviour, and effects of inequitable gender norms (in men) on sexual risk behaviour and binge drinking. The model was applied to South Africa and was calibrated using data from randomized controlled trials of alcohol counselling interventions (n = 9) and gender-transformative interventions (n = 4) in sub-Saharan Africa. The model was also calibrated to South African data on alcohol consumption and acceptance of inequitable gender norms. Binge drinking was defined as five or more drinks on a single day, in the last month. Results: Binge drinking is estimated to be highly prevalent in South Africa (54% in men and 35% in women, in 2021), and over the 2000–2021 period 54% (95% CI: 34–74%) of new HIV infections occurred in binge drinkers. Binge drinking accounted for 6.8% of new HIV infections (0.0–32.1%) over the same period, which was mediated mainly by an effect of binge drinking in women on engaging in casual sex. Inequitable gender norms accounted for 17.5% of incident HIV infections (0.0–68.3%), which was mediated mainly by an effect of inequitable gender norms on male partner concurrency. A multi-session alcohol counselling intervention that reaches all binge drinkers would reduce HIV incidence by 1.2% (0.0–2.5%) over a 5-year period, while a community-based gender-transformative intervention would reduce incidence by 3.2% (0.8–7.2%) or by 7.3% (0.6–21.2%) if there was no waning of intervention impact. Conclusions: Although binge drinking and inequitable gender norms contribute substantially to HIV transmission in South Africa, recently-trialled alcohol counselling and gender-transformative interventions are likely to have only modest effects on HIV incidence. Further innovation in developing locally-relevant interventions to address binge drinking and inequitable gender norms is needed. [ABSTRACT FROM AUTHOR]

Background: During the early stages of the COVID-19 pandemic, there was considerable uncertainty surrounding epidemiological and clinical aspects of SARS-CoV-2. Governments around the world, starting from varying levels of pandemic preparedness, needed to make decisions about how to respond to SARS-CoV-2 with only limited information about transmission rates, disease severity and the likely effectiveness of public health interventions. In the face of such uncertainties, formal approaches to quantifying the value of information can help decision makers to prioritise research efforts. Methods: In this study we use Value of Information (VoI) analysis to quantify the likely benefit associated with reducing three key uncertainties present in the early stages of the COVID-19 pandemic: the basic reproduction number ( R 0 ), case severity (CS), and the relative infectiousness of children compared to adults (CI). The specific decision problem we consider is the optimal level of investment in intensive care unit (ICU) beds. Our analysis incorporates mathematical models of disease transmission and clinical pathways in order to estimate ICU demand and disease outcomes across a range of scenarios. Results: We found that VoI analysis enabled us to estimate the relative benefit of resolving different uncertainties about epidemiological and clinical aspects of SARS-CoV-2. Given the initial beliefs of an expert, obtaining more information about case severity had the highest parameter value of information, followed by the basic reproduction number R 0 . Resolving uncertainty about the relative infectiousness of children did not affect the decision about the number of ICU beds to be purchased for any COVID-19 outbreak scenarios defined by these three parameters. Conclusion: For the scenarios where the value of information was high enough to justify monitoring, if CS and R 0 are known, management actions will not change when we learn about child infectiousness. VoI is an important tool for understanding the importance of each disease factor during outbreak preparedness and can help to prioritise the allocation of resources for relevant information. [ABSTRACT FROM AUTHOR]

Background: To reduce the burden from the COVID-19 pandemic in the United States, federal and state local governments implemented restrictions such as limitations on gatherings, restaurant dining, and travel, and recommended non-pharmaceutical interventions including physical distancing, mask-wearing, surface disinfection, and increased hand hygiene. Resulting behavioral changes impacted other infectious diseases including enteropathogens such as norovirus and rotavirus, which had fairly regular seasonal patterns prior to the COVID-19 pandemic. The study objective was to project future incidence of norovirus and rotavirus gastroenteritis as contacts resumed and other NPIs are relaxed. Methods: We fitted compartmental mathematical models to pre-pandemic U.S. surveillance data (2012–2019) for norovirus and rotavirus using maximum likelihood estimation. Then, we projected incidence for 2022–2030 under scenarios where the number of contacts a person has per day varies from70%, 80%, 90%, and full resumption (100%) of pre-pandemic levels. Results: We found that the population susceptibility to both viruses increased between March 2020 and November 2021. The 70–90% contact resumption scenarios led to lower incidence than observed pre-pandemic for both viruses. However, we found a greater than two-fold increase in community incidence relative to the pre-pandemic period under the 100% contact scenarios for both viruses. With rotavirus, for which population immunity is driven partially by vaccination, patterns settled into a new steady state quickly in 2022 under the 70–90% scenarios. For norovirus, for which immunity is relatively short-lasting and only acquired through infection, surged under the 100% contact scenario projection. Conclusions: These results, which quantify the consequences of population susceptibility build-up, can help public health agencies prepare for potential resurgence of enteric viruses. [ABSTRACT FROM AUTHOR]

Background: In countries with mature generalized HIV epidemics such as Uganda, there are still groups of individuals that are disproportionately affected. Among the key populations in Uganda are fishing communities, which make up about 10% of the population. Compared to the general population, HIV prevalence and incidence among individuals living in these communities is high. This high HIV burden has been attributed to several factors including limited access to prevention and treatment services as well as ongoing high-risk sexual behaviour. Methods: We investigated the impact of combined HIV prevention interventions on HIV transmission dynamics in high-risk fishing communities in Uganda using a deterministic compartmental model. The model was calibrated to seroprevalence data from a census performed in 2014. To account for remaining uncertainty in the calibrated model parameters, 50 000 simulated scenarios were modelled to investigate the impact of combined prevention interventions. Results: The projected HIV incidence decreased from 1.87 per 100 PY without intervention scale-up to 0.25 per 100 PY after 15 years (2014–2029) of intervention scale-up. A potential combination achieving this 87% reduction in incidence over 15 years in Ugandan FCs included condom use in about 60% of sexual acts, 23% of susceptible men circumcised, 87% of people living with HIV aware of their status, 75% of those on ART, and about 3% of susceptible individuals on oral PrEP. Uncertainty analysis revealed relative reductions in incidence ranging from 30.9 to 86.8%. Sensitivity analyses suggested that condom use and early ART were the most important interventions. Conclusion: Reducing HIV incidence, as well as prevalence and AIDS-related mortality, in these high-risk fishing communities in Uganda is attainable over 15 years with a combination prevention package. Our projected intervention coverage levels are well within the national targets set by the Uganda government and enable coming close to reaching the UNAIDS 95-95-95 targets to end AIDS as a public health threat by 2030. [ABSTRACT FROM AUTHOR]

Background: Understanding the overall effectiveness of non-pharmaceutical interventions to control the COVID-19 pandemic and reduce the burden of disease is crucial for future pandemic planning. However, quantifying the effectiveness of specific control measures and the extent of missed infections, in the absence of early large-scale serological surveys or random community testing, has remained challenging. Methods: Combining data on notified local COVID-19 cases with known and unknown sources of infections in Singapore with a branching process model, we reconstructed the incidence of missed infections during the early phase of the wild-type SARS-CoV-2 and Delta variant transmission. We then estimated the relative effectiveness of border control measures, case finding and contact tracing when there was no or low vaccine coverage in the population. We compared the risk of ICU admission and death between the wild-type SARS-CoV-2 and the Delta variant in notified cases and all infections. Results: We estimated strict border control measures were associated with 0.2 (95% credible intervals, CrI 0.04–0.8) missed imported infections per notified case between July and December 2020, a decline from around 1 missed imported infection per notified case in the early phases of the pandemic. Contact tracing was estimated to identify 78% (95% CrI 62–93%) of the secondary infections generated by notified cases before the partial lockdown in Apr 2020, but this declined to 63% (95% CrI 56–71%) during the lockdown and rebounded to 78% (95% CrI 58–94%) during reopening in Jul 2020. The contribution of contact tracing towards overall outbreak control also hinges on ability to find cases with unknown sources of infection: 42% (95% CrI 12–84%) of such cases were found prior to the lockdown; 10% (95% CrI 7–15%) during the lockdown; 47% (95% CrI 17–85%) during reopening, due to increased testing capacity and health-seeking behaviour. We estimated around 63% (95% CrI 49–78%) of the wild-type SARS-CoV-2 infections were undetected during 2020 and around 70% (95% CrI 49–91%) for the Delta variant in 2021. Conclusions: Combining models with case linkage data enables evaluation of the effectiveness of different components of outbreak control measures, and provides more reliable situational awareness when some cases are missed. Using such approaches for early identification of the weakest link in containment efforts could help policy makers to better redirect limited resources to strengthen outbreak control. [ABSTRACT FROM AUTHOR]

Background: Epithelial-mesenchymal plasticity (EMP) involves bidirectional transitions between epithelial, mesenchymal and multiple intermediary hybrid epithelial/mesenchymal phenotypes. While the process of epithelial-mesenchymal transition (EMT) and its associated transcription factors are well-characterised, the transcription factors that promote mesenchymal-epithelial transition (MET) and stabilise hybrid E/M phenotypes are less well understood. Results: Here, we analyse multiple publicly-available transcriptomic datasets at bulk and single-cell level and pinpoint ELF3 as a factor that is strongly associated with an epithelial phenotype and is inhibited during EMT. Using mechanism-based mathematical modelling, we also show that ELF3 inhibits the progression of EMT. This behaviour was also observed in the presence of an EMT inducing factor WT1. Our model predicts that the MET induction capacity of ELF3 is stronger than that of KLF4, but weaker than that of GRHL2. Finally, we show that ELF3 levels correlates with worse patient survival in a subset of solid tumour types. Conclusion: ELF3 is shown to be inhibited during EMT progression and is also found to inhibit the progression of complete EMT suggesting that ELF3 may be able to counteract EMT induction, including in the presence of EMT-inducing factors, such as WT1. The analysis of patient survival data indicates that the prognostic capacity of ELF3 is specific to cell-of-origin or lineage. [ABSTRACT FROM AUTHOR]

Background: The triglyceride-glucose index (TyG) has been proposed as a surrogate marker of insulin resistance, which is a typical trait of pregnancy. However, very few studies analyzed TyG performance as marker of insulin resistance in pregnancy, and they were limited to insulin resistance assessment at fasting rather than in dynamic conditions, i.e., during an oral glucose tolerance test (OGTT), which allows more reliable assessment of the actual insulin sensitivity impairment. Thus, first aim of the study was exploring in pregnancy the relationships between TyG and OGTT-derived insulin sensitivity. In addition, we developed a new version of TyG, for improved performance as marker of insulin resistance in pregnancy. Methods: At early pregnancy, a cohort of 109 women underwent assessment of maternal biometry and blood tests at fasting, for measurements of several variables (visit 1). Subsequently (26 weeks of gestation) all visit 1 analyses were repeated (visit 2), and a subgroup of women (84 selected) received a 2 h-75 g OGTT (30, 60, 90, and 120 min sampling) with measurement of blood glucose, insulin and C-peptide for reliable assessment of insulin sensitivity (PREDIM index) and insulin secretion/beta-cell function. The dataset was randomly split into 70% training set and 30% test set, and by machine learning approach we identified the optimal model, with TyG included, showing the best relationship with PREDIM. For inclusion in the model, we considered only fasting variables, in agreement with TyG definition. Results: The relationship of TyG with PREDIM was weak. Conversely, the improved TyG, called TyGIS, (linear function of TyG, body weight, lean body mass percentage and fasting insulin) resulted much strongly related to PREDIM, in both training and test sets (R2 > 0.64, p < 0.0001). Bland–Altman analysis and equivalence test confirmed the good performance of TyGIS in terms of association with PREDIM. Different further analyses confirmed TyGIS superiority over TyG. Conclusions: We developed an improved version of TyG, as new surrogate marker of insulin sensitivity in pregnancy (TyGIS). Similarly to TyG, TyGIS relies only on fasting variables, but its performances are remarkably improved than those of TyG. [ABSTRACT FROM AUTHOR]

Background: SARS-CoV-2 is known to transmit in hospital settings, but the contribution of infections acquired in hospitals to the epidemic at a national scale is unknown.Methods: We used comprehensive national English datasets to determine the number of COVID-19 patients with identified hospital-acquired infections (with symptom onset > 7 days after admission and before discharge) in acute English hospitals up to August 2020. As patients may leave the hospital prior to detection of infection or have rapid symptom onset, we combined measures of the length of stay and the incubation period distribution to estimate how many hospital-acquired infections may have been missed. We used simulations to estimate the total number (identified and unidentified) of symptomatic hospital-acquired infections, as well as infections due to onward community transmission from missed hospital-acquired infections, to 31st July 2020.Results: In our dataset of hospitalised COVID-19 patients in acute English hospitals with a recorded symptom onset date (n = 65,028), 7% were classified as hospital-acquired. We estimated that only 30% (range across weeks and 200 simulations: 20-41%) of symptomatic hospital-acquired infections would be identified, with up to 15% (mean, 95% range over 200 simulations: 14.1-15.8%) of cases currently classified as community-acquired COVID-19 potentially linked to hospital transmission. We estimated that 26,600 (25,900 to 27,700) individuals acquired a symptomatic SARS-CoV-2 infection in an acute Trust in England before 31st July 2020, resulting in 15,900 (15,200-16,400) or 20.1% (19.2-20.7%) of all identified hospitalised COVID-19 cases.Conclusions: Transmission of SARS-CoV-2 to hospitalised patients likely caused approximately a fifth of identified cases of hospitalised COVID-19 in the "first wave" in England, but less than 1% of all infections in England. Using time to symptom onset from admission for inpatients as a detection method likely misses a substantial proportion (> 60%) of hospital-acquired infections. [ABSTRACT FROM AUTHOR]

Background: Substantial research is underway to develop next-generation interventions that address current malaria control challenges. As there is limited testing in their early development, it is difficult to predefine intervention properties such as efficacy that achieve target health goals, and therefore challenging to prioritize selection of novel candidate interventions. Here, we present a quantitative approach to guide intervention development using mathematical models of malaria dynamics coupled with machine learning. Our analysis identifies requirements of efficacy, coverage, and duration of effect for five novel malaria interventions to achieve targeted reductions in malaria prevalence. Methods: A mathematical model of malaria transmission dynamics is used to simulate deployment and predict potential impact of new malaria interventions by considering operational, health-system, population, and disease characteristics. Our method relies on consultation with product development stakeholders to define the putative space of novel intervention specifications. We couple the disease model with machine learning to search this multi-dimensional space and efficiently identify optimal intervention properties that achieve specified health goals. Results: We apply our approach to five malaria interventions under development. Aiming for malaria prevalence reduction, we identify and quantify key determinants of intervention impact along with their minimal properties required to achieve the desired health goals. While coverage is generally identified as the largest driver of impact, higher efficacy, longer protection duration or multiple deployments per year are needed to increase prevalence reduction. We show that interventions on multiple parasite or vector targets, as well as combinations the new interventions with drug treatment, lead to significant burden reductions and lower efficacy or duration requirements. Conclusions: Our approach uses disease dynamic models and machine learning to support decision-making and resource investment, facilitating development of new malaria interventions. By evaluating the intervention capabilities in relation to the targeted health goal, our analysis allows prioritization of interventions and of their specifications from an early stage in development, and subsequent investments to be channeled cost-effectively towards impact maximization. This study highlights the role of mathematical models to support intervention development. Although we focus on five malaria interventions, the analysis is generalizable to other new malaria interventions. [ABSTRACT FROM AUTHOR]

Background: Recent years have seen important improvements in available preventive treatment regimens for tuberculosis (TB), and research is ongoing to develop these further. To assist with the formulation of target product profiles for future regimens, we examined which regimen properties would be most influential in the epidemiological impact of preventive treatment. Methods: Following expert consultation, we identified 5 regimen properties relevant to the incidence-reducing impact of a future preventive treatment regimen: regimen duration, efficacy, ease-of-adherence (treatment completion rates in programmatic conditions), forgiveness to non-completion and the barrier to developing rifampicin resistance during treatment. For each regimen property, we elicited expert input for minimally acceptable and optimal (ideal-but-feasible) performance scenarios for future regimens. Using mathematical modelling, we then examined how each regimen property would influence the TB incidence reduction arising from full uptake of future regimens according to current WHO guidelines, in four countries: South Africa, Kenya, India and Brazil. Results: Of all regimen properties, efficacy is the single most important predictor of epidemiological impact, while ease-of-adherence plays an important secondary role. These results are qualitatively consistent across country settings; sensitivity analyses show that these results are also qualitatively robust to a range of model assumptions, including the mechanism of action of future preventive regimens. Conclusions: As preventive treatment regimens against TB continue to improve, understanding the key drivers of epidemiological impact can assist in guiding further development. By meeting these key targets, future preventive treatment regimens could play a critical role in global efforts to end TB. [ABSTRACT FROM AUTHOR]

Background: Sub-Saharan Africa has seen substantial reductions in cases and deaths due to malaria over the past two decades. While this reduction is primarily due to an increasing expansion of interventions, urbanisation has played its part as urban areas typically experience substantially less malaria transmission than rural areas. However, this may be partially lost with the invasion and establishment of Anopheles stephensi. A. stephensi, the primary urban malaria vector in Asia, was first detected in Africa in 2012 in Djibouti and was subsequently identified in Ethiopia in 2016, and later in Sudan and Somalia. In Djibouti, malaria cases have increased 30-fold from 2012 to 2019 though the impact in the wider region remains unclear.Methods: Here, we have adapted an existing model of mechanistic malaria transmission to estimate the increase in vector density required to explain the trends in malaria cases seen in Djibouti. To account for the observed plasticity in An. stephensi behaviour, and the unknowns of how it will establish in a novel environment, we sample behavioural parameters in order to account for a wide range of uncertainty. This quantification is then applied to Ethiopia, considering temperature-dependent extrinsic incubation periods, pre-existing vector-control interventions and Plasmodium falciparum prevalence in order to assess the potential impact of An. stephensi establishment on P. falciparum transmission. Following this, we estimate the potential impact of scaling up ITN (insecticide-treated nets)/IRS (indoor residual spraying) and implementing piperonyl butoxide (PBO) ITNs and larval source management, as well as their economic costs.Results: We estimate that annual P. falciparum malaria cases could increase by 50% (95% CI 14-90) if no additional interventions are implemented. The implementation of sufficient control measures to reduce malaria transmission to pre-stephensi levels will cost hundreds of millions of USD.Conclusions: Substantial heterogeneity across the country is predicted and large increases in vector control interventions could be needed to prevent a major public health emergency. [ABSTRACT FROM AUTHOR]

Background: To accelerate progress against malaria in high burden countries, a strategic reorientation of resources at the sub-national level is needed. This paper describes how mathematical modelling was used in mainland Tanzania to support the strategic revision that followed the mid-term review of the 2015–2020 national malaria strategic plan (NMSP) and the epidemiological risk stratification at the council level in 2018. Methods: Intervention mixes, selected by the National Malaria Control Programme, were simulated for each malaria risk strata per council. Intervention mixes included combinations of insecticide-treated bed nets (ITN), indoor residual spraying, larval source management, and intermittent preventive therapies for school children (IPTsc). Effective case management was either based on estimates from the malaria indicator survey in 2016 or set to a hypothetical target of 85%. A previously calibrated mathematical model in OpenMalaria was used to compare intervention impact predictions for prevalence and incidence between 2016 and 2020, or 2022. Results: For each malaria risk stratum four to ten intervention mixes were explored. In the low-risk and urban strata, the scenario without a ITN mass campaign in 2019, predicted high increase in prevalence by 2020 and 2022, while in the very-low strata the target prevalence of less than 1% was maintained at low pre-intervention transmission intensity and high case management. In the moderate and high strata, IPTsc in addition to existing vector control was predicted to reduce the incidence by an additional 15% and prevalence by 22%. In the high-risk strata, all interventions together reached a maximum reduction of 76%, with around 70% of that reduction attributable to high case management and ITNs. Overall, the simulated revised NMSP was predicted to achieve a slightly lower prevalence in 2020 compared to the 2015–2020 NMSP (5.3% vs 6.3%). Conclusion: Modelling supported the choice of intervention per malaria risk strata by providing impact comparisons of various alternative intervention mixes to address specific questions relevant to the country. The use of a council-calibrated model, that reproduces local malaria trends, represents a useful tool for compiling available evidence into a single analytical platform, that complement other evidence, to aid national programmes with decision-making processes. [ABSTRACT FROM AUTHOR]

Background: With the evolving growth of the COVID-19 epidemic, travel restriction policies would need to be adjusted accordingly. Prohibition of mass event may be relaxed for social and economic benefits when virus transmission stops but could bear the risk of epidemic rebound. Against the background of the varied SARS-CoV-2 prevalence internationally, we modelled the potential impacts of pre-event interventions on epidemic risk of holding a mass event when COVID-19 is under control.Methods: We developed a mathematical model of SARS-CoV-2 transmission in Guangdong Province, China, where local virus transmission ceased to occur. A large-scale international trade fair was assumed to be held, with influx of people from overseas and rest of China over a short period of time, who participated for 2-week. Scenarios of pre-event intervention (none, quarantine arrangement and polymerase chain reaction (PCR) testing for participants) were compared. The influence of contact pattern, SARS-CoV-2 prevalence outside the province and China, and testing coverage were examined in sensitivity analyses.Results: In basecase scenario (no event), the epidemic has been under control since March 2020. The event would lead to the detection of 1% more confirmed cases by 31 July when community contact rate increases to pre-epidemic level. In event scenario without additional interventions, there would be 599 (93%) more new infections comparing with basecase scenario. To avert new infections, quarantining all participants before the event would be the most effective strategy, followed by quarantining all overseas participants and testing all other participants, and testing all participants before the event and on day 7. However, testing strategy is likely to be affected by the SARS-CoV-2 prevalence outside the event province.Conclusions: Pre-event interventions are effective for reducing the risk of epidemic rebound caused by an international large-scale event. Universal testing for participants is likely to be an effective and feasible intervention. [ABSTRACT FROM AUTHOR]

Background: The ability to rapidly adapt to adverse environmental conditions represents the key of success of many pathogens and, in particular, of Mycobacterium tuberculosis. Upon exposition to heat shock, antibiotics or other sources of stress, appropriate responses in terms of genes transcription and proteins activity are activated leading part of a genetically identical bacterial population to express a different phenotype, namely to develop persistence. When the stress response network is mathematically described by an ordinary differential equations model, development of persistence in the bacterial population is associated with bistability of the model, since different emerging phenotypes are represented by different stable steady states. Results: In this work, we develop a mathematical model of SigE stress response network that incorporates interactions not considered in mathematical models currently available in the literature. We provide, through involved analytical computations, accurate approximations of the system's nullclines, and exploit the obtained expressions to determine, in a reliable though computationally efficient way, the number of equilibrium points of the system. Conclusions: Theoretical analysis and perturbation experiments point out the crucial role played by the degradation pathway involving RseA, the anti-sigma factor of SigE, for coexistence of two stable equilibria and the emergence of bistability. Our results also indicate that a fine control on RseA concentration is a necessary requirement in order for the system to exhibit bistability. [ABSTRACT FROM AUTHOR]

Background: To reduce the coronavirus disease burden in England, along with many other countries, the government implemented a package of non-pharmaceutical interventions (NPIs) that have also impacted other transmissible infectious diseases such as norovirus. It is unclear what future norovirus disease incidence is likely to look like upon lifting these restrictions.Methods: Here we use a mathematical model of norovirus fitted to community incidence data in England to project forward expected incidence based on contact surveys that have been collected throughout 2020-2021.Results: We report that susceptibility to norovirus infection has likely increased between March 2020 and mid-2021. Depending upon assumptions of future contact patterns incidence of norovirus that is similar to pre-pandemic levels or an increase beyond what has been previously reported is likely to occur once restrictions are lifted. Should adult contact patterns return to 80% of pre-pandemic levels, the incidence of norovirus will be similar to previous years. If contact patterns return to pre-pandemic levels, there is a potential for the expected annual incidence to be up to 2-fold larger than in a typical year. The age-specific incidence is similar across all ages.Conclusions: Continued national surveillance for endemic diseases such as norovirus will be essential after NPIs are lifted to allow healthcare services to adequately prepare for a potential increase in cases and hospital pressures beyond what is typically experienced. [ABSTRACT FROM AUTHOR]

Background: Cyprus addressed the first wave of SARS CoV-2 (COVID-19) by implementing non-pharmaceutical interventions (NPIs). The aims of this study were: a) to estimate epidemiological parameters of this wave including infection attack ratio, infection fatality ratio, and case ascertainment ratio, b) to assess the impact of public health interventions and examine what would have happened if those interventions had not been implemented.Methods: A dynamic, stochastic, individual-based Susceptible-Exposed-Infected-Recovered (SEIR) model was developed to simulate COVID-19 transmission and progression in the population of the Republic of Cyprus. The model was fitted to the observed trends in COVID-19 deaths and intensive care unit (ICU) bed use.Results: By May 8th, 2020, the infection attack ratio was 0.31% (95% Credible Interval [CrI]: 0.15, 0.54%), the infection fatality ratio was 0.71% (95% CrI: 0.44, 1.61%), and the case ascertainment ratio was 33.2% (95% CrI: 19.7, 68.7%). If Cyprus had not implemented any public health measure, the healthcare system would have been overwhelmed by April 14th. The interventions averted 715 (95% CrI: 339, 1235) deaths. If Cyprus had only increased ICU beds, without any social distancing measure, the healthcare system would have been overwhelmed by April 19th.Conclusions: The decision of the Cypriot authorities to launch early NPIs limited the burden of the first wave of COVID-19. The findings of these analyses could help address the next waves of COVID-19 in Cyprus and other similar settings. [ABSTRACT FROM AUTHOR]

Background: SARS-CoV-2 dynamics are driven by human behaviour. Social contact data are of utmost importance in the context of transmission models of close-contact infections.Methods: Using online representative panels of adults reporting on their own behaviour as well as parents reporting on the behaviour of one of their children, we collect contact mixing (CoMix) behaviour in various phases of the COVID-19 pandemic in over 20 European countries. We provide these timely, repeated observations using an online platform: SOCRATES-CoMix. In addition to providing cleaned datasets to researchers, the platform allows users to extract contact matrices that can be stratified by age, type of day, intensity of the contact and gender. These observations provide insights on the relative impact of recommended or imposed social distance measures on contacts and can inform mathematical models on epidemic spread.Conclusion: These data provide essential information for policymakers to balance non-pharmaceutical interventions, economic activity, mental health and wellbeing, during vaccine rollout. [ABSTRACT FROM AUTHOR]

Background: Bluetongue is a serious disease of ruminants caused by the bluetongue virus (BTV). BTV is transmitted by biting midges (Culicoides spp.). Serological evidence from livestock and the presence of at least one competent vector species of Culicoides suggests that transmission of BTV is possible and may have occurred in Kazakhstan. Methods: We estimated the risk of transmission using a mathematical model of the reproduction number R0 for bluetongue. This model depends on livestock density and climatic factors which affect vector density. Data on climate and livestock numbers from the 2466 local communities were used. This, together with previously published model parameters, was used to estimate R0 for each month of the year. We plotted the results on isopleth maps of Kazakhstan using interpolation to smooth the irregular data. We also mapped the estimated proportion of the population requiring vaccination to prevent outbreaks of bluetongue. Results: The results suggest that transmission of bluetongue in Kazakhstan is not possible in the winter from October to March. Assuming there are vector-competent species of Culicoides endemic in Kazakhstan, then low levels of risk first appear in the south of Kazakhstan in April before spreading north and intensifying, reaching maximum levels in northern Kazakhstan in July. The risk declined in September and had disappeared by October. Conclusion: These results should aid in surveillance efforts for the detection and control of bluetongue in Kazakhstan by indicating where and when outbreaks of bluetongue are most likely to occur. The results also indicate where vaccination efforts should be focussed to prevent outbreaks of disease. [ABSTRACT FROM AUTHOR]

The management of infectious diseases in wildlife reservoirs is challenging and faces several limitations. However, detailed knowledge of host–pathogen systems often reveal heterogeneity among the hosts' contribution to transmission. Management strategies targeting specific classes of individuals and/or areas, having a particular role in transmission, could be more effective and more acceptable than population-wide interventions. In the wild population of Alpine ibex (Capra ibex—a protected species) of the Bargy massif (French Alps), females transmit brucellosis (Brucella melitensis) infection in ~90% of cases, and most transmissions occur in the central spatial units ("core area"). Therefore, we expanded an individual-based model, developed in a previous study, to test whether strategies targeting females or the core area, or both, would be more effective. We simulated the relative efficacy of realistic strategies for the studied population, combining test-and-remove (euthanasia of captured animals with seropositive test results) and partial culling of unmarked animals. Targeting females or the core area was more effective than untargeted management options, and strategies targeting both were even more effective. Interestingly, the number of ibex euthanized and culled in targeted strategies were lower than in untargeted ones, thus decreasing the conservation costs while increasing the sanitary benefits. Although there was no silver bullet for the management of brucellosis in the studied population, targeted strategies offered a wide range of promising refinements to classical sanitary measures. We therefore encourage to look for heterogeneity in other wildlife diseases and to evaluate potential strategies for improving management in terms of efficacy but also acceptability. [ABSTRACT FROM AUTHOR]

Background: Larviciding against malaria vectors in Africa has been limited compared to indoor residual spraying and insecticide treated nets but is increasingly being considered by some countries as a complementary strategy. However, despite progress towards improved larvicides and new tools for mapping or treating mosquito-breeding sites, little is known about the optimal deployment strategies for larviciding in different transmission and seasonality settings. Methods: A malaria transmission model, OpenMalaria, was used to simulate varying larviciding strategies and their impact on host-seeking mosquito densities, entomological inoculation rate (EIR) and malaria prevalence. Variations in coverage, duration, frequency, and timing of larviciding were simulated for three transmission intensities and four transmission seasonality profiles. Malaria transmission was assumed to follow rainfall with a lag of one month. Theoretical sub-Saharan African settings with Anopheles gambiae as the dominant vector were chosen to explore impact. Relative reduction compared to no larviciding was predicted for each indicator during the simulated larviciding period. Results: Larviciding immediately reduced the predicted host-seeking mosquito densities and EIRs to a maximum that approached or exceeded the simulated coverage. Reduction in prevalence was delayed by approximately one month. The relative reduction in prevalence was up to four times higher at low than high transmission. Reducing larviciding frequency (i.e., from every 5 to 10 days) resulted in substantial loss in effectiveness (54, 45 and 53% loss of impact for host-seeking mosquito densities, EIR and prevalence, respectively). In seasonal settings the most effective timing of larviciding was during or at the beginning of the rainy season and least impactful during the dry season, assuming larviciding deployment for four months. Conclusion: The results highlight the critical role of deployment strategies on the impact of larviciding. Overall, larviciding would be more effective in settings with low and seasonal transmission, and at the beginning and during the peak densities of the target species populations. For maximum impact, implementers should consider the practical ranges of coverage, duration, frequency, and timing of larviciding in their respective contexts. More operational data and improved calibration would enable models to become a practical tool to support malaria control programmes in developing larviciding strategies that account for the diversity of contexts. [ABSTRACT FROM AUTHOR]

Background: Malaria elimination is the goal for Bioko Island, Equatorial Guinea. Intensive interventions implemented since 2004 have reduced prevalence, but progress has stalled in recent years. A challenge for elimination has been malaria infections in residents acquired during travel to mainland Equatorial Guinea. The present article quantifies how off-island contributes to remaining malaria prevalence on Bioko Island, and investigates the potential role of a pre-erythrocytic vaccine in making further progress towards elimination. Methods: Malaria transmission on Bioko Island was simulated using a model calibrated based on data from the Malaria Indicator Surveys (MIS) from 2015 to 2018, including detailed travel histories and malaria positivity by rapid-diagnostic tests (RDTs), as well as geospatial estimates of malaria prevalence. Mosquito population density was adjusted to fit local transmission, conditional on importation rates under current levels of control and within-island mobility. The simulations were then used to evaluate the impact of two pre-erythrocytic vaccine distribution strategies: mass treat and vaccinate, and prophylactic vaccination for off-island travellers. Lastly, a sensitivity analysis was performed through an ensemble of simulations fit to the Bayesian joint posterior probability distribution of the geospatial prevalence estimates. Results: The simulations suggest that in Malabo, an urban city containing 80% of the population, there are some pockets of residual transmission, but a large proportion of infections are acquired off-island by travellers to the mainland. Outside of Malabo, prevalence was mainly attributable to local transmission. The uncertainty in the local transmission vs. importation is lowest within Malabo and highest outside. Using a pre-erythrocytic vaccine to protect travellers would have larger benefits than using the vaccine to protect residents of Bioko Island from local transmission. In simulations, mass treatment and vaccination had short-lived benefits, as malaria prevalence returned to current levels as the vaccine's efficacy waned. Prophylactic vaccination of travellers resulted in longer-lasting reductions in prevalence. These projections were robust to underlying uncertainty in prevalence estimates. Conclusions: The modelled outcomes suggest that the volume of malaria cases imported from the mainland is a partial driver of continued endemic malaria on Bioko Island, and that continued elimination efforts on must account for human travel activity. [ABSTRACT FROM AUTHOR]

Background: Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired. Results and discussion: We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase. Conclusion: With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes. [ABSTRACT FROM AUTHOR]

Background: Many public health interventions lead to disruption or decrease of transmission, providing a beneficial effect for people in the population regardless of whether or not they individually participate in the intervention. This protective benefit has been referred to as a herd or community effect and is dependent on sufficient population participation. In practice, public health interventions are implemented at different spatial scales (i.e., at the village, district, or provincial level). Populations, however defined (i.e., neighbourhoods, villages, districts) are frequently connected to other populations through human movement or travel, and this connectedness can influence potential herd effects. Methods: The impact of a public health intervention (mass drug administration for malaria) was modelled, for different levels of connectedness between populations that have similar disease epidemiology (e.g., two nearby villages which have similar baseline malaria incidences and similar malaria intervention measures), or between populations of varying disease epidemiology (e.g., two nearby villages which have different baseline malaria incidences and/or malaria intervention measures). Results: The overall impact of the interventions deployed could be influenced either positively (adding value to the intervention) or negatively (reducing the impact of the intervention) by how much the intervention units are connected with each other (e.g., how frequent people go to the other village or town) and how different the disease intensity between them are. This phenomenon is termed the "assembly effect", and it is a meta-population version of the more commonly understood "herd effect". Conclusions: The connectedness of intervention units or populations is an important factor to be considered to achieve success in public health interventions that could provide herd effects. Appreciating the assembly effect can improve the cost-effective strategies for global disease elimination projects. [ABSTRACT FROM AUTHOR]

Background: Cancer stem cells are important for the development of many solid tumors. These cells receive promoting and inhibitory signals that depend on the nature of their environment (their niche) and determine cell dynamics. Mechanical stresses are crucial to the initiation and interpretation of these signals.Methods: A two-population mathematical model of tumorsphere growth is used to interpret the results of a series of experiments recently carried out in Tianjin, China, and extract information about the intraspecific and interspecific interactions between cancer stem cell and differentiated cancer cell populations.Results: The model allows us to reconstruct the time evolution of the cancer stem cell fraction, which was not directly measured. We find that, in the presence of stem cell growth factors, the interspecific cooperation between cancer stem cells and differentiated cancer cells induces a positive feedback loop that determines growth, independently of substrate hardness. In a frustrated attempt to reconstitute the stem cell niche, the number of cancer stem cells increases continuously with a reproduction rate that is enhanced by a hard substrate. For growth on soft agar, intraspecific interactions are always inhibitory, but on hard agar the interactions between stem cells are collaborative while those between differentiated cells are strongly inhibitory. Evidence also suggests that a hard substrate brings about a large fraction of asymmetric stem cell divisions. In the absence of stem cell growth factors, the barrier to differentiation is broken and overall growth is faster, even if the stem cell number is conserved.Conclusions: Our interpretation of the experimental results validates the centrality of the concept of stem cell niche when tumor growth is fueled by cancer stem cells. Niche memory is found to be responsible for the characteristic population dynamics observed in tumorspheres. The model also shows why substratum stiffness has a deep influence on the behavior of cancer stem cells, stiffer substrates leading to a larger proportion of asymmetric doublings. A specific condition for the growth of the cancer stem cell number is also obtained. [ABSTRACT FROM AUTHOR]

Background: Zika virus (ZIKV) emerged as a global epidemic in 2015-2016 from Latin America with its true geographical extent remaining unclear due to widely presumed underreporting. The identification of locations with potential and unknown spread of ZIKV is a key yet understudied component for outbreak preparedness. Here, we aim to identify locations at a high risk of cryptic ZIKV spread during 2015-2016 to further the understanding of the global ZIKV epidemiology, which is critical for the mitigation of the risk of future epidemics.Methods: We developed an importation simulation model to estimate the weekly number of ZIKV infections imported in each susceptible spatial unit (i.e. location that did not report any autochthonous Zika cases during 2015-2016), integrating epidemiological, demographic, and travel data as model inputs. Thereafter, a global risk model was applied to estimate the weekly ZIKV transmissibility during 2015-2016 for each location. Finally, we assessed the risk of onward ZIKV spread following importation in each susceptible spatial unit to identify locations with a high potential for cryptic ZIKV spread during 2015-2016.Results: We have found 24 susceptible spatial units that were likely to have experienced cryptic ZIKV spread during 2015-2016, of which 10 continue to have a high risk estimate within a highly conservative scenario, namely, Luanda in Angola, Banten in Indonesia, Maharashtra in India, Lagos in Nigeria, Taiwan and Guangdong in China, Dakar in Senegal, Maputo in Mozambique, Kinshasa in Congo DRC, and Pool in Congo. Notably, among the 24 susceptible spatial units identified, some have reported their first ZIKV outbreaks since 2017, thus adding to the credibility of our results (derived using 2015-2016 data only).Conclusion: Our study has provided valuable insights into the potentially high-risk locations for cryptic ZIKV circulation during the 2015-2016 pandemic and has also laid a foundation for future studies that attempt to further narrow this key knowledge gap. Our modelling framework can be adapted to identify areas with likely unknown spread of other emerging vector-borne diseases, which has important implications for public health readiness especially in resource-limited settings. [ABSTRACT FROM AUTHOR]

Background: Long-term care facilities (LTCFs) are vulnerable to outbreaks of coronavirus disease 2019 (COVID-19). Timely epidemiological surveillance is essential for outbreak response, but is complicated by a high proportion of silent (non-symptomatic) infections and limited testing resources.Methods: We used a stochastic, individual-based model to simulate transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) along detailed inter-individual contact networks describing patient-staff interactions in a real LTCF setting. We simulated distribution of nasopharyngeal swabs and reverse transcriptase polymerase chain reaction (RT-PCR) tests using clinical and demographic indications and evaluated the efficacy and resource-efficiency of a range of surveillance strategies, including group testing (sample pooling) and testing cascades, which couple (i) testing for multiple indications (symptoms, admission) with (ii) random daily testing.Results: In the baseline scenario, randomly introducing a silent SARS-CoV-2 infection into a 170-bed LTCF led to large outbreaks, with a cumulative 86 (95% uncertainty interval 6-224) infections after 3 weeks of unmitigated transmission. Efficacy of symptom-based screening was limited by lags to symptom onset and silent asymptomatic and pre-symptomatic transmission. Across scenarios, testing upon admission detected just 34-66% of patients infected upon LTCF entry, and also missed potential introductions from staff. Random daily testing was more effective when targeting patients than staff, but was overall an inefficient use of limited resources. At high testing capacity (> 10 tests/100 beds/day), cascades were most effective, with a 19-36% probability of detecting outbreaks prior to any nosocomial transmission, and 26-46% prior to first onset of COVID-19 symptoms. Conversely, at low capacity (< 2 tests/100 beds/day), group testing strategies detected outbreaks earliest. Pooling randomly selected patients in a daily group test was most likely to detect outbreaks prior to first symptom onset (16-27%), while pooling patients and staff expressing any COVID-like symptoms was the most efficient means to improve surveillance given resource limitations, compared to the reference requiring only 6-9 additional tests and 11-28 additional swabs to detect outbreaks 1-6 days earlier, prior to an additional 11-22 infections.Conclusions: COVID-19 surveillance is challenged by delayed or absent clinical symptoms and imperfect diagnostic sensitivity of standard RT-PCR tests. In our analysis, group testing was the most effective and efficient COVID-19 surveillance strategy for resource-limited LTCFs. Testing cascades were even more effective given ample testing resources. Increasing testing capacity and updating surveillance protocols accordingly could facilitate earlier detection of emerging outbreaks, informing a need for urgent intervention in settings with ongoing nosocomial transmission. [ABSTRACT FROM AUTHOR]

Background: The emergence of a novel coronavirus (SARS-CoV-2) in Wuhan, China, at the end of 2019 has caused widespread transmission around the world. As new epicentres in Europe and America have arisen, of particular concern is the increased number of imported coronavirus disease 2019 (COVID-19) cases in Africa, where the impact of the pandemic could be more severe. We aim to estimate the number of COVID-19 cases imported from 12 major epicentres in Europe and America to each African country, as well as the probability of reaching 10,000 cases in total by the end of March, April, May, and June following viral introduction.Methods: We used the reported number of cases imported from the 12 major epicentres in Europe and America to Singapore, as well as flight data, to estimate the number of imported cases in each African country. Under the assumption that Singapore has detected all the imported cases, the estimates for Africa were thus conservative. We then propagated the uncertainty in the imported case count estimates to simulate the onward spread of the virus, until 10,000 cases are reached or the end of June, whichever is earlier. Specifically, 1,000 simulations were run separately under four different combinations of parameter values to test the sensitivity of our results.Results: We estimated Morocco, Algeria, South Africa, Egypt, Tunisia, and Nigeria as having the largest number of COVID-19 cases imported from the 12 major epicentres. Based on our 1,000 simulation runs, Morocco and Algeria's estimated probability of reaching 10,000 cases by end of March was close to 100% under all scenarios. In particular, we identified countries with less than 1,000 cases in total reported by end of June whilst the estimated probability of reaching 10,000 cases by then was higher than 50% even under the most optimistic scenario.Conclusions: Our study highlights particular countries that are likely to reach (or have reached) 10,000 cases far earlier than the reported data suggest, calling for the prioritization of resources to mitigate the further spread of the epidemic. [ABSTRACT FROM AUTHOR]

Background: Seasonal influenza poses a significant public health and economic burden, associated with the outcome of infection and resulting complications. The true burden of the disease is difficult to capture due to the wide range of presentation, from asymptomatic cases to non-respiratory complications such as cardiovascular events, and its seasonal variability. An understanding of the magnitude of the true annual incidence of influenza is important to support prevention and control policy development and to evaluate the impact of preventative measures such as vaccination. Methods: We use a dynamic disease transmission model, laboratory-confirmed influenza surveillance data, and randomized-controlled trial (RCT) data to quantify the underestimation factor, expansion factor, and symptomatic influenza illnesses in the US and Canada during the 2011-2012 and 2012-2013 influenza seasons. Results: Based on 2 case definitions, we estimate between 0.42−3.2% and 0.33−1.2% of symptomatic influenza illnesses were laboratory-confirmed in Canada during the 2011-2012 and 2012-2013 seasons, respectively. In the US, we estimate between 0.08−0.61% and 0.07−0.33% of symptomatic influenza illnesses were laboratory-confirmed in the 2011-2012 and 2012-2013 seasons, respectively. We estimated the symptomatic influenza illnesses in Canada to be 0.32−2.4 million in 2011-2012 and 1.8−8.2 million in 2012-2013. In the US, we estimate the number of symptomatic influenza illnesses to be 4.4−34 million in 2011-2012 and 23−102 million in 2012-2013. Conclusions: We illustrate that monitoring a representative group within a population may aid in effectively modelling the transmission of infectious diseases such as influenza. In particular, the utilization of RCTs in models may enhance the accuracy of epidemiological parameter estimation. [ABSTRACT FROM AUTHOR]

Background: Ability to adapt to temperature changes trough the Heat Shock Response (HSR) pathways is one of the most fundamental and clinically relevant cellular response systems. Heat Shock (HS) affects the signalling and gene expression responses of the Nuclear Factor κB (NF-κB) transcription factor, a critical regulator of proliferation and inflammation, however, our quantitative understanding of how cells sense and adapt to temperature changes is limited. Methods: We used live-cell time-lapse microscopy and mathematical modelling to understand the signalling of the NF-κB system in the human MCF7 breast adenocarcinoma cells in response to pro-inflammatory Interleukin 1β (IL1β) and Tumour Necrosis Factor α (TNFα) cytokines, following exposure to a 37–43 °C range of physiological and clinical temperatures. Results: We show that exposure to 43 °C 1 h HS inhibits the immediate NF-κB signalling response to TNFα and IL1β stimulation although uptake of cytokines is not impaired. Within 4 h after HS treatment IL1β-induced NF-κB responses return to normal levels, but the recovery of the TNFα-induced responses is still affected. Using siRNA knock-down of Heat Shock Factor 1 (HSF1) we show that this stimulus-specificity is conferred via the Inhibitory κB kinase (IKK) signalosome where HSF1-dependent feedback regulates TNFα, but not IL1β-mediated IKK recovery post HS. Furthermore, we demonstrate that through the temperature-dependent denaturation and recovery of IKK, TNFα and IL1β-mediated signalling exhibit different temperature sensitivity and adaptation to repeated HS when exposed to a 37–43 °C temperature range. Specifically, IL1β-mediated NF-κB responses are more robust to temperature changes in comparison to those induced by TNFα treatment. Conclusions: We demonstrate that the kinetics of the NF-κB system following temperature stress is cytokine specific and exhibit differential adaptation to temperature changes. We propose that this differential temperature sensitivity is mediated via the IKK signalosome, which acts as a bona fide temperature sensor trough the HSR cross-talk. This novel quantitative understanding of NF-κB and HSR interactions is fundamentally important for the potential optimization of therapeutic hyperthermia protocols. Video Abstract [ABSTRACT FROM AUTHOR]

Keywords: COVID-19; Novel coronavirus; SARS-CoV-2; Mathematical modelling; Compartmental models; Forecasting; Non-pharmaceutical interventions; Disease control; Flatten the curve; Lockdown EN COVID-19 Novel coronavirus SARS-CoV-2 Mathematical modelling Compartmental models Forecasting Non-pharmaceutical interventions Disease control Flatten the curve Lockdown 1 4 4 05/27/20 20200525 NES 200525 Background The coronavirus disease 2019 (COVID-19) pandemic has been responsible for over three million reported cases worldwide, including more than 200,000 deaths (as of 1 May 2020). We discuss the authors' results about strategies for intervening to maintain intensive care unit (ICU) bed demand below the number of ICU beds in the UK. [Extracted from the article]

Background: More than ever, it is crucial to make the best use of existing country data, and analytical tools for developing malaria control strategies as the heterogeneity in malaria risk within countries is increasing, and the available malaria control tools are expanding while large funding gaps exist. Global and local policymakers, as well as funders, increasingly recognize the value of mathematical modelling as a strategic tool to support decision making. This case study article describes the long-term use of modelling in close collaboration with the National Malaria Control Programme (NMCP) in Tanzania, the challenges encountered and lessons learned. Case description: In Tanzania, a recent rebound in prevalence led to the revision of the national malaria strategic plan with interventions targeted to the malaria risk at the sub-regional level. As part of the revision, a mathematical malaria modelling framework for setting specific predictions was developed and used between 2016 and 2019 to (1) reproduce setting specific historical malaria trends, and (2) to simulate in silico the impact of future interventions. Throughout the project, multiple stakeholder workshops were attended and the use of mathematical modelling interactively discussed. Evaluation: In Tanzania, the model application created an interdisciplinary and multisectoral dialogue platform between modellers, NMCP and partners and contributed to the revision of the national malaria strategic plan by simulating strategies suggested by the NMCP. The uptake of the modelling outputs and sustained interest by the NMCP were critically associated with following factors: (1) effective sensitization to the NMCP, (2) regular and intense communication, (3) invitation for the modellers to participate in the strategic plan process, and (4) model application tailored to the local context. Conclusion: Empirical data analysis and its use for strategic thinking remain the cornerstone for evidence-based decision-making. Mathematical impact modelling can support the process both by unifying all stakeholders in one strategic process and by adding new key evidence required for optimized decision-making. However, without a long-standing partnership, it will be much more challenging to sensibilize programmes to the usefulness and sustained use of modelling and local resources within the programme or collaborating research institutions need to be mobilized. [ABSTRACT FROM AUTHOR]

Background: Individualization and patient-specific optimization of treatment is a major goal of modern health care. One way to achieve this goal is the application of high-resolution diagnostics together with the application of targeted therapies. However, the rising number of different treatment modalities also induces new challenges: Whereas randomized clinical trials focus on proving average treatment effects in specific groups of patients, direct conclusions at the individual patient level are problematic. Thus, the identification of the best patient-specific treatment options remains an open question. Systems medicine, specifically mechanistic mathematical models, can substantially support individual treatment optimization. In addition to providing a better general understanding of disease mechanisms and treatment effects, these models allow for an identification of patient-specific parameterizations and, therefore, provide individualized predictions for the effect of different treatment modalities.Results: In the following we describe a software framework that facilitates the integration of mathematical models and computer simulations into routine clinical processes to support decision-making. This is achieved by combining standard data management and data exploration tools, with the generation and visualization of mathematical model predictions for treatment options at an individual patient level.Conclusions: By integrating model results in an audit trail compatible manner into established clinical workflows, our framework has the potential to foster the use of systems-medical approaches in clinical practice. We illustrate the framework application by two use cases from the field of haematological oncology. [ABSTRACT FROM AUTHOR]

Background: Nutrition plays a crucial role in regulating reproductive hormones and follicular development in cattle. This is visible particularly during the time of negative energy balance at the onset of milk production after calving. Here, elongated periods of anovulation have been observed, resulting from alterations in luteinizing hormone concentrations, likely caused by lower glucose and insulin concentrations in the blood. The mechanisms that result in a reduced fertility are not completely understood, although a close relationship to the glucose-insulin metabolism is widely supported. Results: Following this idea, we developed a mathematical model of the hormonal network combining reproductive hormones and hormones that are coupled to the glucose compartments within the body of the cow. The model is built on ordinary differential equations and relies on previously introduced models on the bovine estrous cycle and the glucose-insulin dynamics. Necessary modifications and coupling mechanisms are thoroughly discussed. Depending on the composition and the amount of feed, in particular the glucose content in the dry matter, the model quantifies reproductive hormones and follicular development over time. Simulation results for different nutritional regimes in lactating and non-lactating dairy cows are examined and compared with experimental studies. The simulations describe realistically the effects of nutritional glucose supply on the ovulatory cycle of dairy cattle. Conclusions: The mathematical model enables the user to explore the relationship between nutrition and reproduction by running simulations and performing parameter studies. Regarding its applicability, this work is an early attempt towards developing in silico feeding strategies and may eventually help to refine and reduce animal experiments. Reviewers: This article was reviewed by John McNamara and Tin Pang (nominated by Martin Lercher). [ABSTRACT FROM AUTHOR]