Genetic constraints and the adaptive evolution of rabies virus in nature

Phylogeographic dispersion and diversification of rabies virus lineages associated with dogs and crab-eating foxes (Cerdocyon thous) in Brazil

Genetic lineages of dog-associated RABV still circulate in some areas of the North and Northeast of Brazil. In parallel, another RABV lineage circulates among wild canids in the Northeast, particularly the crab-eating fox (Cerdocyon thous). Although previous studies and phylogenetic analyses have been carried out, the way in which these lineages are dispersed temporally and spatially remained to be elucidated. In this study, RABV N gene sequences isolated from canids in North and Northeast Brazil were analyzed by the Bayesian Markov Chain Monte Carlo Method, and the results were then used in a phylogeographic study. It was inferred from the findings that the most recent common ancestor became established at the end of the nineteenth century on the border of the Brazilian states of Paraíba and Pernambuco and diversified into the lineages associated with dogs and C. thous. Around 1910, the original C. thous lineage diversified into two main sublineages in the same area, one of which migrated to the south and the other to the north. The dog-associated lineage diversified around 1945 and moved toward the north and south. From the phylogeographic analysis it was possible to infer not only the movement of the virus lineages but also the probable location where dispersion and diversification occurred. The methodology used here enabled the phylogeographic history of RABV in the region to be reconstructed, and the dispersion pattern of the virus can be used to predict its movements, making it easier to stop the advance of a rabies epidemic.

Characterization of a wild rabies virus isolate of porcine origin in China

Rabies virus (RABV) that circulates worldwide in a variety of mammals can cause fatal encephalomyelitis. GD-SH-01, a street rabies virus, was isolated from a rabid pig in China. We investigated the pathogenicity of GD-SH-01 in suckling and adult mice, and compared the susceptibility of NA and BHK-21 cells in the culture to infection by GD-SH-01 and CVS-24. The complete GD-SH-01 genome sequence was determined and compared with known RABV wild strains to understand the mutations and genetic diversity that allow RABV to spread and adapt in new hosts, such as pigs. Our results suggest that GD-SH-01 possesses the characteristics of a virulent strain in Southern China and shows higher pathogenicity index than that of CVS-24 regardless of its lower level of replication in mouse brain. Up to 47 unique nucleotide substitutions were found in the genome, including five missense mutations. These data provide useful information for further understanding the transmission mechanism and the genetic variation of RABV in dissimilar hosts.

Complete genome and molecular epidemiological data infer the maintenance of rabies among kudu (Tragelaphus strepsiceros) in Namibia

Rabies in kudu is unique to Namibia and two major peaks in the epizootic have occurred since it was first noted in 1977. Due to the large numbers of kudu that were affected, it was suspected that horizontal transmission of rabies occurs among kudu and that rabies was being maintained independently within the Namibian kudu population – separate from canid cycles, despite geographic overlap. In this study, it was our aim to show, through phylogenetic analyses, that rabies was being maintained independently within the Namibian kudu population. We also tested, through complete genome sequencing of four rabies virus isolates from jackal and kudu, whether specific mutations occurred in the virus genome due to host adaptation. We found the separate grouping of all rabies isolates from kudu to those of any other canid species in Namibia, suggesting that rabies was being maintained independently in kudu. Additionally, we noted several mutations unique to isolates from kudu, suggesting that these mutations may be due to the adaptation of rabies to a new host. In conclusion, we show clear evidence that rabies is being maintained independently in the Namibian kudu population – a unique phenomenon with ecological and economic impacts.

[Rabies]

Rabies virus, a neurotropic lyssavirus responsible for unavoidable fatal encephalitis, is transmitted by saliva of infected animals through bite, scratch or licking of broken skin or a mucous membrane. Infection can be prevented by timely prevention (wash for several minutes, antisepsis and vaccination completed by antirabies immunoglobulins [Ig] according to the severity of exposure). The 55,000 human deaths estimated annually worldwide result mainly from uncontrolled canine rabies in enzootic countries (particularly in Africa and in Asia), attributable to a lack of resources or interest for this disease. Bat rabies, henceforth first cause of human's rabies in many countries in America, affects a very small number of individuals but seems more difficult to control. Shortened vaccine protocols, rationalized use of Ig and development of products of substitution should enhance access of exposed patients to prevention. Finally, research on the biological cycle, the pathogeny and on escape of virus-induced mechanisms from the immune system should continue to pave the way for presently unknown treatments of clinical rabies.

Antigenic and genetic characterization of rabies virus isolates from Uruguay

After 25 years without any reported cases of rabies in Uruguay, the northern region of the country experienced an epizootic of bovine paralytic rabies in October 2007. The outbreak affected bovines and equines, and the main source of infection was the bat Desmodus rotundus, the only hematophagous species in the country. From October 2007 to July 2008, 42 bovine, 3 equine and 120 chiropteran samples were submitted to the National Veterinary Diagnostic Laboratory for rabies testing. A total of 12 samples (7 bovine, 2 equine and 3 from D. rotundus) were positive by the fluorescent antibody test, and viruses were isolated by the mouse inoculation test. The objective of this study was to compare the antigenic and genetic characteristics of these isolates and three isolates from insectivorous bats from other regions. Antigenic typing using a panel of eight monoclonal antibodies identified all 12 viruses as variant 3 (AgV3), a variant associated with D. rotundus. Two isolates from insectivorous bats (Tadarida brasiliensis and Molossus sp.) were characterized as antigenic variant 4 (AgV4) while the third, from Myotis sp., could not be characterized using this panel as its reactivity pattern did not match that of any of the known antigenic variants. Partial N-gene sequences (nt 149-1420) of these isolates were aligned with homologous sequences derived from GenBank by the CLUSTAL/W method and used to build a neighbor-joining distance tree with the Kimura 2-parameter model. All 12 isolates were genetically grouped into the D. rotundus cluster as they shared 100% identity. In the phylogenetic analysis, the three isolates from insectivorous bats segregated into three clusters: one related to T. brasiliensis, one to Myotis sp. and the other to Lasiurus sp., although the isolate associated with the latter came from a Molossus sp. specimen. These results indicate that AgV3 was associated with the outbreak of bovine paralytic rabies in Uruguay. This is the first report of rabies virus having been detected in non-hematophagous bats in this country.

