Phylodynamics and human-mediated dispersal of a zoonotic virus

Improved effectiveness of vaccination campaigns against rabies by reducing spatial heterogeneity in coverage

Vaccination programs are the mainstay of control for many infectious diseases. Heterogeneous coverage is hypothesized to reduce vaccination program effectiveness, but this impact has not been quantified in real systems. We address this gap using fine-scale data from two decades of rabies contact tracing and dog vaccination campaigns in Serengeti district, Tanzania. We also aimed to identify drivers of the continued circulation of rabies in the district despite annual vaccination campaigns. Using generalized linear mixed models, we find that current focal (village-level) dog rabies incidence decreases with increasing recent focal vaccination coverage. However, current focal incidence depends most on recent incidence, both focally and in the wider district, consistent with high population connectivity. Removing the masking effects of prior non-focal incidence shows that, for the same average prior non-focal (wider-district) vaccination coverage, increased heterogeneity in coverage among the non-focal villages leads to increased focal incidence. These effects led to outbreaks following years when vaccination campaigns missed many villages, whereas when heterogeneity in coverage was reduced, incidence declined to low levels (<0.4 cases/1,000 dogs annually and no human deaths) and short vaccination lapses thereafter did not lead to resurgence. Through transmission-tree reconstruction, we inferred frequent incursions into the district each year (mean of 7). Inferred incursions substantially increased as a percentage of all cases in recent years, reaching 50% in 2022, suggesting regional connectivity is driving residual transmission. Overall, we empirically demonstrate how population connectivity and spatial heterogeneity in vaccination can impact disease outcomes, highlighting the importance of fine-scale monitoring in managing vaccination programs.

How does date-rounding affect phylodynamic inference for public health?

Phylodynamic analyses infer epidemiological parameters from pathogen genome sequences for enhanced genomic surveillance in public health. Pathogen genome sequences and their associated sampling dates are the essential data in every analysis. However, sampling dates are usually associated with hospitalisation or testing and can sometimes be used to identify individual patients, posing a threat to patient confidentiality. To lower this risk, sampling dates are often given with reduced date-resolution to the month or year, which can potentially bias inference. Here, we introduce a practical guideline on when date-rounding biases the inference of epidemiologically important parameters across a diverse range of empirical and simulated datasets. We show that the direction of bias varies for different parameters, datasets, and tree priors, while compounding with lower date-resolution and higher substitution rates. We also find that bias decreases for datasets with longer sampling intervals, implying that our guideline is most applicable to emerging datasets. We conclude by discussing future solutions that prioritise patient confidentiality and propose a method for safer sharing of sampling dates that translates them them uniformly by a random number.

The evolution and epidemiology of H3N2 canine influenza virus after 20 years in dogs

The H3N2 canine influenza virus (CIV) emerged from an avian reservoir in Asia to circulate entirely among dogs for the last 20 years. The virus was first seen circulating outside Asian dog populations in 2015, in North America. Utilizing viral genomic data in addition to clinical reports and diagnostic testing data, we provide an updated analysis of the evolution and epidemiology of the virus in its canine host. CIV in dogs in North America is marked by a complex life history - including local outbreaks, regional lineage die-outs, and repeated reintroductions of the virus (with diverse genotypes) from different regions of Asia. Phylogenetic and Bayesian analysis reveal multiple CIV clades, and viruses from China have seeded recent North American outbreaks, with 2 or 3 introductions in the past 3 years. Genomic epidemiology confirms that within North America the virus spreads very rapidly among dogs in kennels and shelters in different regions - but then dies out locally. The overall epidemic therefore requires longer-distance dispersal of virus to maintain outbreaks over the long term. With a constant evolutionary rate over 20 years, CIV still appears best adapted to transmission in dense populations and has not gained properties for prolonged circulation among dogs.

Reducing spatial heterogeneity in coverage improves the effectiveness of dog vaccination campaigns against rabies

Vaccination programs are the mainstay of control for many infectious diseases. Heterogeneous coverage is hypothesised to reduce vaccination programme effectiveness, but this impact has not been quantified in real systems. We address this gap using fine-scale data from two decades of rabies contact tracing and dog vaccination campaigns in Serengeti district, Tanzania. We also aimed to identify drivers of continued circulation of rabies in the district despite annual vaccination campaigns. Using generalised linear mixed models, we find that current focal (village-level) dog rabies incidence decreases with increasing recent focal vaccination coverage. However, current focal incidence depends most on recent incidence, both focally and in the wider district, consistent with high population connectivity. Removing the masking effects of prior non-focal incidence shows that, for the same average prior non-focal (wider-district) vaccination coverage, increased heterogeneity in coverage among the non-focal villages leads to increased focal incidence. These effects led to outbreaks following years when vaccination campaigns missed many villages, whereas when heterogeneity in coverage was reduced, incidence declined to low levels (<0.4 cases/1,000 dogs annually and no human deaths) and short vaccination lapses thereafter did not lead to resurgence. Through transmission-tree reconstruction, we inferred frequent incursions into the district each year (mean of 7). Inferred incursions substantially increased as a percentage of all cases in recent years, reaching 50% in 2022, suggesting regional connectivity is driving residual transmission. Overall, we empirically demonstrate how population connectivity and spatial heterogeneity in vaccination can impact disease outcomes, highlighting the importance of fine-scale monitoring in managing vaccination programs.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics that can be informative of the dispersal dynamics and the capacity to spread among hosts. Heterogeneous sampling efforts of genomic sequences can however impact the accuracy of phylogeographic dispersal metrics. While the impact of spatial sampling bias on the outcomes of continuous phylogeographic inference has previously been explored, the impact of sampling intensity (i.e., sampling size) when aiming to characterise dispersal patterns through continuous phylogeographic reconstructions has not yet been thoroughly evaluated. In our study, we use simulations to evaluate the robustness of 3 dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance (IBD) signal metric - to the sampling intensity. Our results reveal that both the diffusion coefficient and IBD signal metrics appear to be the most robust to the number of samples considered for the phylogeographic reconstruction. We then use these 2 dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of IBD patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the use of lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Spatiotemporal dynamics of rabies virus detected in rabid dogs in Cameroon, 2010-2021

Rabies is a viral zoonosis that causes an estimated 60,000 human deaths each year, mainly in Africa and Asia. The etiological agent of rabies, the Rabies Lyssavirus or Rabies Virus (RABV) has been characterized in dog populations in Cameroon, in previous studies. However, the dynamics of RABV maintenance and propagation in dogs are still to be documented in Cameroon. This study thus, aimed at investigating the spatial and temporal dynamics of RABV variants in Cameroon. Long genomic sequences of about 4893 nucleotides, encompassing the N, P, M and G genes as well as part of the G-L intergenic region (Ψ), were determined from 56 RABV strains recovered from dog populations in Cameroon from 2010 to 2021. Temporal and spatial dynamics of RABV circulation in Cameroon were investigated by Bayesian analyses with the BEAST 1.10.4 package from extended RABV genomic sequences data combined with their collection dates and the geographical coordinates of their sampling areas. This revealed a genetic evolution rate of 3.14 × 10-4 substitutions/site/year among Africa-1a and Africa-2 clades of RABV from Cameroon. The most recent common ancestor (MRCA) of the studied strains of the Africa-1a lineage was estimated to have emerged between 1880 and 1906 (95 % HPD; mean 1894), while that of the strains of the Africa-2 clade had a slightly later estimated origin between 1907 and 1928 (95 % HPD, mean 1918). Overall, phylogeographic analyses suggested RABV spread in Cameroon between sub-national regions. Our data provides substantial support to previous findings from similar epidemiological settings, indicating human mediated movements of infected dogs between distant cities may be a key factor in the maintenance of the enzootic cycle of rabies among dogs in Cameroon.

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny. Here, we introduce a family of models-the so-called Phylogenetic Integrated Velocity (PIV) models-that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the United States indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny (Lemey et al., 2010). Here, we introduce a new family of models - the so-called "Phylogenetic Integrated Velocity" (PIV) models - that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the U.S.A. indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Molecular characterisation of human rabies in Tanzania and Kenya: a case series report and phylogenetic investigation

