Genetic diversity and molecular evolution of the rabies virus matrix protein gene in China

Identification of the potential matrix protein of invertebrate iridescent virus 6 (IIV6)

Invertebrate iridescent virus 6 (IIV6) is a nucleocytoplasmic virus with a ∼212 kb linear dsDNA genome that encodes 215 putative open reading frames (ORFs). Proteomic analysis has revealed that the IIV6 virion consists of 54 virally encoded proteins. Interactions among the structural proteins were investigated using the yeast two-hybrid system, revealing that the protein of 415R ORF interacts reciprocally with the potential envelope protein 118L and the major capsid protein 274L. This result suggests that 415R might be a matrix protein that plays a role as a bridge between the capsid and the envelope proteins. To elucidate the function of 415R protein, we determined the localization of 415R in IIV6 structure and analyzed the properties of 415R-silenced IIV6. Specific antibodies produced against 415R protein were used to determine the location of the 415R protein in the virion structure. Both western blot hybridization and immunogold electron microscopy analyses showed that the 415R protein was found in virions treated with Triton X-100, which degrades the viral envelope. The 415R gene was silenced by the RNA interference (RNAi) technique. We used gene-specific dsRNA's to target 415R and showed that this treatment resulted in a significant drop in virus titer. Silencing 415R with dsRNA also reduced the transcription levels of other viral genes. These results provide important data on the role and location of IIV6 415R protein in the virion structure. Additionally, these results may also shed light on the identification of the homologs of 415R among the vertebrate iridoviruses.

Generation of Monoclonal Antibodies against Variable Epitopes of the M Protein of Rabies Virus

Rabies virus (RABV), the causative agent of rabies, is highly neurovirulent for warm-blooded animals with a mortality rate of up to 100%. The RABV matrix protein (M) is required for virus particle assembly and budding. However, little is known about antigenic differences in the M protein. In this study, five monoclonal antibodies (mAbs), designated 3B9, 4A1, 2B11, 2C1, and 4B11, against the RABV M protein were generated using a recombinant M protein. All five mAbs reacted with the CVS-11 strain but showed no reactivity against the HEP-Flury strain in indirect immunofluorescence and western blotting. The epitope targeted by these mAbs was further identified by peptide scanning using GST-fused peptides. The ²⁵PPYDDD³⁰ peptide was defined as the minimal linear epitope. Alignment of amino acid sequences and phylogenetic analysis of different RABV strains indicated that the variable epitope ²⁵PPDGDD³⁰ is only present in the HEP-Flury and variant Flury strains of clade III, while the other strains resembling ERA and SRVA9 within the clade had another variable epitope, ²⁵PLDDDD³⁰. A Y27D mutation within the epitope was found among the rest of the RABV strains distributed in different clades. However, a single D28G mutation eliminated the reactivity of these five mAbs. In addition, the mAbs were able to recognize wildtype RABV strain in indirect immunofluorescence and western blotting and detect RABV-infected brain tissue using immunohistochemistry. The newly established mAbs and identified epitope may facilitate future investigations in the structure and function of the M protein and the development of diagnostic methods for the detection of different RABV strains worldwide. Most importantly, the epitope recognized by the mAbs against M protein might serve as a novel target for the development of a vaccine targeting RABV virulent strains.

Modeling the transmission dynamics and control of rabies in China

