Forcing Versus Feedback: Epidemic Malaria and Monsoon Rains in Northwest India

Malaria epidemics in regions with seasonal windows of transmission can vary greatly in size from year to year. A central question has been whether these interannual cycles are driven by climate, are instead generated by the intrinsic dynamics of the disease, or result from the resonance of these two mechanisms. This corresponds to the more general inverse problem of identifying the respective roles of external forcings vs. internal feedbacks from time series for nonlinear and noisy systems. We propose here a quantitative approach to formally compare rival hypotheses on climate vs. disease dynamics, or external forcings vs. internal feedbacks, that combines dynamical models with recently developed, computational inference methods. The interannual patterns of epidemic malaria are investigated here for desert regions of northwest India, with extensive epidemiological records for Plasmodium falciparum malaria for the past two decades. We formulate a dynamical model of malaria transmission that explicitly incorporates rainfall, and we rely on recent advances on parameter estimation for nonlinear and stochastic dynamical systems based on sequential Monte Carlo methods. Results show a significant effect of rainfall in the inter-annual variability of epidemic malaria that involves a threshold in the disease response. The model exhibits high prediction skill for yearly cases in the malaria transmission season following the monsoonal rains. Consideration of a more complex model with clinical immunity demonstrates the robustness of the findings and suggests a role of infected individuals that lack clinical symptoms as a reservoir for transmission. Our results indicate that the nonlinear dynamics of the disease itself play a role at the seasonal, but not the interannual, time scales. They illustrate the feasibility of forecasting malaria epidemics in desert and semi-arid regions of India based on climate variability. This approach should be applicable to malaria in other locations, to other infectious diseases, and to other nonlinear systems under forcing.     

Phylogenetic surveillance of viral genetic diversity and the evolving molecular epidemiology of human immunodeficiency virus type 1

With ongoing generation of viral genetic diversity and increasing levels of migration, the global human immunodeficiency virus type 1 (HIV-1) epidemic is becoming increasingly heterogeneous. In this study, we investigate the epidemiological characteristics of 5,675 HIV-1 pol gene sequences sampled from distinct infections in the United Kingdom. These sequences were phylogenetically analyzed in conjunction with 976 complete-genome and 3,201 pol gene reference sequences sampled globally and representing the broad range of HIV-1 genetic diversity, allowing us to estimate the probable geographic origins of the various strains present in the United Kingdom. A statistical analysis of phylogenetic clustering in this data set identified several independent transmission chains within the United Kingdom involving recently introduced strains and indicated that strains more commonly associated with infections acquired heterosexually in East Africa are spreading among men who have sex with men. Coalescent approaches were also used and indicated that the transmission chains that we identify originated in the late 1980s to early 1990s. Similar changes in the epidemiological structuring of HIV epidemics are likely to be taking in place in other industrialized nations with large immigrant populations. The framework implemented here takes advantage of the vast amount of routinely generated HIV-1 sequence data and can provide epidemiological insights not readily obtainable through standard surveillance methods.     

Quantifying transmission dynamics of acute hepatitis C virus infections in a heterogeneous population using sequence data

Opioid substitution and syringes exchange programs have drastically reduced hepatitis C virus (HCV) spread in France but HCV sexual transmission in men having sex with men (MSM) has recently arisen as a significant public health concern. The fact that the virus is transmitting in a heterogeneous population, with different transmission routes, makes prevalence and incidence rates poorly informative. However, additional insights can be gained by analyzing virus phylogenies inferred from dated genetic sequence data. By combining a phylodynamics approach based on Approximate Bayesian Computation (ABC) and an original transmission model, we estimate key epidemiological parameters of an ongoing HCV epidemic among MSMs in Lyon (France). We show that this new epidemic is largely independent of the previously observed non-MSM HCV epidemics and that its doubling time is ten times lower (0.44 years versus 4.37 years). These results have practical implications for HCV control and illustrate the additional information provided by virus genomics in public health.     

Molecular epidemiology of the foot-and-mouth disease virus outbreak in the United Kingdom in 2001