BACKGROUND

Rabies remains a major public health problem in low- and middle-income countries. However, human rabies deaths are rarely laboratory-confirmed or sequenced, especially in Africa. Five human rabies deaths from Tanzania and Kenya were investigated and the causative rabies viruses sequenced, with the aim of identifying implications for rabies control at individual, healthcare and societal levels.

CASE PRESENTATION

The epidemiological context and care of these cases was contrasting. Four had a clear history of being bitten by dogs, while one had an unclear biting history. Two individuals sought medical attention within a day of being bitten, whereas three sought care only after developing rabies symptoms. Despite seeking medical care, none of the cases received complete post-exposure prophylaxis: one patient received only tetanus vaccination, one did not complete the post-exposure vaccination regimen, one followed an off-label vaccination schedule, and two did not receive any post-exposure vaccinations before the onset of symptoms. These cases highlight serious gaps in health-seeking behaviour, and in health systems providing appropriate care following risky exposures, including in the accessibility and effectiveness of post-exposure prophylaxis as it is administered in the region.

CONCLUSIONS

The viral genomic and epidemiological data confirms dog-mediated rabies as the cause of each of these deaths. The phylogenetic investigation highlights the transboundary circulation of rabies within domestic dog populations, revealing distinct rabies virus clades with evidence of regional spread. These findings underscore the importance of coordinated cross-border control efforts between the two countries. Urgent action is needed to improve awareness around the need for emergency post-exposure vaccines that should be accessible in local communities and administered appropriately, as well as investment in coordinated dog vaccination to control dog-mediated rabies, the underlying cause of these deaths.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics allowing to unveil the dispersal dynamics and evaluate the capacity to spread among hosts. Heterogeneous sampling intensity of genomic sequences can however impact the accuracy of dispersal insights gained through phylogeographic inference. In our study, we implement a simulation framework to evaluate the robustness of three dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance signal metric - to the sampling effort. Our results reveal that both the diffusion coefficient and isolation-by-distance signal metrics appear to be robust to the number of samples considered for the phylogeographic reconstruction. We then use these two dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of isolation-by-distance patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Patterns of Animal Rabies Prevalence in Northern South Africa between 1998 and 2022

Rabies is endemic in South Africa and rabies cycles are maintained in both domestic and wildlife species. The significant number of canine rabies cases reported by the World Organization for Animal Health Reference Laboratory for Rabies at Onderstepoort suggests the need for increased research and mass dog vaccinations on specific targeted foci in the country. This study aimed to investigate the spatiotemporal distribution of animal rabies cases from 1998 to 2017 in northern South Africa and environmental factors associated with highly enzootic municipalities. A descriptive analysis was used to investigate temporal patterns. The Getis-Ord Gi statistical tool was used to exhibit low and high clusters. Logistic regression was used to examine the association between the predictor variables and highly enzootic municipalities. A total of 9580 specimens were submitted for rabies diagnosis between 1998 and 2022. The highest positive case rates were from companion animals (1733 cases, 59.71%), followed by livestock (635 cases, 21.88%) and wildlife (621 cases, 21.39%). Rabies cases were reported throughout the year, with the majority occurring in the mid-dry season. Hot spots were frequently in the northern and eastern parts of Limpopo and Mpumalanga. Thicket bush and grassland were associated with rabies between 1998 and 2002. However, between 2008 and 2012, cultivated commercial crops and waterbodies were associated with rabies occurrence. In the last period, plantations and woodlands were associated with animal rabies. Of the total number of municipalities, five consistently and repeatedly had the highest rabies prevalence rates. These findings suggest that authorities should prioritize resources for those municipalities for rabies elimination and management.

Uncovering the endemic circulation of rabies in Cambodia

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Emerging Genomic Trends on Rabies Virus in Davao Region, Philippines, 2018-2021

Rabies, caused by the rabies virus (RABV), remains a significant public health issue in the Philippines despite efforts to control it. To eliminate rabies by 2030, effective surveillance strategies are crucial. In this study, we examined RABV evolution and phylodynamics in the Davao Region using genome sequences from Davao City and nearby provinces. We adapted the RABV ARTIC Protocol for Oxford Nanopore High-Throughput Sequencing to optimize workflow efficiency under limited resources. Comparing new virus samples collected from June 2019 to June 2021 (n = 38) with baseline samples from June 2018 to May 2019 (n = 49), new sub-clades were observed in the phylogenetic tree, suggesting divergence from older variants that were previously undetected. Most of the new viruses belonged to the Asian SEA4_A1.1.1 lineage, but new (SEA4_B1 and SEA4_B1.1) and emerging (SEA4_B1.1_E1) lineages that have never been reported in the Philippines were also identified. The baseline study reported phylogeographic clustering of RABV isolates from the same areas. However, this pattern was disrupted in the current biosurveillance, with variants detected in areas outside the original cluster. Furthermore, our findings revealed significant transmission routes between Davao City and neighboring provinces, contrasting with the predominantly intra-city transmission observed in the baseline study. These results underscore the need for ongoing and timely genomic surveillance to monitor genetic diversity changes and the emergence of novel strains, as well as to track alterations in transmission pathways. Implementing cost-effective next-generation sequencing workflows will facilitate the integration of genomic surveillance into rabies control programs, particularly in resource-limited settings. Collaborations between different sectors can empower local laboratories and experts in genomic technologies and analysis.

Impact and mitigation of sampling bias to determine viral spread: Evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies, which enables reconstruction of the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (Bayesian structured coalescent approximation [BASTA] and marginal approximation of the structured coalescent [MASCOT]). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on the simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for the CTMC model. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for the CTMC model, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing the sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

Spatial and temporal dynamic analysis of rabies: A review of current methodologies

Rabies continues to be one of the deadliest, high risk diseases worldwide, posing a severe threat to public health. The lack of human-to-human transmission means that the spread of rabies is not significantly affected by the distribution of humans or migra- tion. Thus, the spatiotemporal dynamic of cases in both wild and domestic animals is an important issue that can result in human cases. This paper gives an overview of the methodologies for the spatial and temporal dynamic analysis of this disease. It introduces the most representative research progress of spatial aggregation, dynamic transmission, spatiotemporal distribution, epidemiologi- cal analysis and application of modelling in the study of rabies transmission in recent years. This overview should be useful for investigating the spatial and temporal dynamics of rabies, as it could help understanding the spread of cases as well as contribute to the development of better prevention and control strategies in ecology and epidemiology.

Cross-reactive immunity potentially drives global oscillation and opposed alternation patterns of seasonal influenza A viruses

