Integrating the landscape epidemiology and genetics of RNA viruses: rabies in domestic dogs as a model

Landscape epidemiology and landscape genetics combine advances in molecular techniques, spatial analyses and epidemiological models to generate a more real-world understanding of infectious disease dynamics and provide powerful new tools for the study of RNA viruses. Using dog rabies as a model we have identified how key questions regarding viral spread and persistence can be addressed using a combination of these techniques. In contrast to wildlife rabies, investigations into the landscape epidemiology of domestic dog rabies requires more detailed assessment of the role of humans in disease spread, including the incorporation of anthropogenic landscape features, human movements and socio-cultural factors into spatial models. In particular, identifying and quantifying the influence of anthropogenic features on pathogen spread and measuring the permeability of dispersal barriers are important considerations for planning control strategies, and may differ according to cultural, social and geographical variation across countries or continents. Challenges for dog rabies research include the development of metapopulation models and transmission networks using genetic information to uncover potential source/sink dynamics and identify the main routes of viral dissemination. Information generated from a landscape genetics approach will facilitate spatially strategic control programmes that accommodate for heterogeneities in the landscape and therefore utilise resources in the most cost-effective way. This can include the efficient placement of vaccine barriers, surveillance points and adaptive management for large-scale control programmes

Geography but not alternative host species explain the spread of raccoon rabies virus in Vermont

In North America, the raccoon-associated variant of rabies virus (RRV) is of special concern, given its relatively rapid spread throughout the eastern USA and its potential public health impact due to high raccoon host densities in urban areas. Northward expansion of this epizootic included an outbreak in the Canadian province of Quebec in 2006–2009 due to trans-border spread from the State of Vermont. To inform a more proactive approach to future control efforts, this study uses phylogenetic analyses to explore the role of geography and alternative carnivore hosts in the dynamics of RRV spread within Vermont. Specifically, we sought to examine whether striped skunks, a species frequently infected by RRV, could be part of the maintenance host community. Whole genome sequencing of 160 RRV samples from Vermont and neighbouring US states were used for fine-scale phylogeographic analyses. Results, together with the complete surveillance record of raccoon rabies since its entry into Vermont in 1994, document incursions by two distinct viral lineages and identify topographical features of the landscape which have significantly influenced viral spread, resulting in a complex distribution pattern of viral variants throughout the state. Results of phylogenetic cluster analysis and discrete state reconstruction contained some evidence of skunk-to-skunk and skunk-to-raccoon transmission but overall failed to support a role for skunks as alternative maintenance hosts

Elucidating the phylodynamics of endemic rabies virus in eastern Africa using whole-genome sequencing

Many of the pathogens perceived to pose the greatest risk to humans are viral zoonoses, responsible for a range of emerging and endemic infectious diseases. Phylogeography is a useful tool to understand the processes that give rise to spatial patterns and drive dynamics in virus populations. Increasingly, whole-genome information is being used to uncover these patterns, but the limits of phylogenetic resolution that can be achieved with this are unclear. Here, whole-genome variation was used to uncover fine-scale population structure in endemic canine rabies virus circulating in Tanzania. This is the first whole-genome population study of rabies virus and the first comprehensive phylogenetic analysis of rabies virus in East Africa, providing important insights into rabies transmission in an endemic system. In addition, sub-continental scale patterns of population structure were identified using partial gene data and used to determine population structure at larger spatial scales in Africa. While rabies virus has a defined spatial structure at large scales, increasingly frequent levels of admixture were observed at regional and local levels. Discrete phylogeographic analysis revealed long-distance dispersal within Tanzania, which could be attributed to human-mediated movement, and we found evidence of multiple persistent, co-circulating lineages at a very local scale in a single district, despite on-going mass dog vaccination campaigns. This may reflect the wider endemic circulation of these lineages over several decades alongside increased admixture due to human-mediated introductions. These data indicate that successful rabies control in Tanzania could be established at a national level, since most dispersal appears to be restricted within the confines of country borders but some coordination with neighbouring countries may be required to limit transboundary movements. Evidence of complex patterns of rabies circulation within Tanzania necessitates the use of whole-genome sequencing to delineate finer scale population structure that can that can guide interventions, such as the spatial scale and design of dog vaccination campaigns and dog movement controls to achieve and maintain freedom from disease

Limitations of variable number of tandem repeat typing identified through whole genome sequencing of Mycobacterium avium subsp. paratuberculosis on a national and herd level