Human rabies was first recorded in ancient China in about 556 BC and is still one of the major public-health problems in China. From 1950 to 2015, 130,494 human rabies cases were reported in Mainland China with an average of 1977 cases per year. It is estimated that 95% of these human rabies cases are due to dog bites. The purpose of this article is to provide a review about the models, results, and simulations that we have obtained recently on studying the transmission of rabies in China. We first construct a basic susceptible, exposed, infectious, and recovered (SEIR) type model for the spread of rabies virus among dogs and from dogs to humans and use the model to simulate the human rabies data in China from 1996 to 2010. Then we modify the basic model by including both domestic and stray dogs and apply the model to simulate the human rabies data from Guangdong Province, China. To study the seasonality of rabies, in Section 4 we further propose a SEIR model with periodic transmission rates and employ the model to simulate the monthly data of human rabies cases reported by the Chinese Ministry of Health from January 2004 to December 2010. To understand the spatial spread of rabies, in Section 5 we add diffusion to the dog population in the basic SEIR model to obtain a reaction-diffusion equation model and determine the minimum wave speed connecting the disease-free equilibrium to the endemic equilibrium. Finally, in order to investigate how the movement of dogs affects the geographically inter-provincial spread of rabies in Mainland China, in Section 6 we propose a multi-patch model to describe the transmission dynamics of rabies between dogs and humans and use the two-patch submodel to investigate the rabies virus clades lineages and to simulate the human rabies data from Guizhou and Guangxi, Hebei and Fujian, and Sichuan and Shaanxi, respectively. Some discussions are provided in Section 7.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15^th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

Zoonoses account for most recently emerged infectious diseases of humans, although little is known about the evolutionary mechanisms involved in cross-species virus transmission. Understanding the evolutionary patterns and processes that underpin such cross-species transmission is of importance for predicting the spread of zoonotic infections, and hence to their ultimate control. We present a large-scale and detailed reconstruction of the evolutionary history of rabies virus (RABV) in domestic and wildlife animal species. RABV is of particular interest as it is capable of infecting many mammals but, paradoxically, is only maintained in distinct epidemiological cycles associated with animal species from the orders Carnivora and Chiroptera. We show that bat-related RABV and dog-related RABV have experienced very different evolutionary dynamics, and that host jumps are sometimes characterized by significant increases in evolutionary rate. Among Carnivora, the association between RABV and particular host species most likely arose from a combination of the historical human-mediated spread of the virus and jumps into new primary host species. In addition, we show that changes in host species are associated with multiple evolutionary pathways including the occurrence of host-specific parallel evolution. Overall, our data indicate that the establishment of dog-related RABV in new carnivore hosts may only require subtle adaptive evolution.

Spatio-temporal Analysis of the Genetic Diversity of Arctic Rabies Viruses and Their Reservoir Hosts in Greenland

There has been limited knowledge on spatio-temporal epidemiology of zoonotic arctic fox rabies among countries bordering the Arctic, in particular Greenland. Previous molecular epidemiological studies have suggested the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3), but have been limited by a low number of available samples preventing in-depth high resolution phylogenetic analysis of RABVs at that time. However, an improved knowledge of the evolution, at a molecular level, of the circulating RABVs and a better understanding of the historical perspective of the disease in Greenland is necessary for better direct control measures on the island. These issues have been addressed by investigating the spatio-temporal genetic diversity of arctic RABVs and their reservoir host, the arctic fox, in Greenland using both full and partial genome sequences. Using a unique set of 79 arctic RABV full genome sequences from Greenland, Canada, USA (Alaska) and Russia obtained between 1977 and 2014, a description of the historic context in relation to the genetic diversity of currently circulating RABV in Greenland and neighboring Canadian Northern territories has been provided. The phylogenetic analysis confirmed delineation into four major arctic RABV lineages (arctic 1–4) with viruses from Greenland exclusively grouping into the circumpolar arctic-3 lineage. High resolution analysis enabled distinction of seven geographically distinct subclades (3.I – 3.VII) with two subclades containing viruses from both Greenland and Canada. By combining analysis of full length RABV genome sequences and host derived sequences encoding mitochondrial proteins obtained simultaneously from brain tissues of 49 arctic foxes, the interaction of viruses and their hosts was explored in detail. Such an approach can serve as a blueprint for analysis of infectious disease dynamics and virus-host interdependencies. The results showed a fine-scale spatial population structure in Greenland arctic foxes based on mitochondrial sequences, but provided no evidence for independent isolated evolutionary development of RABV in different arctic fox lineages. These data are invaluable to support future initiatives for arctic fox rabies control and elimination in Greenland.

