We used a spatially explicit agent-based model of SARS-CoV-2 transmission combined with spatially fine-grained COVID-19 observation data from Public Health Scotland to investigate the initial rise of the Omicron (BA.1) variant of concern. We evaluated plausible scenarios for transmission rate advantage and vaccine immune escape relative to the Delta variant based on the data that would have been available at that time. We also explored possible outcomes of different levels of imposed non-pharmaceutical intervention. The initial results of these scenarios were used to inform the Scottish Government in the early outbreak stages of the Omicron variant.Using the model with parameters fit over the Delta variant epidemic, some initial assumptions about Omicron transmission rate advantage and vaccine escape, and a simple growth rate fitting procedure, we were able to capture the initial outbreak dynamics for Omicron. We found that the modelled dynamics hold up to retrospective scrutiny. The modelled imposition of extra non-pharmaceutical interventions planned by the Scottish Government at the time would likely have little effect in light of the transmission rate advantage held by the Omicron variant and the fact that the planned interventions would have occurred too late in the outbreak’s trajectory. Finally, we found that any assumptions made about the projected distribution of vaccines in the model population had little bearing on the outcome, in terms of outbreak size and timing. Instead, it was the landscape of prior immunity that was most important.
Background: Mobility restrictions prevent the spread of infections to disease-free areas, and early in the coronavirus disease 2019 (COVID-19) pandemic, most countries imposed severe restrictions on mobility as soon as it was clear that containment of local outbreaks was insufficient to control spread. These restrictions have adverse impacts on the economy and other aspects of human health, and it is important to quantify their impact for evaluating their future value. Methods: Here we develop Scotland Coronavirus transmission Model (SCoVMod), a model for COVID-19 in Scotland, which presents unusual challenges because of its diverse geography and population conditions. Our fitted model captures spatio-temporal patterns of mortality in the first phase of the epidemic to a fine geographical scale. Results: We find that lockdown restrictions reduced transmission rates down to an estimated 12\% of its pre-lockdown rate. We show that, while the timing of COVID-19 restrictions influences the role of the transmission rate on the number of COVID-related deaths, early reduction in long distance movements does not. However, poor health associated with deprivation has a considerable association with mortality; the Council Area (CA) with the greatest health-related deprivation was found to have a mortality rate 2.45 times greater than the CA with the lowest health-related deprivation considering all deaths occurring outside of carehomes. Conclusions: We find that in even an early epidemic with poor case ascertainment, a useful spatially explicit model can be fit with meaningful parameters based on the spatio-temporal distribution of death counts. Our simple approach is useful to strategically examine trade-offs between travel related restrictions and physical distancing, and the effect of deprivation-related factors on outcomes.
Background: Mobility restrictions prevent the spread of infections to disease-free areas, and early in the coronavirus disease 2019 (COVID-19) pandemic, most countries imposed severe restrictions on mobility as soon as it was clear that containment of local outbreaks was insufficient to control spread. These restrictions have adverse impacts on the economy and other aspects of human health, and it is important to quantify their impact for evaluating their future value. Methods: Here we develop Scotland Coronavirus transmission Model (SCoVMod), a model for COVID-19 in Scotland, which presents unusual challenges because of its diverse geography and population conditions. Our fitted model captures spatio-temporal patterns of mortality in the first phase of the epidemic to a fine geographical scale. Results: We find that lockdown restrictions reduced transmission rates down to an estimated 12\% of its pre-lockdown rate. We show that, while the timing of COVID-19 restrictions influences the role of the transmission rate on the number of COVID-related deaths, early reduction in long distance movements does not. However, poor health associated with deprivation has a considerable association with mortality; the Council Area (CA) with the greatest health-related deprivation was found to have a mortality rate 2.45 times greater than the CA with the lowest health-related deprivation considering all deaths occurring outside of carehomes. Conclusions: We find that in even an early epidemic with poor case ascertainment, a useful spatially explicit model can be fit with meaningful parameters based on the spatio-temporal distribution of death counts. Our simple approach is useful to strategically examine trade-offs between travel related restrictions and physical distancing, and the effect of deprivation-related factors on outcomes.
The second SARS virus, SARS-CoV-2, emerged in December 2019, and within a month was globally distributed. It was first introduced into Scotland in February 2020 associated with returning travellers and visitors. By March it was circulating in communities across the UK, and to control COVID-19 cases, and prevent overwhelming of the National Health Service (NHS), a 'lockdown' was introduced on 23rd March 2020 with a restriction of people's movements. To augment the public health efforts a large-scale genome epidemiology effort (as part of the COVID-19 Genomics UK (COG-UK) consortium) resulted in the sequencing of over 5000 SARS-CoV-2 genomes by 18th August 2020 from Scottish cases, about a quarter of the estimated number of cases at that time. Here we quantify the geographical origins of the first wave introductions into Scotland from abroad and other UK regions, the spread of these SARS-CoV-2 lineages to different regions within Scotland (defined at the level of NHS Health Board) and the effect of lockdown on virus 'success'. We estimate that approximately 300 introductions seeded lineages in Scotland, with around 25% of these lineages composed of more than five viruses, but by June circulating lineages were reduced to low levels, in line with low numbers of recorded positive cases. Lockdown was, thus, associated with a dramatic reduction in infection numbers and the extinguishing of most virus lineages. Unfortunately since the summer cases have been rising in Scotland in a second wave, with >1000 people testing positive on a daily basis, and hospitalisation of COVID-19 cases on the rise again. Examining the available Scottish genome data from the second wave, and comparing it to the first wave, we find that while some UK lineages have persisted through the summer, the majority of lineages responsible for the second wave are new introductions from outside of Scotland and many from outside of the UK. This indicates that, while lockdown in Scotland is directly linked with the first wave case numbers being brought under control, travel-associated imports (mostly from Europe or other parts of the UK) following the easing of lockdown are responsible for seeding the current epidemic population. This demonstrates that the impact of stringent public health measures can be compromised if following this, movements from regions of high to low prevalence are not minimised.
Summary Restrictions on mobility are a key component of infectious disease controls, preventing the spread of infections to as yet unexposed areas, or to regions which have previously eliminated outbreaks. However, even under the most severe restrictions, some travel must inevitably continue, at the very minimum to retain essential services. For COVID-19, most countries imposed severe restrictions on travel at least as soon as it was clear that containment of local outbreaks would not be possible. Such restrictions are known to have had a substantial impact on the economy and other aspects of human health, and so quantifying the impact of such restrictions is an essential part of evaluating the necessity for future implementation of similar measures. In this analysis, we built a simulation model using National statistical data to record patterns of movements to work, and implement levels of mobility recorded in real time via mobile phone apps. This model was fitted to the pattern of deaths due to COVID-19 using approximate Bayesian inference. Our model is able to recapitulate mortality considering the number of deaths and datazones (DZs, which are areas containing approximately 500-1000 residents) with deaths, as measured across 32 individual council areas (CAs) in Scotland. Our model recreates a trajectory consistent with the observed data until 1 st of July. According to the model, most transmission was occurring “locally” (i.e. in the model, 80% of transmission events occurred within spatially defined “communities” of approximately 100 individuals). We show that the net effect of the various restrictions put into place in March can be captured by a reduction in transmission down to 12% of its pre-lockdown rate effective 28 th March. By comparing different approaches to reducing transmission, we show that, while the timing of COVID-19 restrictions influences the role of the transmission rate on the number of COVID-related deaths, early reduction in long distance movements does not reduce death rates significantly. As this movement of individuals from more infected areas to less infected areas has a minimal impact on transmission, this suggests that the fraction of population already immune in infected communities was not a significant factor in these early stages of the national epidemic even when local clustering of infection is taken into account. The best fit model also shows a considerable influence of the health index of deprivation (part of the “index of multiple deprivations”) on mortality. The most likely value has the CA with the highest level of health-related deprivation to have on average, a 2.45 times greater mortality rate due to COVID-19 compared to the CA with the lowest, showing the impact of health-related deprivation even in the early stages of the COVID-19 national epidemic.
Livestock movements are an important mechanism of infectious disease transmission. Where these are well recorded, network analysis tools have been used to successfully identify system properties, highlight vulnerabilities to transmission, and inform targeted surveillance and control. Here we highlight the main uses of network properties in understanding livestock disease epidemiology and discuss statistical approaches to infer network characteristics from biased or fragmented datasets. We use a ‘hurdle model’ approach that predicts (i) the probability of movement and (ii) the number of livestock moved to generate synthetic ‘complete’ networks of movements between administrative wards, exploiting routinely collected government movement permit data from northern Tanzania. We demonstrate that this model captures a significant amount of the observed variation. Combining the cattle movement network with a spatial between-ward contact layer, we create a multiplex, over which we simulated the spread of ‘fast’ (R0 = 3) and ‘slow’ (R0 = 1.5) pathogens, and assess the effects of random versus targeted disease control interventions (vaccination and movement ban). The targeted interventions substantially outperform those randomly implemented for both fast and slow pathogens. Our findings provide motivation to encourage routine collection and centralization of movement data to construct representative networks. This article is part of the theme issue ‘Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control’. This theme issue is linked with the earlier issue ‘Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes’.
G.L. Chaters1,a, P.C.D. Johnson1,a, S. Cleaveland1, J. Crispell2, W.A. de Glanville1, T. Doherty3, L. Matthews1, S. Mohr1, O.M. Nyasebwa4, G. Rossi3, L.C.M. Salvador3, 4, 5, E. Swai6, R.R.Kao3,b 1 Boyd Orr Centre for Population and Ecosystem Health, Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow G12 8QQ, UK 2 School of Veterinary Medicine, University College Dublin, Ireland 3 Royal (Dick) School of Veterinary Studies and Roslin Institute, University of Edinburgh, Easter Bush Campus, Midlothian, Scotland, UK 4 Department of Infectious Diseases, University of Georgia, Athens, Georgia, USA 5 Institute of Bioinformatics, University of Georgia, Athens, USA 6 Department of Veterinary Services, Ministry of Livestock and Fisheries, Tanzania; a These authors contributed equally to this work b Corresponding author
The movements of livestock between premises and markets can be characterised as a dynamic network where the structure of the network itself can critically impact the transmission dynamics of many infectious diseases. As evidenced by the 2001 foot-and-mouth disease (FMD) epidemic in the UK, this can involve transmission over large geographical distances and can result in major economic loss. One consequence of the FMD epidemic was the introduction of mandatory livestock movement restrictions: a 13-day standstill in Scotland for cattle and sheep after moving livestock onto a farm (allowing many exemptions) and a 6-day standstill for cattle and sheep in England and Wales (with minor exemptions, e.g. direct movements to slaughter). Such standstills are known to be effective but commercial considerations result in pressures to relax them. When contemplating legislative changes such as a change in length of movement restrictions we need to consider the consequent effect these could have on the emergent properties of the system, i.e. the network structure itself. In this study, we investigate how disease dynamics change when the local contact structure of the recorded livestock movement network in Scotland is altered through rewiring movements between premises. The network rewiring used here changes the structure of the recorded trade network through a combination of altered movement restrictions and redirection of movements between holdings and markets to avoid nonsensical activity (e.g. movements to markets on days when they are inactive) while conserving other characteristics (e.g. movement date as closely as possible and market sales of the correct animal production type). Rewiring results in networks with higher clustering coefficients and lower network density. The impact of rewiring on a hypothetical foot-and-mouth disease outbreak in Scotland was assessed by stochastic simulation, considering scenarios with and without exemptions to the standstill rules. As expected, rewiring leads to a decrease in outbreak size and - if standstill exemptions are prohibited - higher probability of smaller outbreaks. Without exemptions, a shorter movement standstill is almost as effective as a longer standstill period, indicating that a simpler biosecurity system would offer minimal additional risk for FMD. These results suggest that explicitly manipulating the contact network structure in a sensible way has the potential to significantly impact disease control.
In recent work, a research group led by Prof. Kao identified risk-based surveillance strategies for bovine Tuberculosis (bTB) in Scotland, using a combination of statistical models and a simulation model for freedom from infection, to identify a series of possible strategies to reduce the number of regular herd tests for bTB in Scotland. In a subsequent work on risk-based surveillance (also led by Prof. Kao) (Defra project SE3285), it was shown that that there are opportunities to improve the current surveillance system in areas of low risk of bTB (LRAs). The current proposal will build on previous work to determine if alternative surveillance strategies can improve the current surveillance system by looking at the most recent epidemiological situation in these areas. Different strategies will be considered, contingent on discussion with Defra and APHA colleagues. Broadly speaking reduction strategies can follow several strands: i) Reduction in testing - direct reduction of the amount of testing done, as a strategy aimed at directly reducing costs. Ideally, this will be a 'dominant' strategy, whereby testing is reduced but with minimal risk of an increase in incidence of breakdowns. This strategy would be most similar to the recommendations for Scotland. ii) Onward risk reduction - this aims partly to identify herd breakdowns more quickly by targeting herds that are at a higher risk of infection. iii) Reduction in testing and onward risk – directly reduction of the amount of testing done while identifying herd breakdowns more quickly. We shall adopt a three-step approach to this project, using methodologies and analyses based on existing work, either published or currently undergoing scientific peer review. A). Develop statistical models to identify risk factors for breakdowns within LRAs. B). Define risk-based surveillance strategies based on combinations of the criteria (i) to (iii). C). Based on the best, analytically sound recommendations using A and B we shall consult with Defra and APHA colleagues to determine their feasibility. We shall therefore provide an analysis of the current risk situation in LRAs where the two aims of testing reduction and reduced onward risk will be balanced against the practicality of implementation and relevance to policy requirements.
The centrality dependence of the mean charged-particle multiplicity as a function of pseudorapidity is measured in approximately 1 mu b(-1) of proton-lead collisions at a nucleon-nucleon centre-of-mass energy of root s(NN) = 5.02 TeV using the ATLAS detector at the Large Hadron Collider. Charged particles with absolute pseudorapidity less than 2.7 are reconstructed using the ATLAS pixel detector. The p + Pb collision centrality is characterised by the total transverse energy measured in the Pb-going direction of the forward calorimeter. The charged-particle pseudorapidity distributions are found to vary strongly with centrality, with an increasing asymmetry between the proton-going and Pb-going directions as the collisions become more central. Three different estimations of the number of nucleons participating in the p + Pb collision have been carried out using the Glauber model as well as two Glauber-Gribov inspired extensions to the Glauber model. Charged-particle multiplicities per participant pair are found to vary differently for these three models, highlighting the importance of including colour fluctuations in nucleon-nucleon collisions in the modelling of the initial state of p + Pb collisions.
Background Modelling disease outbreaks often involves integrating the wealth of data that are gathered during modern outbreaks into complex mathematical or computational models of transmission. Incorporating these data into simple compartmental epidemiological models is often challenging, requiring the use of more complex but also more efficient computational models. In this paper we introduce a new framework that allows for a more systematic and user-friendly way of building and running epidemiological models that efficiently handles disease data and reduces much of the boilerplate code that usually associated to these models. We introduce the framework by developing an SIR model on a simple network as an example. Results We develop Broadwick, a modular, object-oriented epidemiological framework that efficiently handles large epidemiological datasets and provides packages for stochastic simulations, parameter inference using Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC) methods. Each algorithm used is fully customisable with sensible defaults that are easily overridden by custom algorithms as required. Conclusion Broadwick is an epidemiological modelling framework developed to increase the productivity of researchers by providing a common framework with which to develop and share complex models. It will appeal to research team leaders as it allows for models to be created prior to a disease outbreak and has the ability to handle large datasets commonly found in epidemiological modelling.
Using a data sample of 106 million $\psi(3686)$ events collected with the BESIII detector operated at the BEPCII storage ring, we study for the first time the decay $\chi_{cJ}\to\phi K^{0}_S K^{\pm}\pi^{\mp}$ and $\chi_{cJ}\to\phi K^{+} K^{-}\pi^{0}$ in the E1 radiative transition $\psi(3686)\to\gamma\chi_{cJ}$. The decays are dominated by the three-body decay $\chi_{cJ}\to \phi K^*(892)\bar{K}$. We measure branching fractions for this reaction via the neutral and charged $K^*(892)$ and find them consistent with each other within the expectation of isospin symmetry. In the $K\bar{K}\pi$ invariant mass distribution a structure near the $K^*(892)\bar{K}$ mass threshold is observed, and the corresponding mass and width are measured to be $1412\pm4(\mathrm{stat.})\pm8(\mathrm{sys.}) \mathrm{MeV}/c^2$ and $\Gamma$ = $84\pm12(\mathrm{stat.})\pm40(\mathrm{sys.}) \mathrm{MeV}$, respectively. The observed state favors an assignment to the $h_1(1380)$, considering its possible $J^{PC}$ and comparing its mass, width and decay mode to those reported in the Particle Data Group.
A search for the decay to a pair of new particles of either the 125 GeV Higgs boson (h) or a second charge parity (CP)-even Higgs boson (H) is presented. The data set corresponds to an integrated luminosity of 20.3 fb−1 of pp collisions at √s=8 TeV recorded by the ATLAS experiment at the LHC in 2012. The search was done in the context of the next-to-minimal supersymmetric standard model, in which the new particles are the lightest neutral pseudoscalar Higgs bosons (a). One of the two a bosons is required to decay to two muons while the other is required to decay to two τ leptons. No significant excess is observed above the expected backgrounds in the dimuon invariant mass range from 3.7 to 50 GeV. Upper limits are placed on the production of h→aa relative to the standard model gg→h production, assuming no coupling of the a boson to quarks. The most stringent limit is placed at 3.5% for ma=3.75 GeV. Upper limits are also placed on the production cross section of H→aa from 2.33 to 0.72 pb, for fixed ma=5 GeV with mH ranging from 100 to 500 GeV.Received 8 May 2015DOI:https://doi.org/10.1103/PhysRevD.92.052002This article is available under the terms of the Creative Commons Attribution 3.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.© 2015 CERN, for the ATLAS Collaboration
The jet energy scale (JES) and its systematic uncertainty are determined for jets measured with the ATLAS detector using proton-proton collision data with a centre-of-mass energy of [Formula: see text] TeV corresponding to an integrated luminosity of [Formula: see text][Formula: see text]. Jets are reconstructed from energy deposits forming topological clusters of calorimeter cells using the anti-[Formula: see text] algorithm with distance parameters [Formula: see text] or [Formula: see text], and are calibrated using MC simulations. A residual JES correction is applied to account for differences between data and MC simulations. This correction and its systematic uncertainty are estimated using a combination of in situ techniques exploiting the transverse momentum balance between a jet and a reference object such as a photon or a [Formula: see text] boson, for [Formula: see text] and pseudorapidities [Formula: see text]. The effect of multiple proton-proton interactions is corrected for, and an uncertainty is evaluated using in situ techniques. The smallest JES uncertainty of less than 1 % is found in the central calorimeter region ([Formula: see text]) for jets with [Formula: see text]. For central jets at lower [Formula: see text], the uncertainty is about 3 %. A consistent JES estimate is found using measurements of the calorimeter response of single hadrons in proton-proton collisions and test-beam data, which also provide the estimate for [Formula: see text] TeV. The calibration of forward jets is derived from dijet [Formula: see text] balance measurements. The resulting uncertainty reaches its largest value of 6 % for low-[Formula: see text] jets at [Formula: see text]. Additional JES uncertainties due to specific event topologies, such as close-by jets or selections of event samples with an enhanced content of jets originating from light quarks or gluons, are also discussed. The magnitude of these uncertainties depends on the event sample used in a given physics analysis, but typically amounts to 0.5-3 %.
Citation Aad, G., B. Abbott, J. Abdallah, S. Abdel Khalek, O. Abdinov, R. Aben, B. Abi, et al. 2015. “Erratum to: Search for production of WW / WZ resonances decaying to a lepton, neutrino and jets in pp collisions at documentclass[12pt]{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$sqrt{s}=8$$end{document}s=8 TeV with the ATLAS detector.” The European Physical Journal. C, Particles and Fields 75 (8): 370. doi:10.1140/epjc/s10052-0153593-4. http://dx.doi.org/10.1140/epjc/s10052-015-3593-4.
The online version of the original article can be found under doi:10.1140/epjc/s10052-015-3425-6. The online version of the original article can be found at http://dx.doi.org/10.1140/epjc/s10052-015-3425-6.
The results of a dedicated search for pair production of scalar partners of charm quarks are reported. The search is based on an integrated luminosity of 20.3 fb^{-1} of pp collisions at sqrt[s]=8 TeV recorded with the ATLAS detector at the LHC. The search is performed using events with large missing transverse momentum and at least two jets, where the two leading jets are each tagged as originating from c quarks. Events containing isolated electrons or muons are vetoed. In an R-parity-conserving minimal supersymmetric scenario in which a single scalar-charm state is kinematically accessible, and where it decays exclusively into a charm quark and a neutralino, 95% confidence-level upper limits are obtained in the scalar-charm-neutralino mass plane such that, for neutralino masses below 200 GeV, scalar-charm masses up to 490 GeV are excluded.