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

Tick-borne viruses and their risk to public health in the Caribbean: Spotlight on bats as reservoirs in Cuba

In recent decades, tick-borne diseases (TBDs) have surged and expanded globally due to factors like changes in human activities, land use patterns, and climate change, and it have been associated with the emergence of zoonotic diseases. Cuba faces the impact of ticks on human health and the economy. Although Cuba has studied TBDs extensively for the past 50 years, focus on tick-borne viral pathogens affecting humans remains scant. Despite TBDs not currently being a major health concern in Cuba, factors like inadequate clinician awareness, climate conditions, global tick emergence, and evidence of zoonotic pathogens in ticks underscore the importance of enhanced TBD surveillance in the country. Here we revised the available information on ticks as vectors of pathogenic viruses to humans, spotlighting bats as potential reservoirs of tick-borne viruses (TBVs). Ticks on bats have gained interest as potential reservoirs of pathogenic viruses to humans in Cuba and worldwide. Understanding their role in maintaining viruses and their potential transmission to humans is crucial for the implementation of surveillance and control programs to reduce the risk of tick-borne viral diseases and public health management.

Phylogeographic Aspects of Bat Lyssaviruses in Europe: A Review

During the last few decades, bat lyssaviruses have become the topic of intensive molecular and epidemiological investigations. Since ancient times, rhabdoviruses have caused fatal encephalitis in humans which has led to research into effective strategies for their eradication. Modelling of potential future cross-species virus transmissions forms a substantial component of the recent infection biology of rabies. In this article, we summarise the available data on the phylogeography of both bats and lyssaviruses in Europe and the adjacent reg ions, especially in the contact zone between the Palearctic and Ethiopian realms. Within these zones, three bat families are present with high potential for cross-species transmission and the spread of lyssaviruses in Phylogroup II to Europe (part of the western Palearctic). The lack of effective therapies for rabies viruses in Phylogroup II and the most divergent lyssaviruses generates impetus for additional phylogenetic and virological research within this geographical region.

Complete Genome Sequence of Matlo Bat Lyssavirus

The genus Lyssavirus includes rabies virus as well as multiple diverse and recently described novel species. Using next-generation sequencing technologies, we have obtained the whole-genome sequence of Matlo bat lyssavirus, which was isolated from a Natal long-fingered bat (Miniopterus natalensis) in South Africa.

The genus Lyssavirus includes rabies virus as well as multiple diverse and recently described novel species. Using next-generation sequencing technologies, we have obtained the whole-genome sequence of Matlo bat lyssavirus, which was isolated from a Natal long-fingered bat (Miniopterus natalensis) in South Africa.

Development of Monoclonal Antibodies for Detection of Conserved and Variable Epitopes of Large Protein of Rabies Virus

Rabies virus (RABV) causes fatal neurological encephalitis and results in approximately 6000 human death cases worldwide every year. The large (L) protein of RABV, possessing conserved domains, is considered as the target for detection. In this study, three monoclonal antibodies (mAbs), designated as 3F3, 3A6 and L-C, against L protein were generated by using the recombinant truncated L protein (aa 1431-1754) and the epitopes were also identified using a series of overlapping truncated polypeptides for testing the reactivity of mAbs with different RABV strains. The ¹⁴⁷⁹EIFSIP¹⁴⁸⁴, ¹⁶⁵⁹RALSK¹⁶⁶³ and ¹⁷²⁴VFNSL¹⁷²⁸ were identified as the minimal linear epitopes recognized by mAbs 3F3, 3A6 and L-C, respectively. Amino acid alignment showed epitope ¹⁷²⁴VFNSL¹⁷²⁸ recognized by mAb L-C is completely conserved among RABV strains, indicating that mAb L-C could be used to detect all of the RABV strains. Epitope ¹⁴⁷⁹EIFSIP¹⁴⁸⁴ is highly conserved among RABV strains except for a P1484S substitution in a China I sub-lineage strain of Asian lineage, which eliminated the reactivity of the epitope with mAb 3F3. However, the epitope ¹⁶⁵⁹RALSK¹⁶⁶³ was only completely conserved in the Africa-2 and Indian lineages, and a single A1660T substitution, mainly appeared in strains of the China I belonging to Asian lineage and a Cosmopolitan lineage strain, still retained the reactivity of the epitope with mAb 3A6. While both A1660T and K1663R substitutions in a China I lineage strain, single K1663R/Q substitution in some China II strains of Asian lineage and some Arctic-like lineage strains and R1659Q mutation in a strain of Africa-3 lineage eliminated the reactivity of the epitope with mAb 3A6, suggesting mAb 3A6 could be used for differentiation of variable epitopes of some strains in different lineages. Thus, variability and conservation of the three epitopes of L protein showed the reactive difference of mAbs among RABV strains of different lineages. These results may facilitate future studies in development of detection methods for RABV infection, the structure and function of RABV L protein.

Whole genome sequence phylogenetic analysis of four Mexican rabies viruses isolated from cattle