Next to dog-mediated rabies, wildlife rabies continues to pose a public health problem, particularly in the northern hemisphere. Control of this zoonosis at the animal source has been proven the most efficient route to reduction of human rabies burden. Successful elimination of red fox-mediated rabies in Western Europe and parts of North America has demonstrated the viability of wildlife rabies control strategies. In some regions, the epidemiology of wildlife rabies is well understood; this is not the case for arctic rabies, particularly in Greenland. Previous molecular epidemiological studies demonstrated the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3) but were limited by low sample numbers and limited sequence length so as to preclude generation of high resolution phylogenetic analysis. Here, a unique set comprised of 79 complete genome sequences of RABVs from Greenland, Canada, USA (Alaska) and Russia collected over the past four decades was analysed. The use of next generation sequencing (NGS) allowed simultaneous determination of host derived sequences encoding mitochondrial proteins from the same brain tissue of 49 arctic foxes. These sequence data combined with geographical and temporal information permit the study of the genetic diversity and evolution of circulating RABVs in Greenland against the background of reservoir host genetics. The results reveal the existence of a single arctic RABV lineage (arctic-3) in Greenland, which has evolved into multiple distinct variants. These analyses provide an improved knowledge of the evolution of the circulating viruses at the molecular level and a better understanding of the historical perspective of the disease in Greenland compared to other parts of the Arctic. This knowledge will support policy on rabies control in mammalian wildlife reservoirs.

Spatial Temporal Dynamics and Molecular Evolution of Re-Emerging Rabies Virus in Taiwan

Taiwan has been recognized by the World Organization for Animal Health as rabies-free since 1961. Surprisingly, rabies virus (RABV) was identified in a dead Formosan ferret badger in July 2013. Later, more infected ferret badgers were reported from different geographic regions of Taiwan. In order to know its evolutionary history and spatial temporal dynamics of this virus, phylogeny was reconstructed by maximum likelihood and Bayesian methods based on the full-length of glycoprotein (G), matrix protein (M), and nucleoprotein (N) genes. The evolutionary rates and phylogeographic were determined using Beast and SPREAD software. Phylogenetic trees showed a monophyletic group containing all of RABV isolates from Taiwan and it further separated into three sub-groups. The estimated nucleotide substitution rates of G, M, and N genes were between 2.49 × 10⁻⁴–4.75 × 10⁻⁴ substitutions/site/year, and the mean ratio of dN/dS was significantly low. The time of the most recent common ancestor was estimated around 75, 89, and 170 years, respectively. Phylogeographic analysis suggested the origin of the epidemic could be in Eastern Taiwan, then the Formosan ferret badger moved across the Central Range of Taiwan to western regions and separated into two branches. In this study, we illustrated the evolution history and phylogeographic of RABV in Formosan ferret badgers.

Rabies and rabies virus in wildlife in mainland China, 1990-2013

The number of wildlife rabies and wildlife-associated human and livestock rabies cases has increased in recent years, particularly in the southeast and northeast regions of mainland China. To better understand wildlife rabies and its role in human and livestock rabies, we reviewed what is known about wildlife rabies from the 1990s to 2013 in mainland China. In addition, the genetic diversity and phylogeny of available wildlife-originated rabies viruses (RABVs) were analyzed. Several wildlife species carry rabies including the bat, Chinese ferret badger, raccoon dog, rat, fox, and wolf. RABVs have been isolated or detected in the bat, Chinese ferret badger, raccoon dog, Apodemus, deer, and vole. Among them, the bat, Chinese ferret badger, and raccoon dog may play a role in the ecology of lyssaviruses in mainland China. All wildlife-originated RABVs were found to belong to genotype 1 RABV except for a bat-originated Irkut virus isolated in 2012. Several substitutions were found between the glycoprotein of wildlife-originated RABVs and vaccine strains. Whether these substitutions could affect the efficacy of currently used vaccines against infections caused by these wildlife-originated RABVs needs to be investigated further. Phylogenetic analysis showed that RABVs in the bat, Chinese ferret badger, and raccoon dog were distinct from local dog-originated RABVs, and almost all collected wildlife-originated isolates were associated with older China clades II to V, suggesting the possibility of wildlife reservoirs in mainland China through the ages.