Your browser doesn't support javascript.
loading
: 20 | 50 | 100
1 - 20 de 47
1.
Nat Commun ; 15(1): 3230, 2024 Apr 22.
Article En | MEDLINE | ID: mdl-38649361

Despite concern that climate change could increase the human risk to malaria in certain areas, the temperature dependency of malaria transmission is poorly characterized. Here, we use a mechanistic model fitted to experimental data to describe how Plasmodium falciparum infection of the African malaria vector, Anopheles gambiae, is modulated by temperature, including its influences on parasite establishment, conversion efficiency through parasite developmental stages, parasite development rate, and overall vector competence. We use these data, together with estimates of the survival of infected blood-fed mosquitoes, to explore the theoretical influence of temperature on transmission in four locations in Kenya, considering recent conditions and future climate change. Results provide insights into factors limiting transmission in cooler environments and indicate that increases in malaria transmission due to climate warming in areas like the Kenyan Highlands, might be less than previously predicted.


Anopheles , Malaria, Falciparum , Mosquito Vectors , Plasmodium falciparum , Temperature , Plasmodium falciparum/physiology , Malaria, Falciparum/transmission , Malaria, Falciparum/parasitology , Malaria, Falciparum/epidemiology , Animals , Anopheles/parasitology , Humans , Kenya/epidemiology , Mosquito Vectors/parasitology , Climate Change , Female
2.
PLoS Pathog ; 20(4): e1011975, 2024 Apr.
Article En | MEDLINE | ID: mdl-38557892

Arboviruses can emerge rapidly and cause explosive epidemics of severe disease. Some of the most epidemiologically important arboviruses, including dengue virus (DENV), Zika virus (ZIKV), Chikungunya (CHIKV) and yellow fever virus (YFV), are transmitted by Aedes mosquitoes, most notably Aedes aegypti and Aedes albopictus. After a mosquito blood feeds on an infected host, virus enters the midgut and infects the midgut epithelium. The virus must then overcome a series of barriers before reaching the mosquito saliva and being transmitted to a new host. The virus must escape from the midgut (known as the midgut escape barrier; MEB), which is thought to be mediated by transient changes in the permeability of the midgut-surrounding basal lamina layer (BL) following blood feeding. Here, we present a mathematical model of the within-mosquito population dynamics of DENV (as a model system for mosquito-borne viruses more generally) that includes the interaction of the midgut and BL which can account for the MEB. Our results indicate a dose-dependency of midgut establishment of infection as well as rate of escape from the midgut: collectively, these suggest that the extrinsic incubation period (EIP)-the time taken for DENV virus to be transmissible after infection-is shortened when mosquitoes imbibe more virus. Additionally, our experimental data indicate that multiple blood feeding events, which more closely mimic mosquito-feeding behavior in the wild, can hasten the course of infections, and our model predicts that this effect is sensitive to the amount of virus imbibed. Our model indicates that mutations to the virus which impact its replication rate in the midgut could lead to even shorter EIPs when double-feeding occurs. Mechanistic models of within-vector viral infection dynamics provide a quantitative understanding of infection dynamics and could be used to evaluate novel interventions that target the mosquito stages of the infection.


Aedes , Dengue Virus , Dengue , Zika Virus Infection , Zika Virus , Animals , Gastrointestinal Tract , Mosquito Vectors
3.
Clocks Sleep ; 6(1): 114-128, 2024 Feb 26.
Article En | MEDLINE | ID: mdl-38534797

In humans, the nocturnal secretion of melatonin by the pineal gland is suppressed by ocular exposure to light. In the laboratory, melatonin suppression is a biomarker for this neuroendocrine pathway. Recent work has found that individuals differ substantially in their melatonin-suppressive response to light, with the most sensitive individuals being up to 60 times more sensitive than the least sensitive individuals. Planning experiments with melatonin suppression as an outcome needs to incorporate these individual differences, particularly in common resource-limited scenarios where running within-subjects studies at multiple light levels is costly and resource-intensive and may not be feasible with respect to participant compliance. Here, we present a novel framework for virtual laboratory melatonin suppression experiments, incorporating a Bayesian statistical model. We provide a Shiny web app for power analyses that allows users to modify various experimental parameters (sample size, individual-level heterogeneity, statistical significance threshold, light levels), and simulate a systematic shift in sensitivity (e.g., due to a pharmacological or other intervention). Our framework helps experimenters to design compelling and robust studies, offering novel insights into the underlying biological variability in melatonin suppression relevant for practical applications.