Molecular characterization of KGH, the first human isolate of rabies virus in Korea

The complete genome sequence of the KGH strain of the first human rabies virus, which was isolated from a skin biopsy of a patient with rabies, whose symptoms developed due to bites from a raccoon dog in 2001. The size of the KGH strain genome was determined to be 11,928 nucleotides (nt) with a leader sequence of 58 nt, nucleoprotein gene of 1,353 nt, phosphoprotein gene of 894 nt, matrix protein gene of 609 nt, glycoprotein gene of 1,575 nt, RNA-dependent RNA polymerase gene of 6,384 nt, and trailer region of 69 nt. Sequence similarity was compared with 39 fully sequenced rabies virus genomes currently available, and the result showed 70.6-91.6 % at the nucleotide level, and 82.8-97.9 % at the amino acid level. The deduced amino acids in the viral protein were compared with those of other rabies viruses, and various functional regions were investigated. As a result, we found that the KGH strain only had a unique amino acid substitution that was identified to be associated either with host immune response and pathogenicity in the N protein, or with a related region regulating STAT1 in the P protein, and related to pathogenicity in G protein. Based on phylogenetic analyses using the complete genome of 39 rabies viruses, the KGH strain was determined to be closely related with the NNV-RAB-H strain and transplant rabies virus serotype 1, which are Indian isolates, and was confirmed to belong to the Arctic-like 2 clade. The KGH strain was most closely related to the SKRRD0204HC and SKRRD0205HC strain when compared with Korean animal isolates, which was separated around the same time and place, and belonged to the Gangwon III subgroup.

Epidemiology and host spectrum of Borna disease virus infections

Borna disease virus (BDV) has gained lot of interest because of its zoonotic potential, ability to introduce cDNA of its RNA transcripts into host genomes, and ability to cause severe neurobehavioural diseases. Classical Borna disease is a progressive meningoencephalomyelitis in horses and sheep, known in central Europe for centuries. According to current knowledge, BDV or a close relative also infects several other species, including humans at least occasionally, in central Europe and elsewhere, but the existence of potential 'human Borna disease' with its suspected neuropsychiatric symptoms is highly controversial. The recent detection of endogenized BDV-like genes in primate and various other vertebrate genomes confirms that at least ancient bornaviruses did infect our ancestors. The epidemiology of BDV is largely unknown, but accumulating evidence indicates vectors and reservoirs among small wild mammals. The aim of this review is to bring together the current knowledge on epidemiology of BDV infections. Specifically, geographical and host distribution are addressed and assessed in the critical light of the detection methods used. We also review some salient clinical aspects.

Recombination and natural selection in hepatitis E virus genotypes

To gain new insights into the evolutionary processes that created the genetic diversity of the hepatitis E virus (HEV), the Recombination Detection Program (RDP) and SimPlot program were employed to detect recombination events in the genome, then the fixed-effects likelihood (FEL) method was used to detect natural selection effects on viral proteins. Recombination analysis provided strong evidence for both intergenotype and intragenotype recombination events in the sequences analyzed. Recombination events were found to be distributed non-randomly, with the highest frequency in the X domain and the helicase. Strain DQ450072 was identified as intergenotype-recombinant. Natural selection analysis revealed that codons under both negative selection and positive selection were distributed non-randomly. ORF1 and ORF2 have experienced strong purifying selection across genotypes. Furthermore, potentially important sites were also found under positive selection in the N-terminal end of ORF2 and the C-terminal end of ORF3. No significant difference was found among the selective pressures on different genotypes.

Variable evolutionary routes to host establishment across repeated rabies virus host shifts among bats

Determining the genetic pathways that viruses traverse to establish in new host species is crucial to predict the outcome of cross-species transmission but poorly understood for most host-virus systems. Using sequences encoding 78% of the rabies virus genome, we explored the extent, repeatability and dynamic outcome of evolution associated with multiple host shifts among New World bats. Episodic bursts of positive selection were detected in several viral proteins, including regions associated with host cell interaction and viral replication. Host shifts involved unique sets of substitutions, and few sites exhibited repeated evolution across adaptation to many bat species, suggesting diverse genetic determinants over host range. Combining these results with genetic reconstructions of the demographic histories of individual viral lineages revealed that although rabies viruses shared consistent three-stage processes of emergence in each new bat species, host shifts involving greater numbers of positively selected substitutions had longer delays between cross-species transmission and enzootic viral establishment. Our results point to multiple evolutionary routes to host establishment in a zoonotic RNA virus that may influence the speed of viral emergence.

Desmodus rotundus and Artibeus spp. bats might present distinct rabies virus lineages

In Brazil, bats have been assigned an increasing importance in public health as they are important rabies reservoirs. Phylogenetic studies have shown that rabies virus (RABV) strains from frugivorous bats Artibeus spp. are closely associated to those from the vampire bat Desmodus rotundus, but little is known about the molecular diversity of RABV in Artibeus spp. The N and G genes of RABV isolated from Artibeus spp. and cattle infected by D. rotundus were sequenced, and phylogenetic trees were constructed. The N gene nucleotides tree showed three clusters: one for D. rotundus and two for Artibeus spp. Regarding putative N amino acid-trees, two clusters were formed, one for D. rotundus and another for Artibeus spp. RABV G gene phylogeny supported the distinction between D. rotundus and Artibeus spp. strains. These results show the intricate host relationship of RABV's evolutionary history, and are invaluable for the determination of RABV infection sources.

Molecular inferences suggest multiple host shifts of rabies viruses from bats to mesocarnivores in Arizona during 2001-2009