Phylogenetic analysis of the rabies virus in molecular epidemiology has been traditionally performed on partial sequences of the genome, such as the N, G, and P genes; however, that approach raises concerns about the discriminatory power compared to whole genome sequencing. In this study we characterized four strains of the rabies virus isolated from cattle in Querétaro, Mexico by comparing the whole genome sequence to that of strains from the American, European and Asian continents. Four cattle brain samples positive to rabies and characterized as AgV11, genotype 1, were used in the study. A cDNA sequence was generated by reverse transcription PCR (RT-PCR) using oligo dT. cDNA samples were sequenced in an Illumina NextSeq 500 platform. The phylogenetic analysis was performed with MEGA 6.0. Minimum evolution phylogenetic trees were constructed with the Neighbor-Joining method and bootstrapped with 1000 replicates. Three large and seven small clusters were formed with the 26 sequences used. The largest cluster grouped strains from different species in South America: Brazil, and the French Guyana. The second cluster grouped five strains from Mexico. A Mexican strain reported in a different study was highly related to our four strains, suggesting common source of infection. The phylogenetic analysis shows that the type of host is different for the different regions in the American Continent; rabies is more related to bats. It was concluded that the rabies virus in central Mexico is genetically stable and that it is transmitted by the vampire bat Desmodus rotundus.

First detection of European bat lyssavirus type 2 (EBLV-2) in Norway

BACKGROUND

In Europe, bat rabies is primarily attributed to European bat lyssavirus type 1 (EBLV-1) and European bat lyssavirus type 2 (EBLV-2) which are both strongly host-specific. Approximately thirty cases of infection with EBLV-2 in Daubenton's bats (Myotis daubentonii) and pond bats (M. dasycneme) have been reported. Two human cases of rabies caused by EBLV-2 have also been confirmed during the last thirty years, while natural spill-over to other non-flying mammals has never been reported. Rabies has never been diagnosed in mainland Norway previously.

CASE PRESENTATION

In late September 2015, a subadult male Daubenton's bat was found in a poor condition 800 m above sea level in the southern part of Norway. The bat was brought to the national Bat Care Centre where it eventually displayed signs of neurological disease and died after two days. EBLV-2 was detected in brain tissues by polymerase chain reaction (PCR) followed by sequencing of a part of the nucleoprotein gene, and lyssavirus was isolated in neuroblastoma cells.

CONCLUSIONS

The detection of EBLV-2 in a bat in Norway broadens the knowledge on the occurrence of this zoonotic agent. Since Norway is considered free of rabies, adequate information to the general public regarding the possibility of human cases of bat-associated rabies should be given. No extensive surveillance of lyssavirus infections in bats has been conducted in the country, and a passive surveillance network to assess rabies prevalence and bat epidemiology is highly desired.

Genetic analysis of a rabies virus host shift event reveals within-host viral dynamics in a new host

Host shift events play an important role in epizootics as adaptation to new hosts can profoundly affect the spread of the disease and the measures needed to control it. During the late 1990s, an epizootic in Turkey resulted in a sustained maintenance of rabies virus (RABV) within the fox population. We used Bayesian inferences to investigate whole genome sequences from fox and dog brain tissues from Turkey to demonstrate that the epizootic occurred in 1997 (±1 year). Furthermore, these data indicated that the epizootic was most likely due to a host shift from locally infected domestic dogs, rather than an incursion of a novel fox or dog RABV. No evidence was observed for genetic adaptation to foxes at consensus sequence level and dN/dS analysis suggested purifying selection. Therefore, the deep sequence data were analysed to investigate the sub-viral population during a host shift event. Viral heterogeneity was measured in all RABV samples; viruses from the early period after the host shift exhibited greater sequence variation in comparison to those from the later stage, and to those not involved in the host shift event, possibly indicating a role in establishing transmission within a new host. The transient increase in variation observed in the new host species may represent virus replication within a new environment, perhaps due to increased replication within the CNS, resulting in a larger population of viruses, or due to the lack of host constraints present in the new host reservoir.

Complete Genome Sequence of Lleida Bat Lyssavirus

All lyssaviruses (family Rhabdoviridae) cause the disease rabies, an acute progressive encephalitis for which, once symptoms occur, there is no effective cure. Using next-generation sequencing, the full-genome sequence for a novel lyssavirus, Lleida bat lyssavirus (LLEBV), from the original brain of a common bent-winged bat has been confirmed.

The Recently Discovered Bokeloh Bat Lyssavirus: Insights Into Its Genetic Heterogeneity and Spatial Distribution in Europe and the Population Genetics of Its Primary Host

In 2010, a novel lyssavirus named Bokeloh bat lyssavirus (BBLV) was isolated from a Natterer's bat (Myotis nattereri) in Germany. Two further viruses were isolated in the same country and in France in recent years, all from the same bat species and all found in moribund or dead bats. Here we report the description and the full-length genome sequence of five additional BBLV isolates from Germany (n=4) and France (n=1). Interestingly, all of them were isolated from the Natterer's bat, except one from Germany, which was found in a common Pipistrelle bat (Pipistrellus pipistrellus), a widespread and abundant bat species in Europe. The latter represents the first case of transmission of BBLV to another bat species. Phylogenetic analysis clearly demonstrated the presence of two different lineages among this lyssavirus species: lineages A and B. The spatial distribution of these two lineages remains puzzling, as both of them comprised isolates from France and Germany; although clustering of isolates was observed on a regional scale, especially in Germany. Phylogenetic analysis based on the mitochondrial cytochrome b (CYTB) gene from positive Natterer's bat did not suggest a circulation of the respective BBLV sublineages in specific Natterer's bat subspecies, as all of them were shown to belong to the M. nattereri sensu stricto clade/subspecies and were closely related (German and French positive bats). At the bat host level, we demonstrated that the distribution of BBLV at the late stage of the disease seems large and massive, as viral RNA was detected in many different organs.