4.
J R Soc Interface ; 21(212): 20230369, 2024 03.
Article En | MEDLINE | ID: mdl-38442857

Most ordinary differential equation (ODE) models used to describe biological or physical systems must be solved approximately using numerical methods. Perniciously, even those solvers that seem sufficiently accurate for the forward problem, i.e. for obtaining an accurate simulation, might not be sufficiently accurate for the inverse problem, i.e. for inferring the model parameters from data. We show that for both fixed step and adaptive step ODE solvers, solving the forward problem with insufficient accuracy can distort likelihood surfaces, which might become jagged, causing inference algorithms to get stuck in local 'phantom' optima. We demonstrate that biases in inference arising from numerical approximation of ODEs are potentially most severe in systems involving low noise and rapid nonlinear dynamics. We reanalyse an ODE change point model previously fit to the COVID-19 outbreak in Germany and show the effect of the step size on simulation and inference results. We then fit a more complicated rainfall run-off model to hydrological data and illustrate the importance of tuning solver tolerances to avoid distorted likelihood surfaces. Our results indicate that, when performing inference for ODE model parameters, adaptive step size solver tolerances must be set cautiously and likelihood surfaces should be inspected for characteristic signs of numerical issues.


Algorithms , COVID-19 , Humans , COVID-19/epidemiology , Computer Simulation , Disease Outbreaks , Germany
5.
Lancet Infect Dis ; 2023 Nov 28.
Article En | MEDLINE | ID: mdl-38040006

The 2023 Marburg virus disease outbreaks in Equatorial Guinea and Tanzania highlighted the importance of better understanding this lethal pathogen. We did a systematic review (PROSPERO CRD42023393345) of peer-reviewed articles reporting historical outbreaks, modelling studies, and epidemiological parameters focused on Marburg virus disease. We searched PubMed and Web of Science from database inception to March 31, 2023. Two reviewers evaluated all titles and abstracts with consensus-based decision making. To ensure agreement, 13 (31%) of 42 studies were double-extracted and a custom-designed quality assessment questionnaire was used for risk of bias assessment. We present detailed information on 478 reported cases and 385 deaths from Marburg virus disease. Analysis of historical outbreaks and seroprevalence estimates suggests the possibility of undetected Marburg virus disease outbreaks, asymptomatic transmission, or cross-reactivity with other pathogens, or a combination of these. Only one study presented a mathematical model of Marburg virus transmission. We estimate an unadjusted, pooled total random effect case fatality ratio of 61·9% (95% CI 38·8-80·6; I2=93%). We identify epidemiological parameters relating to transmission and natural history, for which there are few estimates. This systematic review and the accompanying database provide a comprehensive overview of Marburg virus disease epidemiology and identify key knowledge gaps, contributing crucial information for mathematical models to support future Marburg virus disease epidemic responses.

6.
Front Pharmacol ; 14: 1270443, 2023.
Article En | MEDLINE | ID: mdl-37927586

Treatment response variability across patients is a common phenomenon in clinical practice. For many drugs this inter-individual variability does not require much (if any) individualisation of dosing strategies. However, for some drugs, including chemotherapies and some monoclonal antibody treatments, individualisation of dosages are needed to avoid harmful adverse events. Model-informed precision dosing (MIPD) is an emerging approach to guide the individualisation of dosing regimens of otherwise difficult-to-administer drugs. Several MIPD approaches have been suggested to predict dosing strategies, including regression, reinforcement learning (RL) and pharmacokinetic and pharmacodynamic (PKPD) modelling. A unified framework to study the strengths and limitations of these approaches is missing. We develop a framework to simulate clinical MIPD trials, providing a cost and time efficient way to test different MIPD approaches. Central for our framework is a clinical trial model that emulates the complexities in clinical practice that challenge successful treatment individualisation. We demonstrate this framework using warfarin treatment as a use case and investigate three popular MIPD methods: 1. Neural network regression; 2. Deep RL; and 3. PKPD modelling. We find that the PKPD model individualises warfarin dosing regimens with the highest success rate and the highest efficiency: 75.1% of the individuals display INRs inside the therapeutic range at the end of the simulated trial; and the median time in the therapeutic range (TTR) is 74%. In comparison, the regression model and the deep RL model have success rates of 47.0% and 65.8%, and median TTRs of 45% and 68%. We also find that the MIPD models can attain different degrees of individualisation: the Regression model individualises dosing regimens up to variability explained by covariates; the Deep RL model and the PKPD model individualise dosing regimens accounting also for additional variation using monitoring data. However, the Deep RL model focusses on control of the treatment response, while the PKPD model uses the data also to further the individualisation of predictions.