In nature, rabies virus (RABV; genus Lyssavirus, family Rhabdoviridae) represents an assemblage of phylogenetic lineages, associated with specific mammalian host species. Although it is generally accepted that RABV evolved originally in bats and further shifted to carnivores, mechanisms of such host shifts are poorly understood, and examples are rarely present in surveillance data. Outbreaks in carnivores caused by a RABV variant, associated with big brown bats, occurred repeatedly during 2001–2009 in the Flagstaff area of Arizona. After each outbreak, extensive control campaigns were undertaken, with no reports of further rabies cases in carnivores for the next several years. However, questions remained whether all outbreaks were caused by a single introduction and further perpetuation of bat RABV in carnivore populations, or each outbreak was caused by an independent introduction of a bat virus. Another question of concern was related to adaptive changes in the RABV genome associated with host shifts. To address these questions, we sequenced and analyzed 66 complete and 20 nearly complete RABV genomes, including those from the Flagstaff area and other similar outbreaks in carnivores, caused by bat RABVs, and representatives of the major RABV lineages circulating in North America and worldwide. Phylogenetic analysis demonstrated that each Flagstaff outbreak was caused by an independent introduction of bat RABV into populations of carnivores. Positive selection analysis confirmed the absence of post-shift changes in RABV genes. In contrast, convergent evolution analysis demonstrated several amino acids in the N, P, G and L proteins, which might be significant for pre-adaptation of bat viruses to cause effective infection in carnivores. The substitution S/T₂₄₂ in the viral glycoprotein is of particular merit, as a similar substitution was suggested for pathogenicity of Nishigahara RABV strain. Roles of the amino acid changes, detected in our study, require additional investigations, using reverse genetics and other approaches.

Host shifts of the rabies virus (RABV) from bats to carnivores are important for our understanding of viral evolution and emergence, and have significant public health implications, particularly for the areas where “terrestrial” rabies has been eliminated. In this study we addressed several rabies outbreaks in carnivores that occurred in the Flagstaff area of Arizona during 2001–2009, and caused by the RABV variant associated with big brown bats (Eptesicus fuscus). Based on phylogenetic analysis we demonstrated that each outbreak resulted from a separate introduction of bat RABV into populations of carnivores. No post-shift changes in viral genomes were detected under the positive selection analysis. Trying to answer the question why certain bat RABV variants are capable for host shifts to carnivores and other variants are not, we developed a convergent evolution analysis, and implemented it for multiple RABV lineages circulating worldwide. This analysis identified several amino acids in RABV proteins which may facilitate host shifts from bats to carnivores. Precise roles of these amino acids require additional investigations, using reverse genetics and animal experimentation. In general, our approach and the results obtained can be used for prediction of host shifts and emergence of other zoonotic pathogens.

Evolution and biogeography of an emerging quasispecies: diversity patterns of the fish Viral Hemorrhagic Septicemia virus (VHSv)

Viral Hemorrhagic Septicemia virus (VHSv) is an RNA rhabdovirus that causes one of the most important finfish diseases, affecting over 70 marine and freshwater species. It was discovered in European cultured fish in 1938 and since has been described across the Northern Hemisphere. Four strains and several substrains have been hypothesized, whose phylogenetic relationships and evolutionary radiation are evaluated here in light of a quasispecies model, including an in-depth analysis of the novel and especially virulent new substrain (IVb) that first appeared in the North American Laurentian Great Lakes in 2003. We analyze the evolutionary patterns, genetic diversity, and biogeography of VHSv using all available RNA sequences from the glycoprotein (G), nucleoprotein (N), and non-virion (Nv) genes, with Maximum Likelihood and bayesian approaches. Results indicate that the G gene evolves at an estimated rate of μ=2.58×10(-4) nucleotide substitutions per site per year, the N gene at μ=4.26×10(-4), and Nv fastest at μ=1.25×10(-3). Phylogenetic trees from the three genes largely are congruent, distinguishing strains I-IV as reciprocally monophyletic with high bootstrap and posterior probability support. VHSv appears to have originated from a marine ancestor in the North Atlantic Ocean, diverging into two primary clades: strain IV in North America (the Northwestern Atlantic Ocean), and strains I-III in the Northeastern Atlantic region (Europe). Strain II may comprise the basal group of the latter clade and diverged in Baltic Sea estuarine waters; strains I and III appear to be sister groups (according to the G and Nv genes), with the former mostly in European freshwaters and the latter in North Sea marine/estuarine waters. Strain IV is differentiated into three monophyletic substrains, with IVa infecting Northeastern Pacific salmonids and many marine fishes (with 44 unique G gene haplotypes), IVb endemic to the freshwater Great Lakes (11 haplotypes), and a newly-designated IVc in marine/estuarine North Atlantic waters (five haplotypes). Two separate substrains independently appeared in the Northwestern Pacific region (Asia) in 1996, with Ib originating from the west and IVa from the east. Our results depict an evolutionary history of relatively rapid population diversifications in star-like patterns, following a quasispecies model. This study provides a baseline for future tracking of VHSv spread and interpreting its evolutionary diversification pathways.

Positive selection sites in the surface genes of dengue virus: phylogenetic analysis of the interserotypic branches of the four serotypes

The existence of four dengue serotypes is associated with a phenomenon called "Antibody-Dependent Enhancement" that has been suggested to cause a severe form of dengue hemorrhagic fever and shock syndrome. To study the evolutionary event that drove the serotype separation, we employed the maximum likelihood approach by focusing on the Premembrane (prM) and Envelop (E) genes. We showed that the separation of dengue serotypes had been dominantly under purifying selection. In spite of the strong selective constraint, one codon of prM gene and twelve codons of E gene were detected to be under positive selection. This indicates that the E protein might have been under a stronger positive pressure than the PrM protein. The codons under positive selection were identified along the interserotypic branches, suggesting that changes at these sites were probably associated with the emergence of the four serotypes and/or adaptation to the new transmission environments.

Molecular phylogenetics of the lyssaviruses--insights from a coalescent approach

Technical improvements over the past 2 decades have enormously facilitated the generation of nucleotide sequence data for lyssavirus collections. These databases are amenable to methods of phylogenetic analysis, which attempt to define the taxonomic structure of this genus and predict the evolutionary relationships of current circulating strains. Coupled with a range of mathematical tools to explore the appropriateness of nucleotide substitution models and test for positive selection, the evolutionary process is being explored in detail. Despite the potential for high viral mutation levels, the operation of purifying selection appears to effectively constrain lyssavirus evolution. The recent development of coalescent theory has provided additional approaches to data analysis whereby the time frame of emergence of viral lineages can be most reliably estimated. Such studies suggest that all currently circulating rabies viruses have emerged within the past 1500 years. Moreover, through the capability of analyzing viral population dynamics and determining patterns of population size variation, coalescent approaches can provide insight into the demographics of viral outbreaks. Whereas human-assisted movement of reservoir host species has clearly facilitated transfer of rabies between continents, topographical landscape features significantly influence the rate and extent of contiguous disease spread. Together with empirical studies on virus diversity, the application of coalescent approaches will help to better understand lyssavirus emergence, evolution, and spread. In particular, such methods are presently facilitating exploration of the factors operating to limit the ability of lyssaviruses to establish new persistent virus-host associations and ultimately control the emergence of new species of this genus.

