Detection and Characterisation of SARS-CoV-2 in Eastern Province of Zambia: A Retrospective Genomic Surveillance Study.

Mutations have driven the evolution and development of new variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with potential implications for increased transmissibility, disease severity and vaccine escape among others. Genome sequencing is a technique that allows scientists to read the genetic code of an organism and has become a powerful tool for studying emerging infectious diseases. Here, we conducted a cross-sectional study in selected districts of the Eastern Province of Zambia, from November 2021 to February 2022. We analyzed SARS-CoV-2 samples (n = 76) using high-throughput sequencing. A total of 4097 mutations were identified in 69 SARS-CoV-2 genomes with 47% (1925/4097) of the mutations occurring in the spike protein. We identified 83 unique amino acid mutations in the spike protein of the seven Omicron sublineages (BA.1, BA.1.1, BA.1.14, BA.1.18, BA.1.21, BA.2, BA.2.23 and XT). Of these, 43.4% (36/83) were present in the receptor binding domain, while 14.5% (12/83) were in the receptor binding motif. While we identified a potential recombinant XT strain, the highly transmissible BA.2 sublineage was more predominant (40.8%). We observed the substitution of other variants with the Omicron strain in the Eastern Province. This work shows the importance of pandemic preparedness and the need to monitor disease in the general population.

Current status and molecular epidemiology of rabies virus from different hosts and regions in Malawi.

Although rabies is endemic in Malawi, there have been no studies in which rabies virus was systematically investigated and characterized in multiple animal hosts in that country. In order to provide molecular epidemiological data on rabies virus in Malawi, 683 suspected rabies case reports from 2008 to 2021 were examined, and 46 (dog = 40, cow = 5, and cat = 1) viable rabies-positive brain samples archived at the Central Veterinary Laboratory (CVL), Lilongwe, Malawi, were analyzed genetically. The results showed an increase in the submission of brain samples from 2008 to 2010, with the highest number of submissions observed in 2020. Of the 683 case reports analyzed for the period under review, 38.1% (260/683) (CI: 34.44 - 41.84) were confirmed by direct fluorescent antibody test. Among the confirmed cases, 65.4% (170/260) (CI: 59.23 - 71.09) were canine rabies. Further, phylogenetic analysis revealed that sequences from different animal hosts clustered together within the Africa 1b lineage, suggesting that the strains circulating in livestock are similar to those in domestic dogs. This finding supports the hypothesis that canine rabies is spilling over to livestock and emphasizes the need for further studies to provide data for effective control of rabies in Malawi.

Autochthonous Leishmania infantum in Dogs, Zambia, 2021.

Leishmaniases are neglected tropical diseases of humans and animals. We detected Leishmania infantum in 3 mixed-breed dogs in Zambia that had no travel history outside the country. Our findings suggest presence of and probable emergence of leishmaniasis in Zambia, indicating the need for physicians and veterinarians to consider the disease during diagnosis.

Prevalence and genetic diversity of Shibuyunji virus, a novel tick-borne phlebovirus identified in Zambia.

Tick-borne pathogens are an emerging public health threat worldwide. However, information on tick-borne viruses is scanty in sub-Saharan Africa. Here, by RT-PCR, 363 ticks (Amblyomma, Hyalomma and Rhipicephalus) in the Namwala and Livingstone districts of Zambia were screened for tick-borne phleboviruses (TBPVs). TBPVs (L gene) were detected in 19 (5.2%) Rhipicephalus ticks in Namwala. All the detected TBPVs were Shibuyunji viruses. Phylogenetically, they were closely related to American dog tick phlebovirus. This study highlights the possible role of Rhipicephalus ticks as the main host of Shibuyunji virus and suggests that these viruses may be present outside the area where they were initially discovered.

West Nile Virus in Farmed Crocodiles, Zambia, 2019.

We detected West Nile virus (WNV) nucleic acid in crocodiles (Crocodylus niloticus) in Zambia. Phylogenetically, the virus belonged to lineage 1a, which is predominant in the Northern Hemisphere. These data provide evidence that WNV is circulating in crocodiles in Africa and increases the risk for animal and human transmission.

Molecular detection and characterization of genotype 1 bovine leukemia virus from beef cattle in the traditional sector in Zambia.