7.
bioRxiv ; 2023 Sep 29.
Article En | MEDLINE | ID: mdl-37808804

Flaviviruses are arthropod-borne (arbo)viruses which can emerge rapidly and cause explosive epidemics of severe disease. Some of the most epidemiologically important flaviviruses, including dengue virus (DENV), Zika virus (ZIKV) and yellow fever virus (YFV), are transmitted by Aedes mosquitoes, most notably Aedes aegypti and Aedes albopictus. After a mosquito blood feeds on an infected host, virus enters the midgut and infects the midgut epithelium. The virus must then overcome a series of barriers before reaching the mosquito saliva and being transmitted to a new host. The virus must escape from the midgut (known as the midgut escape barrier; MEB), which is thought to be mediated by transient changes in the permeability of the midgut-surrounding basal lamina layer (BL) following blood feeding. Here, we present a mathematical model of the within-mosquito population dynamics of flaviviruses that includes the interaction of the midgut and BL which can account for the MEB. Our results indicate a dose-dependency of midgut establishment of infection as well as rate of escape from the midgut: collectively, these suggest that the extrinsic incubation period (EIP) - the time taken for DENV virus to be transmissible after infection - is shortened when mosquitoes imbibe more virus. Additionally, our experimental data indicates that multiple blood feeding events, which more closely mimic mosquito-feeding behavior in the wild, can hasten the course of infections, and our model predicts that this effect is sensitive to the amount of virus imbibed. Our model indicates that mutations to the virus which impact its replication rate in the midgut could lead to even shorter EIPs when double-feeding occurs. Mechanistic models of within-vector viral infection dynamics provide a quantitative understanding of infection dynamics and could be used to evaluate novel interventions that target the mosquito stages of the infection. Author summary: Aedes mosquitoes are the main vectors of dengue virus (DENV), Zika virus (ZIKV) and yellow fever virus (YFV), all of which can cause severe disease in humans with dengue alone infecting an estimated 100-400 million people each year. Understanding the processes that affect whether, and at which rate, mosquitoes may transmit such viruses is, hence, paramount. Here, we present a mathematical model of virus dynamics within infected mosquitoes. By combining the model with novel experimental data, we show that the course of infection is sensitive to the initial dose of virus ingested by the mosquito. The data also indicates that mosquitoes which blood feed subsequent to becoming infected may be able to transmit infection earlier, which is reproduced in the model. This is important as many mosquito species feed multiple times during their lifespan and, any reduction in time to dissemination will increase the number of days that a mosquito is infectious and so enhance the risk of transmission. Our study highlights the key and complementary roles played by mathematical models and experimental data for understanding within-mosquito virus dynamics.

8.
Proc Biol Sci ; 290(2007): 20231664, 2023 09 27.
Article En | MEDLINE | ID: mdl-37752839

We introduce the angular reproduction number Ω, which measures time-varying changes in epidemic transmissibility resulting from variations in both the effective reproduction number R, and generation time distribution w. Predominant approaches for tracking pathogen spread infer either R or the epidemic growth rate r. However, R is biased by mismatches between the assumed and true w, while r is difficult to interpret in terms of the individual-level branching process underpinning transmission. R and r may also disagree on the relative transmissibility of epidemics or variants (i.e. rA > rB does not imply RA > RB for variants A and B). We find that Ω responds meaningfully to mismatches and time-variations in w while mostly maintaining the interpretability of R. We prove that Ω > 1 implies R > 1 and that Ω agrees with r on the relative transmissibility of pathogens. Estimating Ω is no more difficult than inferring R, uses existing software, and requires no generation time measurements. These advantages come at the expense of selecting one free parameter. We propose Ω as complementary statistic to R and r that improves transmissibility estimates when w is misspecified or time-varying and better reflects the impact of interventions, when those interventions concurrently change R and w or alter the relative risk of co-circulating pathogens.


