Calculation of disease dynamics in a population of households

Early mathematical representations of infectious disease dynamics assumed a single, large, homogeneously mixing population. Over the past decade there has been growing interest in models consisting of multiple smaller subpopulations (households, workplaces, schools, communities), with the natural assumption of strong homogeneous mixing within each subpopulation, and weaker transmission between subpopulations. Here we consider a model of SIRS (susceptible-infectious-recovered-suscep​tible) infection dynamics in a very large (assumed infinite) population of households, with the simplifying assumption that each household is of the same size (although all methods may be extended to a population with a heterogeneous distribution of household sizes). For this households model we present efficient methods for studying several quantities of epidemiological interest: (i) the threshold for invasion; (ii) the early growth rate; (iii) the household offspring distribution; (iv) the endemic prevalence of infection; and (v) the transient dynamics of the process. We utilize these methods to explore a wide region of parameter space appropriate for human infectious diseases. We then extend these results to consider the effects of more realistic gamma-distributed infectious periods. We discuss how all these results differ from standard homogeneous-mixing models and assess the implications for the invasion, transmission and persistence of infection. The computational efficiency of the methodology presented here will hopefully aid in the parameterisation of structured models and in the evaluation of appropriate responses for future disease outbreaks

Benchmarking multi-rate codon models

CITATION: Delport, W. et al. 2010. Benchmarking multi-rate codon models. PLoS ONE, 5(7): e11587, doi:10.1371/journal.pone.0011587.The original publication is available at http://journals.plos.org/plosoneThe single rate codon model of non-synonymous substitution is ubiquitous in phylogenetic modeling. Indeed, the use of a non-synonymous to synonymous substitution rate ratio parameter has facilitated the interpretation of selection pressure on genomes. Although the single rate model has achieved wide acceptance, we argue that the assumption of a single rate of non-synonymous substitution is biologically unreasonable, given observed differences in substitution rates evident from empirical amino acid models. Some have attempted to incorporate amino acid substitution biases into models of codon evolution and have shown improved model performance versus the single rate model. Here, we show that the single rate model of non-synonymous substitution is easily outperformed by a model with multiple non-synonymous rate classes, yet in which amino acid substitution pairs are assigned randomly to these classes. We argue that, since the single rate model is so easy to improve upon, new codon models should not be validated entirely on the basis of improved model fit over this model. Rather, we should strive to both improve on the single rate model and to approximate the general time-reversible model of codon substitution, with as few parameters as possible, so as to reduce model over-fitting. We hint at how this can be achieved with a Genetic Algorithm approach in which rate classes are assigned on the basis of sequence information content. © 2010 Delport et al.http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0011587Publisher's versio

A Cellular Automata-Based Mathematical Model for Thymocyte Development

Intrathymic T cell development is an important process necessary for the normal formation of cell-mediated immune responses. Importantly, such a process depends on interactions of developing thymocytes with cellular and extracellular elements of the thymic microenvironment. Additionally, it includes a series of oriented and tunely regulated migration events, ultimately allowing mature cells to cross endothelial barriers and leave the organ. Herein we built a cellular automata-based mathematical model for thymocyte migration and development. The rules comprised in this model take into account the main stages of thymocyte development, two-dimensional sections of the normal thymic microenvironmental network, as well as the chemokines involved in intrathymic cell migration. Parameters of our computer simulations with further adjusted to results derived from previous experimental data using sub-lethally irradiated mice, in which thymus recovery can be evaluated. The model fitted with the increasing numbers of each CD4/CD8-defined thymocyte subset. It was further validated since it fitted with the times of permanence experimentally ascertained in each CD4/CD8-defined differentiation stage. Importantly, correlations using the whole mean volume of young normal adult mice revealed that the numbers of cells generated in silico with the mathematical model fall within the range of total thymocyte numbers seen in these animals. Furthermore, simulations made with a human thymic epithelial network using the same mathematical model generated similar profiles for temporal evolution of thymocyte developmental stages. Lastly, we provided in silico evidence that the thymus architecture is important in the thymocyte development, since changes in the epithelial network result in different theoretical profiles for T cell development/migration. This model likely can be used to predict thymocyte evolution following therapeutic strategies designed for recovery of the thymus in diseases coursing with thymus involution, such as some primary immunodeficiencies, acute infections, and malnutrition

COVID-19 Infection Risk amongst 14,104 Vaccinated Care Home Residents: A national observational longitudinal cohort study in Wales, United Kingdom, December 2020 to March 2021