Whilst bovine leukemia virus (BLV) causes considerable economic losses to the dairy industry worldwide, information on its molecular epidemiology and economic impact in beef cattle is limited. Here, blood from 880 animals from Zambia's major cattle-rearing provinces was screened for BLV by nested PCR. Positive pools were sequenced and phylogenetically analyzed. The estimated pooled prevalence was 2.1%. All strains belonged to genotype 1 and formed a distinct phylogenetic cluster. The study suggests circulation of genotype 1 BLV in beef cattle in these regions. This is the first report on molecular detection and characterization of BLV from beef cattle in Africa.

Genetic diversity of rabies virus in different host species and geographic regions of Zambia and Zimbabwe.

Rabies is endemic in Zambia and Zimbabwe. The previously investigated strains of rabies virus in central Zambia belong to the Africa 1b lineage, with similar circulating virus strains found in the various tested hosts and regions. However, prior work assessed only limited regions and host species. Thus, this study aimed to more comprehensively determine the genetic diversity of rabies virus across regions of Zambia and Zimbabwe. RNA (n = 76) was extracted from positive direct fluorescent antibody test brain tissues from dog, cow, goat, cat, pig, human, and jackal collected from Zambia and Zimbabwe. The amplicons of the nucleoprotein and glycoprotein genes were obtained from all examined samples by nested RT-PCR and subsequently sequenced. A phylogenetic analysis of the N gene confirmed that all the endemic strains of rabies virus in Zambia and Zimbabwe belong to the Africa 1b lineage. The obtained viral gene sequences were phylogenetically divided into two clusters. Cluster II comprised only Zambian strains. In contrast, cluster I comprised both Zambia and Zimbabwe strains, with strains from Zimbabwe forming a distinct lineage from Zambian strains, implying viral genetic divergence due to geographical barriers. However, no evidence of clustering based on host or region was observed, implying the circulation of similar virus strains occurs in different hosts and regions of Zambia and Zimbabwe. The clustering of rabies virus strains from jackals with those from domestic animals provides evidence of similar virus strains circulating in both wildlife and domestic animals, and that the jackal might be one of the potential reservoirs of rabies virus infection. In this study, no strains circulating in Zimbabwe were detected in Zambia.

Seroprevalence of Rift Valley fever in cattle of smallholder farmers in Kwilu Province in the Democratic Republic of Congo.

Rift Valley fever (RVF) is a zoonotic mosquito-borne disease caused by RVF virus (RVFV) that causes abortions and high mortalities in livestock and is also associated with acute and fatal disease in humans. In the Democratic Republic of Congo (DRC), information on the epidemiology of RVF is limited, particularly among cattle reared by smallholder farmers. This cross-sectional study was conducted to investigate the seroprevalence of RVF in cattle raised by smallholder farmers in Kwilu Province of DRC, which has not yet reported an RVF epidemic. A total of 677 cattle sera were collected from four territories and tested for anti-RVFV antibodies using immunofluorescent assay and enzyme-linked immunosorbent assay. The overall seroprevalence of anti-RVFV IgG was 6.5% (44/677) (95% CI 4.81-8.7). There was a statistically significant difference in the seroprevalence among the territories (χ2 = 28.79, p < 0.001). Territory seroprevalences were as follows: Idiofa 14.08% (95% CI 9.78-19.76), Bulungu 4.14% (95% CI 1.83-8.68), Gungu 3.21% (95% CI 1.41-6.78), and Masi-Manimba 1.19% (95% CI 0.06-7.37). Seroprevalence differed significantly among age categories (p = 0.0017) and ecosystem (p < 0.001). The seroprevalence of animals aged between 1 and 2 years was 20.0% (95% CI 8.4-39.13) and was higher than group aged <1 year, between 2 and 3 years, and > 3 years. Forest area (18.92% (95% CI 12.35-27.7)) had higher seropositivity than savannah area (4.06% (95% CI 2.65-6.12)). Sex difference was not significant (χ2 = 0.14, p = 0.704). These findings indicate that cattle in Kwilu Province had been exposed to RVFV, which represents a significant risk for both livestock and human health.

Pathological and molecular diagnosis of the 2013 African swine fever outbreak in Lusaka, Zambia.

African swine fever (ASF) is a highly contagious and fatal hemorrhagic viral disease of domestic pigs. The disease is widespread in sub-Saharan Africa and has repeatedly been introduced into other continents. The current study describes the diagnostic investigations of a hemorrhagic disease that was reported in pigs in Lusaka (October 2013), Zambia. Necropsy, histopathology, and molecular diagnosis using polymerase chain reaction and sequence analysis confirmed the disease to be ASF. The sequences obtained showed high similarity to previously isolated ASF viruses. Consistent surveillance and rapid diagnosis of the disease is recommended to prevent future outbreaks and economic losses as there is currently no vaccine against the disease.

