The spatial Principal Component Analysis (sPCA, Jombart (Heredity 101:92-103, 2008) is designed to investigate non-random spatial distributions of genetic variation. Unfortunately, the associated ...tests used for assessing the existence of spatial patterns (global and local test; (Heredity 101:92-103, 2008) lack statistical power and may fail to reveal existing spatial patterns. Here, we present a non-parametric test for the significance of specific patterns recovered by sPCA.
We compared the performance of this new test to the original global and local tests using datasets simulated under classical population genetic models. Results show that our test outperforms the original global and local tests, exhibiting improved statistical power while retaining similar, and reliable type I errors. Moreover, by allowing to test various sets of axes, it can be used to guide the selection of retained sPCA components.
As such, our test represents a valuable complement to the original analysis, and should prove useful for the investigation of spatial genetic patterns.
•Real-time estimation of reproduction numbers during outbreaks can guide control.•Using up-to-date serial interval data and accounting for imported cases is vital.•We develop a framework for ...estimating pathogen transmissibility appropriately.•We demonstrate it using data from outbreaks of influenza, Ebola and MERS.•Our approach is implemented in R package EpiEstim and online application EpiEstim App.
Accurate estimation of the parameters characterising infectious disease transmission is vital for optimising control interventions during epidemics. A valuable metric for assessing the current threat posed by an outbreak is the time-dependent reproduction number, i.e. the expected number of secondary cases caused by each infected individual. This quantity can be estimated using data on the numbers of observed new cases at successive times during an epidemic and the distribution of the serial interval (the time between symptomatic cases in a transmission chain). Some methods for estimating the reproduction number rely on pre-existing estimates of the serial interval distribution and assume that the entire outbreak is driven by local transmission. Here we show that accurate inference of current transmissibility, and the uncertainty associated with this estimate, requires: (i) up-to-date observations of the serial interval to be included, and; (ii) cases arising from local transmission to be distinguished from those imported from elsewhere. We demonstrate how pathogen transmissibility can be inferred appropriately using datasets from outbreaks of H1N1 influenza, Ebola virus disease and Middle-East Respiratory Syndrome. We present a tool for estimating the reproduction number in real-time during infectious disease outbreaks accurately, which is available as an R software package (EpiEstim 2.2). It is also accessible as an interactive, user-friendly online interface (EpiEstim App), permitting its use by non-specialists. Our tool is easy to apply for assessing the transmission potential, and hence informing control, during future outbreaks of a wide range of invading pathogens.
Increasing attention is being devoted to taking landscape information into account in genetic studies. Among landscape variables, space is often considered as one of the most important. To reveal ...spatial patterns, a statistical method should be spatially explicit, that is, it should directly take spatial information into account as a component of the adjusted model or of the optimized criterion. In this paper we propose a new spatially explicit multivariate method, spatial principal component analysis (sPCA), to investigate the spatial pattern of genetic variability using allelic frequency data of individuals or populations. This analysis does not require data to meet Hardy-Weinberg expectations or linkage equilibrium to exist between loci. The sPCA yields scores summarizing both the genetic variability and the spatial structure among individuals (or populations). Global structures (patches, clines and intermediates) are disentangled from local ones (strong genetic differences between neighbors) and from random noise. Two statistical tests are proposed to detect the existence of both types of patterns. As an illustration, the results of principal component analysis (PCA) and sPCA are compared using simulated datasets and real georeferenced microsatellite data of Scandinavian brown bear individuals (Ursus arctos). sPCA performed better than PCA to reveal spatial genetic patterns. The proposed methodology is implemented in the adegenet package of the free software R.
Multivariate analyses such as principal component analysis were among the first statistical methods employed to extract information from genetic markers. From their early applications to current ...innovations, these approaches have proven to be efficient for the analysis of the genetic variability in various contexts such as human genetics, conservation and adaptation studies. However, because multivariate analysis is a wide and diversified area of statistics, choosing a method appropriate to both the data and to the question being asked can be difficult. Moreover, some particularities of genetic markers need to be taken into account when using multivariate methods. As a consequence, multivariate analyses are often used as black boxes, which results in frequent mistakes in the literature. In this review, we provide a critical analysis of the application of multivariate methods to genetic markers, using a general framework that unifies all these methods for the sake of clarity. First, we focus on some common mistakes in these applications and ways to avoid these pitfalls. We then detail the most critical particularities of allele frequencies that demand adaptations of multivariate methods, and we propose solutions to the subsequent problems. Finally, we tackle several questions of interest in which multivariate analysis has a great role to play, such as the study of the typological coherence of different genetic markers, or the investigation of spatial genetic patterns.
Species spatial distributions are the result of population demography, behavioral traits, and species interactions in spatially heterogeneous environmental conditions. Hence the composition of ...species assemblages is an integrative response variable, and its variability can be explained by the complex interplay among several structuring factors. The thorough analysis of spatial variation in species assemblages may help infer processes shaping ecological communities. We suggest that ecological studies would benefit from the combined use of the classical statistical models of community composition data, such as constrained or unconstrained multivariate analyses of site-by-species abundance tables, with rapidly emerging and diversifying methods of spatial pattern analysis. Doing so allows one to deal with spatially explicit ecological models of beta diversity in a biogeographic context through the multiscale analysis of spatial patterns in original species data tables, including spatial characterization of fitted or residual variation from environmental models. We summarize here the recent progress for specifying spatial features through spatial weighting matrices and spatial eigenfunctions in order to define spatially constrained or scale-explicit multivariate analyses. Through a worked example on tropical tree communities, we also show the potential of the overall approach to identify significant residual spatial patterns that could arise from the omission of important unmeasured explanatory variables or processes.
Abstract Important policy questions during infections disease outbreaks include: i) How effective are particular interventions?; ii) When can resource-intensive interventions be removed? We used ...mathematical modelling to address these questions during the 2017 Ebola outbreak in Likati Health Zone, Democratic Republic of the Congo (DRC). Eight cases occurred before 15 May 2017, when the Ebola Response Team (ERT; co-ordinated by the World Health Organisation and DRC Ministry of Health) was deployed to reduce transmission. We used a branching process model to estimate that, pre-ERT arrival, the reproduction number was $$R=1.49$$ R = 1.49 (95% credible interval $$({{{\mathrm{0.67,2.81}}})$$ ( 0.67, 2.81 ) ). The risk of further cases occurring without the ERT was estimated to be 0.97 (97%). However, no cases materialised, suggesting that the ERT’s measures were effective. We also estimated the risk of withdrawing the ERT in real-time. By the actual ERT withdrawal date (2 July 2017), the risk of future cases without the ERT was only 0.01, indicating that the ERT withdrawal decision was safe. We evaluated the sensitivity of our results to the estimated $$R$$ R value and considered different criteria for determining the ERT withdrawal date. This research provides an extensible modelling framework that can be used to guide decisions about when to relax interventions during future outbreaks.
Epidemiology and public health planning will increasingly rely on the analysis of genetic sequence data. In particular, genetic data coupled with dates and locations of sampled isolates can be used ...to reconstruct the spatiotemporal dynamics of pathogens during outbreaks. Thus far, phylogenetic methods have been used to tackle this issue. Although these approaches have proved useful for informing on the spread of pathogens, they do not aim at directly reconstructing the underlying transmission tree. Instead, phylogenetic models infer most recent common ancestors between pairs of isolates, which can be inadequate for densely sampled recent outbreaks, where the sample includes ancestral and descendent isolates. In this paper, we introduce a novel method based on a graph approach to reconstruct transmission trees directly from genetic data. Using simulated data, we show that our approach can efficiently reconstruct genealogies of isolates in situations where classical phylogenetic approaches fail to do so. We then illustrate our method by analyzing data from the early stages of the swine-origin A/H1N1 influenza pandemic. Using 433 isolates sequenced at both the hemagglutinin and neuraminidase genes, we reconstruct the likely history of the worldwide spread of this new influenza strain. The presented methodology opens new perspectives for the analysis of genetic data in the context of disease outbreaks.
In Mekong Delta farms (Vietnam), antimicrobials are extensively used, but limited data are available on levels of antimicrobial resistance (AMR) among Escherichia coli isolates. We performed a ...structured survey of AMR in E. coli isolates (n = 434) from 90 pig, chicken, and duck farms. The results were compared with AMR among E. coli isolates (n = 234) from 66 small wild animals (rats and shrews) trapped on farms and in forests and rice fields. The isolates were susceptibility tested against eight antimicrobials. E. coli isolates from farmed animals were resistant to a median of 4 (interquartile range IQR, 3 to 6) antimicrobials versus 1 (IQR, 1 to 2) among wild mammal isolates (P < 0.001). The prevalences of AMR among farmed species isolates (versus wild animals) were as follows: tetracycline, 84.7% (versus 25.6%); ampicillin, 78.9% (versus 85.9%); trimethoprim-sulfamethoxazole, 52.1% (versus 18.8%); chloramphenicol, 39.9% (versus 22.5%); amoxicillin-clavulanic acid, 36.6% (versus 34.5%); and ciprofloxacin, 24.9% (versus 7.3%). The prevalence of multidrug resistance (MDR) (resistance against three or more antimicrobial classes) among pig isolates was 86.7% compared to 66.9 to 72.7% among poultry isolates. After adjusting for host species, MDR was ∼8 times greater among isolates from wild mammals trapped on farms than among those trapped in forests/rice fields (P < 0.001). Isolates were assigned to unique profiles representing their combinations of susceptibility results. Multivariable analysis of variance indicated that AMR profiles from wild mammals trapped on farms and those from domestic animals were more alike (R(2) range, 0.14 to 0.30) than E. coli isolates from domestic animals and mammals trapped in the wild (R(2) range, 0.25 to 0.45). The results strongly suggest that AMR on farms is a key driver of environmental AMR in the Mekong Delta.
Since April 2015, whole genome sequencing (WGS) has been the routine test for Salmonella identification, surveillance and outbreak investigation at the national reference laboratory in England and ...Wales. In May 2015, an outbreak of Salmonella Enteritidis cases was detected using WGS data and investigated. UK cases were interviewed to obtain a food history and links between suppliers were mapped to produce a food chain network for chicken eggs. The association between the food chain network and the phylogeny was explored using a network comparison approach. Food and environmental samples were taken from premises linked to cases and tested for Salmonella. Within the outbreak single nucleotide polymorphism defined cluster, 136 cases were identified in the UK and 18 in Spain. One isolate from a food containing chicken eggs was within the outbreak cluster. There was a significant association between the chicken egg food chain of UK cases and phylogeny of outbreak isolates. This is the first published Salmonella outbreak to be prospectively detected using WGS. This outbreak in the UK was linked with contemporaneous cases in Spain by WGS. We conclude that UK and Spanish cases were exposed to a common source of Salmonella-contaminated chicken eggs.
Nosocomial outbreaks of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are frequent despite implementation of conventional infection control measures. An outbreak investigation was ...undertaken using advanced genomic and statistical techniques to reconstruct likely transmission chains and assess the role of healthcare workers (HCWs) in SARS-CoV-2 transmission.
A nosocomial SARS-CoV-2 outbreak in a university-affiliated rehabilitation clinic was investigated, involving patients and HCWs, with high coverage of pathogen whole-genome sequences (WGS). The time-varying reproduction number from epidemiological data (Rt) was estimated, and maximum likelihood phylogeny was used to assess genetic diversity of the pathogen. Genomic and epidemiological data were combined into a Bayesian framework to model the directionality of transmission, and a case–control study was performed to investigate risk factors for nosocomial SARS-CoV-2 acquisition in patients.
The outbreak lasted from 14th March to 12th April 2020, and involved 37 patients (31 with WGS) and 39 employees (31 with WGS), 37 of whom were HCWs. Peak Rt was estimated to be between 2.2 and 3.6. The phylogenetic tree showed very limited genetic diversity, with 60 of 62 (96.7%) isolates forming one large cluster of identical genomes. Despite the resulting uncertainty in reconstructed transmission events, the analyses suggest that HCWs (one of whom was the index case) played an essential role in cross-transmission, with a significantly greater fraction of infections (P<2.2e-16) attributable to HCWs (70.7%) than expected given the number of HCW cases (46.7%). The excess of transmission from HCWs was higher when considering infection of patients 79.0%; 95% confidence interval (CI) 78.5–79.5% and frail patients (Clinical Frailty Scale score >5; 82.3%; 95% CI 81.8–83.4%). Furthermore, frail patients were found to be at greater risk for nosocomial COVID-19 than other patients (adjusted odds ratio 6.94, 95% CI 2.13–22.57).
This outbreak report highlights the essential role of HCWs in SARS-CoV-2 transmission dynamics in healthcare settings. Limited genetic diversity in pathogen genomes hampered the reconstruction of individual transmission events, resulting in substantial uncertainty in who infected whom. However, this study shows that despite such uncertainty, significant transmission patterns can be observed.
PDF