Genetic diversity and positive selection analysis of classical swine fever virus isolates in south China

Classical swine fever virus (CSFV) causes a highly contagious disease that leads to significant economic losses in the pig industry worldwide. However, there is a paucity of knowledge on the accurate genotyping of CSFV isolates in south China. This study genotyped the E2 gene of 14 CSFV strains isolated during 2008-2010 from domestic pigs in different districts of south China. Phylogenetic analyses revealed that all of the 14 CSFV isolates were clustered into genetic subgroup 1.1. This contrasts with most parts of China, where group 2 isolates are predominant. Furthermore, the positive selection pressures acting on the E(rns) and E2 envelope protein genes of CSFV were assessed and a site-by-site analysis of the dN/dS ratio was performed to identify specific codons that undergo diversification under positive selection. While no significant evidence for positive selection was observed in E(rns), two positively selected sites at amino acid residues 49 and 72 in the E2 encoding region were identified. Our results revealed that a predominance of subgroup 1.1 CSFV isolates is currently circulating in some districts of south China, which appear to be unrelated to the Chinese C-strain vaccine. Moreover, the envelope protein gene, E2, has undergone positive selection in 14 CSFV strains and two positively selected sites have been identified in this study. Understanding the molecular epidemiology and functional importance of these positively selected amino acid positions could help to predict possible changes in virulence, the development of vaccines and disease control.

Computational analysis suggests that lyssavirus glycoprotein gene plays a minor role in viral adaptation

The Lyssavirus glycoprotein (G) is a membrane protein responsible for virus entry and protective immune responses. To explore possible roles of the glycoprotein in host shift or adaptation of Lyssavirus, we retrieved 53 full-length glycoprotein gene sequences from NCBI GenBank. The sequences were from different host isolates over a period of 70 years in 21 countries. Computational analyses detected 1 recombinant (AY987478, a dog isolate of CHAND03, genotype 1 in India) with incongruent phylogenetic support. No recombination was detected when AY98748 was excluded in the analyses. We applied different selection models to identify selection pressure on the glycoprotein gene. One codon at amino acid residual 483 was found to be under weak positive selection with marginal probability of 95% by using the maximum likelihood method. We found no significant evidence of positive selection on any site of the glycoprotein gene when the putative recombinant AY987478 was excluded. The computational analyses suggest that the G gene has been under purifying selection and that the evolution of the G gene may not play a significant role in Lyssavirus adaptation.

Sequence comparison between three geographically distinct Spodoptera frugiperda multiple nucleopolyhedrovirus isolates: Detecting positively selected genes

The complete genomic sequence of a Nicaraguan plaque purified Spodoptera frugiperda nucleopolyhedrovirus (SfMNPV) genotype SfMNPV-B was determined and compared to previously sequenced isolates from United States (SfMNPV-3AP2) and Brazil (SfMNPV-19). The genome of SfMNPV-B (132,954bp) was 1623bp and 389bp larger than that of SfMNPV-3AP2 and SfMNPV-19, respectively. Genome size differences were mainly due to a deletion located in the SfMNPV-3AP2 egt region and small deletions and point mutations in SfMNPV-19. Nucleotide sequences were strongly conserved (99.35% identity) and a high degree of predicted amino acid sequence identity was observed. A total of 145 open reading frames (ORFs) were identified in SfMNPV-B, two of them (sf39a and sf110a) had not been previously identified in the SfMNPV-3AP2 and SfMNPV-19 genomes and one (sf57a) was absent in both these genomes. In addition, sf6 was not previously identified in the SfMNPV-19 genome. In contrast, SfMNPV-B and SfMNPV-19 both lacked sf129 that had been reported in SfMNPV-3AP2. In an effort to identify genes potentially involved in virulence or in determining population adaptations, selection pressure analysis was performed. Three ORFs were identified undergoing positive selection: sf49 (pif-3), sf57 (odv-e66b) and sf122 (unknown function). Strong selection for ODV envelope protein genes indicates that the initial infection process in the insect midgut is one critical point at which adaptation acts during the transmission of these viruses in geographically distant populations. The function of ORF sf122 is being examined.

Quantifying antigenic relationships among the lyssaviruses

All lyssaviruses cause fatal encephalitis in mammals. There is sufficient antigenic variation within the genus to cause variable vaccine efficacy, but this variation is difficult to characterize quantitatively: sequence analysis cannot yet provide detailed antigenic information, and antigenic neutralization data have been refractory to high-resolution robust interpretation. Here, we address these issues by using state-of-the-art antigenic analyses to generate a high-resolution antigenic map of a global panel of 25 lyssaviruses. We compared the calculated antigenic distances with viral glycoprotein ectodomain sequence data. Although 67% of antigenic variation was predictable from the glycoprotein amino acid sequence, there are in some cases substantial differences between genetic and antigenic distances, thus highlighting the risk of inferring antigenic relationships solely from sequence data at this time. These differences included epidemiologically important antigenic differences between vaccine strains and wild-type rabies viruses. Further, we quantitatively assessed the antigenic relationships measured by using rabbit, mouse, and human sera, validating the use of nonhuman experimental animals as a model for determining antigenic variation in humans. The use of passive immune globulin is a crucial component of rabies postexposure prophylaxis, and here we also show that it is possible to predict the reactivity of immune globulin against divergent lyssaviruses.

Contrasting landscape epidemiology of two sympatric rabies virus strains