Trends, patterns and relationship of antimicrobial use and resistance in bacterial isolates tested between 2015-2020 in a national referral hospital of Zambia.

Increased antimicrobial resistance (AMR) among bacteria underscores the need to strengthen AMR surveillance and promote data-based prescribing. To evaluate trends and associations between antimicrobial usage (AMU) and AMR, we explored a dataset of 34,672 bacterial isolates collected between 2015 and 2020 from clinical samples at the University Teaching Hospital (UTH) in Lusaka, Zambia. The most frequently isolated species were Escherichia coli (4,986/34,672; 14.4%), Staphylococcus aureus (3,941/34,672; 11.4%), and Klebsiella pneumoniae (3,796/34,672; 10.9%). Of the 16 drugs (eight classes) tested, only amikacin and imipenem showed good (> 50%) antimicrobial activity against both E. coli and K. pneumoniae, while nitrofurantoin was effective only in E. coli. Furthermore, 38.8% (1,934/4,980) of E. coli and 52.4% (2,079/3,791) of K. pneumoniae isolates displayed multidrug resistance (MDR) patterns on antimicrobial susceptibility tests. Among S. aureus isolates, 44.6% (973/2,181) were classified as methicillin-resistant (MRSA). Notably, all the MRSA exhibited MDR patterns. The annual hospital AMR rates varied over time, while there was a weak positive relationship (r = 0.38, 95% CI = 0.11-0.60) between the monthly use of third-generation cephalosporins (3GCs) and 3GC resistance among Enterobacterales. Overall, the results revealed high AMR rates that fluctuated over time, with a weak positive relationship between 3GC use and resistance. To our knowledge, this is the first report to evaluate the association between AMU and AMR in Zambia. Our results highlight the need to strengthen antimicrobial stewardship programs and optimize AMU in hospital settings.

Antimicrobial Use Survey and Detection of ESBL-Escherichia coli in Commercial and Medium-/Small-Scale Poultry Farms in Selected Districts of Zambia.

Antimicrobial resistance (AMR) among Escherichia coli from food animals is a rising problem, and heavy antimicrobial use in poultry is a contributing factor. In Zambia, studies linking poultry-associated AMR and antibiotic use (AMU) are rare. This study aimed to investigate commercial and medium-/small-scale poultry farmers' usage of antimicrobials based on a questionnaire survey in ten districts of Zambia. In addition, the study characterized extended-spectrum ß-lactamase (ESBL)-producing E. coli isolates obtained from poultry in the same districts. Data regarding knowledge and usage of antimicrobials were collected from commercial and medium-/small-scale poultry farmers using a pre-tested structured questionnaire. At the same time, cloacal samples were collected and analyzed. One hundred and fifty E. coli isolates were tested for antimicrobial susceptibility using eight antibiotic classes. The isolates were further screened for ESBL production by streaking them on cefotaxime (CTX)-supplemented MacConkey agar, then subjecting them to sequencing on a NextSeq. The questionnaire survey showed that more medium-/small-scale than commercial poultry farmers used antimicrobials (OR = 7.70, 95% CI = 2.88-20.61) but less prescriptions (OR = 0.02, 95% CI = 0.00-0.08). Susceptibility testing revealed that resistance was highest to ampicillin (128/148, 86.5%) and tetracycline (101/136, 74.3%) and that the prevalence of multidrug resistance (MDR) (28/30, 93.3%) was high. Whole-genome sequencing (WGS) of eight (8/30, 26.7%) isolates with CTX Minimum Inhibitory Concentration (MIC) ≥ 4 µg/mL revealed the presence of ESBL-encoding genes blaCTX-M-14, blaCTX-M-55, and blaTEM. WGS also detected other AMR genes for quinolones, aminoglycosides, phenicols, tetracycline, macrolides, and folate-pathway antagonists. Altogether, the questionnaire survey results showed a higher proportion of AMU and lower prescription usage among medium-/small-scale farmers. In addition, our results emphasize the circulation of ESBL-producing E. coli strains with associated MDR. It is critical to educate farmers about AMR risks and to encourage responsible usage of antimicrobials. Furthermore, there is a need to strengthen regulations limiting access to antimicrobials. Finally, there is a need to establish a one health system to guide public health response.