Disease Outbreaks , Epidemics , Basic Reproduction Number , Software
9.
Article En | MEDLINE | ID: mdl-37456558

Insecticide resistance is a growing problem that risks harming the progress made by vector control tools in reducing the malaria burden globally. New methods for quantifying the extent of resistance in wild populations are urgently needed to guide deployment of interventions to improve disease control. Intensity bioassays measure mosquito mortality at a range of insecticide doses and characterise phenotypic resistance in regions where resistance is already detected. These data are increasingly being collected but tend to exhibit high measurement error and there is a lack of formal guidelines on how they should be analysed or compared. This paper introduces a novel Bayesian framework for analysing intensity bioassay data, which uses a flexible statistical model able to capture a wide variety of relationships between mortality and insecticide dose. By accounting for background mortality of mosquitoes, our approach minimises the impact of this source of measurement noise resulting in more precise quantification of resistance. It outputs a range of metrics for describing the intensity and variability in resistance within the sample and quantifies the level of measurement error in the assay. The functionality is illustrated with data from laboratory-reared mosquitoes to show how the lethal dose varies within and between different strains. The framework can also be used to formally test hypotheses by explicitly considering the high heterogeneity seen in these types of data in field samples. Here we show that the intensity of resistance (as measured by the median lethal dose (LC50) of insecticide) increases over 7 years in mosquitoes from one village in Burkina Faso but remains constant in another. This work showcases the benefits of statistically rigorous analysis of insecticide bioassay data and highlights the additional information available from this and other dose-response data.

10.
Science ; 381(6655): 336-343, 2023 07 21.
Article En | MEDLINE | ID: mdl-37471538

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.


COVID-19 , SARS-CoV-2 , Humans , Africa, Southern , COVID-19/transmission , COVID-19/virology , Genomics , SARS-CoV-2/classification , SARS-CoV-2/genetics , SARS-CoV-2/pathogenicity , Phylogeny
11.
Elife ; 122023 Jun 09.
Article En | MEDLINE | ID: mdl-37294299

The factors leading to the global emergence of Enterovirus D68 (EV-D68) in 2014 as a cause of acute flaccid myelitis (AFM) in children are unknown. To investigate potential changes in virus transmissibility or population susceptibility, we measured the seroprevalence of EV-D68-specific neutralising antibodies in serum samples collected in England in 2006, 2011, and 2017. Using catalytic mathematical models, we estimate an approximately 50% increase in the annual probability of infection over the 10-year study period, coinciding with the emergence of clade B around 2009. Despite such increase in transmission, seroprevalence data suggest that the virus was already widely circulating before the AFM outbreaks and the increase of infections by age cannot explain the observed number of AFM cases. Therefore, the acquisition of or an increase in neuropathogenicity would be additionally required to explain the emergence of outbreaks of AFM. Our results provide evidence that changes in enterovirus phenotypes cause major changes in disease epidemiology.


Enterovirus D, Human , Enterovirus Infections , Neuromuscular Diseases , Humans , Seroepidemiologic Studies , Disease Outbreaks , Enterovirus Infections/epidemiology , England/epidemiology
12.
PLoS Comput Biol ; 19(5): e1011135, 2023 05.
Article En | MEDLINE | ID: mdl-37216399

Variability is an intrinsic property of biological systems and is often at the heart of their complex behaviour. Examples range from cell-to-cell variability in cell signalling pathways to variability in the response to treatment across patients. A popular approach to model and understand this variability is nonlinear mixed effects (NLME) modelling. However, estimating the parameters of NLME models from measurements quickly becomes computationally expensive as the number of measured individuals grows, making NLME inference intractable for datasets with thousands of measured individuals. This shortcoming is particularly limiting for snapshot datasets, common e.g. in cell biology, where high-throughput measurement techniques provide large numbers of single cell measurements. We introduce a novel approach for the estimation of NLME model parameters from snapshot measurements, which we call filter inference. Filter inference uses measurements of simulated individuals to define an approximate likelihood for the model parameters, avoiding the computational limitations of traditional NLME inference approaches and making efficient inferences from snapshot measurements possible. Filter inference also scales well with the number of model parameters, using state-of-the-art gradient-based MCMC algorithms such as the No-U-Turn Sampler (NUTS). We demonstrate the properties of filter inference using examples from early cancer growth modelling and from epidermal growth factor signalling pathway modelling.


