Phylogeography of foot-and-mouth disease virus serotype O in Ecuador

Ecological and Anthropogenic Spatial Gradients Shape Patterns of Dispersal of Foot-and-Mouth Disease Virus in Uganda

Using georeferenced phylogenetic trees, phylogeography allows researchers to elucidate interactions between environmental heterogeneities and patterns of infectious disease spread. Concordant with the increasing availability of pathogen genetic sequence data, there is a growing need for tools to test epidemiological hypotheses in this field. In this study, we apply tools traditionally used in ecology to elucidate the epidemiology of foot-and-mouth disease virus (FMDV) in Uganda. We analyze FMDV serotype O genetic sequences and their corresponding spatiotemporal metadata from a cross-sectional study of cattle. We apply step selection function (SSF) models, typically used to study wildlife habitat selection, to viral phylogenies to show that FMDV is more likely to be found in areas of low rainfall. Next, we use a novel approach, a resource gradient function (RGF) model, to elucidate characteristics of viral source and sink areas. An RGF model applied to our data reveals that areas of high cattle density and areas near livestock markets may serve as sources of FMDV dissemination in Uganda, and areas of low rainfall serve as viral sinks that experience frequent reintroductions. Our results may help to inform risk-based FMDV control strategies in Uganda. More broadly, these tools advance the phylogenetic toolkit, as they may help to uncover patterns of spread of other organisms for which genetic sequences and corresponding spatiotemporal metadata exist.

Phylogeographical and cross-species transmission dynamics of SAT1 and SAT2 foot-and-mouth disease virus in Eastern Africa

Understanding the dynamics of foot-and-mouth disease virus (FMDV), an endemic and economically constraining disease, is critical in designing control programmes in Africa. This study investigates the evolutionary epidemiology of SAT1 and SAT2 FMDV in Eastern Africa, as well as between cattle and wild African buffalo. Bayesian phylodynamic models were used to analyse SAT1 and SAT2 VP1 gene segments collected between 1975 and 2016, focusing on the SAT1 and SAT2 viruses currently circulating in Eastern Africa. The root state posterior probabilities inferred from our analyses suggest Zimbabwe as the ancestral location for SAT1 currently circulating in Eastern Africa (p = 0.67). For the SAT2 clade, Kenya is inferred to be the ancestral location for introduction of the virus into other countries in Eastern Africa (p = 0.72). Salient (Bayes factor >10) viral dispersal routes were inferred from Tanzania to Kenya, and from Kenya to Uganda for SAT1 and SAT2, respectively. Results suggest that cattle are the source of the SAT1 and SAT2 clades currently circulating in Eastern Africa. In addition, our results suggest that the majority of SAT1 and SAT2 in livestock come from other livestock rather than wildlife, with limited evidence that buffalo serve as reservoirs for cattle. Insights from the present study highlight the role of cattle movements and anthropogenic activities in shaping the evolutionary history of SAT1 and SAT2 in Eastern Africa. While the results may be affected by inherent limitations of imperfect surveillance, our analysis elucidates the dynamics between host species in this region, which is key to guiding disease intervention activities.

Phylodynamics and evolutionary epidemiology of African swine fever p72-CVR genes in Eurasia and Africa

African swine fever (ASF) is a complex infectious disease of swine that constitutes devastating impacts on animal health and the world economy. Here, we investigated the evolutionary epidemiology of ASF virus (ASFV) in Eurasia and Africa using the concatenated gene sequences of the viral protein 72 and the central variable region of isolates collected between 1960 and 2015. We used Bayesian phylodynamic models to reconstruct the evolutionary history of the virus, to identify virus population demographics and to quantify dispersal patterns between host species. Results suggest that ASFV exhibited a significantly high evolutionary rate and population growth through time since its divergence in the 18th century from East Africa, with no signs of decline till recent years. This increase corresponds to the growing pig trade activities between continents during the 19th century, and may be attributed to an evolutionary drift that resulted from either continuous circulation or maintenance of the virus within Africa and Eurasia. Furthermore, results implicate wild suids as the ancestral host species (root state posterior probability = 0.87) for ASFV in the early 1700s in Africa. Moreover, results indicate the transmission cycle between wild suids and pigs is an important cycle for ASFV spread and maintenance in pig populations, while ticks are an important natural reservoir that can facilitate ASFV spread and maintenance in wild swine populations. We illustrated the prospects of phylodynamic methods in improving risk-based surveillance, support of effective animal health policies, and epidemic preparedness in countries at high risk of ASFV incursion.

Phylogeography of the current rabies viruses in Indonesia

