Reconciling phylodynamics with epidemiology: the case of dengue virus in southern Vietnam

Leveraging insect-specific viruses to elucidate mosquito population structure and dynamics

Several aspects of mosquito ecology that are important for vectored disease transmission and control have been difficult to measure at epidemiologically important scales in the field. In particular, the ability to describe mosquito population structure and movement rates has been hindered by difficulty in quantifying fine-scale genetic variation among populations. The mosquito virome represents a possible avenue for quantifying population structure and movement rates across multiple spatial scales. Mosquito viromes contain a diversity of viruses, including several insect-specific viruses (ISVs) and "core" viruses that have high prevalence across populations. To date, virome studies have focused on viral discovery and have only recently begun examining viral ecology. While nonpathogenic ISVs may be of little public health relevance themselves, they provide a possible route for quantifying mosquito population structure and dynamics. For example, vertically transmitted viruses could behave as a rapidly evolving extension of the host's genome. It should be possible to apply established analytical methods to appropriate viral phylogenies and incidence data to generate novel approaches for estimating mosquito population structure and dispersal over epidemiologically relevant timescales. By studying the virome through the lens of spatial and genomic epidemiology, it may be possible to investigate otherwise cryptic aspects of mosquito ecology. A better understanding of mosquito population structure and dynamics are key for understanding mosquito-borne disease ecology and methods based on ISVs could provide a powerful tool for informing mosquito control programs.

Epidemiological inference for emerging viruses using segregating sites

Epidemiological models are commonly fit to case and pathogen sequence data to estimate parameters and to infer unobserved disease dynamics. Here, we present an inference approach based on sequence data that is well suited for model fitting early on during the expansion of a viral lineage. Our approach relies on a trajectory of segregating sites to infer epidemiological parameters within a Sequential Monte Carlo framework. Using simulated data, we first show that our approach accurately recovers key epidemiological quantities under a single-introduction scenario. We then apply our approach to SARS-CoV-2 sequence data from France, estimating a basic reproduction number of approximately 2.3-2.7 under an epidemiological model that allows for multiple introductions. Our approach presented here indicates that inference approaches that rely on simple population genetic summary statistics can be informative of epidemiological parameters and can be used for reconstructing infectious disease dynamics during the early expansion of a viral lineage.

Epidemiological inference from pathogen genomes: A review of phylodynamic models and applications

Phylodynamics requires an interdisciplinary understanding of phylogenetics, epidemiology, and statistical inference. It has also experienced more intense application than ever before amid the SARS-CoV-2 pandemic. In light of this, we present a review of phylodynamic models beginning with foundational models and assumptions. Our target audience is public health researchers, epidemiologists, and biologists seeking a working knowledge of the links between epidemiology, evolutionary models, and resulting epidemiological inference. We discuss the assumptions linking evolutionary models of pathogen population size to epidemiological models of the infected population size. We then describe statistical inference for phylodynamic models and list how output parameters can be rearranged for epidemiological interpretation. We go on to cover more sophisticated models and finish by highlighting future directions.

Markov genealogy processes

We construct a family of genealogy-valued Markov processes that are induced by a continuous-time Markov population process. We derive exact expressions for the likelihood of a given genealogy conditional on the history of the underlying population process. These lead to a nonlinear filtering equation which can be used to design efficient Monte Carlo inference algorithms. We demonstrate these calculations with several examples. Existing full-information approaches for phylodynamic inference are special cases of the theory.

Phylodynamic analysis in the understanding of the current COVID-19 pandemic and its utility in vaccine and antiviral design and assessment

Over the last decades, the use of phylogenetic methods in the study of emerging infectious diseases has gained considerable traction in public health. Particularly, the integration of phylogenetic analyses with the understanding of the pathogen dynamics at the population level has provided powerful tools for epidemiological surveillance systems. In the same way, the development of statistical methods and theory, as well as improvement of computational efficiency for evolutionary analysis, has expanded the use of these tools for vaccine and antiviral development. Today with the Coronavirus Disease 2019 (COVID-19), this seems to be critical. In this article, we discuss how the application of phylodynamic analysis can improve the understanding of current pandemic dynamics as well as the design, selection, and evaluation of vaccine candidates and antivirals.

Linked surveillance and genetic data uncovers programmatically relevant geographic scale of Guinea worm transmission in Chad

Background

Guinea worm (Dracunculus medinensis) was detected in Chad in 2010 after a supposed ten-year absence, posing a challenge to the global eradication effort. Initiation of a village-based surveillance system in 2012 revealed a substantial number of dogs infected with Guinea worm, raising questions about paratenic hosts and cross-species transmission.