Viral strain evolution and disease emergence are influenced by anthropogenic change to the environment. We investigated viral characteristics, host ecology, and landscape features in the rabies-striped skunk disease system of the central Great Plains to determine how these factors interact to influence disease emergence. We amplified portions of the N and G genes of rabies viral RNA from 269 samples extracted from striped skunk brains throughout the distribution of two different rabies strains for which striped skunks were the reservoir. Because the distribution of these two strains overlapped on the landscape and were present in the same host population, we could evaluate how viral properties influenced epidemiological patterns in the area of sympatry. We found that South Central Skunk rabies (SCSK) exhibited intense purifying selection and high infectivity, which are both characteristics of an epizootic virus. Conversely, North Central Skunk rabies (NCSK) exhibited relaxed purifying selection and comparatively lower infectivity, suggesting the presence of an enzootic virus. The host population in the area of sympatry was highly admixed, and skunks among allopatric and sympatric areas had similar effective population sizes. Spatial analysis indicated that landscape features had minimal influence on NCSK movement across the landscape, but those same features were partial barriers to the spread of SCSK. We conclude that NCSK and SCSK have different epidemiological properties that interact differently with both host and landscape features to influence rabies spread in the central Great Plains. We suggest a holistic approach for future studies of emerging infectious diseases that includes studies of viral properties, host characteristics, and spatial features.

Evolutionary history of African mongoose rabies

Two biotypes or variants of rabies virus (RABV) occur in southern Africa. These variants are respectively adapted to hosts belonging to the Canidae family (the canid variant) and hosts belonging to the Herpestidae family (the mongoose variant). Due to the distinct host adaptation and differences in epidemiology and pathogenesis, it has been hypothesized that the two variants were introduced into Africa at different times. The objective of this study was to investigate the molecular phylogeny of representative RABV isolates of the mongoose variant towards a better understanding of the origins of this group. The study was based on an analysis of the full nucleoprotein and glycoprotein gene sequences of a panel of 27 viruses. Phylogenetic analysis of this dataset confirmed extended evolutionary adaptation of isolates in specific geographic areas. The evolutionary dynamics of this virus variant was investigated using Bayesian methodology, allowing for rate variation among viral lineages. Molecular clock analysis estimated the age of the African mongoose RABV to be approximately 200 years old, which is in concurrence with literature describing rabies in mongooses since the early 1800 s.

Low genetic diversities of rabies virus populations within different hosts in Brazil

The low rates of nonsynonymous evolution observed in natural rabies virus (RABV) isolates are suggested to have arisen in association with the structural and functional constraints operating on the virus protein and the infection strategies employed by RABV within infected hosts to avoid strong selection by the immune response. In order to investigate the relationship between the genetic characteristics of RABV populations within hosts and the virus evolution, the present study examined the genetic heterogeneities of RABV populations within naturally infected dogs and foxes in Brazil, as well as those of bat RABV populations that were passaged once in suckling mice. Sequence analyses of complete RABV glycoprotein (G) genes showed that RABV populations within infected hosts were genetically highly homogeneous whether they were infected naturally or experimentally (nucleotide diversities of 0-0.95x10(-3)). In addition, amino acid mutations were randomly distributed over the entire region of the G protein, and the nonsynonymous/synonymous rate ratios (d(N)/d(S)) for the G protein gene were less than 1. These findings suggest that the low genetic diversities of RABV populations within hosts reflect the stabilizing selection operating on the virus, the infection strategies of the virus, and eventually, the evolutionary patterns of the virus.

A history estimate and evolutionary analysis of rabies virus variants in China

To investigate the evolutionary dynamics of rabies virus (RABV) in China, we collected and sequenced 55 isolates sampled from 14 Chinese provinces over the last 40 years and performed a coalescent-based analysis of the G gene. This revealed that the RABV currently circulating in China is composed of three main groups. Bayesian coalescent analysis estimated the date of the most recent common ancestor for the current RABV Chinese strains to be 1412 (with a 95 % confidence interval of 1006-1736). The estimated mean substitution rate for the G gene sequences (3.961x10(-4) substitutions per site per year) was in accordance with previous reports for RABV.

Evolutionary history and dynamics of dog rabies virus in western and central Africa

The burden of rabies in Africa is estimated at 24,000 human deaths year(-1), almost all of which result from infection with dog rabies viruses (RABV). To investigate the evolutionary dynamics of RABV in western and central Africa, 92 isolates sampled from 27 African countries over 29 years were collected and sequenced. This revealed that RABV currently circulating in dogs in this region fell into a single lineage designated 'Africa 2'. A detailed analysis of the phylogeographical structure of this Africa 2 lineage revealed strong population subdivision at the country level, with only limited movement of virus among localities, including a possible east-to-west spread across Africa. In addition, Bayesian coalescent analysis suggested that the Africa 2 lineage was introduced into this region of Africa only recently (probably <200 years ago), in accordance with the timescale of expanding European colonial influence and urbanization, and then spread relatively slowly, perhaps occupying the entire region in a 100 year period.

The origin and phylogeography of dog rabies virus

Rabies is a progressively fatal and incurable viral encephalitis caused by a lyssavirus infection. Almost all of the 55 000 annual rabies deaths in humans result from infection with dog rabies viruses (RABV). Despite the importance of rabies for human health, little is known about the spread of RABV in dog populations, and patterns of biodiversity have only been studied in limited geographical space. To address these questions on a global scale, we sequenced 62 new isolates and performed an extensive comparative analysis of RABV gene sequence data, representing 192 isolates sampled from 55 countries. From this, we identified six clades of RABV in non-flying mammals, each of which has a distinct geographical distribution, most likely reflecting major physical barriers to gene flow. Indeed, a detailed analysis of phylogeographic structure revealed only limited viral movement among geographical localities. Using Bayesian coalescent methods we also reveal that the sampled lineages of canid RABV derive from a common ancestor that originated within the past 1500 years. Additionally, we found no evidence for either positive selection or widespread population bottlenecks during the global expansion of canid RABV. Overall, our study reveals that the stochastic processes of genetic drift and population subdivision are the most important factors shaping the global phylogeography of canid RABV.

A reassessment of the emergence time of European bat lyssavirus type 1