Several human pathogens exhibit distinct patterns of seasonality and circulate as pairs. For instance, influenza A virus subtypes oscillate and peak during winter seasons of the world’s temperate climate zones. Alternation of dominant strains in successive influenza seasons makes epidemic forecasting a major challenge. From the start of the 2009 influenza pandemic we enrolled influenza A virus infected patients (n = 2980) in a global prospective clinical study. Complete hemagglutinin sequences were obtained from 1078 A/H1N1 and 1033 A/H3N2 viruses. We used phylodynamics to construct high resolution spatio-temporal phylogenetic hemagglutinin trees and estimated global influenza A effective reproductive numbers (R) over time (2009–2013). We demonstrate that R oscillates around R = 1 with a clear opposed alternation pattern between phases of the A/H1N1 and A/H3N2 subtypes. Moreover, we find a similar alternation pattern for the number of global viral spread between the sampled geographical locations. Both observations suggest a between-strain competition for susceptible hosts on a global level. Extrinsic factors that affect person-to-person transmission are a major driver of influenza seasonality. The data presented here indicate that cross-reactive host immunity is also a key intrinsic driver of influenza seasonality, which determines the influenza A virus strain at the onset of each epidemic season.

Making genomic surveillance deliver: A lineage classification and nomenclature system to inform rabies elimination

The availability of pathogen sequence data and use of genomic surveillance is rapidly increasing. Genomic tools and classification systems need updating to reflect this. Here, rabies virus is used as an example to showcase the potential value of updated genomic tools to enhance surveillance to better understand epidemiological dynamics and improve disease control. Previous studies have described the evolutionary history of rabies virus, however the resulting taxonomy lacks the definition necessary to identify incursions, lineage turnover and transmission routes at high resolution. Here we propose a lineage classification system based on the dynamic nomenclature used for SARS-CoV-2, defining a lineage by phylogenetic methods for tracking virus spread and comparing sequences across geographic areas. We demonstrate this system through application to the globally distributed Cosmopolitan clade of rabies virus, defining 96 total lineages within the clade, beyond the 22 previously reported. We further show how integration of this tool with a new rabies virus sequence data resource (RABV-GLUE) enables rapid application, for example, highlighting lineage dynamics relevant to control and elimination programmes, such as identifying importations and their sources, as well as areas of persistence and routes of virus movement, including transboundary incursions. This system and the tools developed should be useful for coordinating and targeting control programmes and monitoring progress as countries work towards eliminating dog-mediated rabies, as well as having potential for broader application to the surveillance of other viruses.

Evaluating the Impact of Anthropogenic Factors on the Dissemination of Contemporary Cosmopolitan, Arctic, and Arctic-like Rabies Viruses

Rabies is a globally prevalent viral zoonosis that causes 59,000 deaths per year and has important economic consequences. Most virus spread is associated with the migration of its primary hosts. Anthropogenic dissemination, mainly via the transportation of rabid dogs, shaped virus ecology a few hundred years ago and is responsible for several current outbreaks. A systematic analysis of aberrant long-distance events in the steppe and Arctic-like groups of rabies virus was performed using statistical (Bayesian) phylogeography and plots of genetic vs. geographic distances. The two approaches produced similar results but had some significant differences and complemented each other. No phylogeographic analysis could be performed for the Arctic group because polar foxes transfer the virus across the whole circumpolar region at high velocity, and there was no correlation between genetic and geographic distances in this virus group. In the Arctic-like group and the steppe subgroup of the cosmopolitan group, a significant number of known sequences (15-20%) was associated with rapid long-distance transfers, which mainly occurred within Eurasia. Some of these events have been described previously, while others have not been documented. Most of the recent long-distance transfers apparently did not result in establishing the introduced virus, but a few had important implications for the phylogeographic history of rabies. Thus, human-mediated long-distance transmission of the rabies virus remains a significant threat that needs to be addressed.

Molecular Epidemiology of Rabies in Wild Canidae in Tunisia

Rabies is a viral zoonosis that is transmissible to humans via domestic and wild animals. There are two epidemiological cycles for rabies, the urban and the sylvatic cycles. In an attempt to study the epidemiological role of wild canidae in rabies transmission, the present study aimed to analyze the genetic characteristics of virus isolates and confirm prior suggestions that rabies is maintained through a dog reservoir in Tunisia. Virus strains isolated from wild canidae were subject to viral sequencing, and Bayesian phylogenetic analysis was performed using Beast2 software. Essentially, the virus strains isolated from wild canidae belonged to the Africa-1 clade, which clearly diverges from fox-related strains. Our study also demonstrated that genetic characteristics of the virus isolates were not as distinct as could be expected if a wild reservoir had already existed. On the contrary, the geographic landscape is responsible for the genetic diversity of the virus. The landscape itself could have also acted as a natural barrier to the spread of the virus.

Use of partial N-gene sequences as a tool to monitor progress on rabies control and elimination efforts in Ethiopia

Ethiopia is one of the African countries most affected by rabies. A coarse catalog of rabies viruses (RABV) was created as a benchmark to assess the impact of control and elimination activities. We evaluated a 726 bp amplicon at the end of the N-gene to infer viral lineages in circulation using maximum likelihood and Bayesian methods for phylogenetic reconstruction. We sequenced 228 brain samples from wild and domestic animals collected in five Ethiopian regions during 2010-2017. Results identified co-circulating RABV lineages that are causing recurrent spillover infections into wildlife and domestic animals. We found no evidence of importation of RABVs from other African countries or vaccine-induced cases in the area studied. A divergent RABV lineage might be involved in an independent rabies cycle in jackals. This investigation provides a feasible approach to assess rabies control and elimination efforts in resource-limited countries.

On the Use of Phylogeographic Inference to Infer the Dispersal History of Rabies Virus: A Review Study

Rabies is a neglected zoonotic disease which is caused by negative strand RNA-viruses belonging to the genus Lyssavirus. Within this genus, rabies viruses circulate in a diverse set of mammalian reservoir hosts, is present worldwide, and is almost always fatal in non-vaccinated humans. Approximately 59,000 people are still estimated to die from rabies each year, leading to a global initiative to work towards the goal of zero human deaths from dog-mediated rabies by 2030, requiring scientific efforts from different research fields. The past decade has seen a much increased use of phylogeographic and phylodynamic analyses to study the evolution and spread of rabies virus. We here review published studies in these research areas, making a distinction between the geographic resolution associated with the available sequence data. We pay special attention to environmental factors that these studies found to be relevant to the spread of rabies virus. Importantly, we highlight a knowledge gap in terms of applying these methods when all required data were available but not fully exploited. We conclude with an overview of recent methodological developments that have yet to be applied in phylogeographic and phylodynamic analyses of rabies virus.

Rabies in Costa Rica - Next Steps Towards Controlling Bat-Borne Rabies After its Elimination in Dogs