Methodology/principal findings

We coupled genomic and surveillance case data from 2012-2018 to investigate the modes of transmission between dog and human hosts and the geographic connectivity of worms. Eighty-six variants across four genes in the mitochondrial genome identified 41 genetically distinct worm genotypes. Spatiotemporal modeling revealed worms with the same genotype (‘genetically identical’) were within a median range of 18.6 kilometers of each other, but largely within approximately 50 kilometers. Genetically identical worms varied in their degree of spatial clustering, suggesting there may be different factors that favor or constrain transmission. Each worm was surrounded by five to ten genetically distinct worms within a 50 kilometer radius. As expected, we observed a change in the genetic similarity distribution between pairs of worms using variants across the complete mitochondrial genome in an independent population.

Conclusions/significance

In the largest study linking genetic and surveillance data to date of Guinea worm cases in Chad, we show genetic identity and modeling can facilitate the understanding of local transmission. The co-occurrence of genetically non-identical worms in quantitatively identified transmission ranges highlights the necessity for genomic tools to link cases. The improved discrimination between pairs of worms from variants identified across the complete mitochondrial genome suggests that expanding the number of genomic markers could link cases at a finer scale. These results suggest that scaling up genomic surveillance for Guinea worm may provide additional value for programmatic decision-making critical for monitoring cases and intervention efficacy to achieve elimination.

The global eradication effort for Guinea worm disease has dramatically decreased the global burden of the disease and enabled 187 countries to be certified by the World Health Organization to be free of endemic transmission. Despite this progress, several countries continue to have endemic transmission. In Chad, a long absence of reported cases was interrupted with the identification of new Guinea worm cases, prompting a substantial scale up of surveillance and intervention efforts. Here, we study the value of increasing genomic surveillance as a tool for programmatic evaluation of surveillance and intervention efforts in Chad. Linking surveillance and genomic samples, parsimonious spatial models help reveal a consistent geographic clustering of similar genetic sequences across Chad. We also demonstrate that expanding the sequencing can offer better resolution for distinguishing Guinea worm samples. In this retrospective study, we found evidence that scaling up genomic surveillance can be an important monitoring and evaluation tool for the eradication program in Chad.

What can phylodynamics bring to animal health research?

Infectious diseases are a major burden to global economies, and public and animal health. To date, quantifying the spread of infectious diseases to inform policy making has traditionally relied on epidemiological data collected during epidemics. However, interest has grown in recent phylodynamic techniques to infer pathogen transmission dynamics from genetic data. Here, we provide examples of where this new discipline has enhanced disease management in public health and illustrate how it could be further applied in animal health. In particular, we describe how phylodynamics can address fundamental epidemiological questions, such as inferring key transmission parameters in animal populations and quantifying spillover events at the wildlife-livestock interface, and generate important insights for the design of more effective control strategies.

Accounting for population structure reveals ambiguity in the Zaire Ebolavirus reservoir dynamics

Ebolaviruses pose a substantial threat to wildlife populations and to public health in Africa. Evolutionary analyses of virus genome sequences can contribute significantly to elucidate the origin of new outbreaks, which can help guide surveillance efforts. The reconstructed between-outbreak evolutionary history of Zaire ebolavirus so far has been highly consistent. By removing the confounding impact of population growth bursts during local outbreaks on the free mixing assumption that underlies coalescent-based demographic reconstructions, we find-contrary to what previous results indicated-that the circulation dynamics of Ebola virus in its animal reservoir are highly uncertain. Our findings also accentuate the need for a more fine-grained picture of the Ebola virus diversity in its reservoir to reliably infer the reservoir origin of outbreak lineages. In addition, the recent appearance of slower-evolving variants is in line with latency as a survival mechanism and with bats as the natural reservoir host.

Choices and trade-offs in inference with infectious disease models

Inference using mathematical models of infectious disease dynamics can be an invaluable tool for the interpretation and analysis of epidemiological data. However, researchers wishing to use this tool are faced with a choice of models and model types, simulation methods, inference methods and software packages. Given the multitude of options, it can be challenging to decide on the best approach. Here, we delineate the choices and trade-offs involved in deciding on an approach for inference, and discuss aspects that might inform this decision. We provide examples of inference with a dataset of influenza cases using the R packages pomp and rbi.

Dengue in Santander State, Colombia: fluctuations in the prevalence of virus serotypes are linked to dengue incidence and genetic diversity of the circulating viruses

OBJECTIVE

