A Bayesian approach for inferring the dynamics of partially observed endemic infectious diseases from space-time-genetic data

We describe a statistical framework for reconstructing the sequence of transmission events between observed cases of an endemic infectious disease using genetic, temporal and spatial information. Previous approaches to reconstructing transmission trees have assumed all infections in the study area originated from a single introduction and that a large fraction of cases were observed. There are as yet no approaches appropriate for endemic situations in which a disease is already well established in a host population and in which there may be multiple origins of infection, or that can enumerate unobserved infections missing from the sample. Our proposed framework addresses these shortcomings, enabling reconstruction of partially observed transmission trees and estimating the number of cases missing from the sample. Analyses of simulated datasets show the method to be accurate in identifying direct transmissions, while introductions and transmissions via one or more unsampled intermediate cases could be identified at high to moderate levels of case detection. When applied to partial genome sequences of rabies virus sampled from an endemic region of South Africa, our method reveals several distinct transmission cycles with little contact between them, and direct transmission over long distances suggesting significant anthropogenic influence in the movement of infected dogs.     

Population structure of two rabies hosts relative to the known distribution of rabies virus variants in Alaska

For pathogens that infect multiple species, the distinction between reservoir hosts and spillover hosts is often difficult. In Alaska, three variants of the arctic rabies virus exist with distinct spatial distributions. We tested the hypothesis that rabies virus variant distribution corresponds to the population structure of the primary rabies hosts in Alaska, arctic foxes (Vulpes lagopus) and red foxes (Vulpes vulpes) to possibly distinguish reservoir and spillover hosts. We used mitochondrial DNA (mtDNA) sequence and nine microsatellites to assess population structure in those two species. mtDNA structure did not correspond to rabies virus variant structure in either species. Microsatellite analyses gave varying results. Bayesian clustering found two groups of arctic foxes in the coastal tundra region, but for red foxes it identified tundra and boreal types. Spatial Bayesian clustering and spatial principal components analysis identified 3 and 4 groups of arctic foxes, respectively, closely matching the distribution of rabies virus variants in the state. Red foxes, conversely, showed eight clusters comprising two regions (boreal and tundra) with much admixture. These results run contrary to previous beliefs that arctic fox show no fine‐scale spatial population structure. While we cannot rule out that the red fox is part of the maintenance host community for rabies in Alaska, the distribution of virus variants appears to be driven primarily by the arctic fox. Therefore, we show that host population genetics can be utilized to distinguish between maintenance and spillover hosts when used in conjunction with other approaches.     

Rethinking Human–Nonhuman Primate Contact and Pathogenic Disease Spillover

Zoonotic transmissions are a major global health risk, and human–animal contact is frequently raised as an important driver of transmission. A literature examining zooanthroponosis largely agrees that more human–animal contact leads to more risk. Yet the basis of this proposition, the term contact, has not been rigorously analyzed. To understand how contact is used to explain cross-species spillovers, we conducted a multi-disciplinary review of studies addressing human–nonhuman primate (NHP) engagements and pathogenic transmissions and employing the term contact. We find that although contact is frequently invoked, it is employed inconsistently and imprecisely across these studies, overlooking the range of pathogens and their transmission routes and directions. We also examine a related but more expansive approach focusing on human and NHP habitats and their spatial overlap, which can potentially facilitate pathogenic transmission. Contact and spatial overlap investigations cannot, however, explain the processes that bring together people, animals and pathogens. We therefore examine another approach that enhances our understanding of zoonotic spillovers: anthropological studies identifying such historical, social, environmental processes. Comparable to a One Health approach, our ongoing research in Cameroon draws contact, spatial overlap and anthropological–historical approaches into dialog to suggest where, when and how pathogenic transmissions between people and NHPs may occur. In conclusion, we call for zoonotic disease researchers to specify more precisely the human–animal contacts they investigate and to attend to how broader ecologies, societies and histories shape pathogen–human–animal interactions.     

Ecological Drivers of Virus Evolution: Astrovirus as a Case Study