Rabies is a major fatal zoonotic disease in Indonesia. This study was conducted to determine the recent dynamics of rabies virus (RABV) in various areas and animal species throughout Indonesia. A total of 27 brain samples collected from rabid animals of various species in Bali, Sumatra, Kalimantan, Sulawesi, Java, and Flores in 2008 to 2010 were investigated. The cDNA of the nucleoprotein gene from each sample was generated and amplified by one-step reverse transcription-PCR, after which the products were sequenced and analyzed. The symmetric substitution model of a Bayesian stochastic search variable selection extension of the discrete phylogeographic model of the social network was applied in BEAST ver. 1.7.5 software. The spatial dispersal was visualized in Cartographica using Spatial Phylogenetic Reconstruction of Evolutionary Dynamics. We demonstrated inter-island introduction and reintroduction, and dog was found to be the only source of infection of other animals. Ancestors of Indonesian RABVs originated in Java and its descendants were transmitted to Kalimantan, then further to Sumatra, Flores, and Bali. The Flores descendent was subsequently transmitted to Sulawesi and back to Kalimantan. The viruses found in various animal species were transmitted by the dog.

Spatiotemporal dynamics of simian immunodeficiency virus brain infection in CD8+ lymphocyte-depleted rhesus macaques with neuroAIDS

Despite the success of combined antiretroviral therapy in controlling viral replication in human immunodeficiency virus (HIV)-infected individuals, HIV-associated neurocognitive disorders, commonly referred to as neuroAIDS, remain a frequent and poorly understood complication. Infection of CD8(+) lymphocyte-depleted rhesus macaques with the SIVmac251 viral swarm is a well-established rapid disease model of neuroAIDS that has provided critical insight into HIV-1-associated neurocognitive disorder onset and progression. However, no studies so far have characterized in depth the relationship between intra-host viral evolution and pathogenesis in this model. Simian immunodeficiency virus (SIV) env gp120 sequences were obtained from six infected animals. Sequences were sampled longitudinally from several lymphoid and non-lymphoid tissues, including individual lobes within the brain at necropsy, for four macaques; two animals were sacrificed at 21 days post-infection (p.i.) to evaluate early viral seeding of the brain. Bayesian phylodynamic and phylogeographic analyses of the sequence data were used to ascertain viral population dynamics and gene flow between peripheral and brain tissues, respectively. A steady increase in viral effective population size, with a peak occurring at ~50-80 days p.i., was observed across all longitudinally monitored macaques. Phylogeographic analysis indicated continual viral seeding of the brain from several peripheral tissues throughout infection, with the last migration event before terminal illness occurring in all macaques from cells within the bone marrow. The results strongly supported the role of infected bone marrow cells in HIV/SIV neuropathogenesis. In addition, our work demonstrated the applicability of Bayesian phylogeography to intra-host studies in order to assess the interplay between viral evolution and pathogenesis.

Spatiotemporal Phylogenetic Analysis and Molecular Characterisation of Infectious Bursal Disease Viruses Based on the VP2 Hyper-Variable Region

Background

Infectious bursal disease is a highly contagious and acute viral disease caused by the infectious bursal disease virus (IBDV); it affects all major poultry producing areas of the world. The current study was designed to rigorously measure the global phylogeographic dynamics of IBDV strains to gain insight into viral population expansion as well as the emergence, spread and pattern of the geographical structure of very virulent IBDV (vvIBDV) strains.

Methodology/Principal Findings

Sequences of the hyper-variable region of the VP2 (HVR-VP2) gene from IBDV strains isolated from diverse geographic locations were obtained from the GenBank database; Cuban sequences were obtained in the current work. All sequences were analysed by Bayesian phylogeographic analysis, implemented in the Bayesian Evolutionary Analysis Sampling Trees (BEAST), Bayesian Tip-association Significance testing (BaTS) and Spatial Phylogenetic Reconstruction of Evolutionary Dynamics (SPREAD) software packages. Selection pressure on the HVR-VP2 was also assessed. The phylogeographic association-trait analysis showed that viruses sampled from individual countries tend to cluster together, suggesting a geographic pattern for IBDV strains. Spatial analysis from this study revealed that strains carrying sequences that were linked to increased virulence of IBDV appeared in Iran in 1981 and spread to Western Europe (Belgium) in 1987, Africa (Egypt) around 1990, East Asia (China and Japan) in 1993, the Caribbean Region (Cuba) by 1995 and South America (Brazil) around 2000. Selection pressure analysis showed that several codons in the HVR-VP2 region were under purifying selection.

Conclusions/Significance

To our knowledge, this work is the first study applying the Bayesian phylogeographic reconstruction approach to analyse the emergence and spread of vvIBDV strains worldwide.