To investigate the link between fluctuations in the prevalence of dengue virus (DENV) serotypes and the number of dengue cases in the metropolitan area of Bucaramanga, Santander State, Colombia, in the 2007-2010 and 2014-2017 periods.

METHOD

Viruses were isolated from febrile patient samples by direct application to C6/36-HT cells and typed using monoclonal antibodies. We performed autocorrelation and cross-correlation analyses to determine whether fluctuations in the prevalence of DENV serotypes and dengue cases were correlated. Full envelope (E) gene sequences were employed to examine the genetic diversity of serotypes circulating by using a phylogenetic approach.

RESULTS

All four dengue virus serotypes were detected. DENV-1 was the dominant serotype in both periods followed by DENV-3 or DENV-2 depending on the period; DENV-4 was the least prevalent virus in both periods. Cross-correlation analyses suggest a temporal relation between the fluctuations in the prevalence of DENV serotypes, which were almost simultaneous (lag = 0) or related to recent past fluctuations (lag > 1.0) in the number of dengue cases. Data suggest that a sustained predominance of DENV-1, an increase of the DENV-4 prevalence, and a switch from DENV-3 to DENV-2 could be linked to an outbreak. Circulating viruses were grouped into Genotype V, Asia/American III and II for DENV-1, -2, -3 and -4, respectively; intragenotypic diversity was detected.

CONCLUSIONS

The present work highlights the need of comprehensive studies on dynamics of DENV in Colombia to understand transmission of dengue and evaluate the effectiveness of a vaccination programme.

Evolutionary Virology at 40

RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus-host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

Highlights on the Application of Genomics and Bioinformatics in the Fight Against Infectious Diseases: Challenges and Opportunities in Africa

Genomics and bioinformatics are increasingly contributing to our understanding of infectious diseases caused by bacterial pathogens such as Mycobacterium tuberculosis and parasites such as Plasmodium falciparum. This ranges from investigations of disease outbreaks and pathogenesis, host and pathogen genomic variation, and host immune evasion mechanisms to identification of potential diagnostic markers and vaccine targets. High throughput genomics data generated from pathogens and animal models can be combined with host genomics and patients’ health records to give advice on treatment options as well as potential drug and vaccine interactions. However, despite accounting for the highest burden of infectious diseases, Africa has the lowest research output on infectious disease genomics. Here we review the contributions of genomics and bioinformatics to the management of infectious diseases of serious public health concern in Africa including tuberculosis (TB), dengue fever, malaria and filariasis. Furthermore, we discuss how genomics and bioinformatics can be applied to identify drug and vaccine targets. We conclude by identifying challenges to genomics research in Africa and highlighting how these can be overcome where possible.

Bayesian phylodynamic inference with complex models

Population genetic modeling can enhance Bayesian phylogenetic inference by providing a realistic prior on the distribution of branch lengths and times of common ancestry. The parameters of a population genetic model may also have intrinsic importance, and simultaneous estimation of a phylogeny and model parameters has enabled phylodynamic inference of population growth rates, reproduction numbers, and effective population size through time. Phylodynamic inference based on pathogen genetic sequence data has emerged as useful supplement to epidemic surveillance, however commonly-used mechanistic models that are typically fitted to non-genetic surveillance data are rarely fitted to pathogen genetic data due to a dearth of software tools, and the theory required to conduct such inference has been developed only recently. We present a framework for coalescent-based phylogenetic and phylodynamic inference which enables highly-flexible modeling of demographic and epidemiological processes. This approach builds upon previous structured coalescent approaches and includes enhancements for computational speed, accuracy, and stability. A flexible markup language is described for translating parametric demographic or epidemiological models into a structured coalescent model enabling simultaneous estimation of demographic or epidemiological parameters and time-scaled phylogenies. We demonstrate the utility of these approaches by fitting compartmental epidemiological models to Ebola virus and Influenza A virus sequence data, demonstrating how important features of these epidemics, such as the reproduction number and epidemic curves, can be gleaned from genetic data. These approaches are provided as an open-source package PhyDyn for the BEAST2 phylogenetics platform.

Identifying spatio-temporal dynamics of Ebola in Sierra Leone using virus genomes