Rabies is an acute, progressive encephalitis caused by a lyssavirus, with the highest case fatality of any conventional infectious disease. More than 17 different lyssaviruses have been described, but rabies virus is the most widely distributed and important member of the genus. Globally, tens of thousands of human fatalities still occur each year. Although all mammals are susceptible, most human fatalities are caused by the bites of rabid dogs, within lesser developed countries. A global plan envisions the elimination of human rabies cases caused via dogs by the year 2030. The combination of prophylaxis of exposed humans and mass vaccination of dogs is an essential strategy for such success. Regionally, the Americas are well on the way to meet this goal. As one example of achievement, Costa Rica, a small country within Central America, reported the last autochthonous case of human rabies transmitted by a dog at the end of the 1970s. Today, rabies virus transmitted by the common vampire bat, Desmodus rotundus, as well as other wildlife, remains a major concern for humans, livestock, and other animals throughout the region. This review summarizes the historical occurrence of dog rabies and its elimination in Costa Rica, describes the current occurrence of the disease with a particular focus upon affected livestock, discusses the ecology of the vampire bat as the primary reservoir relevant to management, details the clinical characteristics of recent human rabies cases, and provides suggestions for resolution of global challenges posed by this zoonosis within a One Health context.

Mathematical modelling and phylodynamics for the study of dog rabies dynamics and control: A scoping review

BACKGROUND

Rabies is a fatal yet vaccine-preventable disease. In the last two decades, domestic dog populations have been shown to constitute the predominant reservoir of rabies in developing countries, causing 99% of human rabies cases. Despite substantial control efforts, dog rabies is still widely endemic and is spreading across previously rabies-free areas. Developing a detailed understanding of dog rabies dynamics and the impact of vaccination is essential to optimize existing control strategies and developing new ones. In this scoping review, we aimed at disentangling the respective contributions of mathematical models and phylodynamic approaches to advancing the understanding of rabies dynamics and control in domestic dog populations. We also addressed the methodological limitations of both approaches and the remaining issues related to studying rabies spread and how this could be applied to rabies control.

METHODOLOGY/PRINCIPAL FINDINGS

We reviewed how mathematical modelling of disease dynamics and phylodynamics have been developed and used to characterize dog rabies dynamics and control. Through a detailed search of the PubMed, Web of Science, and Scopus databases, we identified a total of n = 59 relevant studies using mathematical models (n = 30), phylodynamic inference (n = 22) and interdisciplinary approaches (n = 7). We found that despite often relying on scarce rabies epidemiological data, mathematical models investigated multiple aspects of rabies dynamics and control. These models confirmed the overwhelming efficacy of massive dog vaccination campaigns in all settings and unraveled the role of dog population structure and frequent introductions in dog rabies maintenance. Phylodynamic approaches successfully disentangled the evolutionary and environmental determinants of rabies dispersal and consistently reported support for the role of reintroduction events and human-mediated transportation over long distances in the maintenance of rabies in endemic areas. Potential biases in data collection still need to be properly accounted for in most of these analyses. Finally, interdisciplinary studies were determined to provide the most comprehensive assessments through hypothesis generation and testing. They also represent new avenues, especially concerning the reconstruction of local transmission chains or clusters through data integration.

CONCLUSIONS/SIGNIFICANCE

Despite advances in rabies knowledge, substantial uncertainty remains regarding the mechanisms of local spread, the role of wildlife in dog rabies maintenance, and the impact of community behavior on the efficacy of control strategies including vaccination of dogs. Future integrative approaches that use phylodynamic analyses and mechanistic models within a single framework could take full advantage of not only viral sequences but also additional epidemiological information as well as dog ecology data to refine our understanding of rabies spread and control. This would represent a significant improvement on past studies and a promising opportunity for canine rabies research in the frame of the One Health concept that aims to achieve better public health outcomes through cross-sector collaboration.

Rabies control in Liberia: Joint efforts towards zero by 30

Despite declaration as a national priority disease, dog rabies remains endemic in Liberia, with surveillance systems and disease control activities still developing. The objective of these initial efforts was to establish animal rabies diagnostics, foster collaboration between all rabies control stakeholders, and develop a short-term action plan with estimated costs for rabies control and elimination in Liberia. Four rabies diagnostic tests, the direct fluorescent antibody (DFA) test, the direct immunohistochemical test (dRIT), the reverse transcriptase polymerase chain reaction (RT-PCR) assay and the rapid immunochromatographic diagnostic test (RIDT), were implemented at the Central Veterinary Laboratory (CVL) in Monrovia between July 2017 and February 2018. Seven samples (n=7) out of eight suspected animals were confirmed positive for rabies lyssavirus, and molecular analyses revealed that all isolates belonged to the Africa 2 lineage, subgroup H. During a comprehensive in-country One Health rabies stakeholder meeting in 2018, a practical workplan, a short-term action plan and an accurately costed mass dog vaccination strategy were developed. Liberia is currently at stage 1.5/5 of the Stepwise Approach towards Rabies Elimination (SARE) tool, which corresponds with countries that are scaling up local-level interventions (e.g. dog vaccination campaigns) to the national level. Overall an estimated 5.3 - 8 million USD invested over 13 years is needed to eliminate rabies in Liberia by 2030. Liberia still has a long road to become free from dog-rabies. However, the dialogue between all relevant stakeholders took place, and disease surveillance considerably improved through implementing rabies diagnosis at the CVL. The joint efforts of diverse national and international stakeholders laid important foundations to achieve the goal of zero dog-mediated human rabies deaths by 2030.

Circuitscape in Julia: Empowering Dynamic Approaches to Connectivity Assessment

The conservation field is experiencing a rapid increase in the amount, variety, and quality of spatial data that can help us understand species movement and landscape connectivity patterns. As interest grows in more dynamic representations of movement potential, modelers are often limited by the capacity of their analytic tools to handle these datasets. Technology developments in software and high-performance computing are rapidly emerging in many fields, but uptake within conservation may lag, as our tools or our choice of computing language can constrain our ability to keep pace. We recently updated Circuitscape, a widely used connectivity analysis tool developed by Brad McRae and Viral Shah, by implementing it in Julia, a high-performance computing language. In this initial re-code (Circuitscape 5.0) and later updates, we improved computational efficiency and parallelism, achieving major speed improvements, and enabling assessments across larger extents or with higher resolution data. Here, we reflect on the benefits to conservation of strengthening collaborations with computer scientists, and extract examples from a collection of 572 Circuitscape applications to illustrate how through a decade of repeated investment in the software, applications have been many, varied, and increasingly dynamic. Beyond empowering continued innovations in dynamic connectivity, we expect that faster run times will play an important role in facilitating co-production of connectivity assessments with stakeholders, increasing the likelihood that connectivity science will be incorporated in land use decisions.

Combining viral genetic and animal mobility network data to unravel peste des petits ruminants transmission dynamics in West Africa

Peste des petits ruminants (PPR) is a deadly viral disease that mainly affects small domestic ruminants. This disease threaten global food security and rural economy but its control is complicated notably because of extensive, poorly monitored animal movements in infected regions. Here we combined the largest PPR virus genetic and animal mobility network data ever collected in a single region to improve our understanding of PPR endemic transmission dynamics in West African countries. Phylogenetic analyses identified the presence of multiple PPRV genetic clades that may be considered as part of different transmission networks evolving in parallel in West Africa. A strong correlation was found between virus genetic distance and network-related distances. Viruses sampled within the same mobility communities are significantly more likely to belong to the same genetic clade. These results provide evidence for the importance of animal mobility in PPR transmission in the region. Some nodes of the network were associated with PPRV sequences belonging to different clades, representing potential "hotspots" for PPR circulation. Our results suggest that combining genetic and mobility network data could help identifying sites that are key for virus entrance and spread in specific areas. Such information could enhance our capacity to develop locally adapted control and surveillance strategies, using among other risk factors, information on animal mobility.