The previous study of the evolutionary rates of European bat lyssavirus type 1 (EBLV-1) used a strict molecular clock to estimate substitution rates of the nucleoprotein gene and in turn times of the most recent common ancestor (tMRCA) of the entire genotype and the two major EBLV-1 lineages (EBLV-1A and EBLV-1B). The results of that study suggested that the evolutionary rate of EBLV-1 was one of the lowest recorded for RNA viruses and that genetic diversity of EBLV-1 arose 500-750 years ago. Here I have shown that the use of a relaxed molecular clock (allowing branch rates to vary within a phylogeny) shows that these previous estimates should be revised. The relaxed clock provides a significantly better fit to all datasets. The substitution rate of EBLV-1B is compatible to that expected given previous estimates for the N gene of rabies virus whilst rate estimations for EBLV-1A appear to be confounded by substantial rate variation within the phylogeny. The relaxed clock substitution rate for EBLV-1 (1.1 x 10(-4)) is higher than had been estimated previously, and closer to that expected for the N gene. Moreover, tMRCA estimates for EBLV-1 are substantially reduced using the relaxed molecular clock (70-300 years) although the differing dynamics of EBLV-1A and EBLV-1B confound the confidence in this estimate. Current diversity of both EBLV-1A and EBLV-1B appears to have emerged within the last 100 years. Reconstruction of the population histories suggests that EBLV-1B may be emerging whilst the signal derived from the EBLV-1A phylogeny may be dampened by clade-specific dynamics.

Molecular diversity and phylogeny of Hantaan virus in Guizhou, China: evidence for Guizhou as a radiation center of the present Hantaan virus

To gain further insight into the molecular epidemiology of Hantaan virus (HTNV) in Guizhou, China, rodents were captured in this region in 2004 and 2005. In addition, serum samples were collected from four patients. Ten hantaviruses were isolated successfully in cell culture from four humans, two Apodemus agrarius, three Rattus norvegicus and one Rattus nitidus. The nucleotide sequences for their small (S), medium (M) and partial large (L) segments were determined. Phylogenetic analysis of the S and M segment sequences revealed that all of these isolates belong to the species HTNV, suggesting a spillover of HTNV from A. agrarius to Rattus rats. All available isolates from Guizhou were divided into four distinct groups either in the S segment tree or in the M segment tree. The clustering pattern of these isolates in the S segment tree was not in agreement with that in the M or L segment tree, showing that genetic reassortment between HTNV had occurred naturally. Analysis of the S segment sequences from available HTNV strains indicated that they formed three clades. The first clade, which comprised only viruses from Guizhou, was the outgroup of clades II and III. The viruses in the second clade were found in Guizhou and mainly in the far-east Asian region, including China. However, the viruses in the third clade were found in most areas of China, including Guizhou, in which haemorrhagic fever with renal syndrome (HFRS) is endemic. Our results reveal that the highest genetic diversity of HTNV is in a limited geographical region of Guizhou, and suggest that Guizhou might be a radiation centre of the present form of HTNV.

Detecting molecular adaptation at individual codons in the glycoprotein gene of the geographically diversified infectious hematopoietic necrosis virus, a fish rhabdovirus

Salmonid fishes, the principal hosts of the infectious hematopoietic necrosis virus (IHNV), are a candidate species for aquaculture in many countries. IHNV causes an acute disease resulting in severe economic loss in salmonid fish farming. Previous phylogenetic analyses revealed the existence of multiple genogroups of this virus throughout the geographical range of its host. Here, we report the importance of natural selection in shaping the evolution of certain codons at the surface glycoprotein (G-protein) gene of this virus. Maximum likelihood (ML)-based codon substitution analyses revealed that approximately 2.8% of the codons for the entire G-protein are shown to have higher nonsynonymous substitution per nonsynonymous site (dn) than the synonymous substitutions per synonymous site (ds) (dn/ds=omega>4.335). Thus, the data suggest that positive selection (omega>1) is the major driving force in the evolution of certain codons. However, majority of these positively selected sites cannot be mapped to the regions of antigenic determinants of IHNV. Based on the reports of previous studies, epitopes with positively selected sites are immunodominant and viruses can escape from immune responses by producing antigenic variation at positively selected sites, therefore, vaccines directed against these neutralizing epitopes of IHNV that consist of no positively selected sites will be more effective. Some of the positively selected sites showed radical change in amino acids with respect to their charge and polarity; however, it is unclear how these changes affect the fitness of the virus.

Complete genome sequence of a raccoon rabies virus isolate

The entire genome of a mid-Atlantic raccoon strain rabies virus (RRV) isolated in Canada was sequenced; this is the second North American wildlife rabies virus isolate to be fully characterized. The overall organization and length of the genome was similar to that of other lyssaviruses. The nucleotide sequence identity of the raccoon strain ranged between 32.7% and 85.0% when compared to other lyssaviruses, while the deduced amino acid sequence identity ranged between 22.9% and 94.2% with the nucleoprotein and polymerase being the most conserved. Notable features of RRV include the phosphoprotein's four amino acid extension compared to most other rabies viruses, and a nucleotide substitution immediately prior to the normal start codon that results in an additional methionine at the beginning of the L protein. This is the first report of the RRV L gene sequence and its 2128 amino acid product. Rates of non-synonymous and synonymous nucleotide changes within the lyssavirus L gene identified the conserved blocks II, III and IV as being most constrained. Analysis of L gene codon substitution patterns favoured models that supported positive selection, but only one site, corresponding to Leu62 of the RRV L protein, was identified as being under weak positive selection.

Complete genomes of Aravan, Khujand, Irkut and West Caucasian bat viruses, with special attention to the polymerase gene and non-coding regions

The purpose of this study was to generate complete genome sequences of Aravan (ARAV), Khujand (KHUV), Irkut (IRKV) and West Caucasian bat (WCBV) viruses, and to compare them with genomes of other lyssaviruses. We focused on RNA-dependent RNA-polymerase (L) and non-coding regions, because other genes of these viruses have been described previously. The L protein is organized into six conserved blocks (I-VI), previously detected in all Mononegavirales. Furthermore, lyssaviruses have two additional conserved regions, L1 and L2, located in the COOH part of the L. L1 may be responsible for methylation of viral mRNA cap structures, whereas the significance of L2 is unclear. Phylogenetic patterns based on the L are similar to those described for the nucleoprotein. The WCBV is the most divergent member of the genus. Besides phylogeny, it has a short trailer region (57 nucleotides versus 69-70 nucleotides in other lyssaviruses) and different intergenic region lengths, including an exceptionally long non-coding region of the glycoprotein (697 nucleotides) containing a potential open reading frame of 180 nucleotides. The absence of a flanking transcription initiation signal, as well as Northern and Western blot data, suggests that this region is not independently transcribed but is a part of G mRNA.

