Abstract
Rabies is a fatal zoonotic disease that poses a threat to public health. Rabies virus (RABV) is excreted in the saliva of infected animals, and is primarily transmitted by bite. The role of the salivary glands in virus propagation is significant, but has been less studied in the pathogenic mechanisms of RABV. To identify functionally important genes in the salivary glands, we used RNA sequencing (RNA-seq) to establish and analyze mRNA expression profiles in parotid tissue infected with two RABV strains, CVS-11 and PB4. The biological functions of differentially expressed genes (DEGs) were determined by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, which revealed 3,764 DEGs (678 up-regulated and 3,086 down-regulated) in the CVS-11 infected group and 4,557 DEGs (874 up-regulated and 3,683 down-regulated) in the PB4 infected group. Various biological processes are involved, including the salivary secretion pathway and the phosphatidylinositol 3-kinase-Akt (PI3K-Akt) signaling pathway. This study provides the first mapping of the transcriptome changes in response to RABV infection in parotid tissue, offering new insights into the study of RABV-affected salivary gland function and RABV pathogenic mechanisms in parotid tissue. The salivary gland-enriched transcripts may be potential targets of interest for rabies disease control.