Background: Mycobacterium avium subsp. paratuberculosis (MAP), the causative bacterium of Johne’s disease in dairy cattle, is widespread in the Canadian dairy industry and has significant economic and animal welfare implications. An understanding of the population dynamics of MAP can be used to identify introduction events, improve control efforts and target transmission pathways, although this requires an adequate understanding of MAP diversity and distribution between herds and across the country. Whole genome sequencing (WGS) offers a detailed assessment of the SNP-level diversity and genetic relationship of isolates, whereas several molecular typing techniques used to investigate the molecular epidemiology of MAP, such as variable number of tandem repeat (VNTR) typing, target relatively unstable repetitive elements in the genome that may be too unpredictable to draw accurate conclusions. The objective of this study was to evaluate the diversity of bovine MAP isolates in Canadian dairy herds using WGS and then determine if VNTR typing can distinguish truly related and unrelated isolates.<p></p> Results: Phylogenetic analysis based on 3,039 SNPs identified through WGS of 124 MAP isolates identified eight genetically distinct subtypes in dairy herds from seven Canadian provinces, with the dominant type including over 80% of MAP isolates. VNTR typing of 527 MAP isolates identified 12 types, including “bison type” isolates, from seven different herds. At a national level, MAP isolates differed from each other by 1–2 to 239–240 SNPs, regardless of whether they belonged to the same or different VNTR types. A herd-level analysis of MAP isolates demonstrated that VNTR typing may both over-estimate and under-estimate the relatedness of MAP isolates found within a single herd.<p></p> Conclusions: The presence of multiple MAP subtypes in Canada suggests multiple introductions into the country including what has now become one dominant type, an important finding for Johne’s disease control. VNTR typing often failed to identify closely and distantly related isolates, limiting the applicability of using this typing scheme to study the molecular epidemiology of MAP at a national and herd-level.<p></p&gt

Genomic and immunogenic protein diversity of Erysipelothrix rhusiopathiae isolated from pigs in Great Britain: implications for vaccine protection

Erysipelas, caused by the bacterium Erysipelothrix rhusiopathiae, is re-emerging in swine and poultry production systems worldwide. While the global genomic diversity of this species has been characterized, how much of this genomic and functional diversity is maintained at smaller scales is unclear. Specifically, while several key immunogenic surface proteins have been identified for E. rhusiopathiae, little is known about their presence among field strains and their divergence from vaccines, which could result in vaccine failure. Here, a comparative genomics approach was taken to determine the diversity of E. rhusiopathiae strains in pigs in Great Britain over nearly three decades, as well as to assess the field strains’ divergence from the vaccine strain most commonly used in British pigs. In addition, the presence/absence and variability of 13 previously described immunogenic surface proteins was determined, including SpaA which is considered a key immunogen. We found a high diversity of E. rhusiopathiae strains in British pigs, similar to the situation described in European poultry but in contrast to swine production systems in Asia. Of the four clades of E. rhusiopathiae found globally, three were represented among British pig isolates, with Clade 2 being the most common. All British pig isolates had one amino acid difference in the immunoprotective domain of the SpaA protein compared to the vaccine strain. However, we were able to confirm using in silico structural protein analyses that this difference is unlikely to compromise vaccine protection. Of 12 other known immunogenic surface proteins of E. rhusiopathiae examined, 11 were found to be present in all British pig isolates and the vaccine strain, but with highly variable degrees of conservation at the amino acid sequence level, ranging from 0.3 to 27% variant positions. Moreover, the phylogenetic incongruence of these proteins suggests that horizontal transfer of genes encoding for antigens is commonplace for this bacterium. We hypothesize that the sequence variants in these proteins could be responsible for differences in the efficacy of the immune response. Our results provide the necessary basis for testing this hypothesis through in vitro and in vivo studies

Bacterial genomics reveal the complex epidemiology of an emerging pathogen in Arctic and boreal ungulates

Northern ecosystems are currently experiencing unprecedented ecological change, largely driven by a rapidly changing climate. Pathogen range expansion, and emergence and altered patterns of infectious disease, are increasingly reported in wildlife at high latitudes. Understanding the causes and consequences of shifting pathogen diversity and host-pathogen interactions in these ecosystems is important for wildlife conservation, and for indigenous populations that depend on wildlife. Among the key questions are whether disease events are associated with endemic or recently introduced pathogens, and whether emerging strains are spreading throughout the region. In this study, we used a phylogenomic approach to address these questions of pathogen endemicity and spread for Erysipelothrix rhusiopathiae, an opportunistic multi-host bacterial pathogen associated with recent mortalities in arctic and boreal ungulate populations in North America. We isolated E. rhusiopathiae from carcasses associated with large-scale die-offs of muskoxen in the Canadian Arctic Archipelago, and from contemporaneous mortality events and/or population declines among muskoxen in northwestern Alaska and caribou and moose in western Canada. Bacterial genomic diversity differed markedly among these locations; minimal divergence was present among isolates from muskoxen in the Canadian Arctic, while in caribou and moose populations, strains from highly divergent clades were isolated from the same location, or even from within a single carcass. These results indicate that mortalities among northern ungulates are not associated with a single emerging strain of E. rhusiopathiae, and that alternate hypotheses need to be explored. Our study illustrates the value and limitations of bacterial genomic data for discriminating between ecological hypotheses of disease emergence, and highlights the importance of studying emerging pathogens within the broader context of environmental and host factors