Containing the recent West African outbreak of Ebola virus (EBOV) required the deployment of substantial global resources. Despite recent progress in analysing and modelling EBOV epidemiological data, a complete characterization of the spatio-temporal spread of Ebola cases remains a challenge. In this work, we offer a novel perspective on the EBOV epidemic in Sierra Leone that uses individual virus genome sequences to inform population-level, spatial models. Calibrated to phylogenetic linkages of virus genomes, these spatial models provide unique insight into the disease mobility of EBOV in Sierra Leone without the need for human mobility data. Consistent with other investigations, our results show that the spread of EBOV during the beginning and middle portions of the epidemic strongly depended on the size of and distance between populations. Our phylodynamic analysis also revealed a change in model preference towards a spatial model with power-law characteristics in the latter portion of the epidemic, correlated with the timing of major intervention campaigns. More generally, we believe this framework, pairing molecular diagnostics with a dynamic model selection procedure, has the potential to be a powerful forecasting tool along with offering operationally relevant guidance for surveillance and sampling strategies during an epidemic.

Tracking external introductions of HIV using phylodynamics reveals a major source of infections in rural KwaZulu-Natal, South Africa

Despite increasing access to antiretrovirals, HIV incidence in rural KwaZulu-Natal remains among the highest ever reported in Africa. While many epidemiological factors have been invoked to explain such high incidence, widespread human mobility and viral movement suggest that transmission between communities may be a major source of new infections. High cross-community transmission rates call into question how effective increasing the coverage of antiretroviral therapy locally will be at preventing new infections, especially if many new cases arise from external introductions. To help address this question, we use a phylodynamic model to reconstruct epidemic dynamics and estimate the relative contribution of local transmission versus external introductions to overall incidence in KwaZulu-Natal from HIV-1 phylogenies. By comparing our results with population-based surveillance data, we show that we can reliably estimate incidence from viral phylogenies once viral movement in and out of the local population is accounted for. Our analysis reveals that early epidemic dynamics were largely driven by external introductions. More recently, we estimate that 35 per cent (95% confidence interval: 20-60%) of new infections arise from external introductions. These results highlight the growing need to consider larger-scale regional transmission dynamics when designing and testing prevention strategies.

Understanding dengue virus evolution to support epidemic surveillance and counter-measure development

Dengue virus (DENV) causes a profound burden of morbidity and mortality, and its global burden is rising due to the co-circulation of four divergent DENV serotypes in the ecological context of globalization, travel, climate change, urbanization, and expansion of the geographic range of the Ae.aegypti and Ae.albopictus vectors. Understanding DENV evolution offers valuable opportunities to enhance surveillance and response to DENV epidemics via advances in RNA virus sequencing, bioinformatics, phylogenetic and other computational biology methods. Here we provide a scoping overview of the evolution and molecular epidemiology of DENV and the range of ways that evolutionary analyses can be applied as a public health tool against this arboviral pathogen.

Challenges in dengue research: A computational perspective

The dengue virus is now the most widespread arbovirus affecting human populations, causing significant economic and social impact in South America and South-East Asia. Increasing urbanization and globalization, coupled with insufficient resources for control, misguided policies or lack of political will, and expansion of its mosquito vectors are some of the reasons why interventions have so far failed to curb this major public health problem. Computational approaches have elucidated on dengue's population dynamics with the aim to provide not only a better understanding of the evolution and epidemiology of the virus but also robust intervention strategies. It is clear, however, that these have been insufficient to address key aspects of dengue's biology, many of which will play a crucial role for the success of future control programmes, including vaccination. Within a multiscale perspective on this biological system, with the aim of linking evolutionary, ecological and epidemiological thinking, as well as to expand on classic modelling assumptions, we here propose, discuss and exemplify a few major computational avenues-real-time computational analysis of genetic data, phylodynamic modelling frameworks, within-host model frameworks and GPU-accelerated computing. We argue that these emerging approaches should offer valuable research opportunities over the coming years, as previously applied and demonstrated in the context of other pathogens.

Genomic and epidemiological characterisation of a dengue virus outbreak among blood donors in Brazil

Outbreaks caused by Dengue, Zika and Chikungunya viruses can spread rapidly in immunologically naïve populations. By analysing 92 newly generated viral genome sequences from blood donors and recipients, we assess the dynamics of dengue virus serotype 4 during the 2012 outbreak in Rio de Janeiro. Phylogenetic analysis indicates that the outbreak was caused by genotype II, although two isolates of genotype I were also detected for the first time in Rio de Janeiro. Evolutionary analysis and modelling estimates are congruent, indicating a reproduction number above 1 between January and June, and at least two thirds of infections being unnoticed. Modelling analysis suggests that viral transmission started in early January, which is consistent with multiple introductions, most likely from the northern states of Brazil, and with an increase in within-country air travel to Rio de Janeiro. The combination of genetic and epidemiological data from blood donor banks may be useful to anticipate epidemic spread of arboviruses.

