Although RABV is widely studied, there are still gaps in understanding its genetic evolution. One of the difficulties in conducting more comprehensive evolutionary studies on this virus is the lack of consensus on the target gene to be analyzed. The P and G genes are the most variable in the RABV genome. Analyzes of the P gene are rare, while the same does not hold true for the G gene, probably due to its antigenicity.
In this study, a partial analysis of the RABV genome was performed, including the P and G genes. Studies by Wu (2002), Conselheiro et al. (2022), Oliveira (2010) and Kobayashi (2007, 2010) analyzed part of the RABV genome, focusing on the N and G genes, alone or in association. Despite advances in genetic sequencing tools, studies like Bonnaud et al. (2019), which analyzed the complete RABV genome of samples of wild virus, comparing RABV variation in transmission between dog and fox species, are still a minority.
Although the mutations identified in this study are not synonymous, they did not persist in subsequent passages. This may be associated with the high rate of negative selection of RABV and the high adaptation of CVS-31 to the experimental model used (mice). It is worth noting that the CVS-31 sample is a standard virus sample, therefore more adapted than wild virus samples, where mutations can occur more frequently to achieve more efficient replication.
For a more comprehensive understanding of how RABV evolves, it is important to understand the various factors that can contribute to greater or lesser variation of the virus. Heterogeneity in RABV, as in other viruses, is associated with various factors such as the infection period, transmission route, viral load, host immune system, and interaction between virus and host (Batista et al., 2015).
In this study, the intracerebral route of inoculation was used, which is the most effective in experimental infection. This strategy, however, may have limited virus modulation, as there is no need for the virus to be transported from one tissue to the nervous tissue. Different virus entry points can result in different variations or mutations, depending on its path in the host's body and how the immune system reacts, possibly exerting greater selective pressure and, consequently, greater variation (Conselheiro et al., 2022).
The rate of variation in identity of the RABV among different genetic lineages is approximately 2% across the five RABV genes, a low variation mainly when compared to other RNA genome viruses (Sato et al., 2004, 2006; Kobayash, 2007). When analyzing mutations over successive passages, 148 mutations were identified after five successive passages of RABV isolated from foxes and dogs (Bonnaud et al., 2019).
It is known that host adaptation to the virus is an important factor for RABV alteration. Samples used by Bonnaud et al. were wild viruses, which may explain the high number of mutations, contrasting the results obtained in this study with the standard virus sample CVS-31, where only 3 mutations were identified despite a higher number of passages (10).
In the research conducted by Kissi et al. (1999), successive passages of RABV were carried out in different animal species. Inoculations in mice were not enough to identify virus variations after three passages, while in dogs and cats, 12 mutations were identified in the N, P, and mainly G genes after three successive passages.
A study conducted by our team analyzed the genetic stability of RABV samples with different genetic lineages already identified in Brazil. In three of the RABV isolates, no genetic alterations were identified after ten successive passages. In only one of the analyzed genetic lineages, a non-synonymous mutation in the P gene was identified, indicating that RABV could be less adapted to the non-hematophagous bat (Eptesicus furinalis) in comparison with the other samples analyzed (Batista et al., 2015)..
Mice are widely used as an experimental model for RABV (Fagundes & Taha, 2004). In the present study, three mutations in two RABV genes were identified over the ten successive passages performed. It is possible that a greater number of mutations could have occurred if another animal had been used for these passages. Additionally, the complete genome of the virus was not analyzed, so mutations that were not identified may have occurred; however, the most variable genes of the genome were analyzed in this study.
It is worth noting that the mutations identified were non-synonymous, indicating that the protein structure could be altered. Depending on the position of the amino acid in a protein, it can be involved in one or more functions within its biological activity.
The first 19 amino acids represent the signal peptides (SP), indicating the beginning of the protein and assisting in its entry into the rough endoplasmic reticulum. The region between amino acids 20 and 458 is called the ectodomain (EC), primarily responsible for the pathogenicity of RABV (Carnieli et al. 2009). Mutations identified in the G gene in our study occurred at position G20, in the SP region, and G249, in the EC region.
Certain amino acids, such as 36, 37, and 38, are involved in protein glycosylation. Amino acids between 189 and 214 are involved in RABV internalization, while amino acids between 102 and 179 are involved in virus G binding to receptors at a low pH. Considering the positions of the identified alterations in this study, there was likely no significant change in protein activity, as the mutations were not maintained in subsequent passages, and the lethality remained at 100% in all passages. This lack of change may be due to the fact that CVS-31 is already efficient in infecting and killing mice; therefore, if any alteration occurred, it would likely be a reduction in pathogenicity, which did not happen over the passages.
The study by Okada et al. (2016) was used to understand the activity regions of P. It comprises 297 amino acids, and its translation can produce both the full P gene (P1) and the tPs (P2, P3, P4, and P5), which are important inhibitors of cellular antiviral actions, recognized by interferons and involved in RABV neuroinvasiveness.
These tPs are translations of parts of the P gene. The mutation found in this study occurred in amino acid P223, which plays a role in tP production and acts in virus neuroinvasiveness, but this alteration did not persist in subsequent passages.
The results of this study confirm the genetic stability of the standard virus CVS-31. This result should be considered for the development of laboratory techniques and vaccines against rabies, as its immunogenic and pathogenic capacity was not altered after several successive passages.
Specific mutations in the P and G genes of RABV CVS-31 were identified over ten successive passages in mice, indicating an evolutionary dynamic of the virus in response to selective pressures imposed by the host and environment. However, it is important to note that these specific mutations were not maintained in subsequent passages, demonstrating the high rate of negative selection and efficient adaptation of CVS-31 to the mouse host.
The analysis revealed that despite intracerebral inoculation being an effective strategy to infect mice, allowing RABV multiplication in the central nervous system, it may have limited the virus modulation due to the absence of neurotropism.
Although there are gaps in the comprehensive understanding of RABV evolution, this study contributed to the investigation of specific genetic mutations, as mutations were found in less studied regions of the proteins, providing significant insights into the evolutionary dynamics of the virus in an animal model. However, additional studies are needed to understand the relationship between the identified genetic mutations and the functional changes in the P and G proteins, especially regarding the virus biological activity.