Although RNA viruses exhibit a high frequency of host jumps, major differences exist among the different virus families. Astroviruses infect a wide range of hosts, affecting both public health systems and economic production chains. Here we delineate the ecological and adaptive processes that drive the cross-species transmission of astroviruses. We observe that distinct transmission zones determine the prevailing astrovirus host and virus diversity, which in turn suggests that no single host group (e.g., bats) can be the natural reservoir, as illustrated through our phylogenetic analysis.     

Rates of Viral Evolution Are Linked to Host Geography in Bat Rabies

Rates of evolution span orders of magnitude among RNA viruses with important implications for viral transmission and emergence. Although the tempo of viral evolution is often ascribed to viral features such as mutation rates and transmission mode, these factors alone cannot explain variation among closely related viruses, where host biology might operate more strongly on viral evolution. Here, we analyzed sequence data from hundreds of rabies viruses collected from bats throughout the Americas to describe dramatic variation in the speed of rabies virus evolution when circulating in ecologically distinct reservoir species. Integration of ecological and genetic data through a comparative Bayesian analysis revealed that viral evolutionary rates were labile following historical jumps between bat species and nearly four times faster in tropical and subtropical bats compared to temperate species. The association between geography and viral evolution could not be explained by host metabolism, phylogeny or variable selection pressures, and instead appeared to be a consequence of reduced seasonality in bat activity and virus transmission associated with climate. Our results demonstrate a key role for host ecology in shaping the tempo of evolution in multi-host viruses and highlight the power of comparative phylogenetic methods to identify the host and environmental features that influence transmission dynamics.     

Frequent transmission of immunodeficiency viruses among bobcats and pumas

With the exception of human immunodeficiency virus (HIV), which emerged in humans after cross-species transmissions of simian immunodeficiency viruses from nonhuman primates, immunodeficiency viruses of the family Lentiviridae represent species-specific viruses that rarely cross species barriers to infect new hosts. Among the Felidae, numerous immunodeficiency-like lentiviruses have been documented, but only a few cross-species transmissions have been recorded, and these have not been perpetuated in the recipient species. Lentivirus seroprevalence was determined for 79 bobcats (Lynx rufus) and 31 pumas (Puma concolor) from well-defined populations in Southern California. Partial genomic sequences were subsequently obtained from 18 and 12 seropositive bobcats and pumas, respectively. Genotypes were analyzed for phylogenic relatedness and genotypic composition among the study set and archived feline lentivirus sequences. This investigation of feline immunodeficiency virus infection in bobcats and pumas of Southern California provides evidence that cross-species infection has occurred frequently among these animals. The data suggest that transmission has occurred in multiple locations and are most consistent with the spread of the virus from bobcats to pumas. Although the ultimate causes remain unknown, these transmission events may occur as a result of puma predation on bobcats, a situation similar to that which fostered transmission of HIV to humans, and likely represent the emergence of a lentivirus with relaxed barriers to cross-species transmission. This unusual observation provides a valuable opportunity to evaluate the ecological, behavioral, and molecular conditions that favor repeated transmissions and persistence of lentivirus between species.     

Phylodynamics and human-mediated dispersal of a zoonotic virus

Understanding the role of humans in the dispersal of predominantly animal pathogens is essential for their control. We used newly developed Bayesian phylogeographic methods to unravel the dynamics and determinants of the spread of dog rabies virus (RABV) in North Africa. Each of the countries studied exhibited largely disconnected spatial dynamics with major geopolitical boundaries acting as barriers to gene flow. Road distances proved to be better predictors of the movement of dog RABV than accessibility or raw geographical distance, with occasional long distance and rapid spread within each of these countries. Using simulations that bridge phylodynamics and spatial epidemiology, we demonstrate that the contemporary viral distribution extends beyond that expected for RABV transmission in African dog populations. These results are strongly supportive of human-mediated dispersal, and demonstrate how an integrated phylogeographic approach will turn viral genetic data into a powerful asset for characterizing, predicting, and potentially controlling the spatial spread of pathogens.     

Influenza as a model system for studying the cross-species transfer and evolution of the SARS coronavirus