Algorithms , Nonlinear Dynamics , Humans , Time Factors , Probability
13.
J Surg Case Rep ; 2023(3): rjad189, 2023 Mar.
Article En | MEDLINE | ID: mdl-36998254

We present a unique case of bowel obstruction with a hiatus hernia causing atypical chest pain with dynamic ST-segment elevation in a regional Australian emergency department. The ST elevation only resolved after nasogastric decompression of the bowel obstruction. Early thrombolysis of presumed myocardial infarction led to upper gastrointestinal tract bleeding that could have been avoided with timely diagnosis. An extensive review of literature, in addition to our case report, suggests bowel obstruction is a differential diagnosis for patients who have inferior pattern ST elevation but normal troponin presenting with atypical chest pain, nausea, vomiting and previous abdominal surgery.

14.
Parasit Vectors ; 16(1): 21, 2023 Jan 20.
Article En | MEDLINE | ID: mdl-36670470

BACKGROUND: The continued spread of insecticide resistance in mosquito vectors of malaria and arboviral diseases may lead to operational failure of insecticide-based interventions if resistance is not monitored and managed efficiently. This study aimed to develop and validate a new WHO glass bottle bioassay method as an alternative to the WHO standard insecticide tube test to monitor mosquito susceptibility to new public health insecticides with particular modes of action, physical properties or both. METHODS: A multi-centre study involving 21 laboratories worldwide generated data on the susceptibility of seven mosquito species (Aedes aegypti, Aedes albopictus, Anopheles gambiae sensu stricto [An. gambiae s.s.], Anopheles funestus, Anopheles stephensi, Anopheles minimus and Anopheles albimanus) to seven public health insecticides in five classes, including pyrethroids (metofluthrin, prallethrin and transfluthrin), neonicotinoids (clothianidin), pyrroles (chlorfenapyr), juvenile hormone mimics (pyriproxyfen) and butenolides (flupyradifurone), in glass bottle assays. The data were analysed using a Bayesian binomial model to determine the concentration-response curves for each insecticide-species combination and to assess the within-bioassay variability in the susceptibility endpoints, namely the concentration that kills 50% and 99% of the test population (LC50 and LC99, respectively) and the concentration that inhibits oviposition of the test population by 50% and 99% (OI50 and OI99), to measure mortality and the sterilizing effect, respectively. RESULTS: Overall, about 200,000 mosquitoes were tested with the new bottle bioassay, and LC50/LC99 or OI50/OI99 values were determined for all insecticides. Variation was seen between laboratories in estimates for some mosquito species-insecticide combinations, while other test results were consistent. The variation was generally greater with transfluthrin and flupyradifurone than with the other compounds tested, especially against Anopheles species. Overall, the mean within-bioassay variability in mortality and oviposition inhibition were < 10% for most mosquito species-insecticide combinations. CONCLUSION: Our findings, based on the largest susceptibility dataset ever produced on mosquitoes, showed that the new WHO bottle bioassay is adequate for evaluating mosquito susceptibility to new and promising public health insecticides currently deployed for vector control. The datasets presented in this study have been used recently by the WHO to establish 17 new insecticide discriminating concentrations (DCs) for either Aedes spp. or Anopheles spp. The bottle bioassay and DCs can now be widely used to monitor baseline insecticide susceptibility of wild populations of vectors of malaria and Aedes-borne diseases worldwide.


Anopheles , Insecticides , Malaria , Pyrethrins , Animals , Female , Insecticides/pharmacology , Mosquito Vectors , Public Health , Bayes Theorem , Mosquito Control/methods , Pyrethrins/pharmacology , Insecticide Resistance , Biological Assay , World Health Organization
15.
J Theor Biol ; 558: 111351, 2023 02 07.
Article En | MEDLINE | ID: mdl-36379231