Using geostatistics to better understand the epidemiology of animal rabies in Morocco: what is the contribution of the predictive value?

This study aims to characterize the spatial distribution of animal rabies in Morocco in order to provide appropriate control approaches. Descriptive analyses of the epidemiological data show that the number of reported canine rabies cases greatly underestimates the true incidence of the disease. Underreporting subsequently affects the coherence of its spatial distribution. To perform accurate geographic distribution mapping of the disease based on interpolation methods, a data set was created using data between 2000 and 2018 to compare the derived disease cases with known true values in order to identify disease clusters. The subsequent interpolation was conducted using Ordinary Kriging regression methods and the semi variogram to focus on short distances and reduce uncertainty. The estimated clusters of rabies were evaluated using a cross validation step which revealed predicted cases close to the true values. To improve the precision of analysis, the authors displayed georeferenced dog and human rabies cases reported during the last three years, demonstrating reliable results that correspond to the estimated cluster areas similar to the true disease incidence on the field. This work highlights a strong correlation between infrastructure projects (i.e. railways, roads, facilities) and rabies epizootics for several specific locations. This study is the first attempt to use geostatistics to build upon the understanding of animal rabies in Morocco and shed light on the most appropriate strategies to sustainably reduce and mitigate the risk of rabies. There has been little literature on the use of kriging methods in animal health research. Thus, this study also aimed to explore a novel method in the veterinary sciences to establish kriging as a valid and coherent analysis tool to identify the extent to which the geostatistic area can objectively support understanding on animal rabies and saw it as being highly instrumental in coping with gaps in the data.

Inference of Nipah virus evolution, 1999-2015

Despite near-annual human outbreaks of Nipah virus (NiV) disease in Bangladesh, typically due to individual spillover events from the local bat population, only twenty whole-genome NiV sequences exist from humans and ten from bats. NiV whole-genome sequences from annual outbreaks have been challenging to generate, primarily due to the low viral load in human throat swab and serum specimens. Here, we used targeted enrichment with custom NiV-specific probes and generated thirty-five additional unique full-length genomic sequences directly from human specimens and viral isolates. We inferred the temporal and geographic evolutionary history of NiV in Bangladesh and expanded a tool to visualize NiV spatio-temporal spread from a Bayesian continuous diffusion analysis. We observed that strains from Bangladesh segregated into two distinct clades that have intermingled geographically in Bangladesh over time and space. As these clades expanded geographically and temporally, we did not observe evidence for significant branch and site-specific selection, except for a single site in the Henipavirus L polymerase. However, the Bangladesh 1 and 2 clades are differentiated by mutations initially occurring in the polymerase, with additional mutations accumulating in the N, G, F, P, and L genes on external branches. Modeling the historic geographical and temporal spread demonstrates that while widespread, NiV does not exhibit significant genetic variation in Bangladesh. Thus, future public health measures should address whether NiV within in the bat population also exhibits comparable genetic variation, if zoonotic transmission results in a genetic bottleneck and if surveillance techniques are detecting only a subset of NiV.

[New aspects of rabies control]

Rabies still causes about 60,000 human deaths per year, mainly in poor populations in Africa and Asia. However, since Louis Pasteur developed the first vaccine 130 years ago, prophylactic measures have been considerably improved and simplified. They now consist of the vaccine combined with purified rabies immunoglobulins of equine or human origin. In general, however, post-exposure prophylaxis protocols are long and expensive. Furthermore, the immunoglobulins used for associated serotherapy are costly and not widely available in developing countries. Approaches have been developed to deal with these two issues that offer hope for a paradigm shift for the benefit of exposed populations. Finally, mass rabies vaccination in dogs, which are the most cost-effective measure for preventing rabies in humans, are difficult to implement and sometimes have moderate effectiveness. The identification and analysis of the epidemiological drivers conditioning the circulation of the virus in dog populations allow a better understanding of the key control points that need to be associated with these campaigns for a better efficacy.

Spatio-temporal evolution of canine rabies in Tunisia, 2011-2016

Tunisia is an endemic country for dog mediated rabies. An increase in canine rabies cases during the last decade has been suspected. Since no studies have been conducted on rabies spatial distribution, the present work was focused on spatiotemporal evolution of rabies in Tunisia during the 2011-2016 period with a special focus on the reservoir species. Data collected concerned suspected dogs that originate from the whole country. Surveillance indicators such as positive fractions and number of suspected dogs received at the laboratory have been calculated. Spatiotemporal hotspots were then mapped, spatial and spatio-temporal analysis were carried out using discrete Poisson spatial model and space-time permutation models available in SaTScan9 software. The study revealed that an actual increase in canine rabies incidence occurred in Tunisia since 2012. Spatial and spatio-temporal analysis identified clusters centered in the North and in the Center East of the country. Spatio-temporal clusters were non overlapping, indicating that this spatial distribution is not fixed through time. A large heterogeneity in surveillance indicators such as number of suspected dogs was associated to the distance to the laboratory or to insufficient coordination between governorates.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Long-range movements coupled with heterogeneous incubation period sustain dog rabies at the national scale in Africa

Dog-transmitted rabies is responsible for more than 98% of human cases worldwide, remaining a persistent problem in developing countries. Mass vaccination targets predominantly major cities, often compromising disease control due to re-introductions. Previous work suggested that areas neighboring cities may behave as the source of these re-introductions. To evaluate this hypothesis, we introduce a spatially explicit metapopulation model for rabies diffusion in Central African Republic. Calibrated on epidemiological data for the capital city, Bangui, the model predicts that long-range movements are essential for continuous re-introductions of rabies-exposed dogs across settlements, eased by the large fluctuations of the incubation period. Bangui's neighborhood, instead, would not be enough to self-sustain the epidemic, contrary to previous expectations. Our findings suggest that restricting long-range travels may be very efficient in limiting rabies persistence in a large and fragmented dog population. Our framework can be applied to other geographical contexts where dog rabies is endemic.

Movement patterns of free-roaming dogs on heterogeneous urban landscapes: Implications for rabies control

In 2015, a case of canine rabies in Arequipa, Peru indicated the re-emergence of rabies virus in the city. Despite mass dog vaccination campaigns across the city and reactive ring vaccination and other control activities around positive cases (e.g. elimination of unowned dogs), the outbreak has spread. Here we explore how the urban landscape of Arequipa affects the movement patterns of free-roaming dogs, the main reservoirs of the rabies virus in the area. We tracked 23 free-roaming dogs using Global Positioning System (GPS) collars. We analyzed the spatio-temporal GPS data using the time- local convex hull method. Dog movement patterns varied across local environments. We found that water channels, an urban feature of Arequipa that are dry most of the year, promote movement. Dogs that used the water channels extensively move on average 7 times further (p = 0.002) and 1.2 times more directionally (p = 0.027) than dogs that do not use the water channels at all. They were also 1.3 times faster on average, but this difference was not statistically significant (p = 0.197). Our findings suggest that water channels can be used by dogs as 'highways' to transverse the city and have the potential to spread disease far beyond the radius of control practices. Control efforts should focus on a robust vaccination campaign attuned to the geography of the city, and not limited to small-scale rings surrounding cases.

In Search of Covariates of HIV-1 Subtype B Spread in the United States-A Cautionary Tale of Large-Scale Bayesian Phylogeography