Severe acute respiratory syndrome coronavirus (SARS-CoV) moved into humans from a reservoir species and subsequently caused an epidemic in its new host. We know little about the processes that allowed the cross-species transfer of this previously unknown virus. I discuss what we have learned about the movement of viruses into humans from studies of influenza A, both how it crossed from birds to humans and how it subsequently evolved within the human population. Starting with a brief review of severe acute respiratory syndrome to highlight the kinds of problems we face in learning about this viral disease, I then turn to influenza A, focusing on three topics. First, I present a reanalysis of data used to test the hypothesis that swine served as a "mixing vessel" or intermediate host in the transmission of avian influenza to humans during the 1918 "Spanish flu" pandemic. Second, I review studies of archived viruses from the three recent influenza pandemics. Third, I discuss current limitations in using molecular data to study the evolution of infectious disease. Although influenza A and SARS-CoV differ in many ways, our knowledge of influenza A may provide important clues about what limits or favours cross-species transfers and subsequent epidemics of newly emerging pathogens.     

Metapopulation dynamics of rabies and the efficacy of vaccination

Spatial structure in a host population results in heterogeneity in transmission dynamics. We used a Bayesian framework to evaluate competing metapopulation models of rabies transmission among domestic dog populations in villages in Tanzania. A proximate indicator of disease, medical records of animal-bite injuries, is used to infer the occurrence (presence/absence) of suspected rabid dog cases in one month intervals. State-space models were used to explore the implications of different levels of reporting probability on model parameter estimates. We find evidence for a relatively high rate of infection of these populations from neighbouring districts or from other species distributed throughout the study area, rather than from adjacent wildlife protected areas, suggesting wildlife is unlikely to be implicated in the long-term persistence of rabies. Stochastic simulation of our highest ranked models in vaccinated and hypothetical unvaccinated populations indicated that pulsed vaccination campaigns occurring from 2002 to 2007 reduced rabies occurrence by 57.3 per cent in vaccinated villages in the 1 year following each pulse, and that a similar regional campaign would deliver an 80.9 per cent reduction in occurrence. This work demonstrates how a relatively coarse, proximate sentinel of rabies infection is useful for making inferences about spatial disease dynamics and the efficacy of control measures.     

Disease,habitat fragmentation and conservation

Habitat loss and the resultant fragmentation of remaining habitat is the primary cause of loss of biological diversity. How do these processes affect the dynamics of parasites and pathogens? Hess has provided some important insights into this problem using metapopulation models for pathogens that exhibit ''S-I'' dynamics; for example, pathogens such as rabies in which the host population may be divided into susceptible and infected individuals. A major assumption of Hess''s models is that infected patches become extinct, rather than recovering and becoming resistant to future infections. In this paper, we build upon this framework in two different ways: first, we examine the consequences of including patches that are resistant to infection; second, we examine the consequences of including a second species of host that can act as a reservoir for the pathogen. Both of these effects are likely to be important from a conservation perspective. The results of both sets of analysis indicate that the benefits of corridors and other connections that allow species to disperse through the landscape far outweigh the possible risks of increased pathogen transmission. Even in the commonest case, where harmful pathogens are maintained by a common reservoir host, increased landscape connectance still allows greater coexistence and persistence of a threatened or endangered host.     

Molecular epidemiology identifies only a single rabies virus variant circulating in complex carnivore communities of the Serengeti

Understanding the transmission dynamics of generalist pathogens that infect multiple host species is essential for their effective control. Only by identifying those host populations that are critical to the permanent maintenance of the pathogen, as opposed to populations in which outbreaks are the result of 'spillover' infections, can control measures be appropriately directed. Rabies virus is capable of infecting a wide range of host species, but in many ecosystems, particular variants circulate among only a limited range of potential host populations. The Serengeti ecosystem (in northwestern Tanzania) supports a complex community of wild carnivores that are threatened by generalist pathogens that also circulate in domestic dog populations surrounding the park boundaries. While the combined assemblage of host species appears capable of permanently maintaining rabies in the ecosystem, little is known about the patterns of circulation within and between these host populations. Here we use molecular phylogenetics to test whether distinct virus-host associations occur in this species-rich carnivore community. Our analysis identifies a single major variant belonging to the group of southern Africa canid-associated viruses (Africa 1b) to be circulating within this ecosystem, and no evidence for species-specific grouping. A statistical parsimony analysis of nucleoprotein and glycoprotein gene sequence data is consistent with both within- and between-species transmission events. While likely differential sampling effort between host species precludes a definitive inference, the results are most consistent with dogs comprising the reservoir of rabies and emphasize the importance of applying control efforts in dog populations.     