Genomic diversity and evolution of the lyssaviruses

Lyssaviruses are RNA viruses with single-strand, negative-sense genomes responsible for rabies-like diseases in mammals. To date, genomic and evolutionary studies have most often utilized partial genome sequences, particularly of the nucleoprotein and glycoprotein genes, with little consideration of genome-scale evolution. Herein, we report the first genomic and evolutionary analysis using complete genome sequences of all recognised lyssavirus genotypes, including 14 new complete genomes of field isolates from 6 genotypes and one genotype that is completely sequenced for the first time. In doing so we significantly increase the extent of genome sequence data available for these important viruses. Our analysis of these genome sequence data reveals that all lyssaviruses have the same genomic organization. A phylogenetic analysis reveals strong geographical structuring, with the greatest genetic diversity in Africa, and an independent origin for the two known genotypes that infect European bats. We also suggest that multiple genotypes may exist within the diversity of viruses currently classified as 'Lagos Bat'. In sum, we show that rigorous phylogenetic techniques based on full length genome sequence provide the best discriminatory power for genotype classification within the lyssaviruses.

Pathogen evolution and disease emergence in carnivores

Emerging infectious diseases constitute some of the most pressing problems for both human and domestic animal health, and biodiversity conservation. Currently it is not clear whether the removal of past constraints on geographical distribution and transmission possibilities for pathogens alone are sufficient to give rise to novel host-pathogen combinations, or whether pathogen evolution is also generally required for establishment in novel hosts. Canine distemper virus (CDV) is a morbillivirus that is prevalent in the world dog population and poses an important conservation threat to a diverse range of carnivores. We performed an extensive phylogenetic and molecular evolution analysis on complete sequences of all CDV genes to assess the role of selection and recombination in shaping viral genetic diversity and driving the emergence of CDV in non-dog hosts. We tested the specific hypothesis that molecular adaptation at known receptor-binding sites of the haemagglutinin gene is associated with independent instances of the spread of CDV to novel non-dog hosts in the wild. This hypothesis was upheld, providing compelling evidence that repeated evolution at known functional sites (in this case residues 530 and 549 of the haemagglutinin molecule) is associated with multiple independent occurrences of disease emergence in a range of novel host species.

Lyssaviruses: current trends

Various technological developments have revitalized the approaches employed to study the disease of rabies. In particular, reverse genetics has facilitated the generation of novel viruses used to improve our understanding of the fundamental aspects of rabies virus (RABV) biology and pathogenicity and yielded novel constructs potentially useful as vaccines against rabies and other diseases. Other techniques such as high throughput methods to examine the impact of rabies virus infection on host cell gene expression and two hybrid systems to explore detailed protein-protein interactions also contribute substantially to our understanding of virus-host interactions. This review summarizes much of the increased knowledge about rabies that has resulted from such studies but acknowledges that this is still insufficient to allow rational attempts at curing those who present with clinical disease.

Human monoclonal antibody and vaccine approaches to prevent human rabies

Rabies, being a major zoonotic disease, significantly impacts global public health. It is invariably fatal once clinical signs are apparent. The majority of human rabies deaths occur in developing countries. India alone reports more than 50% of the global rabies deaths. Although it is a vaccine-preventable disease, effective rabies prevention in humans with category III bites requires the combined administration of rabies immunoglobulin (RIG) and vaccine. Cell culture rabies vaccines have become widely available in developing countries, virtually replacing the inferior and unsafe nerve tissue vaccines. Limitations inherent to the conventional RIG of either equine or human origin have prompted scientists to look for monoclonal antibody-based human RIG as an alternative. Fully human monoclonal antibodies have been found to be safer and equally efficacious than conventional RIG when tested in mice and hamsters. In this chapter, rabies epidemiology, reservoir control measures, post-exposure prophylaxis of human rabies, and combination therapy for rabies are discussed. Novel human monoclonal antibodies, their production, and the significance of plants as expression platforms are emphasized.

Are all lyssavirus genes equal for phylogenetic analyses?

Individual lyssavirus genes were evaluated for phylogenetic studies from available full genome sequences. The full genome of the ERA rabies virus was sequenced and its accuracy was confirmed through virus recovery by reverse genetics. The full length of the ERA is 11,931 nucleotides (nt), with a leader sequence of 58 nt, the nucleoprotein (N) gene of 1350 nt, phosphoprotein (P) gene of 891 nt, matrix protein (M) gene of 606 nt, glycoprotein (G) gene of 1572 nt, RNA-dependent RNA polymerase (L) gene of 6384 nt, Psi-region (or G-L intergenic region) of 400 nt, and a trailer region of 70 nt. The five mono-cistrons are separated by intergenic regions of 2, 5, 5 and 24 nt, respectively. One obvious difference between the ERA and SAD-B19 rabies virus strains was the putative stop/polyadenylation signal of the G gene, with a poly(A(8)) tract for ERA, and a poly(A(5)) for SAD-B19. The TGpoly(A(8)) sequence tract was identified to be a leaky termination signal in the ERA strain. Through analyses of nt diversity, protein co-variations, structural and functional constraints, and reconstruction of phylogenetic trees from comprehensive datasets, we propose lyssavirus genes probably are of similar value for phylogenetic analyses.

Molecular epidemiology of rabies in Guangxi Province, south of China

BACKGROUND

Surveillance data for rabies in Guangxi Province in China showed that human rabies cases have gradually increased since 1996.

OBJECTIVE

To evaluate the epidemiology of rabies at the molecular level and provide suggestions for effective prevention of rabies in Guangxi.

STUDY DESIGN

Since 2000, 1569 brains from suspected rabid animals were collected from different areas of Guangxi. Rabies virus was isolated from 42 samples. RT-PCR was used to amplify a 455 nucleotide segment of the 3'-terminal of the N gene. The sequencing data from that segment was used for phylogenetic analysis.

RESULTS