Multidrug-Resistant Escherichia coli from Raw Cow Milk in Namwala District, Zambia: Public Health Implications.

Escherichia coli (E. coli), a major foodborne disease-causing pathogen found in raw cow milk, has even far more reaching public health ramifications as it encodes for antimicrobial resistance (AMR). This study aimed to identify multidrug-resistant (MDR) E. coli from raw cow's milk and evaluate their antimicrobial-resistant profiles. In total, 418 pooled raw cow milk samples were collected from milk collection centers and analysed using standard culture methods to isolate E. coli. Antimicrobial Susceptibility Testing (AST) was conducted using the Kirby Bauer disk diffusion method and PCR was used to identify cefotaxime (CTX) resistant genes. Overall isolation of E. coli was 51.2% (214/418) with MDR observed in 21% (45/214) of isolates across different antibiotic combinations. Resistance was observed towards ampicillin (107/214, 50%), tetracycline (86/214, 40.1%), trimethoprim/sulfamethoxazole (61/214, 28.5%), and amoxicillin/clavulanic acid (CTX) (50/214, 23.4%). Notably, 15% (32/214) resistance to CTX was observed, while 12.6% (27/214) exhibited resistance to imipenem. The blaCTX-M and blaTEM genes were detected in CTX-resistant isolates. The findings of MDR E. coli that harbour blaCTX-M and blaTEM genes in raw cow's milk indicate serious public health risks for consumers.

Detection of Old and New World Relapsing Fever Borreliae in Ornithodoros Ticks Collected from Warthog Burrows in Zambia.

Relapsing fever (RF) is an arthropod-borne disease caused by Borrelia spirochete, which is one of the major public health concerns in endemic regions including Africa. However, information on Borrelia spirochetes is limited in Zambia. Here, we investigate the Borrelia spirochetes harbored by Ornithodoros ticks in Zambian National Parks. We analyzed 182 DNA samples pooled from 886 Ornithodoros ticks. Of these, 43 tested positive, and their sequence revealed that the ticks harbored both Old and New World RF borreliae. This research presents the first evidence of Old-World RF borreliae in Zambia. The New World RF borreliae detected herein differed from the Candidatus Borrelia fainii previously reported in Zambia and were closely related to the pathogenic Borrelia sp. VS4 identified in Tanzania. Additionally, Borrelia theileri was recently reported in Zambia. Hence, at least four different Borrelia species occur in Zambia, and the organisms causing relapsing fever there might be more complex than previously thought. We empirically confirmed that real-time PCR with TaqMan minor groove binder probes accurately and simultaneously detected both Old and New World RF. In this manner, they could facilitate quantitative analyses of both types of RF borreliae. Subsequent investigations should endeavor to isolate the aforementioned Borrelia spp. and perform serosurveys on patients with RF.

Prevalence and Genomic Characterization of Rotavirus A from Domestic Pigs in Zambia: Evidence for Possible Porcine-Human Interspecies Transmission.

Rotavirus is a major cause of diarrhea globally in animals and young children under 5 years old. Here, molecular detection and genetic characterization of porcine rotavirus in smallholder and commercial pig farms in the Lusaka Province of Zambia were conducted. Screening of 148 stool samples by RT-PCR targeting the VP6 gene revealed a prevalence of 22.9% (34/148). Further testing of VP6-positive samples with VP7-specific primers produced 12 positives, which were then Sanger-sequenced. BLASTn of the VP7 positives showed sequence similarity to porcine and human rotavirus strains with identities ranging from 87.5% to 97.1%. By next-generation sequencing, the full-length genetic constellation of the representative strains RVA/pig-wt/ZMB/LSK0137 and RVA/pig-wt/ZMB/LSK0147 were determined. Genotyping of these strains revealed a known Wa-like genetic backbone, and their genetic constellations were G4-P[6]-I5-R1-C1-M1-A8-N1-T1-E1-H1 and G9-P[13]-I5-R1-C1-M1-A8-N1-T1-E1-H1, respectively. Phylogenetic analysis revealed that these two viruses might have their ancestral origin from pigs, though some of their gene segments were related to human strains. The study shows evidence of reassortment and possible interspecies transmission between pigs and humans in Zambia. Therefore, the "One Health" surveillance approach for rotavirus A in animals and humans is recommended to inform the design of effective control measures.