Infections with HIV-1 group M subtype B viruses account for the majority of the HIV epidemic in the Western world. Phylogeographic studies have placed the introduction of subtype B in the United States in New York around 1970, where it grew into a major source of spread. Currently, it is estimated that over one million people are living with HIV in the US and that most are infected with subtype B variants. Here, we aim to identify the drivers of HIV-1 subtype B dispersal in the United States by analyzing a collection of 23,588 pol sequences, collected for drug resistance testing from 45 states during 2004-2011. To this end, we introduce a workflow to reduce this large collection of data to more computationally-manageable sample sizes and apply the BEAST framework to test which covariates associate with the spread of HIV-1 across state borders. Our results show that we are able to consistently identify certain predictors of spread under reasonable run times across datasets of up to 10,000 sequences. However, the general lack of phylogenetic structure and the high uncertainty associated with HIV trees make it difficult to interpret the epidemiological relevance of the drivers of spread we are able to identify. While the workflow we present here could be applied to other virus datasets of a similar scale, the characteristic star-like shape of HIV-1 phylogenies poses a serious obstacle to reconstructing a detailed evolutionary and spatial history for HIV-1 subtype B in the US.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Distinct rates and patterns of spread of the major HIV-1 subtypes in Central and East Africa

Since the ignition of the HIV-1 group M pandemic in the beginning of the 20th century, group M lineages have spread heterogeneously throughout the world. Subtype C spread rapidly through sub-Saharan Africa and is currently the dominant HIV lineage worldwide. Yet the epidemiological and evolutionary circumstances that contributed to its epidemiological expansion remain poorly understood. Here, we analyse 346 novel pol sequences from the DRC to compare the evolutionary dynamics of the main HIV-1 lineages, subtypes A1, C and D. Our results place the origins of subtype C in the 1950s in Mbuji-Mayi, the mining city of southern DRC, while subtypes A1 and D emerged in the capital city of Kinshasa, and subtypes H and J in the less accessible port city of Matadi. Following a 15-year period of local transmission in southern DRC, we find that subtype C spread at least three-fold faster than other subtypes circulating in Central and East Africa. In conclusion, our results shed light on the origins of HIV-1 main lineages and suggest that socio-historical rather than evolutionary factors may have determined the epidemiological fate of subtype C in sub-Saharan Africa.

Using phylogeographic approaches to analyse the dispersal history, velocity and direction of viral lineages - Application to rabies virus spread in Iran

Recent years have seen the extensive use of phylogeographic approaches to unveil the dispersal history of virus epidemics. Spatially explicit reconstructions of viral spread represent valuable sources of lineage movement data that can be exploited to investigate the impact of underlying environmental layers on the dispersal of pathogens. Here, we performed phylogeographic inference and applied different post hoc approaches to analyse a new and comprehensive data set of viral genomes to elucidate the dispersal history and dynamics of rabies virus (RABV) in Iran, which have remained largely unknown. We first analysed the association between environmental factors and variations in dispersal velocity among lineages. Second, we present, test and apply a new approach to study the link between environmental conditions and the dispersal direction of lineages. The statistical performance (power of detection, false-positive rate) of this new method was assessed using simulations. We performed phylogeographic analyses of RABV genomes, allowing us to describe the large diversity of RABV in Iran and to confirm the cocirculation of several clades in the country. Overall, we estimate a relatively high lineage dispersal velocity, similar to previous estimates for dog rabies virus spread in northern Africa. Finally, we highlight a tendency for RABV lineages to spread in accessible areas associated with high human population density. Our analytical workflow illustrates how phylogeographic approaches can be used to investigate the impact of environmental factors on several aspects of viral dispersal dynamics.

High-resolution contact networks of free-ranging domestic dogs Canis familiaris and implications for transmission of infection

Contact patterns strongly influence the dynamics of disease transmission in both human and non-human animal populations. Domestic dogs Canis familiaris are a social species and are a reservoir for several zoonotic infections, yet few studies have empirically determined contact patterns within dog populations. Using high-resolution proximity logging technology, we characterised the contact networks of free-ranging domestic dogs from two settlements (n = 108 dogs, covering >80% of the population in each settlement) in rural Chad. We used these data to simulate the transmission of an infection comparable to rabies and investigated the effects of including observed contact heterogeneities on epidemic outcomes. We found that dog contact networks displayed considerable heterogeneity, particularly in the duration of contacts and that the network had communities that were highly correlated with household membership. Simulations using observed contact networks had smaller epidemic sizes than those that assumed random mixing, demonstrating the unsuitability of homogenous mixing models in predicting epidemic outcomes. When contact heterogeneities were included in simulations, the network position of the individual initially infected had an important effect on epidemic outcomes. The risk of an epidemic occurring was best predicted by the initially infected individual's ranked degree, while epidemic size was best predicted by the individual's ranked eigenvector centrality. For dogs in one settlement, we found that ranked eigenvector centrality was correlated with range size. Our results demonstrate that observed heterogeneities in contacts are important for the prediction of epidemiological outcomes in free-ranging domestic dogs. We show that individuals presenting a higher risk for disease transmission can be identified by their network position and provide evidence that observable traits hold potential for informing targeted disease management strategies.

Transmission ecology of canine parvovirus in a multi-host, multi-pathogen system

Understanding multi-host pathogen maintenance and transmission dynamics is critical for disease control. However, transmission dynamics remain enigmatic largely because they are difficult to observe directly, particularly in wildlife. Here, we investigate the transmission dynamics of canine parvovirus (CPV) using state-space modelling of 20 years of CPV serology data from domestic dogs and African lions in the Serengeti ecosystem. We show that, although vaccination reduces the probability of infection in dogs, and despite indirect enhancement of population seropositivity as a result of vaccine shedding, the vaccination coverage achieved has been insufficient to prevent CPV from becoming widespread. CPV is maintained by the dog population and has become endemic with approximately 3.5-year cycles and prevalence reaching approximately 80%. While the estimated prevalence in lions is lower, peaks of infection consistently follow those in dogs. Dogs exposed to CPV are also more likely to become infected with a second multi-host pathogen, canine distemper virus. However, vaccination can weaken this coupling, raising questions about the value of monovalent versus polyvalent vaccines against these two pathogens. Our findings highlight the need to consider both pathogen- and host-level community interactions when seeking to understand the dynamics of multi-host pathogens and their implications for conservation, disease surveillance and control programmes.

Transmission dynamics of re-emerging rabies in domestic dogs of rural China

Despite ongoing efforts to control transmission, rabies prevention remains a challenge in many developing countries, especially in rural areas of China where re-emerging rabies is under-reported due to a lack of sustained animal surveillance. By taking advantage of detailed genomic and epidemiological data for the re-emerging rabies outbreak in Yunnan Province, China, collected between 1999 and 2015, we reconstruct the demographic and dispersal history of domestic dog rabies virus (RABV) as well as the dynamics of dog-to-dog and dog-to-human transmission. Phylogeographic analyses reveal a lower diffusion coefficient than previously estimated for dog RABV dissemination in northern Africa. Furthermore, epidemiological analyses reveal transmission rates between dogs, as well as between dogs and humans, lower than estimates for Africa. Finally, we show that reconstructed epidemic history of RABV among dogs and the dynamics of rabid dogs are consistent with the recorded human rabies cases. This work illustrates the benefits of combining phylogeographic and epidemic modelling approaches for uncovering the spatiotemporal dynamics of zoonotic diseases, with both approaches providing estimates of key epidemiological parameters.