Backgroundvaccinations for COVID-19 have been prioritised for older people living in care homes. However, vaccination trials included limited numbers of older people.Aimwe aimed to study infection rates of SARS-CoV-2 for older care home residents following vaccination and identify factors associated with increased risk of infection.Study Design and Settingwe conducted an observational data-linkage study including 14,104 vaccinated older care home residents in Wales (UK) using anonymised electronic health records and administrative data.Methodswe used Cox proportional hazards models to estimate hazard ratios (HRs) for the risk of testing positive for SARS-CoV-2 infection following vaccination, after landmark times of either 7 or 21 days post-vaccination. We adjusted HRs for age, sex, frailty, prior SARS-CoV-2 infections and vaccination type.Resultswe observed a small proportion of care home residents with positive polymerase chain reaction (tests following vaccination 1.05% (N = 148), with 90% of infections occurring within 28 days. For the 7-day landmark analysis we found a reduced risk of SARS-CoV-2 infection for vaccinated individuals who had a previous infection; HR (95% confidence interval) 0.54 (0.30, 0.95). For the 21-day landmark analysis, we observed high HRs for individuals with low and intermediate frailty compared with those without; 4.59 (1.23, 17.12) and 4.85 (1.68, 14.04), respectively.Conclusionsincreased risk of infection after 21 days was associated with frailty. We found most infections occurred within 28 days of vaccination, suggesting extra precautions to reduce transmission risk should be taken in this time frame

Disease control across urban–rural gradients

Controlling the regional re-emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) after its initial spread in ever-changing personal contact networks and disease landscapes is a challenging task. In a landscape context, contact opportunities within and between populations are changing rapidly as lockdown measures are relaxed and a number of social activities re-activated. Using an individual-based metapopulation model, we explored the efficacy of different control strategies across an urban–rural gradient in Wales, UK. Our model shows that isolation of symptomatic cases or regional lockdowns in response to local outbreaks have limited efficacy unless the overall transmission rate is kept persistently low. Additional isolation of non-symptomatic infected individuals, who may be detected by effective test-and-trace strategies, is pivotal to reducing the overall epidemic size over a wider range of transmission scenarios. We define an ‘urban–rural gradient in epidemic size' as a correlation between regional epidemic size and connectivity within the region, with more highly connected urban populations experiencing relatively larger outbreaks. For interventions focused on regional lockdowns, the strength of such gradients in epidemic size increased with higher travel frequencies, indicating a reduced efficacy of the control measure in the urban regions under these conditions. When both non-symptomatic and symptomatic individuals are isolated or regional lockdown strategies are enforced, we further found the strongest urban–rural epidemic gradients at high transmission rates. This effect was reversed for strategies targeted at symptomatic individuals only. Our results emphasize the importance of test-and-trace strategies and maintaining low transmission rates for efficiently controlling SARS-CoV-2 spread, both at landscape scale and in urban areas

Graduate Entry Medicine: Selection Criteria and Student Performance

Background: Graduate entry medicine raises new questions about the suitability of students with different backgrounds. We examine this, and the broader issue of effectiveness of selection and assessment procedures. Methods: The data included background characteristics, academic record, interview score and performance in pre-clinical modular assessment for two years intake of graduate entry medical students. Exploratory factor analysis is a powerful method for reducing a large number of measures to a smaller group of underlying factors. It was used here to identify patterns within and between the selection and performance data. Principal Findings: Basic background characteristics were of little importance in predicting exam success. However, easily interpreted components were detected within variables comprising the ‘selection ’ and ‘assessment ’ criteria. Three selection components were identified (‘Academic’, ‘GAMSAT’, ‘Interview’) and four assessment components (‘General Exam’, ‘Oncology’, ‘OSCE’, ‘Family Case Study’). There was a striking lack of relationships between most selection and performance factors. Only ‘General Exam ’ and ‘Academic ’ showed a correlation (Pearson’s r = 0.55, p,0.001). Conclusions: This study raises questions about methods of student selection and their effectiveness in predicting performance and assessing suitability for a medical career. Admissions tests and most exams only confirmed previous academic achievement, while interview scores were not correlated with any consequent assessment

Classification of accelerometer wear and non-wear events in seconds for monitoring free-living physical activity

_____________________________________________________________ This article is brought to you by Swansea University. Any person downloading material is agreeing to abide by the terms of the repository licence. Authors are personally responsible for adhering to publisher restrictions or conditions. When uploading content they are required to comply with their publisher agreement and the SHERPA RoMEO database to judge whether or not it is copyright safe to add this version of the paper to this repository. Design: A bi-moving-window-based approach was used to combine acceleration and skin temperature data to identify wear and non-wear time events in triaxial accelerometer data that monitor physical activity. Setting: Local residents in Swansea, Wales, UK. Participants: 50 participants aged under 16 years (n=23) and over 17 years (n=27) were recruited in two phases: phase 1: design of the wear/non-wear algorithm (n=20) and phase 2: validation of the algorithm (n=30). Methods: Participants wore a triaxial accelerometer (GeneActiv) against the skin surface on the wrist (adults) or ankle (children). Participants kept a diary to record the timings of wear and non-wear and were asked to ensure that events of wear/non-wear last for a minimum of 15 min. Results: The overall sensitivity of the proposed method was 0.94 (95% CI 0.90 to 0.98) and specificity 0.91 (95% CI 0.88 to 0.94). It performed equally well for children compared with adults, and females compared with males. Using surface skin temperature data in combination with acceleration data significantly improved the classification of wear/non-wear time when compared with methods that used acceleration data only ( p<0.01). Conclusions: Using either accelerometer seismic information or temperature information alone is prone to considerable error. Combining both sources of data can give accurate estimates of non-wear periods thus giving better classification of sedentary behaviour. This method can be used in population studies of physical activity in free-living environments