Phylodynamic Inference across Epidemic Scales

Within-host genetic diversity and large transmission bottlenecks confound phylodynamic inference of epidemiological dynamics. Conventional phylodynamic approaches assume that nodes in a time-scaled pathogen phylogeny correspond closely to the time of transmission between hosts that are ancestral to the sample. However, when hosts harbor diverse pathogen populations, node times can substantially pre-date infection times. Imperfect bottlenecks can cause lineages sampled in different individuals to coalesce in unexpected patterns. To address realistic violations of standard phylodynamic assumptions we developed a new inference approach based on a multi-scale coalescent model, accounting for nonlinear epidemiological dynamics, heterogeneous sampling through time, non-negligible genetic diversity of pathogens within hosts, and imperfect transmission bottlenecks. We apply this method to HIV-1 and Ebola virus (EBOV) outbreak sequence data, illustrating how and when conventional phylodynamic inference may give misleading results. Within-host diversity of HIV-1 causes substantial upwards bias in the number of infected hosts using conventional coalescent models, but estimates using the multi-scale model have greater consistency with reported number of diagnoses through time. In contrast, we find that within-host diversity of EBOV has little influence on estimated numbers of infected hosts or reproduction numbers, and estimates are highly consistent with the reported number of diagnoses through time. The multi-scale coalescent also enables estimation of within-host effective population size using single sequences from a random sample of patients. We find within-host population genetic diversity of HIV-1 p17 to be 2Nμ=0.012 (95% CI 0.0066-0.023), which is lower than estimates based on HIV envelope serial sequencing of individual patients.

Optimal point process filtering and estimation of the coalescent process

The coalescent process is a widely used approach for inferring the demographic history of a population, from samples of its genetic diversity. Several parametric and non-parametric coalescent inference methods, involving Markov chain Monte Carlo, Gaussian processes, and other algorithms, already exist. However, these techniques are not always easy to adapt and apply, thus creating a need for alternative methodologies. We introduce the Bayesian Snyder filter as an easily implementable and flexible minimum mean square error estimator for parametric demographic functions on fixed genealogies. By reinterpreting the coalescent as a self-exciting Markov process, we show that the Snyder filter can be applied to both isochronously and heterochronously sampled datasets. We analytically solve the filter equations for the constant population size Kingman coalescent, derive expressions for its mean squared estimation error, and estimate its robustness to prior distribution specification. For populations with deterministically time-varying size we numerically solve the Snyder equations, and test this solution on common demographic models. We find that the Snyder filter accurately recovers the true demographic history for these models. We also apply the filter to a well-studied, dataset of hepatitis C virus sequences and show that the filter compares well to a popular phylodynamic inference method. The Snyder filter is an exact (given discretised priors, it does not approximate the posterior) and direct Bayesian estimation method that has the potential to become a useful alternative tool for coalescent inference.

Emerging Concepts of Data Integration in Pathogen Phylodynamics