Demographic studies of owned dogs in the Northern Peninsula Area, Australia, to inform population and disease management strategies

OBJECTIVE

To generate domestic dog demographic information to aid population and disease management in the Aboriginal and Torres Strait Islander communities of the Northern Peninsula Area, Queensland, Australia.

METHODS

Sight-resight surveys using standard and modified methods were conducted to estimate the free-roaming dog population size. A cross-sectional questionnaire of dog owners was used to gather dog demographic information and investigate owners' dog management behaviours. A survey was also conducted to estimate the total dog population size.

RESULTS

The mean total dog population size was estimated to be 813 (range, 770-868). The roaming dog population was 430 or 542 (95% confidence interval (CI) 254-608; 95% CI 405-680, standard and modified methods, respectively). Therefore, the roaming population represents 52.8% or 66.7% of the total population based on the sight-resight methodology. We surveyed 65 dog owners who owned 165 dogs (1 : 1 ratio of male : female dogs). Only 14% (95% CI 9-19) of dogs were sterilised and significantly more males were entire (P = 0.02). Although most dogs were pets (65%), hunting dogs were significantly more likely to be taken outside of the resident community (P < 0.001). The birth rate was 2.4 puppies/dog-owning house/year, which was higher than the death rate (1.7 dogs/dog-owning house/year). In the previous 12 months, 90% of the 109 deaths were dogs aged 0-2 years old.

CONCLUSION

This study demonstrated that most of the dog population in the NPA is free-roaming and that the population has increased, likely because of a lack of population management strategies such as sterilisation. This information will be used to develop population and disease management strategies in the NPA.

Vaccination of dogs in an African city interrupts rabies transmission and reduces human exposure

Despite the existence of effective rabies vaccines for dogs, dog-transmitted human rabies persists and has reemerged in Africa. Two consecutive dog vaccination campaigns took place in Chad in 2012 and 2013 (coverage of 71% in both years) in the capital city of N'Djaména, as previously published. We developed a deterministic model of dog-human rabies transmission fitted to weekly incidence data of rabid dogs and exposed human cases in N'Djaména. Our analysis showed that the effective reproductive number, that is, the number of new dogs infected by a rabid dog, fell to below one through November 2014. The modeled incidence of human rabies exposure fell to less than one person per million people per year. A phylodynamic estimation of the effective reproductive number from 29 canine rabies virus genetic sequences of the viral N-protein confirmed the results of the deterministic transmission model, implying that rabies transmission between dogs was interrupted for 9 months. However, new dog rabies cases appeared earlier than the transmission and phylodynamic models predicted. This may have been due to the continuous movement of rabies-exposed dogs into N'Djaména from outside the city. Our results show that canine rabies transmission to humans can be interrupted in an African city with currently available dog rabies vaccines, provided that the vaccination area includes larger adjacent regions, and local communities are informed and engaged.

Rabies elimination research: juxtaposing optimism, pragmatism and realism

More than 100 years of research has now been conducted into the prevention, control and elimination of rabies with safe and highly efficacious vaccines developed for use in human and animal populations. Domestic dogs are a major reservoir for rabies, and although considerable advances have been made towards the elimination and control of canine rabies in many parts of the world, the disease continues to kill tens of thousands of people every year in Africa and Asia. Policy efforts are now being directed towards a global target of zero human deaths from dog-mediated rabies by 2030 and the global elimination of canine rabies. Here we demonstrate how research provides a cause for optimism as to the feasibility of these goals through strategies based around mass dog vaccination. We summarize some of the pragmatic insights generated from rabies epidemiology and dog ecology research that can improve the design of dog vaccination strategies in low- and middle-income countries and which should encourage implementation without further delay. We also highlight the need for realism in reaching the feasible, although technically more difficult and longer-term goal of global elimination of canine rabies. Finally, we discuss how research on rabies has broader relevance to the control and elimination of a suite of diseases of current concern to human and animal health, providing an exemplar of the value of a ‘One Health’ approach.

The spread and evolution of rabies virus: conquering new frontiers

Rabies is a lethal zoonotic disease that is caused by lyssaviruses, most often rabies virus. Despite control efforts, sporadic outbreaks in wildlife populations are largely unpredictable, underscoring our incomplete knowledge of what governs viral transmission and spread in reservoir hosts. Furthermore, the evolutionary history of rabies virus and related lyssaviruses remains largely unclear. Robust surveillance efforts combined with diagnostics and disease modelling are now providing insights into the epidemiology and evolution of rabies virus. The immune status of the host, the nature of exposure and strain differences all clearly influence infection and transmission dynamics. In this Review, we focus on rabies virus infections in the wildlife and synthesize current knowledge in the rapidly advancing fields of rabies virus epidemiology and evolution, and advocate for multidisciplinary approaches to advance our understanding of this disease.

Human-Mediated Gene Flow Contributes to Metapopulation Genetic Structure of the Pathogenic Fungus Alternaria alternata from Potato

Metapopulation structure generated by recurrent extinctions and recolonizations plays an important role in the evolution of species but is rarely considered in agricultural systems. In this study, generation and mechanism of metapopulation structure were investigated by microsatellite assaying 725 isolates of Alternaria alternata sampled from potato hosts at 16 locations across China. We found a single major cluster, no isolate-geography associations and no bottlenecks in the A. alternata isolates, suggesting a metapopulation genetic structure of the pathogen. We also found weak isolation-by-distance, lower among than within cropping region population differentiation, concordant moving directions of potato products and net gene flow and the highest gene diversity in the region with the most potato imports. These results indicate that in addition to natural dispersal, human-mediated gene flow also contributes to the generation and dynamics of the metapopulation genetic structure of A. alternata in China. Metapopulation structure increases the adaptive capacity of the plant pathogen as a result of enhanced genetic variation and reduced population fragmentation. Consequently, rigid quarantine regulations may be required to reduce population connectivity and the evolutionary potential of A. alternata and other pathogens with a similar population dynamics for a sustainable plant disease management.

Landscape attributes governing local transmission of an endemic zoonosis: Rabies virus in domestic dogs

Landscape heterogeneity plays an important role in disease spread and persistence, but quantifying landscape influences and their scale dependence is challenging. Studies have focused on how environmental features or global transport networks influence pathogen invasion and spread, but their influence on local transmission dynamics that underpin the persistence of endemic diseases remains unexplored. Bayesian phylogeographic frameworks that incorporate spatial heterogeneities are promising tools for analysing linked epidemiological, environmental and genetic data. Here, we extend these methodological approaches to decipher the relative contribution and scale-dependent effects of landscape influences on the transmission of endemic rabies virus in Serengeti district, Tanzania (area ~4,900 km2 ). Utilizing detailed epidemiological data and 152 complete viral genomes collected between 2004 and 2013, we show that the localized presence of dogs but not their density is the most important determinant of diffusion, implying that culling will be ineffective for rabies control. Rivers and roads acted as barriers and facilitators to viral spread, respectively, and vaccination impeded diffusion despite variable annual coverage. Notably, we found that landscape effects were scale-dependent: rivers were barriers and roads facilitators on larger scales, whereas the distribution of dogs was important for rabies dispersal across multiple scales. This nuanced understanding of the spatial processes that underpin rabies transmission can be exploited for targeted control at the scale where it will have the greatest impact. Moreover, this research demonstrates how current phylogeographic frameworks can be adapted to improve our understanding of endemic disease dynamics at different spatial scales.