Full genome sequence and sfRNA interferon antagonist activity of Zika virus from Recife, Brazil

Background: The outbreak of Zika virus (ZIKV) in the Americas has transformed a previously obscure mosquito-transmitted arbovirus of the Flaviviridae family into a major public health concern. Little is currently known about the evolution and biology of ZIKV and the factors that contribute to the associated pathogenesis. Determining genomic sequences of clinical viral isolates and characterization of elements within these are an important prerequisite to advance our understanding of viral replicative processes and virus-host interactions. Methodology/Principal findings: We obtained a ZIKV isolate from a patient who presented with classical ZIKV-associated symptoms, and used high throughput sequencing and other molecular biology approaches to determine its full genome sequence, including non-coding regions. Genome regions were characterized and compared to the sequences of other isolates where available. Furthermore, we identified a subgenomic flavivirus RNA (sfRNA) in ZIKV-infected cells that has antagonist activity against RIG-I induced type I interferon induction, with a lesser effect on MDA-5 mediated action. Conclusions/Significance: The full-length genome sequence including non-coding regions of a South American ZIKV isolate from a patient with classical symptoms will support efforts to develop genetic tools for this virus. Detection of sfRNA that counteracts interferon responses is likely to be important for further understanding of pathogenesis and virus-host interactions

Practical and effective diagnosis of animal anthrax in endemic low-resource settings

Anthrax threatens human and animal health, and people’s livelihoods in many rural communities in Africa and Asia. In these areas, anthrax surveillance is challenged by a lack of tools for on-site detection. Furthermore, cultural practices and infrastructure may affect sample availability and quality. Practical yet accurate diagnostic solutions are greatly needed to quantify anthrax impacts. We validated microscopic and molecular methods for the detection of Bacillus anthracis in field-collected blood smears and identified alternative samples suitable for anthrax confirmation in the absence of blood smears. We investigated livestock mortalities suspected to be caused by anthrax in northern Tanzania. Field-prepared blood smears (n = 152) were tested by microscopy using four staining techniques as well as polymerase chain reaction (PCR) followed by Bayesian latent class analysis. Median sensitivity (91%, CI 95% [84–96%]) and specificity (99%, CI 95% [96–100%]) of microscopy using azure B were comparable to those of the recommended standard, polychrome methylene blue, PMB (92%, CI 95% [84–97%] and 98%, CI 95% [95–100%], respectively), but azure B is more available and convenient. Other commonly-used stains performed poorly. Blood smears could be obtained for <50% of suspected anthrax cases due to local customs and conditions. However, PCR on DNA extracts from skin, which was almost always available, had high sensitivity and specificity (95%, CI 95% [90–98%] and 95%, CI 95% [87–99%], respectively), even after extended storage at ambient temperature. Azure B microscopy represents an accurate diagnostic test for animal anthrax that can be performed with basic laboratory infrastructure and in the field. When blood smears are unavailable, PCR using skin tissues provides a valuable alternative for confirmation. Our findings lead to a practical diagnostic approach for anthrax in low-resource settings that can support surveillance and control efforts for anthrax-endemic countries globally

Demographic and environmental drivers of metagenomic viral diversity in vampire bats

Indexación: Scopus.Viruses infect all forms of life and play critical roles as agents of disease, drivers of biochemical cycles and sources of genetic diversity for their hosts. Our understanding of viral diversity derives primarily from comparisons among host species, precluding insight into how intraspecific variation in host ecology affects viral communities or how predictable viral communities are across populations. Here we test spatial, demographic and environmental hypotheses explaining viral richness and community composition across populations of common vampire bats, which occur in diverse habitats of North, Central and South America. We demonstrate marked variation in viral communities that was not consistently predicted by a null model of declining community similarity with increasing spatial or genetic distances separating populations. We also find no evidence that larger bat colonies host greater viral diversity. Instead, viral diversity follows an elevational gradient, is enriched by juvenile-biased age structure, and declines with local anthropogenic food resources as measured by livestock density. Our results establish the value of linking the modern influx of metagenomic sequence data with comparative ecology, reveal that snapshot views of viral diversity are unlikely to be representative at the species level, and affirm existing ecological theories that link host ecology not only to single pathogen dynamics but also to viral communities.https://onlinelibrary.wiley.com/doi/10.1111/mec.1525