Phylodynamics has become an increasingly popular statistical framework to extract evolutionary and epidemiological information from pathogen genomes. By harnessing such information, epidemiologists aim to shed light on the spatio-temporal patterns of spread and to test hypotheses about the underlying interaction of evolutionary and ecological dynamics in pathogen populations. Although the field has witnessed a rich development of statistical inference tools with increasing levels of sophistication, these tools initially focused on sequences as their sole primary data source. Integrating various sources of information, however, promises to deliver more precise insights in infectious diseases and to increase opportunities for statistical hypothesis testing. Here, we review how the emerging concept of data integration is stimulating new advances in Bayesian evolutionary inference methodology which formalize a marriage of statistical thinking and evolutionary biology. These approaches include connecting sequence to trait evolution, such as for host, phenotypic and geographic sampling information, but also the incorporation of covariates of evolutionary and epidemic processes in the reconstruction procedures. We highlight how a full Bayesian approach to covariate modeling and testing can generate further insights into sequence evolution, trait evolution, and population dynamics in pathogen populations. Specific examples demonstrate how such approaches can be used to test the impact of host on rabies and HIV evolutionary rates, to identify the drivers of influenza dispersal as well as the determinants of rabies cross-species transmissions, and to quantify the evolutionary dynamics of influenza antigenicity. Finally, we briefly discuss how data integration is now also permeating through the inference of transmission dynamics, leading to novel insights into tree-generative processes and detailed reconstructions of transmission trees. [Bayesian inference; birth–death models; coalescent models; continuous trait evolution; covariates; data integration; discrete trait evolution; pathogen phylodynamics.

Synchrony of Dengue Incidence in Ho Chi Minh City and Bangkok

Background

Ho Chi Minh City and Bangkok are highly dengue endemic. The extent to which disease patterns are attributable to local versus regional dynamics remains unclear. To address this gap we compared key transmission parameters across the locations.

Methods and Principal Findings

We used 2003–2009 age-stratified case data to inform catalytic transmission models. Further, we compared the spatial clustering of serotypes within each city. We found that annual case numbers were highly consistent across the two cities (correlation of 0.77, 95% CI: 0.74–0.79) as was the annual force of infection (correlation of 0.57, 95% CI: 0.46–0.68). Serotypes were less similar with serotype-specific correlations ranging from 0.65 for DENV1 to -0.14 for DENV4. Significant spatial clustering of serotypes was observed in HCMC at distances <500m, similar to previous observations from Bangkok.

Discussions

Dengue dynamics are comparable across these two hubs. Low correlation in serotype distribution suggests that similar built environments, vector populations and climate, rather than viral flow drives these observations.

All four serotypes of dengue have circulated endemically throughout Southeast Asia for decades. However, despite the enormous burden of disease, there remains poor understanding of the similarity in disease patterns across the region. We analyzed data from over 100,000 cases of dengue from two of the largest cities in the region, Bangkok and Ho Chi Minh City between 2001 and 2009. We use basic statistical methods to reconstruct the annual probability of infection in the two cities during this time period using methods that are robust to differences in reporting mechanisms. We find that both the epidemic curves and annual probabilities of infection were highly correlated across the cities, however, serotype-specific correlations were far more variable. Finally, we used geocoded case homes from Ho Chi Minh to demonstrate that cases in the city clustered at spatial scales (<500m) similar to that previously observed in Bangkok. These findings show that dengue dynamics are highly comparable across these two urban hubs; however, the low correlation in serotype distribution suggests that similar built environments and climate, rather than viral flow drives these observations.

Identifying predictors of time-inhomogeneous viral evolutionary processes

Various factors determine the rate at which mutations are generated and fixed in viral genomes. Viral evolutionary rates may vary over the course of a single persistent infection and can reflect changes in replication rates and selective dynamics. Dedicated statistical inference approaches are required to understand how the complex interplay of these processes shapes the genetic diversity and divergence in viral populations. Although evolutionary models accommodating a high degree of complexity can now be formalized, adequately informing these models by potentially sparse data, and assessing the association of the resulting estimates with external predictors, remains a major challenge. In this article, we present a novel Bayesian evolutionary inference method, which integrates multiple potential predictors and tests their association with variation in the absolute rates of synonymous and non-synonymous substitutions along the evolutionary history. We consider clinical and virological measures as predictors, but also changes in population size trajectories that are simultaneously inferred using coalescent modelling. We demonstrate the potential of our method in an application to within-host HIV-1 sequence data sampled throughout the infection of multiple patients. While analyses of individual patient populations lack statistical power, we detect significant evidence for an abrupt drop in non-synonymous rates in late stage infection and a more gradual increase in synonymous rates over the course of infection in a joint analysis across all patients. The former is predicted by the immune relaxation hypothesis while the latter may be in line with increasing replicative fitness during the asymptomatic stage.

Trends in the Mechanistic and Dynamic Modeling of Infectious Diseases

The dynamics of infectious disease epidemics are driven by interactions between individuals with differing disease status (e.g., susceptible, infected, immune). Mechanistic models that capture the dynamics of such “dependent happenings” are a fundamental tool of infectious disease epidemiology. Recent methodological advances combined with access to new data sources and computational power have resulted in an explosion in the use of dynamic models in the analysis of emerging and established infectious diseases. Increasing use of models to inform practical public health decision making has challenged the field to develop new methods to exploit available data and appropriately characterize the uncertainty in the results. Here, we discuss recent advances and areas of active research in the mechanistic and dynamic modeling of infectious disease. We highlight how a growing emphasis on data and inference, novel forecasting methods, and increasing access to “big data” are changing the field of infectious disease dynamics. We showcase the application of these methods in phylodynamic research, which combines mechanistic models with rich sources of molecular data to tie genetic data to population-level disease dynamics. As dynamics and mechanistic modeling methods mature and are increasingly tied to principled statistical approaches, the historic separation between the infectious disease dynamics and “traditional” epidemiologic methods is beginning to erode; this presents new opportunities for cross pollination between fields and novel applications.

Epidemiological dynamics of an urban Dengue 4 outbreak in São Paulo, Brazil

Background: Dengue studies at the urban scale are scarce and required for guiding control efforts. In Brazil, the burden of dengue is high and challenges city public health administrations with limited resources. Here we studied the dynamics of a dengue epidemic in a single city. Methods: Serum samples from dengue suspected cases were collected and tested, from December 2012 and July 2013 in Guarujá, Brazil. We use incidence series analysis to provide a detailed view of the reproduction number dynamics and a Bayesian analysis to infer the spread of the serotype using geographic and temporal data. Results: We obtained nucleotide sequences from 354 envelope genes and georeferenced 286 samples during the course of the outbreak. Serotype 4 was responsible for the epidemic. We identified at least two major lineages that overlapped in distribution. We observed high reproduction numbers and high cladogenesis prior to the escalation of clinical case notifications. Three densely populated non-adjacent neighborhoods played a pivotal role during the onset and/or course of the epidemic. Discussion: Our findings point to high dengue virus transmission with a substantial proportion of unapparent cases that led to a late recognition of an outbreak. Usually source reductions initiatives tend to be insufficient once an epidemic has been established. Nevertheless, health authorities in Guarujá prioritized vector control on specific places with clusters of georeferenced viremic patients, which appear to have diminished the epidemic impact.

Getting to the root of epidemic spread with phylodynamic analysis of genomic data

When epidemiological and evolutionary dynamics occur on similar timescales, pathogen genomes sampled from infected hosts carry a signature of the dynamics of epidemic spread. Phylodynamic inference methods aim to extract this signature from genetic data. We discuss the contribution of phylodynamics toward understanding the 2014 West African Ebola virus epidemic.

Eight challenges in phylodynamic inference

The field of phylodynamics, which attempts to enhance our understanding of infectious disease dynamics using pathogen phylogenies, has made great strides in the past decade. Basic epidemiological and evolutionary models are now well characterized with inferential frameworks in place. However, significant challenges remain in extending phylodynamic inference to more complex systems. These challenges include accounting for evolutionary complexities such as changing mutation rates, selection, reassortment, and recombination, as well as epidemiological complexities such as stochastic population dynamics, host population structure, and different patterns at the within-host and between-host scales. An additional challenge exists in making efficient inferences from an ever increasing corpus of sequence data.

Molecular epidemiology, phylogeny and evolution of the filarial nematode Wuchereria bancrofti

Wuchereria bancrofti (Wb) is the most widely distributed of the three nematodes known to cause lymphatic filariasis (LF), the other two being Brugia malayi and Brugia timori. Current tools available to monitor LF are limited to diagnostic tests targeting DNA repeats, filarial antigens, and anti-filarial antibodies. While these tools are useful for detection and surveillance, elimination programs have yet to take full advantage of molecular typing for inferring infection history, strain fingerprinting, and evolution. To date, molecular typing approaches have included whole mitochondrial genomes, genotyping, targeted sequencing, and random amplified polymorphic DNA (RAPDs). These studies have revealed much about Wb biology. For example, in one study in Papua New Guinea researchers identified 5 major strains that were widespread and many minor strains some of which exhibit geographic stratification. Genome data, while rare, has been utilized to reconstruct evolutionary relationships among taxa of the Onchocercidae (the clade of filarial nematodes) and identify gene synteny. Their phylogeny reveals that speciation from the common ancestor of both B. malayi and Wb occurred around 5-6 millions years ago with shared ancestry to other filarial nematodes as recent as 15 million years ago. These discoveries hold promise for gene discovery and identifying drug targets in species that are more amenable to in vivo experiments. Continued technological developments in whole genome sequencing and data analysis will likely replace many other forms of molecular typing, multiplying the amount of data available on population structure, genetic diversity, and phylogenetics. Once widely available, the addition of population genetic data from genomic studies should hasten the elimination of LF parasites like Wb. Infectious disease control programs have benefited greatly from population genetics data and recently from population genomics data. However, while there is currently a surplus of data for diseases like malaria and HIV, there is a scarcity of this data for filarial nematodes. With the falling cost of genome sequencing, research on filarial nematodes could benefit from the addition of population genetics statistics and phylogenetics especially in dealing with elimination programs. A comprehensive review focusing on population genetics of filarial nematode does not yet exist. Here our goal is to provide a current overview of the molecular epidemiology of W. bancrofti (Wb) the primary causative agent of LF. We begin by reviewing studies utilizing molecular typing techniques with specific focus on genomic and population datasets. Next, we used whole mitochondrial genome data to construct a phylogeny and examine the evolutionary history of the Onchocercidae. Then, we provide a perspective to aid in understanding how population genetic techniques translate to modern epidemiology. Finally, we introduce the concept of genomic epidemiology and provide some examples that will aid in future studies of Wb.

Genomic analysis of emerging pathogens: methods, application and future trends

The number of emerging infectious diseases is increasing. Characterizing novel or re-emerging infections is aided by the availability of pathogen genomes. In this review, we evaluate methods that exploit pathogen sequences and the contribution of genomic analysis to understand the epidemiology of recently emerged infectious diseases.

Epidemiologic data and pathogen genome sequences: a powerful synergy for public health

Epidemiologists aim to inform the design of public health interventions with evidence on the evolution, emergence and spread of infectious diseases. Sequencing of pathogen genomes, together with date, location, clinical manifestation and other relevant data about sample origins, can contribute to describing nearly every aspect of transmission dynamics, including local transmission and global spread. The analyses of these data have implications for all levels of clinical and public health practice, from institutional infection control to policies for surveillance, prevention and treatment. This review highlights the range of epidemiological questions that can be addressed from the combination of genome sequence and traditional ‘line lists’ (tables of epidemiological data where each line includes demographic and clinical features of infected individuals). We identify opportunities for these data to inform interventions that reduce disease incidence and prevalence. By considering current limitations of, and challenges to, interpreting these data, we aim to outline a research agenda to accelerate the genomics-driven transformation in public health microbiology.

Phylodynamic inference for structured epidemiological models

Coalescent theory is routinely used to estimate past population dynamics and demographic parameters from genealogies. While early work in coalescent theory only considered simple demographic models, advances in theory have allowed for increasingly complex demographic scenarios to be considered. The success of this approach has lead to coalescent-based inference methods being applied to populations with rapidly changing population dynamics, including pathogens like RNA viruses. However, fitting epidemiological models to genealogies via coalescent models remains a challenging task, because pathogen populations often exhibit complex, nonlinear dynamics and are structured by multiple factors. Moreover, it often becomes necessary to consider stochastic variation in population dynamics when fitting such complex models to real data. Using recently developed structured coalescent models that accommodate complex population dynamics and population structure, we develop a statistical framework for fitting stochastic epidemiological models to genealogies. By combining particle filtering methods with Bayesian Markov chain Monte Carlo methods, we are able to fit a wide class of stochastic, nonlinear epidemiological models with different forms of population structure to genealogies. We demonstrate our framework using two structured epidemiological models: a model with disease progression between multiple stages of infection and a two-population model reflecting spatial structure. We apply the multi-stage model to HIV genealogies and show that the proposed method can be used to estimate the stage-specific transmission rates and prevalence of HIV. Finally, using the two-population model we explore how much information about population structure is contained in genealogies and what sample sizes are necessary to reliably infer parameters like migration rates.

Mathematical models play an important role in our understanding of what processes drive the complex population dynamics of infectious pathogens. Yet developing statistical methods for fitting models to epidemiological data is difficult. Epidemiological data is often noisy, incomplete, aggregated across different scales and generally provides only a partial picture of the underlying disease dynamics. Using nontraditional sources of data, like molecular sequences of pathogens, can provide additional information about epidemiological dynamics. But current “phylodynamic” inference methods for fitting models to genealogies reconstructed from sequence data have a number of major limitations. We present a statistical framework that builds upon earlier work to address two of these limitations: population structure and stochasticity. By incorporating population structure, our framework can be applied in cases where the host population is divided into different subpopulations, such as by spatial isolation. Our framework also takes into consideration stochastic noise and can therefore capture the inherent variability of epidemiological dynamics. These advances allow for a much wider class of epidemiological models to be fit to genealogies in order to estimate key epidemiological parameters and to reconstruct past disease dynamics.

Global spread of dengue virus types: mapping the 70 year history

Since the first isolation of dengue virus (DENV) in 1943, four types have been identified. Global phenomena such as urbanization and international travel are key factors in facilitating the spread of dengue. Documenting the type-specific record of DENV spread has important implications for understanding patterns in dengue hyperendemicity and disease severity as well as vaccine design and deployment strategies. Existing studies have examined the spread of DENV types at regional or local scales, or described phylogeographic relationships within a single type. Here we summarize the global distribution of confirmed instances of each DENV type from 1943 to 2013 in a series of global maps. These show the worldwide expansion of the types, the expansion of disease hyperendemicity, and the establishment of an increasingly important infectious disease of global public health significance.