A Bayesian approach for inferring the dynamics of partially observed endemic infectious diseases from space-time-genetic data

We describe a statistical framework for reconstructing the sequence of transmission events between observed cases of an endemic infectious disease using genetic, temporal and spatial information. Previous approaches to reconstructing transmission trees have assumed all infections in the study area originated from a single introduction and that a large fraction of cases were observed. There are as yet no approaches appropriate for endemic situations in which a disease is already well established in a host population and in which there may be multiple origins of infection, or that can enumerate unobserved infections missing from the sample. Our proposed framework addresses these shortcomings, enabling reconstruction of partially observed transmission trees and estimating the number of cases missing from the sample. Analyses of simulated datasets show the method to be accurate in identifying direct transmissions, while introductions and transmissions via one or more unsampled intermediate cases could be identified at high to moderate levels of case detection. When applied to partial genome sequences of rabies virus sampled from an endemic region of South Africa, our method reveals several distinct transmission cycles with little contact between them, and direct transmission over long distances suggesting significant anthropogenic influence in the movement of infected dogs.     

Reconstructing disease outbreaks from genetic data: a graph approach

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.     

A computationally tractable birth-death model that combines phylogenetic and epidemiological data

Inferring the dynamics of pathogen transmission during an outbreak is an important problem in infectious disease epidemiology. In mathematical epidemiology, estimates are often informed by time series of confirmed cases, while in phylodynamics genetic sequences of the pathogen, sampled through time, are the primary data source. Each type of data provides different, and potentially complementary, insight. Recent studies have recognised that combining data sources can improve estimates of the transmission rate and the number of infected individuals. However, inference methods are typically highly specialised and field-specific and are either computationally prohibitive or require intensive simulation, limiting their real-time utility. We present a novel birth-death phylogenetic model and derive a tractable analytic approximation of its likelihood, the computational complexity of which is linear in the size of the dataset. This approach combines epidemiological and phylodynamic data to produce estimates of key parameters of transmission dynamics and the unobserved prevalence. Using simulated data, we show (a) that the approximation agrees well with existing methods, (b) validate the claim of linear complexity and (c) explore robustness to model misspecification. This approximation facilitates inference on large datasets, which is increasingly important as large genomic sequence datasets become commonplace.     

Characterizing the Transmission Potential of Zoonotic Infections from Minor Outbreaks

The transmission potential of a novel infection depends on both the inherent transmissibility of a pathogen, and the level of susceptibility in the host population. However, distinguishing between these pathogen- and population-specific properties typically requires detailed serological studies, which are rarely available in the early stages of an outbreak. Using a simple transmission model that incorporates age-stratified social mixing patterns, we present a novel method for characterizing the transmission potential of subcritical infections, which have effective reproduction number R<1, from readily available data on the size of outbreaks. We show that the model can identify the extent to which outbreaks are driven by inherent pathogen transmissibility and pre-existing population immunity, and can generate unbiased estimates of the effective reproduction number. Applying the method to real-life infections, we obtained accurate estimates for the degree of age-specific immunity against monkeypox, influenza A(H5N1) and A(H7N9), and refined existing estimates of the reproduction number. Our results also suggest minimal pre-existing immunity to MERS-CoV in humans. The approach we describe can therefore provide crucial information about novel infections before serological surveys and other detailed analyses are available. The methods would also be applicable to data stratified by factors such as profession or location, which would make it possible to measure the transmission potential of emerging infections in a wide range of settings.     

Total infectome characterization of respiratory infections in pre-COVID-19 Wuhan,China

At the end of 2019 Wuhan witnessed an outbreak of “atypical pneumonia” that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their “total infectome”, including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen–Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016–2017.     

Comparison of Salmonella enterica serotype Infantis isolates from a veterinary teaching hospital

AIMS: To compare Salmonella enterica serotype Infantis isolates obtained from patients or the environment of a veterinary teaching hospital over a period of 9 years following a nosocomial outbreak to determine whether isolates were epidemiologically related or represented unrelated introductions into the hospital environment. METHODS AND RESULTS: Fifty-six S. Infantis isolates were compared based on their phenotypic (antimicrobial drug [AMD] susceptibility pattern) and genotypic (pulsed-field gel electrophoresis [PFGE] pattern and presence of integrons) characteristics. Epidemiologically unrelated S. Infantis isolates clustered separately from all but two of the hospital isolates, and several isolates from different years and various sources were indistinguishable from each other in cluster analysis of two-enzyme PFGE results. A high percentage of isolates (80.3%) were resistant to at least one AMD, with 67.8% showing resistance to >5 AMD. The majority (74.1%) of isolates tested contained type 1 integrons. CONCLUSION: Results strongly suggest that there was nosocomial transmission of S. Infantis during the initial outbreak, and that contamination arising from this outbreak persisted across years despite rigorous hygiene and biosecurity precautions and may have led to subsequent nosocomial infections. SIGNIFICANCE AND IMPACT OF THE STUDY: Evidence of persistence and transmission of Salmonella clones across years, even in the face of rigorous preventive measures, has important implications for other facilities that have experienced outbreaks of Salmonella infections.     