The objective of this study was to quantify the extent to which the genetic diversity of foot-and-mouth disease virus (FMDV) arising over the course of infection of an individual animal becomes fixed, is transmitted to other animals, and thereby accumulates over the course of an outbreak. Complete consensus sequences of 23 genomes (each of 8,200 nucleotides) of FMDV were recovered directly from epithelium tissue acquired from 21 farms infected over a nearly 7-month period during the 2001 FMDV outbreak in the United Kingdom. An analysis of these consensus sequences revealed very few apparently ambiguous sites but clear evidence of 197 nucleotide substitutions at 191 different sites. We estimated the rate of nucleotide substitution to be 2.26 x 10(-5) per site per day (95% confidence interval [CI], 1.75 x 10(-5) to 2.80 x 10(-5)) and nucleotide substitutions to accrue in the consensus sequence at an average rate of 1.5 substitutions per farm infection. This is a sufficiently high rate showing that detailed histories of the transmission pathways can be reliably reconstructed. Coalescent methods indicated that the date at which FMDV first infected livestock in the United Kingdom was 7 February 2001 (95% CI, 20 January to 19 February 2001), which was identical to estimates obtained on the basis of purely clinical evidence. Nucleotide changes appeared to have occurred evenly across the genome, and within the open reading frame, the ratio of nonsynonymous-to-synonymous change was 0.09. The ability to recover particular transmission pathways of acutely acting RNA pathogens from genetic data will help resolve uncertainties about how virus is spread and could help in the control of future epidemics.     

The epidemiology and iatrogenic transmission of hepatitis C virus in Egypt: a Bayesian coalescent approach

Hepatitis C virus (HCV) is a leading cause of liver cancer and cirrhosis, and Egypt has possibly the highest HCV prevalence worldwide. In this article we use a newly developed Bayesian inference framework to estimate the transmission dynamics of HCV in Egypt from sampled viral gene sequences, and to predict the public health impact of the virus. Our results indicate that the effective number of HCV infections in Egypt underwent rapid exponential growth between 1930 and 1955. The timing and speed of this spread provides quantitative genetic evidence that the Egyptian HCV epidemic was initiated and propagated by extensive antischistosomiasis injection campaigns. Although our results show that HCV transmission has since decreased, we conclude that HCV is likely to remain prevalent in Egypt for several decades. Our combined population genetic and epidemiological analysis provides detailed estimates of historical changes in Egyptian HCV prevalence. Because our results are consistent with a demographic scenario specified a priori, they also provide an objective test of inference methods based on the coalescent process.     

Transmission parameters of the 2001 foot and mouth epidemic in Great Britain

Despite intensive ongoing research, key aspects of the spatial-temporal evolution of the 2001 foot and mouth disease (FMD) epidemic in Great Britain (GB) remain unexplained. Here we develop a Markov Chain Monte Carlo (MCMC) method for estimating epidemiological parameters of the 2001 outbreak for a range of simple transmission models. We make the simplifying assumption that infectious farms were completely observed in 2001, equivalent to assuming that farms that were proactively culled but not diagnosed with FMD were not infectious, even if some were infected. We estimate how transmission parameters varied through time, highlighting the impact of the control measures on the progression of the epidemic. We demonstrate statistically significant evidence for assortative contact patterns between animals of the same species. Predictive risk maps of the transmission potential in different geographic areas of GB are presented for the fitted models.     

Regression-Based Ranking of Pathogen Strains with Respect to Their Contribution to Natural Epidemics

Genetic variation in pathogen populations may be an important factor driving heterogeneity in disease dynamics within their host populations. However, to date, we understand poorly how genetic diversity in diseases impact on epidemiological dynamics because data and tools required to answer this questions are lacking. Here, we combine pathogen genetic data with epidemiological monitoring of disease progression, and introduce a statistical exploratory method to investigate differences among pathogen strains in their performance in the field. The method exploits epidemiological data providing a measure of disease progress in time and space, and genetic data indicating the relative spatial patterns of the sampled pathogen strains. Applying this method allows to assign ranks to the pathogen strains with respect to their contributions to natural epidemics and to assess the significance of the ranking. This method was first tested on simulated data, including data obtained from an original, stochastic, multi-strain epidemic model. It was then applied to epidemiological and genetic data collected during one natural epidemic of powdery mildew occurring in its wild host population. Based on the simulation study, we conclude that the method can achieve its aim of ranking pathogen strains if the sampling effort is sufficient. For powdery mildew data, the method indicated that one of the sampled strains tends to have a higher fitness than the four other sampled strains, highlighting the importance of strain diversity for disease dynamics. Our approach allowing the comparison of pathogen strains in natural epidemic is complementary to the classical practice of using experimental infections in controlled conditions to estimate fitness of different pathogen strains. Our statistical tool, implemented in the R package StrainRanking, is mainly based on regression and does not rely on mechanistic assumptions on the pathogen dynamics. Thus, the method can be applied to a wide range of pathogens.     