Quantification of Lyssavirus-Neutralizing Antibodies Using Vesicular Stomatitis Virus Pseudotype Particles

Rabies is a highly fatal zoonotic disease which is primarily caused by rabies virus (RABV) although other members of the genus Lyssavirus can cause rabies as well. As yet, 14 serologically and genetically diverse lyssaviruses have been identified, mostly in bats. To assess the quality of rabies vaccines and immunoglobulin preparations, virus neutralization tests with live RABV are performed in accordance with enhanced biosafety standards. In the present work, a novel neutralization test is presented which takes advantage of a modified vesicular stomatitis virus (VSV) from which the glycoprotein G gene has been deleted and replaced by reporter genes. This single-cycle virus was trans-complemented with RABV envelope glycoprotein. Neutralization of this pseudotype virus with RABV reference serum or immune sera from vaccinated mice showed a strong correlation with the rapid fluorescent focus inhibition test (RFFIT). Importantly, pseudotype viruses containing the envelope glycoproteins of other lyssaviruses were neutralized by reference serum to a significantly lesser extent or were not neutralized at all. Taken together, a pseudotype virus system has been successfully developed which allows the safe, fast, and sensitive detection of neutralizing antibodies directed against different lyssaviruses.

Chimeric Filoviruses for Identification and Characterization of Monoclonal Antibodies

Recent experiments suggest that some glycoprotein (GP)-specific monoclonal antibodies (MAbs) can protect experimental animals against the filovirus Ebola virus (EBOV). There is a need for isolation of MAbs capable of neutralizing multiple filoviruses. Antibody neutralization assays for filoviruses frequently use surrogate systems such as the rhabdovirus vesicular stomatitis Indiana virus (VSV), lentiviruses or gammaretroviruses with their envelope proteins replaced with EBOV GP or pseudotyped with EBOV GP. It is optimal for both screening and in-depth characterization of newly identified neutralizing MAbs to generate recombinant filoviruses that express a reporter fluorescent protein in order to more easily monitor and quantify the infection. Our study showed that unlike neutralization-sensitive chimeric VSV, authentic filoviruses are highly resistant to neutralization by MAbs. We used reverse genetics techniques to replace EBOV GP with its counterpart from the heterologous filoviruses Bundibugyo virus (BDBV), Sudan virus, and even Marburg virus and Lloviu virus, which belong to the heterologous genera in the filovirus family. This work resulted in generation of multiple chimeric filoviruses, demonstrating the ability of filoviruses to tolerate swapping of the envelope protein. The sensitivity of chimeric filoviruses to neutralizing MAbs was similar to that of authentic biologically derived filoviruses with the same GP. Moreover, disabling the expression of the secreted GP (sGP) resulted in an increased susceptibility of an engineered virus to the BDBV52 MAb isolated from a BDBV survivor, suggesting a role for sGP in evasion of antibody neutralization in the context of a human filovirus infection.

IMPORTANCE

The study demonstrated that chimeric rhabdoviruses in which G protein is replaced with filovirus GP, widely used as surrogate targets for characterization of filovirus neutralizing antibodies, do not accurately predict the ability of antibodies to neutralize authentic filoviruses, which appeared to be resistant to neutralization. However, a recombinant EBOV expressing a fluorescent protein tolerated swapping of GP with counterparts from heterologous filoviruses, allowing high-throughput screening of B cell lines to isolate MAbs of any filovirus specificity. Human MAb BDBV52, which was isolated from a survivor of BDBV infection, was capable of partially neutralizing a chimeric EBOV carrying BDBV GP in which expression of sGP was disabled. In contrast, the parental virus expressing sGP was resistant to the MAb. Thus, the ability of filoviruses to tolerate swapping of GP can be used for identification of neutralizing MAbs specific to any filovirus and for the characterization of MAb specificity and mechanism of action.

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.

Evolution of genome size and complexity in the rhabdoviridae

RNA viruses exhibit substantial structural, ecological and genomic diversity. However, genome size in RNA viruses is likely limited by a high mutation rate, resulting in the evolution of various mechanisms to increase complexity while minimising genome expansion. Here we conduct a large-scale analysis of the genome sequences of 99 animal rhabdoviruses, including 45 genomes which we determined de novo, to identify patterns of genome expansion and the evolution of genome complexity. All but seven of the rhabdoviruses clustered into 17 well-supported monophyletic groups, of which eight corresponded to established genera, seven were assigned as new genera, and two were taxonomically ambiguous. We show that the acquisition and loss of new genes appears to have been a central theme of rhabdovirus evolution, and has been associated with the appearance of alternative, overlapping and consecutive ORFs within the major structural protein genes, and the insertion and loss of additional ORFs in each gene junction in a clade-specific manner. Changes in the lengths of gene junctions accounted for as much as 48.5% of the variation in genome size from the smallest to the largest genome, and the frequency with which new ORFs were observed increased in the 3' to 5' direction along the genome. We also identify several new families of accessory genes encoded in these regions, and show that non-canonical expression strategies involving TURBS-like termination-reinitiation, ribosomal frame-shifts and leaky ribosomal scanning appear to be common. We conclude that rhabdoviruses have an unusual capacity for genomic plasticity that may be linked to their discontinuous transcription strategy from the negative-sense single-stranded RNA genome, and propose a model that accounts for the regular occurrence of genome expansion and contraction throughout the evolution of the Rhabdoviridae.

Analysis of adaptive evolution in Lyssavirus genomes reveals pervasive diversifying selection during species diversification