'SEEDY' (Simulation of Evolutionary and Epidemiological Dynamics): An R Package to Follow Accumulation of Within-Host Mutation in Pathogens

Genome sequencing is an increasingly common component of infectious disease outbreak investigations. However, the relationship between pathogen transmission and observed genetic data is complex, and dependent on several uncertain factors. As such, simulation of pathogen dynamics is an important tool for interpreting observed genomic data in an infectious disease outbreak setting, in order to test hypotheses and to explore the range of outcomes consistent with a given set of parameters. We introduce ‘seedy’, an R package for the simulation of evolutionary and epidemiological dynamics (http://cran.r-project.org/web/packages/seedy/). Our software implements stochastic models for the accumulation of mutations within hosts, as well as individual-level disease transmission. By allowing variables such as the transmission bottleneck size, within-host effective population size and population mixing rates to be specified by the user, our package offers a flexible framework to investigate evolutionary dynamics during disease outbreaks. Furthermore, our software provides theoretical pairwise genetic distance distributions to provide a likelihood of person-to-person transmission based on genomic observations, and using this framework, implements transmission route assessment for genomic data collected during an outbreak. Our open source software provides an accessible platform for users to explore pathogen evolution and outbreak dynamics via simulation, and offers tools to assess observed genomic data in this context.     

Distinguishing epidemic waves from disease spillover in a wildlife population

Serengeti lions frequently experience viral outbreaks. In 1994, one-third of Serengeti lions died from canine distemper virus (CDV). Based on the limited epidemiological data available from this period, it has been unclear whether the 1994 outbreak was propagated by lion-to-lion transmission alone or involved multiple introductions from other sympatric carnivore species. More broadly, we do not know whether contacts between lions allow any pathogen with a relatively short infectious period to percolate through the population (i.e. reach epidemic proportions). We built one of the most realistic contact network models for a wildlife population to date, based on detailed behavioural and movement data from a long-term lion study population. The model allowed us to identify previously unrecognized biases in the sparse data from the 1994 outbreak and develop methods for judiciously inferring disease dynamics from typical wildlife samples. Our analysis of the model in light of the 1994 outbreak data strongly suggest that, although lions are sufficiently well connected to sustain epidemics of CDV-like diseases, the 1994 epidemic was fuelled by multiple spillovers from other carnivore species, such as jackals and hyenas.     

Validation of the Early Warning and Response System (EWARS) for dengue outbreaks: Evidence from the national vector control program in Mexico

BackgroundDuring 2017, twenty health districts (locations) implemented a dengue outbreak Early Warning and Response System (EWARS) in Mexico, which processes epidemiological, meteorological and entomological alarm indicators to predict dengue outbreaks and triggers early response activities.Out of the 20 priority districts where more than one fifth of all national disease transmission in Mexico occur, eleven districts were purposely selected and analyzed. Nine districts presented outbreak alarms by EWARS but without subsequent outbreaks (“non-outbreak districts”) and two presented alarms with subsequent dengue outbreaks (“outbreak districts”). This evaluation study assesses and compares the impact of alarm-informed response activities and the consequences of failing a timely and adequate response across the outbreak groups.MethodsFive indicators of dengue outbreak response (larval control, entomological studies with water container interventions, focal spraying and indoor residual spraying) were quantitatively analyzed across two groups (”outbreak districts” and “non-outbreak districts”). However, for quality control purposes, only qualitative concluding remarks were derived from the fifth response indicator (fogging).ResultsThe average coverage of vector control responses was significantly higher in non-outbreak districts and across all four indicators. In the “outbreak districts” the response activities started late and were of much lower intensity compared to “non-outbreak districts”. Vector control teams at districts-level demonstrated diverse levels of compliance with local guidelines for ‘initial’, ‘early’ and ‘late’ responses to outbreak alarms, which could potentially explain the different outcomes observed following the outbreak alarms.ConclusionFailing timely and adequate response of alarm signals generated by EWARS showed to negatively impact the disease outbreak control process. On the other hand, districts with adequate and timely response guided by alarm signals demonstrated successful records of outbreak prevention. This study presents important operational scenarios when failing or successding EWARS but warrants investigating the effectiveness and cost-effectiveness of EWARS using a more robust designs.     

A new phylodynamic model of Mycobacterium bovis transmission in a multi-host system uncovers the role of the unobserved reservoir

Multi-host pathogens are particularly difficult to control, especially when at least one of the hosts acts as a hidden reservoir. Deep sequencing of densely sampled pathogens has the potential to transform this understanding, but requires analytical approaches that jointly consider epidemiological and genetic data to best address this problem. While there has been considerable success in analyses of single species systems, the hidden reservoir problem is relatively under-studied. A well-known exemplar of this problem is bovine Tuberculosis, a disease found in British and Irish cattle caused by Mycobacterium bovis, where the Eurasian badger has long been believed to act as a reservoir but remains of poorly quantified importance except in very specific locations. As a result, the effort that should be directed at controlling disease in badgers is unclear. Here, we analyse densely collected epidemiological and genetic data from a cattle population but do not explicitly consider any data from badgers. We use a simulation modelling approach to show that, in our system, a model that exploits available cattle demographic and herd-to-herd movement data, but only considers the ability of a hidden reservoir to generate pathogen diversity, can be used to choose between different epidemiological scenarios. In our analysis, a model where the reservoir does not generate any diversity but contributes to new infections at a local farm scale are significantly preferred over models which generate diversity and/or spread disease at broader spatial scales. While we cannot directly attribute the role of the reservoir to badgers based on this analysis alone, the result supports the hypothesis that under current cattle control regimes, infected cattle alone cannot sustain M. bovis circulation. Given the observed close phylogenetic relationship for the bacteria taken from cattle and badgers sampled near to each other, the most parsimonious hypothesis is that the reservoir is the infected badger population. More broadly, our approach demonstrates that carefully constructed bespoke models can exploit the combination of genetic and epidemiological data to overcome issues of extreme data bias, and uncover important general characteristics of transmission in multi-host pathogen systems.     

Evolutionary History and Phylogeography of Rabies Viruses Associated with Outbreaks in Trinidad

Bat rabies is an emerging disease of public health significance in the Americas. The Caribbean island of Trinidad experiences periodic outbreaks within the livestock population. We performed molecular characterisation of Trinidad rabies virus (RABV) and used a Bayesian phylogeographic approach to investigate the extent to which outbreaks are a result of in situ evolution versus importation of virus from the nearby South American mainland. Trinidadian RABV sequences were confirmed as bat variant and clustered with Desmodus rotundus (vampire bat) related sequences. They fell into two largely temporally defined lineages designated Trinidad I and II. The Trinidad I lineage which included sequences from 1997–2000 (all but two of which were from the northeast of the island) was most closely related to RABV from Ecuador (2005, 2007), French Guiana (1990) and Venezuela (1993, 1994). Trinidad II comprised sequences from the southwest of the island, which clustered into two groups: Trinidad IIa, which included one sequence each from 2000 and 2007, and Trinidad IIb including all 2010 sequences. The Trinidad II sequences were most closely related to sequences from Brazil (1999, 2004) and Uruguay (2007, 2008). Phylogeographic analyses support three separate RABV introductions from the mainland from which each of the three Trinidadian lineages arose. The estimated dates for the introductions and subsequent lineage expansions suggest periods of in situ evolution within Trinidad following each introduction. These data also indicate co-circulation of Trinidad lineage I and IIa during 2000. In light of these findings and the likely vampire bat origin of Trinidadian RABV, further studies should be conducted to investigate the relationship between RABV spatiotemporal dynamics and vampire bat population ecology, in particular any movement between the mainland and Trinidad.     

There are no equal opportunity infectors: Epidemiological modelers must rethink our approach to inequality in infection risk

Mathematical models have come to play a key role in global pandemic preparedness and outbreak response: helping to plan for disease burden, hospital capacity, and inform nonpharmaceutical interventions. Such models have played a pivotal role in the COVID-19 pandemic, with transmission models—and, by consequence, modelers—guiding global, national, and local responses to SARS-CoV-2. However, these models have largely not accounted for the social and structural factors, which lead to socioeconomic, racial, and geographic health disparities. In this piece, we raise and attempt to clarify several questions relating to this important gap in the research and practice of infectious disease modeling: Why do epidemiologic models of emerging infections typically ignore known structural drivers of disparate health outcomes? What have been the consequences of a framework focused primarily on aggregate outcomes on infection equity? What should be done to develop a more holistic approach to modeling-based decision-making during pandemics? In this review, we evaluate potential historical and political explanations for the exclusion of drivers of disparity in infectious disease models for emerging infections, which have often been characterized as “equal opportunity infectors” despite ample evidence to the contrary. We look to examples from other disease systems (HIV, STIs) and successes in including social inequity in models of acute infection transmission as a blueprint for how social connections, environmental, and structural factors can be integrated into a coherent, rigorous, and interpretable modeling framework. We conclude by outlining principles to guide modeling of emerging infections in ways that represent the causes of inequity in infection as central rather than peripheral mechanisms.     

Characteristics of SARS-CoV-2 transmission in a medium-sized city with traditional communities during the early COVID-19 epidemic in China

The nationwide COVID-19 epidemic ended in 2020, a few months after its outbreak in Wuhan, China at the end of 2019. Most COVID-19 cases occurred in Hubei Province, with a few local outbreaks in other provinces of China. A few studies have reported the early SARS-CoV-2 epidemics in several large cities or provinces of China. However, information regarding the early epidemics in small and medium-sized cities, where there are still traditionally large families and community culture is more strongly maintained and thus, transmission profiles may differ, is limited. In this study, we characterized 60 newly sequenced SARS-CoV-2 genomes from Anyang as a representative of small and medium-sized Chinese cities, compared them with more than 400 reference genomes from the early outbreak, and studied the SARS-CoV-2 transmission profiles. Genomic epidemiology revealed multiple SARS-CoV-2 introductions in Anyang and a large-scale expansion of the epidemic because of the large family size. Moreover, our study revealed two transmission patterns in a single outbreak, which were attributed to different social activities. We observed the complete dynamic process of single-nucleotide polymorphism development during community transmission and found that intrahost variant analysis was an effective approach to studying cluster infections. In summary, our study provided new SARS-CoV-2 transmission profiles representative of small and medium-sized Chinese cities as well as information on the evolution of SARS-CoV-2 strains during the early COVID-19 epidemic in China.     

Near real-time surveillance of the SARS-CoV-2 epidemic with incomplete data

When responding to infectious disease outbreaks, rapid and accurate estimation of the epidemic trajectory is critical. However, two common data collection problems affect the reliability of the epidemiological data in real time: missing information on the time of first symptoms, and retrospective revision of historical information, including right censoring. Here, we propose an approach to construct epidemic curves in near real time that addresses these two challenges by 1) imputation of dates of symptom onset for reported cases using a dynamically-estimated “backward” reporting delay conditional distribution, and 2) adjustment for right censoring using the NobBS software package to nowcast cases by date of symptom onset. This process allows us to obtain an approximation of the time-varying reproduction number (R_t) in real time. We apply this approach to characterize the early SARS-CoV-2 outbreak in two Spanish regions between March and April 2020. We evaluate how these real-time estimates compare with more complete epidemiological data that became available later. We explore the impact of the different assumptions on the estimates, and compare our estimates with those obtained from commonly used surveillance approaches. Our framework can help improve accuracy, quantify uncertainty, and evaluate frequently unstated assumptions when recovering the epidemic curves from limited data obtained from public health systems in other locations.     