Phylodynamics of the Emergence of Influenza Viruses after Cross-Species Transmission

Human populations are constantly exposed to emerging pathogens such as influenza A viruses that result from cross-species transmissions. Generally these sporadic events are evolutionary dead-ends, but occasionally, viruses establish themselves in a new host that offers a novel genomic context to which the virus must adjust to avoid attenuation. However, the dynamics of this process are unknown. Here we present a novel method to characterize the time it takes to G+C composition at third codon positions (GC3 content) of influenza viruses to adjust to that of a new host. We compare the inferred dynamics in two subtypes, H1N1 and H3N2, based on complete genomes of viruses circulating in humans, swine and birds between 1900–2009. Our results suggest that both subtypes have the same fast-adjusting genes, which are not necessarily those with the highest absolute rates of evolution, but those with the most relaxed selective pressures. Our analyses reveal that NA and NS2 genes adjust the fastest to a new host and that selective pressures of H3N2 viruses are relaxed faster than for H1N1. The asymmetric nature of these processes suggests that viruses with the greatest adjustment potential to humans are coming from both birds and swine for H3N2, but only from birds for H1N1.     

The Genealogical Population Dynamics of HIV-1 in a Large Transmission Chain: Bridging within and among Host Evolutionary Rates

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the ‘store and retrieve’ hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.     

Within-host genetic diversity of endemic and emerging parvoviruses of dogs and cats

Viral emergence can result from the adaptation of endemic pathogens to new or altered host environments, a process that is strongly influenced by the underlying sequence diversity. To determine the extent and structure of intrahost genetic diversity in a recently emerged single-stranded DNA virus, we analyzed viral population structures during natural infections of animals with canine parvovirus (CPV) or its ancestor, feline panleukopenia virus (FPV). We compared infections that occurred shortly after CPV emerged with more recent infections and examined the population structure of CPV after experimental cross-species transmission to cats. Infections with CPV and FPV showed limited genetic diversity regardless of the analyzed host tissue or year of isolation. Coinfections with genetically distinct viral strains were detected in some cases, and rearranged genomes were seen in both FPV and CPV. The sporadic presence of some sequences with multiple mutations suggested the occurrence of either particularly error-prone viral replication or coinfection by more distantly related strains. Finally, some potentially organ-specific host effects were seen during experimental cross-species transmission, with many of the mutations located in the nonstructural protein NS2. These included residues with evidence of positive selection at the population level, which is compatible with a role of this protein in host adaptation.     

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 km²). 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.     

Feature Selection Methods for Identifying Genetic Determinants of Host Species in RNA Viruses

Despite environmental, social and ecological dependencies, emergence of zoonotic viruses in human populations is clearly also affected by genetic factors which determine cross-species transmission potential. RNA viruses pose an interesting case study given their mutation rates are orders of magnitude higher than any other pathogen – as reflected by the recent emergence of SARS and Influenza for example. Here, we show how feature selection techniques can be used to reliably classify viral sequences by host species, and to identify the crucial minority of host-specific sites in pathogen genomic data. The variability in alleles at those sites can be translated into prediction probabilities that a particular pathogen isolate is adapted to a given host. We illustrate the power of these methods by: 1) identifying the sites explaining SARS coronavirus differences between human, bat and palm civet samples; 2) showing how cross species jumps of rabies virus among bat populations can be readily identified; and 3) de novo identification of likely functional influenza host discriminant markers.     

A Multi-scale Analysis of Influenza A Virus Fitness Trade-offs due to Temperature-dependent Virus Persistence