Lyssavirus is a diverse genus of viruses that infect a variety of mammalian hosts, typically causing encephalitis. The evolution of this lineage, particularly the rabies virus, has been a focus of research because of the extensive occurrence of cross-species transmission, and the distinctive geographical patterns present throughout the diversification of these viruses. Although numerous studies have examined pattern-related questions concerning Lyssavirus evolution, analyses of the evolutionary processes acting on Lyssavirus diversification are scarce. To clarify the relevance of positive natural selection in Lyssavirus diversification, we conducted a comprehensive scan for episodic diversifying selection across all lineages and codon sites of the five coding regions in lyssavirus genomes. Although the genomes of these viruses are generally conserved, the glycoprotein (G), RNA-dependent RNA polymerase (L) and polymerase (P) genes were frequently targets of adaptive evolution during the diversification of the genus. Adaptive evolution is particularly manifest in the glycoprotein gene, which was inferred to have experienced the highest density of positively selected codon sites along branches. Substitutions in the L gene were found to be associated with the early diversification of phylogroups. A comparison between the number of positively selected sites inferred along the branches of RABV population branches and Lyssavirus intespecies branches suggested that the occurrence of positive selection was similar on the five coding regions of the genome in both groups.

Molecular characterization of cryptically circulating rabies virus from ferret badgers, Taiwan

After the last reported cases of rabies in a human in 1959 and a nonhuman animal in 1961, Taiwan was considered free from rabies. However, during 2012-2013, an outbreak occurred among ferret badgers in Taiwan. To examine the origin of this virus strain, we sequenced 3 complete genomes and acquired multiple rabies virus (RABV) nucleoprotein and glycoprotein sequences. Phylogeographic analyses demonstrated that the RABV affecting the Taiwan ferret badgers (RABV-TWFB) is a distinct lineage within the group of lineages from Asia and that it has been differentiated from its closest lineages, China I (including isolates from Chinese ferret badgers) and the Philippines, 158-210 years ago. The most recent common ancestor of RABV-TWFB originated 91-113 years ago. Our findings indicate that RABV could be cryptically circulating in the environment. An understanding of the underlying mechanism might shed light on the complex interaction between RABV and its host.

Phylogenetic and epidemiologic evidence of multiyear incubation in human rabies

Eight years after emigrating from Brazil, an otherwise healthy man developed rabies. An exposure prior to immigration was reported. Genetic analysis revealed a canine rabies virus variant found only in the patient's home country, and the patient had not traveled internationally since immigrating to the United States. We describe how epidemiological, phylogenetic, and viral sequencing data provided confirmation that rabies encephalomyelitis may present after a long, multiyear incubation period, a consideration that previously has been hypothesized without the ability to exclude a more recent exposure. Accordingly, rabies should be considered in the diagnosis of any acute encephalitis, myelitis, or encephalomyelitis.

Novel lyssavirus in bat, Spain

A new tentative lyssavirus, Lleida bat lyssavirus, was found in a bent-winged bat (Miniopterus schreibersii) in Spain. It does not belong to phylogroups I or II, and it seems to be more closely related to the West Causasian bat virus, and especially to the Ikoma lyssavirus.

Analysis of the complete genome of the first Irkut virus isolate from China: comparison across the Lyssavirus genus

The genome of Irkut virus, isolate IRKV-THChina12, the first non-rabies lyssavirus from China (of bat origin), has been completely sequenced. In general, coding and non-coding regions of this viral genome are similar to those of other lyssaviruses. However, alignment of the deduced amino acid sequences of the structural proteins of IRKV-THChina12 with those of other lyssavirus representatives revealed significant variability between viral species. The nucleoprotein and matrix protein were found to be the most conserved, followed by the large protein, glycoprotein and phosphoprotein. Differences in the antigenic sites in glycoprotein may result in only partial protection of the available rabies biologics against Irkut virus, which is of particular concern for pre- and post-exposure rabies prophylaxis.

Next generation sequencing of viral RNA genomes

Background

With the advent of Next Generation Sequencing (NGS) technologies, the ability to generate large amounts of sequence data has revolutionized the genomics field. Most RNA viruses have relatively small genomes in comparison to other organisms and as such, would appear to be an obvious success story for the use of NGS technologies. However, due to the relatively low abundance of viral RNA in relation to host RNA, RNA viruses have proved relatively difficult to sequence using NGS technologies. Here we detail a simple, robust methodology, without the use of ultra-centrifugation, filtration or viral enrichment protocols, to prepare RNA from diagnostic clinical tissue samples, cell monolayers and tissue culture supernatant, for subsequent sequencing on the Roche 454 platform.

Results

As representative RNA viruses, full genome sequence was successfully obtained from known lyssaviruses belonging to recognized species and a novel lyssavirus species using these protocols and assembling the reads using de novo algorithms. Furthermore, genome sequences were generated from considerably less than 200 ng RNA, indicating that manufacturers’ minimum template guidance is conservative. In addition to obtaining genome consensus sequence, a high proportion of SNPs (Single Nucleotide Polymorphisms) were identified in the majority of samples analyzed.

Conclusions

The approaches reported clearly facilitate successful full genome lyssavirus sequencing and can be universally applied to discovering and obtaining consensus genome sequences of RNA viruses from a variety of sources.

Complete genome sequence of rabies virus CVS-24 from China