Relating Phylogenetic Trees to Transmission Trees of Infectious Disease Outbreaks

Transmission events are the fundamental building blocks of the dynamics of any infectious disease. Much about the epidemiology of a disease can be learned when these individual transmission events are known or can be estimated. Such estimations are difficult and generally feasible only when detailed epidemiological data are available. The genealogy estimated from genetic sequences of sampled pathogens is another rich source of information on transmission history. Optimal inference of transmission events calls for the combination of genetic data and epidemiological data into one joint analysis. A key difficulty is that the transmission tree, which describes the transmission events between infected hosts, differs from the phylogenetic tree, which describes the ancestral relationships between pathogens sampled from these hosts. The trees differ both in timing of the internal nodes and in topology. These differences become more pronounced when a higher fraction of infected hosts is sampled. We show how the phylogenetic tree of sampled pathogens is related to the transmission tree of an outbreak of an infectious disease, by the within-host dynamics of pathogens. We provide a statistical framework to infer key epidemiological and mutational parameters by simultaneously estimating the phylogenetic tree and the transmission tree. We test the approach using simulations and illustrate its use on an outbreak of foot-and-mouth disease. The approach unifies existing methods in the emerging field of phylodynamics with transmission tree reconstruction methods that are used in infectious disease epidemiology.     