Defining objective clusters for rabies virus sequences using affinity propagation clustering

Rabies is caused by lyssaviruses, and is one of the oldest known zoonoses. In recent years, more than 21,000 nucleotide sequences of rabies viruses (RABV), from the prototype species rabies lyssavirus, have been deposited in public databases. Subsequent phylogenetic analyses in combination with metadata suggest geographic distributions of RABV. However, these analyses somewhat experience technical difficulties in defining verifiable criteria for cluster allocations in phylogenetic trees inviting for a more rational approach. Therefore, we applied a relatively new mathematical clustering algorythm named ‘affinity propagation clustering’ (AP) to propose a standardized sub-species classification utilizing full-genome RABV sequences. Because AP has the advantage that it is computationally fast and works for any meaningful measure of similarity between data samples, it has previously been applied successfully in bioinformatics, for analysis of microarray and gene expression data, however, cluster analysis of sequences is still in its infancy. Existing (516) and original (46) full genome RABV sequences were used to demonstrate the application of AP for RABV clustering. On a global scale, AP proposed four clusters, i.e. New World cluster, Arctic/Arctic-like, Cosmopolitan, and Asian as previously assigned by phylogenetic studies. By combining AP with established phylogenetic analyses, it is possible to resolve phylogenetic relationships between verifiably determined clusters and sequences. This workflow will be useful in confirming cluster distributions in a uniform transparent manner, not only for RABV, but also for other comparative sequence analyses.

Rabies is one of the oldest known zoonoses, caused by lyssaviruses. In recent years, more than 21,000 nucleotide sequences for rabies viruses (RABV) have been deposited in public databases. In this study, a novel mathematical approach called affinity propagation (AP) clustering, a highly powerful tool, to verifiably divide full genome RABV sequences into genetic clusters, was used. A panel of existing and novel RABV full genome sequences was used to demonstrate the application of AP for RABV clustering. Using a combination of AP with established phylogenetic analyses is useful in resolving phylogenetic relationships between more objectively determined clusters and sequences. This workflow will help to substantiate a transparent cluster distribution, not only for RABV, but also for other comparative sequence analyses.

Using Viral Gene Sequences to Compare and Explain the Heterogeneous Spatial Dynamics of Virus Epidemics

Rabies is an important zoonotic disease distributed worldwide. A key question in rabies epidemiology is the identification of factors that impact virus dispersion. Here we apply new analytical methods, based on phylogeographic reconstructions of viral lineage movement, to undertake a comparative evolutionary-epidemiological study of the spatial dynamics of rabies virus (RABV) epidemics in different hosts and habitats. We compiled RABV data sets from skunk, raccoon, bat and domestic dog populations in order to investigate the viral diffusivity of different RABV epidemics, and to detect and compare the environmental factors that impact the velocity of viral spread in continuous spatial landscapes. We build on a recently developed statistical framework that uses spatially- and temporally-referenced phylogenies. We estimate several spatial statistics of virus spread, which reveal a higher diffusivity of RABV in domestic dogs compared with RABV in other mammals. This finding is explained by subsequent analyses of environmental heterogeneity, which indicate that factors relating to human geography play a significant role in RABV dispersion in domestic dogs. More generally, our results suggest that human-related factors are important worldwide in explaining RABV dispersion in terrestrial host species. Our study shows that phylogenetically informed viral movements can be used to elucidate the factors that impact virus dispersal, opening new opportunities for a better understanding of the impact of host species and environmental conditions on the spatial dynamics of rapidly evolving populations.

Cross-border spread, lineage displacement and evolutionary rate estimation of rabies virus in Yunnan Province, China

Background

Rabies is an important but underestimated threat to public health, with most cases reported in Asia. Since 2000, a new epidemic wave of rabies has emerged in Yunnan Province, southwestern China, which borders three countries in Southeast Asia.

Method

We estimated gene-specific evolutionary rates for rabies virus using available data in GenBank, then used this information to calibrate the timescale of rabies virus (RABV) spread in Asia. We used 452 publicly available geo-referenced complete nucleoprotein (N) gene sequences, including 52 RABV sequences that were recently generated from samples collected in Yunnan between 2008 and 2012.

Results

The RABV N gene evolutionary rate was estimated to be 1.88 × 10⁻⁴ (1.37–2.41 × 10⁻⁴, 95% Bayesian credible interval, BCI) substitutions per site per year. Phylogenetic reconstructions show that the currently circulating RABV lineages in Yunnan result from at least seven independent introductions (95% BCI: 6–9 introductions) and represent each of the three main Asian RABV lineages, SEA-1, -2 and -3. We find that Yunnan is a sink location for the domestic spread of RABV and connects RABV epidemics in North China, South China, and Southeast Asia. Cross-border spread from southeast Asia (SEA) into South China, and intermixing of the North and South China epidemics is also well supported. The influx of RABV into Yunnan from SEA was not well-supported, likely due to the poor sampling of SEA RABV diversity. We found evidence for a lineage displacement of the Yunnan SEA-2 and -3 lineages by Yunnan SEA-1 strains, and considered whether this could be attributed to fitness differences.

Conclusion

Overall, our study contributes to a better understanding of the spread of RABV that could facilitate future rabies virus control and prevention efforts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0769-6) contains supplementary material, which is available to authorized users.

A Century Spent Combating Rabies in Morocco (1911-2015): How Much Longer?

Rabies has no known beginning in Morocco and to date, government control efforts and plans fail to eradicate the disease. A review and analysis of available epidemiological data are crucial to learn lessons from the past and to propose effective actions. Legally, animal rabies is a notifiable disease since 1913 and legislation has been updated periodically since. Dogs have always been considered as both the disease's vector and reservoir, while cattle, other herbivores, and humans are victims. Animal rabies cases evolution from 1942 to 2015 is characterized by ascending phase then decreasing one following structured rabies control plan implementation in 1980s. Indeed, from 1986 to 2010, three rabies control plans have been conducted based on free of charge rabies vaccination of owned dogs through mass campaigns. The geographical distribution of rabies is stable over the years with highest cases number in rich rural areas and around cities. Human rabies cases are decreasing over the time (1976-2015) thanks to the opening of new antirabic treatment centers in the last decade which permit the administration of more PEPs. After a century of rabies control, Morocco registered an average of 301 animal cases and 21 human cases annually for the last decade (2005-2015). Few reasons led to those limited results. The lack in law enforcement and, moreover, the fact that the law do not take into account responsible dog ownership aspect are of importance. Lack of dog population knowledge and management and intersectoral coordination deficiency are additional failure reasons. The gathered data will help to build a new strategy with a focus on a "One Health" approach. Dog population ecology parameters' study is of primary importance. We estimated dog population to be 2.8 million dogs based on human:dog ratio. Enhancing vaccination coverage of dog population is feasible by combining parenteral vaccination and complementary oral vaccination. Updating legislation by inclusion of responsible dog ownership and law enforcement are crucial. Over the last century, Morocco registered a slow decreasing tendency in the number of animal and human rabies cases. Urgent strategy need to be implemented because rabies elimination is an achievable goal in Morocco.