The entire genome of the mouse-adapted rabies virus strain CVS-24 (challenge virus standard 24), was sequenced. The overall length of the genome was 11,927 nucleotide (nt), comprising a leader sequence of 58 nt, a nucleoprotein (N) gene of 1353 nt, phosphoprotein (P) gene of 894 nt, a matrix protein (M) gene of 609 nt, a glycoprotein (G) gene of 1575 nt, an RNA-dependent RNA polymerase (RdRp, L) gene of 6384 nt and a trailer region of 70 nt. There was a TGAAAAAAA (TG7) consensus sequence at the end of each gene, except the G gene which had an AGAAAAAAA sequence at the end, and the L/trailer region had the sequence CGAAAAAAA. Three were AACAYYYCT consensus start signals close to TG7. The five cistrons were separated by intergenic regions (IGRS) of 2, 5, 5, 24 nt, respectively. Residue 333 of the mature G protein, which is considered to be associated with pathogenicity, was Ala in CVS-24. The topology of the phylogenetic trees generated using N protein sequences suggested that CVS-11 and CVS-N2C have a close relationship to CVS-24.

Isolation of Irkut virus from a Murina leucogaster bat in China

Background and objectives

Bats are recognized as a major reservoir of lyssaviruses; however, no bat lyssavirus has been isolated in Asia except for Aravan and Khujand virus in Central Asia. All Chinese lyssavirus isolates in previous reports have been of species rabies virus, mainly from dogs. Following at least two recent bat-associated human rabies-like cases in northeast China, we have initiated a study of the prevalence of lyssaviruses in bats in Jilin province and their public health implications. A bat lyssavirus has been isolated and its pathogenicity in mice and genomic alignment have been determined.

Results

We report the first isolation of a bat lyssavirus in China, from the brain of a northeastern bat, Murina leucogaster. Its nucleoprotein gene shared 92.4%/98.9% (nucleotide) and 92.2%/98.8% (amino acid) identity with the two known Irkut virus isolates from Russia, and was designated IRKV-THChina12. Following intracranial and intramuscular injection, IRKV-THChina12 produced rabies-like symptoms in adult mice with a short inoculation period and high mortality. Nucleotide sequence analysis showed that IRKV-THChina12 has the same genomic organization as other lyssaviruses and its isolation provides an independent origin for the species IRKV.

Conclusions

We have identified the existence of a bat lyssavirus in a common Chinese bat species. Its high pathogenicity in adult mice suggests that public warnings and medical education regarding bat bites in China should be increased, and that surveillance be extended to provide a better understanding of Irkut virus ecology and its significance for public health.

The Lyssavirus genus presently comprises 12 species and two unapproved species with different antigenic characteristics. Rabies virus is detectable worldwide; Lagos bat virus, Mokola virus, Duvenhage virus, Shimoni bat virus, and Ikoma lyssavirus circulate in Africa; European bat lyssavirus types 1 and 2, Irkut virus, West Caucasian bat virus, and Bokeloh bat lyssavirus are found in Europe; and Australian bat lyssavirus has been isolated in Australia. Only Aravan and Khujand viruses have been identified in central Asia. Bats are recognized as the most important reservoirs of lyssaviruses. In China, all lyssavirus isolates in previous reports have been rabies virus, mainly from dogs; none has been from bats. Recently, however, at least two bat-associated human rabies or rabies-like cases have been reported in northeast China. Therefore, we conducted a search for bat lyssaviruses in Jilin province, close to where the first bat-associated human rabies case was recorded. We isolated a bat lyssavirus, identified as an Irkut virus isolate with high pathogenicity in experimental mice. Our findings suggest that public warnings and medical education regarding bat bites in China should be increased, and that surveillance should be extended to provide a better understanding of Irkut virus ecology and its significance for public health.

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.

Complete genome sequence of Ikoma lyssavirus

Lyssaviruses (family Rhabdoviridae) constitute one of the most important groups of viral zoonoses globally. All lyssaviruses cause the disease rabies, an acute progressive encephalitis for which, once symptoms occur, there is no effective cure. Currently available vaccines are highly protective against the predominantly circulating lyssavirus species. Using next-generation sequencing technologies, we have obtained the whole-genome sequence for a novel lyssavirus, Ikoma lyssavirus (IKOV), isolated from an African civet in Tanzania displaying clinical signs of rabies. Genetically, this virus is the most divergent within the genus Lyssavirus. Characterization of the genome will help to improve our understanding of lyssavirus diversity and enable investigation into vaccine-induced immunity and protection.

Rabies virus vaccines: is there a need for a pan-lyssavirus vaccine?

All members of the lyssavirus genus are capable of causing disease that invariably results in death following the development of clinical symptoms. The recent detection of several novel lyssavirus species across the globe, in different animal species, has demonstrated that the lyssavirus genus contains a greater degree of genetic and antigenic variation than previously suspected. The divergence of species within the genus has led to a differentiation of lyssavirus isolates based on both antigenic and genetic data into two, and potentially a third phylogroup. Critically, from both a human and animal health perspective, current rabies vaccines appear able to protect against lyssaviruses classified within phylogroup I. However no protection is afforded against phylogroup II viruses or other more divergent viruses. Here we review current knowledge regarding the diversity and antigenicity of the lyssavirus glycoprotein. We review the degree of cross protection afforded by rabies vaccines, the genetic and antigenic divergence of the lyssaviruses and potential mechanisms for the development of novel lyssavirus vaccines for use in areas where divergent lyssaviruses are known to circulate, as well as for use by those at occupational risk from these pathogens.

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.