Expanding diversity of bunyaviruses identified in mosquitoes.

Mosquitoes interact with various organisms in the environment, and female mosquitoes in particular serve as vectors that directly transmit a number of microorganisms to humans and animals by blood-sucking. Comprehensive analysis of mosquito-borne viruses has led to the understanding of the existence of diverse viral species and to the identification of zoonotic arboviruses responsible for significant outbreaks and epidemics. In the present study on mosquito-borne bunyaviruses we employed a broad-spectrum RT-PCR approach and identified eighteen different additional species in the Phenuiviridae family and also a number of related but unclassified bunyaviruses in mosquitoes collected in Zambia. The entire RNA genome segments of the newly identified viruses were further analyzed by RNA sequencing with a ribonuclease R (RNase R) treatment to reduce host-derived RNAs and enrich viral RNAs, taking advantage of the dsRNA panhandle structure of the bunyavirus genome. All three or four genome segments were identified in eight bunyavirus species. Furthermore, L segments of three different novel viruses related to the Leishbunyaviridae were found in mosquitoes together with genes from the suspected host, the Crithidia parasite. In summary, our virus detection approach using a combination of broad-spectrum RT-PCR and RNA sequencing analysis with a simple virus enrichment method allowed the discovery of novel bunyaviruses. The diversity of bunyaviruses is still expanding and studies on this will allow a better understanding of the ecology of hematophagous mosquitoes.

Insight into the genetic diversity of Mycobacterium bovis isolated from cattle in Malawi.

We describe the genetic diversity and phylogenetic relationships of Mycobacterium bovis, isolated from cattle in Malawi. Deletion analysis, spoligotyping, and MIRU-VNTR typing were used to genotype the isolates. Combined with a larger dataset from neighboring countries, the overall M. bovis diversity in Southern Africa was contextualized. From the southern and northern regions of Malawi, 24 isolates were confirmed as M. bovis. We pooled data for the central region (60 isolates) from our recent publication to conceptualize the genetic and phylogenetic relationships of M. bovis in Malawi. European 1 was the dominant M. bovis clonal complex, with 10 unique spoligotype patterns, and SB0131 was ubiquitous. High genetic diversity, a low clustering rate, and many singletons, coupled with a low mutation transmission index, infer a low level of recent transmission, and suggest an endemic status of bovine tuberculosis (bTB) in Malawi. M. bovis isolates from Zambia, Mozambique, and South Africa were genetically related to Malawian isolates, whereas Tanzanian isolates were distantly related. The diversity and phylogenetic analysis suggest earlier introductions and maintenance of M. bovis by constant reinfection from reservoir animals. These findings are fundamental to understanding the source and route of infection in order to establish alternative management strategies for bTB.

Genomic Analysis of Mycobacterium tuberculosis Strains Resistant to Second-Line Anti-Tuberculosis Drugs in Lusaka, Zambia.

The emergence of pre-extensively drug-resistant tuberculosis (pre-XDR-TB) is a threat to TB control programs in developing countries such as Zambia. Studies in Zambia have applied molecular techniques to understand drug-resistance-associated mutations, circulating lineages and transmission patterns of multi-drug-resistant (MDR) Mycobacterium tuberculosis. However, none has reported genotypes and mutations associated with pre-XDR TB. This study characterized 63 drug-resistant M. tuberculosis strains from the University Teaching Hospital between 2018 and 2019 using targeted gene sequencing and conveniently selected 50 strains for whole genome sequencing. Sixty strains had resistance mutations associated to MDR, one polyresistant, and two rifampicin resistant. Among MDR strains, seven percent (4/60) had mutations associated with pre-XDR-TB. While four, one and nine strains had mutations associated with ethionamide, para-amino-salicylic acid and streptomycin resistances, respectively. All 50 strains belonged to lineage 4 with the predominant sub-lineage 4.3.4.2.1 (38%). Three of four pre-XDR strains belonged to sub-lineage 4.3.4.2.1. Sub-lineage 4.3.4.2.1 strains were less clustered when compared to sub-lineages L4.9.1 and L4.3.4.1 based on single nucleotide polymorphism differences. The finding that resistances to second-line drugs have emerged among MDR-TB is a threat to TB control. Hence, the study recommends a strengthened routine drug susceptibility testing for second-line TB drugs to stop the progression of pre-XDR to XDR-TB and improve patient treatment outcomes.