Serotype-Specific Transmission and Waning Immunity of Endemic Foot-and-Mouth Disease Virus in Cameroon

Foot-and-mouth disease virus (FMDV) causes morbidity and mortality in a range of animals and threatens local economies by acting as a barrier to international trade. The outbreak in the United Kingdom in 2001 that cost billions to control highlighted the risk that the pathogen poses to agriculture. In response, several mathematical models have been developed to parameterize and predict both transmission dynamics and optimal disease control. However, a lack of understanding of the multi-strain etiology prevents characterization of multi-strain dynamics. Here, we use data from FMDV serology in an endemic setting to probe strain-specific transmission and immunodynamics. Five serotypes of FMDV affect cattle in the Far North Region of Cameroon. We fit both catalytic and reverse catalytic models to serological data to estimate the force of infection and the rate of waning immunity, and to detect periods of sustained transmission. For serotypes SAT2, SAT3, and type A, a model assuming life-long immunity fit better. For serotypes SAT1 and type O, the better-fit model suggests that immunity may wane over time. Our analysis further indicates that type O has the greatest force of infection and the longest duration of immunity. Estimates for the force of infection were time-varying and indicated that serotypes SAT1 and O displayed endemic dynamics, serotype A displayed epidemic dynamics, and SAT2 and SAT3 did not sustain local chains of transmission. Since these results were obtained from the same population at the same time, they highlight important differences in transmission specific to each serotype. They also show that immunity wanes at rates specific to each serotype, which influences patterns of local persistence. Overall, this work shows that viral serotypes can differ significantly in their epidemiological and immunological characteristics. Patterns and processes that drive transmission in endemic settings must consider complex viral dynamics for accurate representation and interpretation.     

Modeling routes of chronic wasting disease transmission: environmental prion persistence promotes deer population decline and extinction

Chronic wasting disease (CWD) is a fatal disease of deer, elk, and moose transmitted through direct, animal-to-animal contact, and indirectly, via environmental contamination. Considerable attention has been paid to modeling direct transmission, but despite the fact that CWD prions can remain infectious in the environment for years, relatively little information exists about the potential effects of indirect transmission on CWD dynamics. In the present study, we use simulation models to demonstrate how indirect transmission and the duration of environmental prion persistence may affect epidemics of CWD and populations of North American deer. Existing data from Colorado, Wyoming, and Wisconsin's CWD epidemics were used to define plausible short-term outcomes and associated parameter spaces. Resulting long-term outcomes range from relatively low disease prevalence and limited host-population decline to host-population collapse and extinction. Our models suggest that disease prevalence and the severity of population decline is driven by the duration that prions remain infectious in the environment. Despite relatively low epidemic growth rates, the basic reproductive number, R(0), may be much larger than expected under the direct-transmission paradigm because the infectious period can vastly exceed the host's life span. High prion persistence is expected to lead to an increasing environmental pool of prions during the early phases (i.e. approximately during the first 50 years) of the epidemic. As a consequence, over this period of time, disease dynamics will become more heavily influenced by indirect transmission, which may explain some of the observed regional differences in age and sex-specific disease patterns. This suggests management interventions, such as culling or vaccination, will become increasingly less effective as CWD epidemics progress.     

Estimating the basic reproductive number from viral sequence data

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.     

Unravelling transmission trees of infectious diseases by combining genetic and epidemiological data