Whether an outbreak of infectious disease is likely to grow or dissipate is determined through the time-varying reproduction number, Rt. Real-time or retrospective identification of changes in Rt following the imposition or relaxation of interventions can thus contribute important evidence about disease transmission dynamics which can inform policymaking. Here, we present a method for estimating shifts in Rt within a renewal model framework. Our method, which we call EpiCluster, is a Bayesian nonparametric model based on the Pitman-Yor process. We assume that Rt is piecewise-constant, and the incidence data and priors determine when or whether Rt should change and how many times it should do so throughout the series. We also introduce a prior which induces sparsity over the number of changepoints. Being Bayesian, our approach yields a measure of uncertainty in Rt and its changepoints. EpiCluster is fast, straightforward to use, and we demonstrate that it provides automated detection of rapid changes in transmission, either in real-time or retrospectively, for synthetic data series where the Rt profile is known. We illustrate the practical utility of our method by fitting it to case data of outbreaks of COVID-19 in Australia and Hong Kong, where it finds changepoints coinciding with the imposition of non-pharmaceutical interventions. Bayesian nonparametric methods, such as ours, allow the volume and complexity of the data to dictate the number of parameters required to approximate the process and should find wide application in epidemiology. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".


COVID-19 , Humans , Bayes Theorem , Retrospective Studies , COVID-19/epidemiology , Pandemics , Disease Outbreaks
16.
Nature ; 610(7930): 154-160, 2022 10.
Article En | MEDLINE | ID: mdl-35952712

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).


COVID-19 , SARS-CoV-2 , COVID-19/epidemiology , COVID-19/prevention & control , COVID-19/transmission , COVID-19/virology , Cities/epidemiology , Contact Tracing , England/epidemiology , Genome, Viral/genetics , Humans , Quarantine/legislation & jurisprudence , SARS-CoV-2/genetics , SARS-CoV-2/growth & development , SARS-CoV-2/isolation & purification , Travel/legislation & jurisprudence
17.
PLoS One ; 17(8): e0272442, 2022.
Article En | MEDLINE | ID: mdl-35981055

A large range of prognostic models for determining the risk of COVID-19 patient mortality exist, but these typically restrict the set of biomarkers considered to measurements available at patient admission. Additionally, many of these models are trained and tested on patient cohorts from a single hospital, raising questions about the generalisability of results. We used a Bayesian Markov model to analyse time series data of biomarker measurements taken throughout the duration of a COVID-19 patient's hospitalisation for n = 1540 patients from two hospitals in New York: State University of New York (SUNY) Downstate Health Sciences University and Maimonides Medical Center. Our main focus was to quantify the mortality risk associated with both static (e.g. demographic and patient history variables) and dynamic factors (e.g. changes in biomarkers) throughout hospitalisation, by so doing, to explain the observed patterns of mortality. By using our model to make predictions across the hospitals, we assessed how predictive factors generalised between the two cohorts. The individual dynamics of the measurements and their associated mortality risk were remarkably consistent across the hospitals. The model accuracy in predicting patient outcome (death or discharge) was 72.3% (predicting SUNY; posterior median accuracy) and 71.3% (predicting Maimonides) respectively. Model sensitivity was higher for detecting patients who would go on to be discharged (78.7%) versus those who died (61.8%). Our results indicate the utility of including dynamic clinical measurements when assessing patient mortality risk but also highlight the difficulty of identifying high risk patients.


COVID-19 , Bayes Theorem , Biomarkers , Hospitalization , Hospitals , Humans , New York/epidemiology , Retrospective Studies , SARS-CoV-2 , Time Factors
18.
Math Biosci ; 349: 108824, 2022 07.
Article En | MEDLINE | ID: mdl-35537550

The COVID-19 epidemic continues to rage in many parts of the world. In the UK alone, an array of mathematical models have played a prominent role in guiding policymaking. Whilst considerable pedagogical material exists for understanding the basics of transmission dynamics modelling, there is a substantial gap between the relatively simple models used for exposition of the theory and those used in practice to model the transmission dynamics of COVID-19. Understanding these models requires considerable prerequisite knowledge and presents challenges to those new to the field of epidemiological modelling. In this paper, we introduce an open-source R package, comomodels, which can be used to understand the complexities of modelling the transmission dynamics of COVID-19 through a series of differential equation models. Alongside the base package, we describe a host of learning resources, including detailed tutorials and an interactive web-based interface allowing dynamic investigation of the model properties. We then use comomodels to illustrate three key lessons in the transmission of COVID-19 within R Markdown vignettes.


COVID-19 , Epidemics , Humans , Learning , Models, Theoretical
19.
Sci Rep ; 12(1): 1402, 2022 01 26.
Article En | MEDLINE | ID: mdl-35082312