Impact of COVID-19 pandemic on community medication dispensing: a national cohort analysis in Wales, UK

BackgroundPopulation-level information on dispensed medication provides insight on the distribution of treated morbidities, particularly if linked to other population-scale data at an individual-level.ObjectiveTo evaluate the impact of COVID-19 on dispensing patterns of medications.MethodsRetrospective observational study using population-scale, individual-level dispensing records in Wales, UK. Total dispensed drug items for the population between 1st January 2016 and 31st December 2019 (3-years, pre-COVID-19) were compared to 2020 with follow up until 27th July 2021 (COVID-19 period). We compared trends across all years and British National Formulary (BNF) chapters and highlighted the trends in three major chapters for 2019-21: 1-Cardiovascular system (CVD); 2-Central Nervous System (CNS); 3-Immunological & Vaccine. We developed an interactive dashboard to enable monitoring of changes as the pandemic evolves.ResultAmongst all BNF chapters, 73,410,543 items were dispensed in 2020 compared to 74,121,180 items in 2019 demonstrating -0.96% relative decrease in 2020. Comparison of monthly patterns showed average difference (D) of -59,220 and average Relative Change (RC) of -0.74% between the number of dispensed items in 2020 and 2019. Maximum RC was observed in March 2020 (D= +1,224,909 and RC= +20.62%), followed by second peak in June 2020 (D= +257,920, RC= +4.50%). A third peak was observed in September 2020 (D= +264,138, RC= +4.35%). Large increases in March 2020 were observed for CVD and CNS medications across all age groups. The Immunological and Vaccine products dropped to very low levels across all age groups and all months (including the March dispensing peak).ConclusionsReconfiguration of routine clinical services during COVID-19 led to substantial changes in community pharmacy drug dispensing. This change may contribute to a long-term burden of COVID-19, raising the importance of a comprehensive and timely monitoring of changes for evaluation of the potential impact on clinical care and outcomes

Long term extension of a randomised controlled trial of probiotics using electronic health records

Most randomised controlled trials (RCTs) are relatively short term and, due to costs and available resources, have limited opportunity to be re-visited or extended. There is no guarantee that effects of treatments remain unchanged beyond the study. Here, we illustrate the feasibility, benefits and cost-effectiveness of enriching standard trial design with electronic follow up. We completed a 5-year electronic follow up of a RCT investigating the impact of probiotics on asthma and eczema in children born 2005-2007, with traditional fieldwork follow up to two years. Participants and trial outcomes were identified and analysed after five years using secure, routine, anonymised, person-based electronic health service databanks. At two years, we identified 93% of participants and compared fieldwork with electronic health records, highlighting areas of agreement and disagreement. Retention of children from lower socio-economic groups was improved, reducing volunteer bias. At 5 years we identified a reduced 82% of participants. These data allowed the trial's first robust analysis of asthma endpoints. We found no indication that probiotic supplementation to pregnant mothers and infants protected against asthma or eczema at 5 years. Continued longer-term follow up is technically straightforward

Calibration Belt for Quality-of-Care Assessment Based on Dichotomous Outcomes

Prognostic models applied in medicine must be validated on independent samples, before their use can be recommended. The assessment of calibration, i.e., the model's ability to provide reliable predictions, is crucial in external validation studies. Besides having several shortcomings, statistical techniques such as the computation of the standardized mortality ratio (SMR) and its confidence intervals, the Hosmer–Lemeshow statistics, and the Cox calibration test, are all non-informative with respect to calibration across risk classes. Accordingly, calibration plots reporting expected versus observed outcomes across risk subsets have been used for many years. Erroneously, the points in the plot (frequently representing deciles of risk) have been connected with lines, generating false calibration curves. Here we propose a methodology to create a confidence band for the calibration curve based on a function that relates expected to observed probabilities across classes of risk. The calibration belt allows the ranges of risk to be spotted where there is a significant deviation from the ideal calibration, and the direction of the deviation to be indicated. This method thus offers a more analytical view in the assessment of quality of care, compared to other approaches