Phylogenetic analysis of partial RNA-polymerase blocks II and III of Rabies virus isolated from the main rabies reservoirs in Brazil

This study describes the results of the sequencing and analysis of segments of Blocks II and III of the RNA polymerase L gene of Rabies virus isolates from different reservoir species of Brazil. The phylogenetic relations of the virus were determined and a variety of species-specific nucleotides were found in the analyzed areas, but the majority of these mutations were found to be synonymous. However, an analysis of the putative amino acid sequences were shown to have some characteristic mutations between some reservoir species of Brazil, indicating that there was positive selection in the RNA polymerase L gene of Rabies virus. On comparing the putative viral sequences obtained from the Brazilian isolates and other Lyssavirus, it was determined that amino acid mutations occurred in low-restriction areas. This study of the L gene of Rabies virus is the first to be conducted with samples of virus isolates from Brazil, and the results obtained will help in the determination of the phylogenetic relations of the virus.

Molecular characterization of China rabies virus vaccine strain

BACKGROUND

Rabies virus (RV), the agent of rabies, can cause a severe encephalomyelitis in several species of mammals, including humans. As a human rabies vaccine strain employed in China, the genetic knowledge of the aG strain has not been fully studied. The main goal of the present study is to amplify the whole genome of aG strain, and genetic relationships between other vaccine strains and wild strains were analyzed.

RESULTS

The entire genome of human rabies virus vaccine strain aG employed in China was sequenced; this is the second rabies virus vaccine strain from China to be fully characterized. The overall organization and the length of the genome were similar to that of other lyssaviruses. The length of aG strain was 11925nt, comprising a leader sequence of 58nt, nucleoprotein (N) gene of 1353nt, phosphoprotein (P) gene of 894 nt, matrix protein (M) gene of 609nt, glycoprotein (G) gene of 1575nt, RNA-dependent RNA polymerase (RdRp,L) gene of 6384nt, and a trailer region of 70 nt. There was TGAAAAAAA (TGA7) consensus sequence in the end of each gene, except AGA7 at the end of G gene. There was AACAYYYCT consensus start signal at the beginning of each gene.

CONCLUSIONS

In this report, we analyzed the full genome of China human rabies vaccine strain aG. Our studies indicated that the genome of aG retained the basic characteristics of RV. At gene level, N was the most conserved among the five coding genes, indicating this gene is the most appropriate for quantitative genotype definition. The phylogenetic analysis of the N indicated the aG strain clustered most closely with Japanese and Russian rabies vaccine strains, suggesting that they may share the same ancestor; also, the aG strain did not share high homology with wild strains isolated from China, making it may not be the best vaccine strain, more research is needed to elucidate the genetic relationship among the RV circulating in China.

Rhabdovirus accessory genes

The Rhabdoviridae is one of the most ecologically diverse families of RNA viruses with members infecting a wide range of organisms including placental mammals, marsupials, birds, reptiles, fish, insects and plants. The availability of complete nucleotide sequences for an increasing number of rhabdoviruses has revealed that their ecological diversity is reflected in the diversity and complexity of their genomes. The five canonical rhabdovirus structural protein genes (N, P, M, G and L) that are shared by all rhabdoviruses are overprinted, overlapped and interspersed with a multitude of novel and diverse accessory genes. Although not essential for replication in cell culture, several of these genes have been shown to have roles associated with pathogenesis and apoptosis in animals, and cell-to-cell movement in plants. Others appear to be secreted or have the characteristics of membrane-anchored glycoproteins or viroporins. However, most encode proteins of unknown function that are unrelated to any other known proteins. Understanding the roles of these accessory genes and the strategies by which rhabdoviruses use them to engage, divert and re-direct cellular processes will not only present opportunities to develop new anti-viral therapies but may also reveal aspects of cellar function that have broader significance in biology, agriculture and medicine.

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.

Rabies virus transcription and replication

Rabies virus (RABV) is a negative-stranded RNA virus. Its genome is tightly encapsidated by the viral nucleoprotein (N) and this RNA-N complex is the template for transcription and replication by the viral RNA-dependent RNA polymerase (L) and its cofactor, the phosphoprotein (P). We present molecular, structural, and cellular aspects of RABV transcription and replication. We first summarize the characteristics and molecular biology of both RNA synthesis processes. We then discuss biochemical and structural data on the viral proteins (N, P, and L) and their interactions with regard to their role in viral transcription and replication. Finally, we review evidence that rabies viral transcription and replication take place in cytoplasmic inclusion bodies formed in RABV-infected cells and discuss the role of this cellular compartmentalization.

Bats and lyssaviruses

Numerous bat species have been identified as important reservoirs of zoonotic viral pathogens. Rabies and rabies-related viruses constitute one of the most important viral zoonoses and pose a significant threat to public health across the globe. Whereas rabies virus (RABV) appears to be restricted to bats of the New World, related lyssavirus species have not been detected in the Americas and have only been detected in bat populations across Africa, Eurasia, and Australia. Currently, 11 distinct species of lyssavirus have been identified, 10 of which have been isolated from bat species and all of which appear to be able to cause encephalitis consistent with that seen with RABV infection of humans. In contrast, whereas lyssaviruses are apparently able to cause clinical disease in bats, it appears that these lyssaviruses may also be able to circulate within bat populations in the absence of clinical disease. This feature of these highly encephalitic viruses, alongside many other aspects of lyssavirus infection in bats, is poorly understood. Here, we review what is known of the complex relationship between bats and lyssaviruses, detailing both natural and experimental infections of these viruses in both chiropteran and nonchiropteran models. We also discuss potential mechanisms of virus excretion, transmission both to conspecifics and spill-over of virus into nonvolant species, and mechanisms of maintenance within bat populations. Importantly, we review the significance of neutralizing antibodies reported within bat populations and discuss the potential mechanisms by which highly neurovirulent viruses such as the lyssaviruses are able to infect bat species in the absence of clinical disease.