Knowledge on the transmission tree of an epidemic can provide valuable insights into disease dynamics. The transmission tree can be reconstructed by analysing either detailed epidemiological data (e.g. contact tracing) or, if sufficient genetic diversity accumulates over the course of the epidemic, genetic data of the pathogen. We present a likelihood-based framework to integrate these two data types, estimating probabilities of infection by taking weighted averages over the set of possible transmission trees. We test the approach by applying it to temporal, geographical and genetic data on the 241 poultry farms infected in an epidemic of avian influenza A (H7N7) in The Netherlands in 2003. We show that the combined approach estimates the transmission tree with higher correctness and resolution than analyses based on genetic or epidemiological data alone. Furthermore, the estimated tree reveals the relative infectiousness of farms of different types and sizes.     

Transmission of Equine Influenza Virus during an Outbreak Is Characterized by Frequent Mixed Infections and Loose Transmission Bottlenecks

The ability of influenza A viruses (IAVs) to cross species barriers and evade host immunity is a major public health concern. Studies on the phylodynamics of IAVs across different scales – from the individual to the population – are essential for devising effective measures to predict, prevent or contain influenza emergence. Understanding how IAVs spread and evolve during outbreaks is critical for the management of epidemics. Reconstructing the transmission network during a single outbreak by sampling viral genetic data in time and space can generate insights about these processes. Here, we obtained intra-host viral sequence data from horses infected with equine influenza virus (EIV) to reconstruct the spread of EIV during a large outbreak. To this end, we analyzed within-host viral populations from sequences covering 90% of the infected yards. By combining gene sequence analyses with epidemiological data, we inferred a plausible transmission network, in turn enabling the comparison of transmission patterns during the course of the outbreak and revealing important epidemiological features that were not apparent using either approach alone. The EIV populations displayed high levels of genetic diversity, and in many cases we observed distinct viral populations containing a dominant variant and a number of related minor variants that were transmitted between infectious horses. In addition, we found evidence of frequent mixed infections and loose transmission bottlenecks in these naturally occurring populations. These frequent mixed infections likely influence the size of epidemics.     

Enhancing Bayesian risk prediction for epidemics using contact tracing

Contact-tracing data (CTD) collected from disease outbreaks has received relatively little attention in the epidemic modeling literature because it is thought to be unreliable: infection sources might be wrongly attributed, or data might be missing due to resource constraints in the questionnaire exercise. Nevertheless, these data might provide a rich source of information on the disease transmission rate. This paper presents a novel methodology for combining CTD with rate-based contact network data to improve posterior precision, and therefore predictive accuracy. We present an advancement in Bayesian inference for epidemics that assimilates these data and is robust to partial contact tracing. Using a simulation study based on the British poultry industry, we show how the presence of CTD improves posterior predictive accuracy and can directly inform a more effective control strategy.     

Reconstruction of the Transmission History of RNA Virus Outbreaks Using Full Genome Sequences: Foot-and-Mouth Disease Virus in Bulgaria in 2011

Improvements to sequencing protocols and the development of computational phylogenetics have opened up opportunities to study the rapid evolution of RNA viruses in real time. In practical terms, these results can be combined with field data in order to reconstruct spatiotemporal scenarios that describe the origin and transmission pathways of viruses during an epidemic. In the case of notifiable diseases, such as foot-and-mouth disease (FMD), these analyses provide important insights into the epidemiology of field outbreaks that can support disease control programmes. This study reconstructs the origin and transmission history of the FMD outbreaks which occurred during 2011 in Burgas Province, Bulgaria, a country that had been previously FMD-free-without-vaccination since 1996. Nineteen full genome sequences (FGS) of FMD virus (FMDV) were generated and analysed, including eight representative viruses from all of the virus-positive outbreaks of the disease in the country and 11 closely-related contemporary viruses from countries in the region where FMD is endemic (Turkey and Israel). All Bulgarian sequences shared a single putative common ancestor which was closely related to the index case identified in wild boar. The closest relative from outside of Bulgaria was a FMDV collected during 2010 in Bursa (Anatolia, Turkey). Within Bulgaria, two discrete genetic clusters were detected that corresponded to two episodes of outbreaks that occurred during January and March-April 2011. The number of nucleotide substitutions that were present between, and within, these separate clusters provided evidence that undetected FMDV infection had occurred. These conclusions are supported by laboratory data that subsequently identified three additional FMDV-infected livestock premises by serosurveillance, as well as a number of antibody positive wild boar on both sides of the border with Turkish Thrace. This study highlights how FGS analysis can be used as an effective on-the-spot tool to support and help direct epidemiological investigations of field outbreaks.     