An unusually long Rift valley fever inter-epizootic period in Zambia: Evidence for enzootic virus circulation and risk for disease outbreak.

Rift valley fever (RVF) is a mosquito-borne disease of animals and humans. Although RVF outbreaks are usually reported at 5-15-year intervals in sub-Saharan Africa, Zambia has experienced an unusually long inter-epizootic/-epidemic period of more than three decades. However, serological evidence of RVF virus (RVFV) infection in domestic ruminants during this period underscores the need for comprehensive investigation of the mechanisms of virus perpetuation and disease emergence. Mosquitoes (n = 16,778) captured from eight of the ten provinces of Zambia between April 2014 and May 2019 were pooled (n = 961) and screened for RVFV genome by a pan-phlebo RT-PCR assay. Aedes mosquito pools (n = 85) were further screened by nested RT-PCR assay. Sera from sheep (n = 13), goats (n = 259) and wild ungulates (n = 285) were screened for RVFV antibodies by ELISA while genome detection in pooled sera (n = 276) from domestic (n = 248) and wild ungulates (n = 37) was performed by real-time RT-PCR assay. To examine the association between the long inter-epizootic period and climatic variables, we examined El Niño-Southern Oscillation indices, precipitation anomalies, and normalized difference vegetation index. We then derived RVF risk maps by exploring climatic variables that would favor emergence of primary RVFV vectors. While no RVFV genome could be detected in pooled mosquito and serum samples, seroprevalence was significantly high (OR = 8.13, 95% CI [4.63-14.25]) in wild ungulates (33.7%; 96/285) compared to domestic ruminants (5.6%; 16/272). Retrospective analysis of RVF epizootics in Zambia showed a positive correlation between anomalous precipitation (La Niña) and disease emergence. On risk mapping, whilst northern and eastern parts of the country were at high risk, domestic ruminant population density was low (< 21 animals/km2) in these areas compared to low risk areas (>21 animals/km2). Besides evidence of silent circulation of RVFV and the risk of disease emergence in some areas, wildlife may play a role in the maintenance of RVFV in Zambia.

Clonal relationship between multidrug-resistant Escherichia coli ST69 from poultry and humans in Lusaka, Zambia.

Multidrug-resistant (MDR) Escherichia coli in food animals such as chickens is an emerging public health concern in Zambia. Additionally, the country's high demand for poultry products necessitates further investigation into the link between poultry and human MDR E. coli. Twenty cefotaxime-resistant E. coli isolates collected from poultry in Lusaka, Zambia, were screened for multidrug resistance and sequenced on MiSeq and MinION platforms. Genomes were assembled de novo and compared with 36 previously reported cefotaxime-resistant E. coli isolates from inpatients at the University Teaching Hospital, Lusaka. All (20/20, 100%) poultry isolates exhibited resistance to ampicillin, chloramphenicol and doxycycline. Phylogenetic analysis and hierarchical clustering showed a high degree of genetic relatedness between E. coli O17:H18-ST69 from poultry and humans. The E. coli O17:H18-ST69 clone accounted for 4/20 (20%) poultry- and 9/36 (25%) human-associated isolates that shared two plasmids harboring 14 antimicrobial resistance (AMR) genes. However, comparison analysis showed that the isolates also had other AMR plasmids distinct for each niche. Our results suggested clonal transmission of MDR E. coli between poultry and humans, with the potential acquisition of niche-specific AMR plasmids. Thus, the control of MDR E. coli requires a One Health approach involving both human and animal health sectors.

Glu333 in rabies virus glycoprotein is involved in virus attenuation through astrocyte infection and interferon responses.

The amino acid residue at position 333 of the rabies virus (RABV) glycoprotein (G333) is a major determinant of RABV pathogenicity. Virulent RABV strains possess Arg333, whereas the attenuated strain HEP-Flury (HEP) possesses Glu333. To investigate the potential attenuation mechanism dependent on a single amino acid at G333, comparative analysis was performed between HEP and HEP333R mutant with Arg333. We examined their respective tropism for astrocytes and the subsequent immune responses in astrocytes. Virus replication and subsequent interferon (IFN) responses in astrocytes infected with HEP were increased compared with HEP333R both in vitro and in vivo. Furthermore, involvement of IFN in the avirulency of HEP was demonstrated in IFN-receptor knockout mice. These results indicate that Glu333 contributes to RABV attenuation by determining the ability of the virus to infect astrocytes and stimulate subsequent IFN responses.