Transcriptomic Analysis of mRNA Expression Profiles in the Microglia of Mouse Brains Infected with Rabies Viruses of Varying Virulence

Rabies is a lethal encephalitis caused by the rabies virus (RABV) with a fatality rate near 100% after the onset of clinical symptoms in humans and animals. Microglia are resident immune cells in the central nervous system. Few studies have been conducted on the functional role of microglia in RABV infection. Here, we performed a transcriptomic analysis of mRNA expression profiles in the microglia of mouse brains intracerebrally infected with RABV. We successfully isolated single microglial cells from the mouse brains. The survival rate of dissociated microglial cells was 81.91%–96.7%, and the purity was 88.3%. Transcriptomic analysis revealed 22,079 differentially expressed mRNAs identified in the microglia of mouse brains infected with RABV strains (rRC-HL, GX074, and CVS-24) of varying virulence at 4 and 7 days post-infection (dpi) compared to the control group. The numbers of DEGs versus the control at 4 and 7 dpi in mice infected with rRC-HL, GX074, and CVS-24 were 3622 and 4590, 265 and 4901, and 4079 and 6337. The GO enrichment analysis showed that response to stress, response to external stimulus, regulation of response to stimulus, and immune system process were abundant during RABV infection. The KEGG analysis indicated that the Tlr, Tnf, RIG-I, NOD, NF-κB, MAPK, and Jak-STAT signaling pathways were involved in RABV infection at both 4 and 7 dpi. However, some phagocytosis and cell signal transduction processes, such as endocytosis, p53, phospholipase D, and oxidative phosphorylation signaling pathways, were only expressed at 7 dpi. The involvement of the Tnf and Tlr signaling pathways prompted us to construct a protein–protein interaction (PPI) network of these pathways. The PPI revealed 8 DEGs, including Mmp9, Jun, Pik3r1, and Mapk12. Notably, Il-1b interacted with Tnf and Il-6 with combined scores of 0.973 and 0.981, respectively. RABV causes significant changes in mRNA expression profiles in the microglia in mice. 22,079 differentially expressed mRNAs were identified in the microglia of mice infected with RABV strains of varying virulence at 4 and 7 dpi. The DEGs were evaluated using GO, KEGG, and PPI network analysis. Many immune pathways were up-regulated in RABV-infected groups. The findings will help elucidate the microglial molecular mechanisms of cellular metabolism dysregulated by RABV and may provide important information for investigating RABV pathogenesis and therapeutic methods.