Nucleotide homology comparisons and phylogenetic tree analysis based on this sequence indicated that all the rabies virus isolates from Guangxi belonged to genotype 1 and could be divided into four groups. Groups I, II and IV included 23, 10 and 8 isolates, respectively. These had nucleotide homologies of 97.1-100%, 98.2-100% and 99.1-99.6%, respectively. Only the GXN119 strain belonged to group III. Group I had two group-specific mutations: T90N and E110D. Group II had one group-specific mutation of T42S.

CONCLUSIONS

This study showed that rabies virus isolates from Guangxi have a close genetic relationship and topographical distribution.

A high-resolution genetic signature of demographic and spatial expansion in epizootic rabies virus

Emerging pathogens potentially undergo rapid evolution while expanding in population size and geographic range during the course of invasion, yet it is generally difficult to demonstrate how these processes interact. Our analysis of a 30-yr data set covering a large-scale rabies virus outbreak among North American raccoons reveals the long lasting effect of the initial infection wave in determining how viral populations are genetically structured in space. We further find that coalescent-based estimates derived from the genetic data yielded an amazingly accurate reconstruction of the known spatial and demographic dynamics of the virus over time. Our study demonstrates the combined evolutionary and population dynamic processes characterizing the spread of pathogen after its introduction into a fully susceptible host population. Furthermore, the results provide important insights regarding the spatial scale of rabies persistence and validate the use of coalescent approaches for uncovering even relatively complex population histories. Such approaches will be of increasing relevance for understanding the epidemiology of emerging zoonotic diseases in a landscape context.

Identification of novel canine rabies virus clades in the Middle East and North Africa

Four novel phylogenetic clades of canine rabies virus (RABV) variants have been identified in the Middle East and North Africa. The three novel Middle Eastern clades comprise RABV isolates from the borders between Israel and neighbouring countries. The North African clade (Africa 4) comprises four RABV isolates from Egypt and one from Israel. We characterized various RABV lineages antigenically by using a panel of monoclonal antibodies to the nucleoprotein (N) and phylogenetically by analysis of entire N gene sequences. The estimated mean substitution rate for the N gene alignment (2.7x10(-4) substitutions per site per year) is comparable with previous estimates for RABV. The application of a molecular clock indicates the emergence of current canine RABV diversity to have occurred at about the same time (approx. 1870) in the Middle East and Europe, following divergence from established lineages in Africa and Asia.

The evolutionary dynamics of canid and mongoose rabies virus in Southern Africa

Two variants of rabies virus (RABV) currently circulate in southern Africa: canid RABV, mainly associated with dogs, jackals, and bat-eared foxes, and mongoose RABV. To investigate the evolutionary dynamics of these variants, we performed coalescent-based analyses of the G-L inter-genic region, allowing for rate variation among viral lineages through the use of a relaxed molecular clock. This revealed that mongoose RABV is evolving more slowly than canid RABV, with mean evolutionary rates of 0.826 and 1.676 x 10(-3) nucleotide substitutions per site, per year, respectively. Additionally, mongoose RABV exhibits older genetic diversity than canid RABV, with common ancestors dating to 73 and 30 years, respectively, and while mongoose RABV has experienced exponential population growth over its evolutionary history in Africa, populations of canid RABV have maintained a constant size. Hence, despite circulating in the same geographic region, these two variants of RABV exhibit striking differences in evolutionary dynamics which are likely to reflect differences in their underlying ecology.

Introduction: conceptualizing and partitioning the emergence process of zoonotic viruses from wildlife to humans

This introduction provides a telegraphic overview of the processes of zoonotic viral emergence, the intricacies of host-virus interactions, and the distinct role of biological transitions and modifying factors. The process of emergence is conceptualized as two transition stages which are common and required for all disease emergence, (1) human contact with the infectious agent and (2) cross-species transmission of the agent, and two transition stages which are not required for emergence and appear unavailable to many zoonotic pathogens, (3) sustained human-to-human transmission and (4) genetic adaptation to the human host. The latter two transitions are presumably prerequisites for the pandemic emergence of a pathogen. The themes introduced herein are amplified and explored in detail by the contributors to this volume. Each author explores the mechanisms and unique circumstances by which evolution, biology, history, and current context have contrived to drive the emergence of different zoonotic agents by a series of related events; although recognizable similarities exist among the events leading to emergence the details and circumstances are never repetitive.

Use of filter paper (FTA) technology for sampling, recovery and molecular characterisation of rabies viruses

This study evaluates the feasibility of the use of the FTA Gene Guard System (a commercial product consisting of filter paper impregnated with patented chemicals supplied by the Whatman company) for the shipment, storage and detection of RNA rabies viruses by a simplified hemi-nested reverse transcriptase polymerase chain reaction. HnRT-PCR of the rabies virus nucleoprotein gene with specific primers showed that viral RNA extracted from crude infected tissues remained stable after fixation on the filter paper under diverse environmental conditions for at least 35 days. The sequence analysis of the products amplified from five out of the seven known genotypes of Lyssaviruses showed the stability of viral RNA viruses after fixation on the filter paper. Furthermore, the sensitivity of the hnRT-PCR following RNA fixation on the filter paper was equivalent to that of standard hnRT-PCR. In conclusion, the stability of viral RNA and the inactivation of infectivity make the FTA technology useful for the storage, transport, collection and subsequent molecular analysis of viral rabies RNA, facilitating epidemiological investigations in the field.

A molecular epidemiological study targeting the glycoprotein gene of rabies virus isolates from China

A group of 31 rabies viruses (RABVs), recovered primarily from dogs, one deer and one human case, were collected from various areas in China between 1989 and 2006. Complete G gene sequences determined for these isolates indicated identities of nucleotide and amino acid sequences of >or=87% and 93.8%, respectively. Phylogenetic analysis of these and some additional Chinese isolates clearly supported the placement of all Chinese viruses in Lyssavirus genotype 1 and divided all Chinese isolates between four distinct groups (I-IV). Several variants identified within the most commonly encountered group I were distributed according to their geographical origins. A comparison of representative Chinese viruses with other isolates retrieved world-wide indicated a close evolutionary relationship between China group I and II viruses and those of Indonesia while China group III viruses formed an outlying branch to variants from Malaysia and Thailand. China group IV viruses were closely related to several vaccine strains. The predicted glycoprotein sequences of these RABVs variants are presented and discussed with respect to the utility of the anti-rabies biologicals currently employed in China.