Epidemiology of human rabies in South Africa, 1983-2007

Rabies remains a global public health problem but increasingly so in the developing world. Given a lack of awareness, priority and diagnostic capability, very few developing countries, especially in Africa, report on laboratory confirmed human rabies cases. Here we present a retrospective study on the epidemiology of human rabies in Republic of South Africa for a 25-year period, 1983-2007, based on laboratory confirmed cases. The study highlights the role of the domestic dog as a reservoir and vector of rabies and contrasts this to the almost negligible contribution of wildlife vectors to the overall burden of human rabies in dog rabies endemic areas. From the collective data set, epidemiological aspects that include various features of these human rabies cases as well as failures in or towards the treatment of exposures are reported. Molecular analysis of virus isolates did not identify any additional cases of rabies attributed to infection with the Duvenhage, Lagos bat or Mokola or any other rabies-related viruses.

Virus neutralising activity of African fruit bat (Eidolon helvum) sera against emerging lyssaviruses

It is likely that phylogroup 2 lyssaviruses circulate within bat reservoirs. We adapted a pseudotype (pt) neutralisation assay (PNA) to a multiplex format enabling serosurveillance for Lagos bat virus (LBV), Mokola virus (MOKV) and West Caucasian bat virus (WCBV) in a potential reservoir, the African straw-coloured fruit bat, Eidolon helvum. Highly correlated titres were observed between single and multiplex PNAs using ptLBV and ptMOKV (r=0.97, p<0.0001), validating its use for bat serosurveillance. Of the bat serum samples screened 56% neutralised ptLBV, 27% ptMOKV and 1% ptWCBV. Mean VNAb titres were 1:266, 1:35 and 1:7 against ptLBV, ptMOKV and ptWCBV respectively. The high seroprevalence estimates suggest that the infection rate of LBV in E. helvum remains high enough to persist in this species. This supports the hypothesis that LBV is endemic in Ghanaian E. helvum and we speculate that LBV may have co-evolved with African megachiroptera.

Characterization of Durham virus, a novel rhabdovirus that encodes both a C and SH protein

The family Rhabdoviridae is a diverse group of non-segmented, negative-sense RNA viruses that are distributed worldwide and infect a wide range of hosts including vertebrates, invertebrates, and plants. Of the 114 currently recognized vertebrate rhabdoviruses, relatively few have been well characterized at both the antigenic and genetic level; hence, the phylogenetic relationships between many of the vertebrate rhabdoviruses remain unknown. The present report describes a novel rhabdovirus isolated from the brain of a moribund American coot (Fulica americana) that exhibited neurological signs when found in Durham County, North Carolina, in 2005. Antigenic characterization of the virus revealed that it was serologically unrelated to 68 other known vertebrate rhabdoviruses. Genomic sequencing of the virus indicated that it shared the highest identity to Tupaia rhabdovirus (TUPV), and as only previously observed in TUPV, the genome encoded a putative C protein in an overlapping open reading frame (ORF) of the phosphoprotein gene and a small hydrophobic (SH) protein located in a novel ORF between the matrix and glycoprotein genes. Phylogenetic analysis of partial amino acid sequences of the nucleoprotein and polymerase protein indicated that, in addition to TUPV, the virus was most closely related to avian and small mammal rhabdoviruses from Africa and North America. In this report, we present the morphological, pathological, antigenic, and genetic characterization of the new virus, tentatively named Durham virus (DURV), and discuss its potential evolutionary relationship to other vertebrate rhabdoviruses.

Generation of recombinant European bat lyssavirus type 1 and inter-genotypic compatibility of lyssavirus genotype 1 and 5 antigenome promoters

Bat lyssaviruses (Fam. Rhabdoviridae) represent a source for the infection of terrestial mammals and the development of rabies disease. Molecular differences in the replication of bat and non-bat lyssaviruses and their contribution to pathogenicity, however, are unknown. One reason for this is the lack of reverse genetics systems for bat-restricted lyssaviruses. To investigate bat lyssavirus replication and host adaptation, we developed a reverse genetics system for European bat lyssavirus type 1 (EBLV-1; genotype 5). This was achieved by co-transfection of HEK-293T cells with a full-length EBLV-1 genome cDNA and expression plasmids for EBLV-1 proteins, resulting in recombinant EBLV-1 (rEBLV-1). Replication of rEBLV-1 was comparable to that of parental virus, showing that rEBLV-1 is a valid tool to investigate EBLV-1 replication functions. In a first approach, we tested whether the terminal promoter sequences of EBLV-1 are genotype-specific. Although genotype 1 (rabies virus) minigenomes were successfully amplified by EBLV-1 helper virus, in the context of the complete virus, only the antigenome promoter (AGP) sequence of EBLV-1 was replaceable, as indicated by comparable replication of rEBLV-1 and the chimeric virus. These analyses demonstrate that the terminal AGPs of genotype 1 and genotype 5 lyssaviruses are compatible with those of the heterologous genotype.