Burkina Faso has one of the highest malaria burdens in sub-Saharan Africa despite the mass deployment of insecticide-treated nets (ITNs) and use of seasonal malaria chemoprevention (SMC) in children aged up to 5 years. Identification of risk factors for Plasmodium falciparum infection in rural Burkina Faso could help to identify and target malaria control measures. A cross-sectional survey of 1,199 children and adults was conducted during the peak malaria transmission season in the Cascades Region of south-west Burkina Faso in 2017. Logistic regression was used to identify risk factors for microscopically confirmed P. falciparum infection. A malaria transmission dynamic model was used to determine the impact on malaria cases averted of administering SMC to children aged 5-15 year old. P. falciparum prevalence was 32.8% in the study population. Children aged 5 to < 10 years old were at 3.74 times the odds (95% CI = 2.68-5.22, P < 0.001) and children aged 10 to 15 years old at 3.14 times the odds (95% CI = 1.20-8.21, P = 0.02) of P. falciparum infection compared to children aged less than 5 years old. Administration of SMC to children aged up to 10 years is predicted to avert an additional 57 malaria cases per 1000 population per year (9.4% reduction) and administration to children aged up to 15 years would avert an additional 89 malaria cases per 1000 population per year (14.6% reduction) in the Cascades Region, assuming current coverage of pyrethroid-piperonyl butoxide ITNs. Malaria infections were high in all age strata, although highest in children aged 5 to 15 years, despite roll out of core malaria control interventions. Given the burden of infection in school-age children, extension of the eligibility criteria for SMC could help reduce the burden of malaria in Burkina Faso and other countries in the region.


Amodiaquine/therapeutic use , Antimalarials/therapeutic use , Malaria, Falciparum/epidemiology , Malaria, Falciparum/prevention & control , Plasmodium falciparum/immunology , Pyrimethamine/therapeutic use , Seasons , Sulfadoxine/therapeutic use , Adolescent , Adult , Antigens, Protozoan/blood , Antigens, Protozoan/immunology , Burkina Faso/epidemiology , Child , Child, Preschool , Cross-Sectional Studies , Drug Combinations , Drug Therapy, Combination/methods , Female , Humans , Insecticide-Treated Bednets , Malaria, Falciparum/blood , Malaria, Falciparum/parasitology , Male , Middle Aged , Plasmodium falciparum/isolation & purification , Prevalence , Risk Factors , Rural Population , Treatment Outcome , Young Adult
20.
Clin Transplant ; 36(2): e14523, 2022 02.
Article En | MEDLINE | ID: mdl-34724254

BACKGROUND: Intravascular ultrasound (IVUS) and optical coherence tomography (OCT) improve sensitivity of cardiac allograft vasculopathy (CAV) detection compared to invasive coronary angiography (ICA), but their ability to predict clinical events is unknown. We determined whether severe CAV detected with ICA, IVUS, or OCT correlates with graft function. METHODS: Comparison of specific vessel parameters between IVUS and OCT on 20 patients attending for angiography 12-24 months post-orthotopic heart transplant. Serial left ventricular ejection fraction (EF) was recorded prospectively. RESULTS: Analyzing 55 coronary arteries, OCT and IVUS correlated well for vessel CAV characteristics. A mean intimal thickness (MIT)OCT  > .25 mm had a sensitivity of 86.7% and specificity of 74.3% at detecting Stanford grade 4 CAV. Those with angiographically evident CAV had significant reduction in graft EF over 7.3 years follow-up (median ΔEF -2% vs +1.5%, P = .03). Patients with MITOCT  > .25 mm in at least one vessel had a lower median EF at time of surveillance (57% vs 62%, P = .014). Two MACEs were noted. CONCLUSION: Imaging with OCT correlates well with IVUS for CAV detection. Combined angiography and OCT to screen for CAV within 12-24 months of transplant predicts concurrent and future deterioration in graft function.


Coronary Artery Disease , Heart Diseases , Heart Transplantation , Allografts , Coronary Angiography/methods , Coronary Artery Disease/diagnostic imaging , Coronary Artery Disease/etiology , Heart Transplantation/adverse effects , Heart Transplantation/methods , Humans , Stroke Volume , Ultrasonography, Interventional , Ventricular Function, Left
...