Successful replication within an infected host and successful transmission between hosts are key to the continued spread of most pathogens. Competing selection pressures exerted at these different scales can lead to evolutionary trade-offs between the determinants of fitness within and between hosts. Here, we examine such a trade-off in the context of influenza A viruses and the differential pressures exerted by temperature-dependent virus persistence. For a panel of avian influenza A virus strains, we find evidence for a trade-off between the persistence at high versus low temperatures. Combining a within-host model of influenza infection dynamics with a between-host transmission model, we study how such a trade-off affects virus fitness on the host population level. We show that conclusions regarding overall fitness are affected by the type of link assumed between the within- and between-host levels and the main route of transmission (direct or environmental). The relative importance of virulence and immune response mediated virus clearance are also found to influence the fitness impacts of virus persistence at low versus high temperatures. Based on our results, we predict that if transmission occurs mainly directly and scales linearly with virus load, and virulence or immune responses are negligible, the evolutionary pressure for influenza viruses to evolve toward good persistence at high within-host temperatures dominates. For all other scenarios, influenza viruses with good environmental persistence at low temperatures seem to be favored.     

Body-size scaling in an SEI model of wildlife diseases

A number of wildlife pathogens are generalist and can affect different host species characterized by a wide range of body sizes. In this work we analyze the role of allometric scaling of host vital and epidemiological rates in a Susceptible-Exposed-Infected (SEI) model. Our analysis shows that the transmission coefficient threshold for the disease to establish in the population scales allometrically (exponent = 0.45) with host size as well as the threshold at which limit cycles occur. In contrast, the threshold of the basic reproduction number for sustained oscillations to occur is independent of the host size and is always greater than 5. In the case of rabies, we show that the oscillation periods predicted by the model match those observed in the field for a wide range of host sizes.The population dynamics of the SEI model is also analyzed in the case of pathogens affecting multiple coexisting hosts with different body sizes. Our analyses show that the basic reproduction number for limit cycles to occur depends on the ratio between host sizes, that the oscillation period in a multihost community is set by the smaller species dynamics, and that intermediate interspecific disease transmission can stabilize the epidemic occurrence in wildlife communities.     

Frontiers in research on biodiversity and disease

Global losses of biodiversity have galvanised efforts to understand how changes to communities affect ecological processes, including transmission of infectious pathogens. Here, we review recent research on diversity–disease relationships and identify future priorities. Growing evidence from experimental, observational and modelling studies indicates that biodiversity changes alter infection for a range of pathogens and through diverse mechanisms. Drawing upon lessons from the community ecology of free‐living organisms, we illustrate how recent advances from biodiversity research generally can provide necessary theoretical foundations, inform experimental designs, and guide future research at the interface between infectious disease risk and changing ecological communities. Dilution effects are expected when ecological communities are nested and interactions between the pathogen and the most competent host group(s) persist or increase as biodiversity declines. To move beyond polarising debates about the generality of diversity effects and develop a predictive framework, we emphasise the need to identify how the effects of diversity vary with temporal and spatial scale, to explore how realistic patterns of community assembly affect transmission, and to use experimental studies to consider mechanisms beyond simple changes in host richness, including shifts in trophic structure, functional diversity and symbiont composition.     

Linking anthropogenic resources to wildlife–pathogen dynamics: a review and meta‐analysis

Urbanisation and agriculture cause declines for many wildlife, but some species benefit from novel resources, especially food, provided in human‐dominated habitats. Resulting shifts in wildlife ecology can alter infectious disease dynamics and create opportunities for cross‐species transmission, yet predicting host–pathogen responses to resource provisioning is challenging. Factors enhancing transmission, such as increased aggregation, could be offset by better host immunity due to improved nutrition. Here, we conduct a review and meta‐analysis to show that food provisioning results in highly heterogeneous infection outcomes that depend on pathogen type and anthropogenic food source. We also find empirical support for behavioural and immune mechanisms through which human‐provided resources alter host exposure and tolerance to pathogens. A review of recent theoretical models of resource provisioning and infection dynamics shows that changes in host contact rates and immunity produce strong non‐linear responses in pathogen invasion and prevalence. By integrating results of our meta‐analysis back into a theoretical framework, we find provisioning amplifies pathogen invasion under increased host aggregation and tolerance, but reduces transmission if provisioned food decreases dietary exposure to parasites. These results carry implications for wildlife disease management and highlight areas for future work, such as how resource shifts might affect virulence evolution.