Probabilistic reconstruction of measles transmission clusters from routinely collected surveillance data

Inference of epidemic dynamics in the COVID-19 era and beyond

The COVID-19 pandemic demonstrated the key role that epidemiology and modelling play in analysing infectious threats and supporting decision making in real-time. Motivated by the unprecedented volume and breadth of data generated during the pandemic, we review modern opportunities for analysis to address questions that emerge during a major modern epidemic. Following the broad chronology of insights required - from understanding initial dynamics to retrospective evaluation of interventions, we describe the theoretical foundations of each approach and the underlying intuition. Through a series of case studies, we illustrate real life applications, and discuss implications for future work.

Quantifying the value of viral genomics when inferring who infected whom in the 2014-16 Ebola virus outbreak in Guinea

Transmission trees can be established through detailed contact histories, statistical or phylogenetic inference, or a combination of methods. Each approach has its limitations, and the extent to which they succeed in revealing a 'true' transmission history remains unclear. In this study, we compared the transmission trees obtained through contact tracing investigations and various inference methods to identify the contribution and value of each approach. We studied eighty-six sequenced cases reported in Guinea between March and November 2015. Contact tracing investigations classified these cases into eight independent transmission chains. We inferred the transmission history from the genetic sequences of the cases (phylogenetic approach), their onset date (epidemiological approach), and a combination of both (combined approach). The inferred transmission trees were then compared to those from the contact tracing investigations. Inference methods using individual data sources (i.e. the phylogenetic analysis and the epidemiological approach) were insufficiently informative to accurately reconstruct the transmission trees and the direction of transmission. The combined approach was able to identify a reduced pool of infectors for each case and highlight likely connections among chains classified as independent by the contact tracing investigations. Overall, the transmissions identified by the contact tracing investigations agreed with the evolutionary history of the viral genomes, even though some cases appeared to be misclassified. Therefore, collecting genetic sequences during outbreak is key to supplement the information contained in contact tracing investigations. Although none of the methods we used could identify one unique infector per case, the combined approach highlighted the added value of mixing epidemiological and genetic information to reconstruct who infected whom.

Understanding the transmission dynamics of a large-scale measles outbreak in Southern Vietnam

OBJECTIVES

During 2018-2020, Southern Vietnam experienced a large measles outbreak of over 26,000 cases. We aimed to understand and quantify the measles spread in space-time dependence and the transmissibility during the outbreak.

METHODS

Measles surveillance reported cases between 1/2018 and 6/2020, vaccination coverage, and population data at provincial level were used. To illustrate the spatiotemporal pattern of disease spread, we employed the endemic-epidemic multivariate time series model decomposing measles risk additively into autoregressive, spatiotemporal, and endemic component. Likelihood-based estimation procedures were performed to determine the time-varying reproductive number R_e of measles.

RESULTS

Our analysis shows that measles incidence was associated with vaccination coverage heterogeneity and spatial interaction between provincial units. The risk of infections was dominated by between-province transmission (36.1% to 78.8%), followed by local endogenous transmission (4.1% to 61.5%) whereas the endemic behavior had a relatively small contribution (2.1% to 33.4%) across provinces. In the exponential phase of the epidemic, R_e was above the threshold with a maximum value of 2.34 (95%CI: 2.20-2.46).

CONCLUSION

Local vaccination coverage and human mobility are important factors contributing to the measles dynamics in Southern Vietnam and the high risk of inter-provincial transmission is of most concern. Strengthening disease surveillance is recommended, and further research is essential to understand the relative contribution of population immunity and control measures in measles epidemics.

Inferring person-to-person networks of Plasmodium falciparum transmission: are analyses of routine surveillance data up to the task?

BACKGROUND

Inference of person-to-person transmission networks using surveillance data is increasingly used to estimate spatiotemporal patterns of pathogen transmission. Several data types can be used to inform transmission network inferences, yet the sensitivity of those inferences to different data types is not routinely evaluated.

METHODS

The influence of different combinations of spatial, temporal, and travel-history data on transmission network inferences for Plasmodium falciparum malaria were evaluated.

RESULTS

The information content of these data types may be limited for inferring person-to-person transmission networks and may lead to an overestimate of transmission. Only when outbreaks were temporally focal or travel histories were accurate was the algorithm able to accurately estimate the reproduction number under control, R_c. Applying this approach to data from Eswatini indicated that inferences of R_c and spatiotemporal patterns therein depend upon the choice of data types and assumptions about travel-history data.

CONCLUSIONS

These results suggest that transmission network inferences made with routine malaria surveillance data should be interpreted with caution.

Basic reproduction number of African swine fever in wild boars (Sus scrofa) and its spatiotemporal heterogeneity in South Korea

Background

African swine fever (ASF) is a hemorrhagic fever occurring in wild boars (Sus scrofa) and domestic pigs. The epidemic situation of ASF in South Korean wild boars has increased the risk of ASF in domestic pig farms. Although basic reproduction number (R₀) can be applied for control policies, it is challenging to estimate the R₀ for ASF in wild boars due to surveillance bias, lack of wild boar population data, and the effect of ASF-positive wild boar carcass on disease dynamics.

Objectives

This study was undertaken to estimate the R₀ of ASF in wild boars in South Korea, and subsequently analyze the spatiotemporal heterogeneity.

Methods

We detected the local transmission clusters using the spatiotemporal clustering algorithm, which was modified to incorporate the effect of ASF-positive wild boar carcass. With the assumption of exponential growth, R₀ was estimated for each cluster. The temporal change of the estimates and its association with the habitat suitability of wild boar were analyzed.

Results

Totally, 22 local transmission clusters were detected, showing seasonal patterns occurring in winter and spring. Mean value of R₀ of each cluster was 1.54. The estimates showed a temporal increasing trend and positive association with habitat suitability of wild boar.

Conclusions

The disease dynamics among wild boars seems to have worsened over time. Thus, in areas with a high elevation and suitable for wild boars, practical methods need to be contrived to ratify the control policies for wild boars.

o2geosocial: Reconstructing who-infected-whom from routinely collected surveillance data

Reconstructing the history of individual transmission events between cases is key to understanding what factors facilitate the spread of an infectious disease. Since conducting extended contact-tracing investigations can be logistically challenging and costly, statistical inference methods have been developed to reconstruct transmission trees from onset dates and genetic sequences. However, these methods are not as effective if the mutation rate of the virus is very slow, or if sequencing data is sparse. We developed the package o2geosocial to combine variables from routinely collected surveillance data with a simple transmission process model. The model reconstructs transmission trees when full genetic sequences are unavailable, or uninformative. Our model incorporates the reported age-group, onset date, location and genotype of infected cases to infer probabilistic transmission trees. The package also includes functions to summarise and visualise the inferred cluster size distribution. The results generated by o2geosocial can highlight regions where importations repeatedly caused large outbreaks, which may indicate a higher regional susceptibility to infections. It can also be used to generate the individual number of secondary transmissions, and show the features associated with individuals involved in high transmission events. The package is available for download from the Comprehensive R Archive Network (CRAN) and GitHub.