A Mathematical Model of Chikungunya Dynamics and Control: The Major Epidemic on Réunion Island

Chikungunya is a re-emerging arboviral disease transmitted by Aedes spp. mosquitoes. Although principally endemic to Africa and Asia, recent outbreaks have occurred in Europe following introductions by returning travellers. A particularly large outbreak occurred on Réunion Island in 2006, the published data from which forms the basis of the current study. A simple, deterministic mathematical model of the transmission of the virus between humans and mosquitoes was constructed and parameterised with the up-to-date literature on infection biology. The model is fitted to the large Réunion epidemic, resulting in an estimate of 4.1 for the type reproduction number of chikungunya. Although simplistic, the model provided a close approximation of both the peak incidence of the outbreak and the final epidemic size. Sensitivity analysis using Monte Carlo simulation demonstrated the strong influence that both the latent period of infection in humans and the pre-patent period have on these two epidemiological outcomes. We show why separating these variables, which are epidemiologically distinct in chikungunya infections, is not only necessary for accurate model fitting but also important in informing control.     

Transmission to Eels, Portals of Entry, and Putative Reservoirs of Vibrio vulnificus Serovar E (Biotype 2)

Vibrio vulnificus serovar E (formerly biotype 2) is the etiologic agent that is responsible for the main infectious disease affecting farmed eels. Although the pathogen can theoretically use water as a vehicle for disease transmission, it has not been isolated from tank water during epizootics to date. In this work, the mode of transmission of the disease to healthy eels, the portals of entry of the pathogen into fish, and their putative reservoirs have been investigated by means of laboratory and field experiments. Results of the experiments of direct and indirect host-to-host transmission, patch contact challenges, and oral-anal intubations suggest that water is the prime vehicle for disease transmission and that gills are the main portals of entry into the eel body. The pathogen mixed with food can also come into the fish through the gastrointestinal tract and develop the disease. These conclusions were supported by field data obtained during a natural outbreak in which we were able to isolate this microorganism from tank water for the first time. The examination of some survivors from experimental infections by indirect immunofluorescence and scanning electron microscopy showed that V. vulnificus serovar E formed a biofilm-like structure on the eel skin surface. In vitro assays demonstrated that the ability of the pathogen to colonize both hydrophilic and hydrophobic surfaces was inhibited by glucose. The capacity to form biofilms on eel surface could constitute a strategy for surviving between epizootics or outbreaks, and coated survivors could act as reservoirs for the disease.     