Coalescent inference for infectious disease: meta-analysis of hepatitis C

Genetic analysis of pathogen genomes is a powerful approach to investigating the population dynamics and epidemic history of infectious diseases. However, the theoretical underpinnings of the most widely used, coalescent methods have been questioned, casting doubt on their interpretation. The aim of this study is to develop robust population genetic inference for compartmental models in epidemiology. Using a general approach based on the theory of metapopulations, we derive coalescent models under susceptible–infectious (SI), susceptible–infectious–susceptible (SIS) and susceptible–infectious–recovered (SIR) dynamics. We show that exponential and logistic growth models are equivalent to SI and SIS models, respectively, when co-infection is negligible. Implementing SI, SIS and SIR models in BEAST, we conduct a meta-analysis of hepatitis C epidemics, and show that we can directly estimate the basic reproductive number (R₀) and prevalence under SIR dynamics. We find that differences in genetic diversity between epidemics can be explained by differences in underlying epidemiology (age of the epidemic and local population density) and viral subtype. Model comparison reveals SIR dynamics in three globally restricted epidemics, but most are better fit by the simpler SI dynamics. In summary, metapopulation models provide a general and practical framework for integrating epidemiology and population genetics for the purposes of joint inference.     

A clonal analysis of Neisseria meningitidis serogroup A

A typing scheme has recently been developed for Neisseria meningitidis serogroup A based on the clonal population structure of these bacteria. An international strain collection consisting of 423 group A strains isolated from 23 epidemics or outbreaks since 1963, as well as from older epidemics and numerous non-epidemic situations was used in the analysis. Strains were first segregated into electrophoretic types, depending on the combined score for the electrophoretic mobilities of 7 cytoplasmic isoenzymes resolved by starch gel electrophoresis and of 2 outer membrane proteins resolved by sodium dodecyl sulphate polyacrylamide gel electrophoresis. The bacteria were subsequently assigned to one of 21 clones after numerical analysis of their electrophoretic types.The epidemiological value of the typing scheme was assessed by examining case and carrier strains isolated during (1982–83) and subsequent to (1984–85) an epidemic in the Gambia, West Africa. The case isolates, all of which were serogroup A, were of a single clonal type. All serogroup A carrier isolates were also of this clone, while carrier strains of other serogroups showed greater clonal diversity. These results indicate that case strains during an epidemic show little clonal diversity and thus that the typing scheme is of value in distinguishing the etiology of epidemics.A retrospective epidemiological analysis of the strains in the international collection showed that most serogroup A epidemics were associated with a single or predominant clone, although some epidemics were of mixed etiology. The survey included 256 isolates from 15 African epidemics since 1963, a period which covers 3 major epidemic waves (1960–63; 1967–73 and 1981–83), thus permitting a detailed epidemiological analysis of serogroup A epidemics in this continent.Epidemiological records indicate that seven clones have been responsible for sets of epidemics throughout the world since 1915 and that at least two of these sets can be considered to represent mutually exclusive pandemics, linking numerous epidemics between 1967–75 and 1973–83, respectively.     

Inferring the Source of Transmission with Phylogenetic Data

Identifying the source of transmission using pathogen genetic data is complicated by numerous biological, immunological, and behavioral factors. A large source of error arises when there is incomplete or sparse sampling of cases. Unsampled cases may act as either a common source of infection or as an intermediary in a transmission chain for hosts infected with genetically similar pathogens. It is difficult to quantify the probability of common source or intermediate transmission events, which has made it difficult to develop statistical tests to either confirm or deny putative transmission pairs with genetic data. We present a method to incorporate additional information about an infectious disease epidemic, such as incidence and prevalence of infection over time, to inform estimates of the probability that one sampled host is the direct source of infection of another host in a pathogen gene genealogy. These methods enable forensic applications, such as source-case attribution, for infectious disease epidemics with incomplete sampling, which is usually the case for high-morbidity community-acquired pathogens like HIV, Influenza and Dengue virus. These methods also enable epidemiological applications such as the identification of factors that increase the risk of transmission. We demonstrate these methods in the context of the HIV epidemic in Detroit, Michigan, and we evaluate the suitability of current sequence databases for forensic and epidemiological investigations. We find that currently available sequences collected for drug resistance testing of HIV are unlikely to be useful in most forensic investigations, but are useful for identifying transmission risk factors.     