[The use of passive rabies immunotherapy: from the past to the future]

Rabies is a fatal disease transmitted by infected animals by bite, scratch, licking on broken skin or contamination of mucosis by saliva. The regimen of post-exposure prophylaxis for people not previously vaccinated, that is currently recommended by WHO, consists of a combination of wound cleaning, active immunization and passive immunization when the exposure is of category 3. Most of the products available on the market, in particular human rabies immunoglobulins, highly purified equine rabies immunoglobulins and the derived F(ab')(2) fragments, are now characterized by high potency and safety. Although the interest of passive anti-rabies immunization was first demonstrated in the first half of the 20th century, there is still an inadequate supply of these products to the target populations mostly in developing countries. Therefore, it is urgent to set-up training and information actions for healthcare personnel on the need to use passive immunotherapy and the lack of adverse effects of the related products. For the future, we hope that a scale up of production and a lower price will improve the accessibility to these products. The development of new products based on monoclonal antibodies and molecular biology, and which may be cheaper, is promising.

Human rabies due to lyssavirus infection of bat origin

Rabies is a fatal viral encephalitis and results from infection with viruses belonging to the genus Lyssavirus. Infection usually results from a bite from a dog infected with classical rabies virus. However, a small number of cases result from contact with bats. It is within bats that most lyssavirus variants, referred to as genotypes, are found. The lyssaviruses found in bats have a distinct geographical distribution and are often restricted to specific bat species. Most have been associated with rabies in humans and in some cases spill-over to domestic animals. Many diagnostic techniques are unable to differentiate rabies virus from other genotypes so it is possible that some human and animal cases go unreported. Furthermore, current vaccines have limited efficacy against some genotypes.

Development of a mouse monoclonal antibody cocktail for post-exposure rabies prophylaxis in humans

As the demand for rabies post-exposure prophylaxis (PEP) treatments has increased exponentially in recent years, the limited supply of human and equine rabies immunoglobulin (HRIG and ERIG) has failed to provide the required passive immune component in PEP in countries where canine rabies is endemic. Replacement of HRIG and ERIG with a potentially cheaper and efficacious alternative biological for treatment of rabies in humans, therefore, remains a high priority. In this study, we set out to assess a mouse monoclonal antibody (MoMAb) cocktail with the ultimate goal to develop a product at the lowest possible cost that can be used in developing countries as a replacement for RIG in PEP. Five MoMAbs, E559.9.14, 1112-1, 62-71-3, M727-5-1, and M777-16-3, were selected from available panels based on stringent criteria, such as biological activity, neutralizing potency, binding specificity, spectrum of neutralization of lyssaviruses, and history of each hybridoma. Four of these MoMAbs recognize epitopes in antigenic site II and one recognizes an epitope in antigenic site III on the rabies virus (RABV) glycoprotein, as determined by nucleotide sequence analysis of the glycoprotein gene of unique MoMAb neutralization-escape mutants. The MoMAbs were produced under Good Laboratory Practice (GLP) conditions. Unique combinations (cocktails) were prepared, using different concentrations of the MoMAbs that were capable of targeting non-overlapping epitopes of antigenic sites II and III. Blind in vitro efficacy studies showed the MoMab cocktails neutralized a broad spectrum of lyssaviruses except for lyssaviruses belonging to phylogroups II and III. In vivo, MoMAb cocktails resulted in protection as a component of PEP that was comparable to HRIG. In conclusion, all three novel combinations of MoMAbs were shown to have equal efficacy to HRIG and therefore could be considered a potentially less expensive alternative biological agent for use in PEP and prevention of rabies in humans.

Human mortality from endemic canine rabies is estimated to be 55,000 deaths per year in Africa and Asia, yet rabies remains a neglected disease throughout most of these countries. More than 99% of human rabies cases are caused by infections resulting from a dog-bite injury. In the vast majority of human exposures to rabies, patients require post-exposure prophylaxis (PEP), which includes both passive (rabies immunoglobulin, RIG) and active immunization (rabies vaccine). The number of victims requiring PEP has increased exponentially in recent years, and human and equine RIG (HRIG and ERIG) were not sufficiently available in countries where canine rabies is endemic. Rabies virus-neutralizing monoclonal antibodies (MAbs) of mouse (Mo) origin have been identified as promising alternatives to HRIG and ERIG. We have developed and assessed both in vitro and in vivo unique mouse monoclonal antibody (MoMAb) cocktails, which are highly efficacious. Three novel combinations were shown to have an equal or superior efficacy to HRIG and therefore could be considered a potentially less expensive alternative for passive prophylactic use to prevent the development of rabies in humans, particularly where needed most in developing countries.

Molecular characterization of the complete genome of a street rabies virus isolated in China

In this study, the complete genomic sequence of a rabies virus isolate HN10, recovered from brain tissue of a rabid patient in China, was determined. This is the first Chinese street isolate that has been fully characterized. The overall organization of this virus is typical of that observed for all other rabies viruses. Alignments of amino acid sequences of the phosphoprotein, glycoprotein and large protein of HN10 with those of other rabies viruses were used to examine the extent of conservation of known functional regions. Phylogenetic analysis using either the complete or partial genomic sequence of HN10 determined that this isolate is most closely associated with viruses previously shown to circulate in Guangxi and Hunan provinces. In addition, of all vaccine strains used for comparison, the attenuated Chinese vaccine strain CTN181 is most closely related to HN10.