Using Neisseria meningitidis genomic diversity to inform outbreak strain identification

Meningococcal disease is a life-threatening illness caused by the human-restricted bacterium Neisseria meningitidis. Outbreaks in the USA involve at least two cases in an organization or community caused by the same serogroup within three months. Genome comparisons, including phylogenetic analysis and quantification of genome distances can provide confirmatory evidence of pathogen transmission during an outbreak. Interpreting genome distances depends on understanding their distribution both among isolates from outbreaks and among those not from outbreaks. Here, we identify outbreak strains based on phylogenetic relationships among 141 N. meningitidis isolates collected from 28 outbreaks in the USA during 2010–2017 and 1516 non-outbreak isolates collected through contemporaneous meningococcal surveillance. We show that genome distance thresholds based on the maximum SNPs and allele distances among isolates in the phylogenetically defined outbreak strains are sufficient to separate most pairs of non-outbreak isolates into separate strains. Non-outbreak isolate pairs that could not be distinguished from each other based on genetic distances were concentrated in the clonal complexes CC11, CC103, and CC32. Within each of these clonal complexes, phylodynamic analysis identified a group of isolates with extremely low diversity, collected over several years and multiple states. Clusters of isolates with low genetic diversity could indicate increased pathogen transmission, potentially resulting in local outbreaks or nationwide clonal expansions.     

Determinants of the Final Size and Case Rate of Nosocomial Outbreaks

Different nosocomial pathogen species have varying infectivity and durations of infectiousness, while the transmission route determines the contact rate between pathogens and susceptible patients. To determine if the pathogen species and transmission route affects the size and spread of outbreaks, we perform a meta-analysis that examines data from 933 outbreaks of hospital-acquired infection representing 14 pathogen species and 8 transmission routes. We find that the mean number of cases in an outbreak is best predicted by the pathogen species and the mean number of cases per day is best predicted by the species-transmission route combination. Our fitted model predicts the largest mean number of cases for Salmonella outbreaks (22.3) and the smallest mean number of cases for Streptococci outbreaks (8.5). The largest mean number of cases per day occurs during Salmonella outbreaks spread via the environment (0.33) and the smallest occurs for Legionella outbreaks spread by multiple transmission routes (0.005). When combined with information on the frequency of outbreaks these findings could inform the design of infection control policies in hospitals.