Molecular characterization of nucleoprotein gene of rabies virus from Maharashtra, India

Large-scale phylogenetic analysis reveals genetic diversity and geographic distribution of rabies virus in South-East and South Asia

South-East Asia (SEA) and South Asia (SA) are two important geographic regions with the most severe enzootic rabies in the world. In these regions, phylogenetic analysis of rabies virus (RABV) has been conducted only at a country level; the results obtained from different countries are scattered and unequal, with a non-uniform system to name RABV genotypes. Therefore, it is difficult to undertake origin-tracking and compare inter-country RABV evolution and transmission. To avoid the confusion in understanding and to generate a panoramic picture of RABV genetic diversity, distribution, and transmission in SEA and SA, the present study conducted a systematic phylogenetic analysis by combining all sequences representing 2368 RABV strains submitted to GenBank by 14 rabies endemic SEA and SA countries. The results showed that RABVs circulating in two regions were classified into four major clades and many subclades: the Asia clade is circulating only in SEA, the Indian subcontinent, and Arctic-like clades only in SA, while the Cosmopolitan clade has been detected in both regions. The results also showed a wide range of hosts were infected by divergent RABV subclades, with dogs being the major transmission source. However, wildlife rabies was also found to be an important issue with 6 wild carnivore species identified as potential sources of spillover risk for sylvatic rabies to humans, domestic animals, and other wild animals. Current findings indicate that the two regions have separate virus clades circulating thus indicating the absence of cross-transmission between the regions. The study emphasizes the importance of phylogenetic analysis in the regions using uniform genotyping and naming systems for rabies surveillance, to coordinate actions of member countries to eliminate dog-mediated human rabies by 2030.

Sequencing and sequence analysis of partial nucleoprotein (N) gene and phylogenetic analysis of rabies virus field isolates from Gujarat state, India

The present study was undertaken with an aim of characterization of rabies virus (genus Lyssavirus of the family Rhabdoviridae under the order Mononegavirales) by sequencing of partial nucleoprotein (N) gene of rabies virus and phylogenetic analysis to know the genotype and lineage of rabies virus present in Gujarat state of India. A total of 32 samples (18 brain samples and 14 saliva samples) were aseptically collected from live and dead animals (viz. dog, buffalo, cow, goat, donkey and hyena) for rabies virus detection. Out of 32 samples, 24 samples were found positive by Reverse Transcriptase Polymerase Chain Reactions and from these 24 positive samples, 20 samples were selected for sequencing having good concentration of gene product. ClustalW alignment of nucleotide sequences and amino acid sequences of field rabies isolates revealed 95.20-100 and 97.95-100% similarity among themselves, respectively. Multiple sequence alignment of field rabies isolates and reference vaccine strains [Pasteur strain and Challenge Virus Strain (CVS)] indicated single nucleotide mutations at total 91 positions and amino acid mutations at total 17 different positions. Phylogenetic analysis of N gene sequences using our 20 field rabies isolates and 21 other reported isolates in Genbank resulted in 3 phylogenetic clusters. All the field rabies isolates showed same genetic lineage among themselves and with other earlier reported Indian rabies isolates placing them in Arctic like lineage of Genotype 1 Rabies virus. However, they were at genetic distance with reference Pasteur and CVS strains, which grouped in different phylogenetic cluster.

Codon usage bias in the N gene of rabies virus

Since its emergence, rabies virus (RABV) has been a major worldwide concern especially in developing countries. The nucleoprotein (N) of RABV is highly conserved and key for genetic typing, thus a better understanding of the N gene evolutionary trajectory can assist the development of control measures. We found that the N gene of RABV has a low codon usage bias with a mean effective number of codons (ENC) value of 56.33 influenced by both mutation pressure and natural selection. However, neutrality analysis indicated that natural selection dominates over mutation pressure. Additionally, we found that dinucleotide bias partly contributed to RABV codon usage bias. On the other hand, based on the clades of phylogenetic tree, we found that the evolutionary rate of the Africa 2 clade was the highest with a mean value of 3.75×10^-3 substitutions per site per year. Above all, our results regarding N gene of RABV codon usage will serve future RABV evolution research.