Evolutionary analysis of inter-farm transmission dynamics in a highly pathogenic avian influenza epidemic

Phylogenetic studies have largely contributed to better understand the emergence, spread and evolution of highly pathogenic avian influenza during epidemics, but sampling of genetic data has never been detailed enough to allow mapping of the spatiotemporal spread of avian influenza viruses during a single epidemic. Here, we present genetic data of H7N7 viruses produced from 72% of the poultry farms infected during the 2003 epidemic in the Netherlands. We use phylogenetic analyses to unravel the pathways of virus transmission between farms and between infected areas. In addition, we investigated the evolutionary processes shaping viral genetic diversity, and assess how they could have affected our phylogenetic analyses. Our results show that the H7N7 virus was characterized by a high level of genetic diversity driven mainly by a high neutral substitution rate, purifying selection and limited positive selection. We also identified potential reassortment in the three genes that we have tested, but they had only a limited effect on the resolution of the inter-farm transmission network. Clonal sequencing analyses performed on six farm samples showed that at least one farm sample presented very complex virus diversity and was probably at the origin of chronological anomalies in the transmission network. However, most virus sequences could be grouped within clearly defined and chronologically sound clusters of infection and some likely transmission events between farms located 0.8-13 Km apart were identified. In addition, three farms were found as most likely source of virus introduction in distantly located new areas. These long distance transmission events were likely facilitated by human-mediated transport, underlining the need for strict enforcement of biosafety measures during outbreaks. This study shows that in-depth genetic analysis of virus outbreaks at multiple scales can provide critical information on virus transmission dynamics and can be used to increase our capacity to efficiently control epidemics.     

Modelling under-reporting in epidemics

Under-reporting of infected cases is crucial for many diseases because of the bias it can introduce when making inference for the model parameters. The objective of this paper is to study the effect of under-reporting in epidemics by considering the stochastic Markovian SIR epidemic in which various reporting processes are incorporated. In particular, we first investigate the effect on the estimation process of ignoring under-reporting when it is present in an epidemic outbreak. We show that such an approach leads to under-estimation of the infection rate and the reproduction number. Secondly, by allowing for the fact that under-reporting is occurring, we develop suitable models for estimation of the epidemic parameters and explore how well the reporting rate and other model parameters can be estimated. We consider the case of a constant reporting probability and also more realistic assumptions which involve the reporting probability depending on time or the source of infection for each infected individual. Due to the incomplete nature of the data and reporting process, the Bayesian approach provides a natural modelling framework and we perform inference using data augmentation and reversible jump Markov chain Monte Carlo techniques.     

Evolution of the HIV-1 V3 region in the Italian epidemic

The epidemic of Human Immunodeficiency Virus type 1 infection in Italy is mostly ascribed to the B subtype, which represents the prevalent subtype in Western Countries. The virus isolates of the B subtype, moreover, show an increasing nucleotide heterogeneity over time, indicating a continuous intra-subtype dynamic evolution, typical of long-lasting epidemics. In recent years, however, the progressive decrease in the transmission rate among the historically defined risk groups (i.e. homosexuals and IDUs) and the parallel increase in heterosexual transmission are slowly introducing variants of non-B subtype into the Italian HIV-1 epidemic. This appears to be strictly linked to the growing number of immigrants from non-Western Countries, where non-B clades and CRFs are prevalent, and consequent inter-racial blending. The distribution of these novel genetic forms needs to be evaluated by continuous molecular monitoring nationwide to verify whether they will overcome the pre-existing B-clade epidemic, which could have significant implications for diagnosis, treatment and vaccine development. Here we review the genetic evolution of HIV-1 spreading within the Italian epidemic.