Transcriptome analysis of salivary glands of rabies-virus-infected mice

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.

1 Introduction

Rabies is a zoonotic disease caused by the rabies virus (RABV), which belongs to the genus Lyssavirus of the family Rhabdoviridae (Fooks et al., 2017; Singh et al., 2017). Human rabies is almost invariably fatal once clinical symptoms occur. With a high global death count per year, rabies continues to pose a severe threat to human and animal health (Fooks et al., 2014). Initial symptoms of rabies include generalized symptoms such as fever, pain, and abnormal tingling or burning at the wound. As the virus spreads to the central nervous system (CNS), the disease progresses to fatal encephalitis and encephalomyelitis, leading to the development of the typical clinical signs of furious rabies or paralytic rabies and ultimately death from asphyxia caused by pharyngospasm or respiratory circulatory failure (Mitrabhakdi et al., 2005; Hemachudha et al., 2013; Jackson, 2016).

It has been reported that up to 99% of human rabies cases are caused by a deep bite or scratch from rabid dogs (Fooks et al., 2017). Cases of RABV infection from organ transplantation from infected patients and aerosol transmission are relatively rare (Srinivasan et al., 2005; Maier et al., 2010). Typically, following transmission by a deep bite from an infected animal, RABV replicates in muscle cells and binds to receptors at the neuromuscular junction, initiating early peripheral infection (Lewis et al., 2000). The virus travels through the peripheral nerve bypassing the spinal cord and reaches the CNS by centripetal spread. Once in the brain, RABV replicates on a large scale and spreads centrifugally to the peripheral nervous system, infecting surrounding tissues such as the heart, salivary glands, adrenal glands, gastrointestinal tract and pancreas (Hemachudha et al., 2002; Jogai et al., 2002; Dietzschold et al., 2008; Tobiume et al., 2009). The salivary gland is the primary outlet for the centrifugal spread of RABV and supports viral replication. Previous studies have shown that infected dogs can excrete viral particles in their saliva for up to 14 days before showing obvious clinical signs of rabies (Fekadu and Shaddock, 1984). The salivary gland epithelium contains a high concentration of RABV antigens and can be used as an alternative to brain and cerebrospinal fluid samples for RABV diagnosis (Beauregard and Casey, 1969; Dierks et al., 1969; Charlton et al., 1983).

Saliva is a unique oral bodily fluid. Approximately 90% of saliva is produced by the salivary glands, which mainly include the submandibular gland (SMG), sublingual gland (SLG), and parotid gland (PG) (Zhang et al., 2016a). The secretion of saliva is controlled by the cerebral cortex, and can be influenced by factors such as diet, environment, age, pathogen invasion and salivary gland lesions. Saliva contains various substances, including amylase, lysozyme, peroxidase, mucin, phospholipids, sodium, potassium, calcium, and magnesium, which contribute to its functions of food digestion, antibacterial activity, and protection of the gastric mucosa. Additionally, saliva has the ability to excrete certain pathogens such as RABV (Li et al., 1995) and enteric viruses like murine norovirus (MNV), rotavirus (RV) and astrovirus (AstV) (Kirby et al., 2010; Pisanic et al., 2019; Anfruns-Estrada et al., 2020; Ghosh et al., 2022). Epstein–Barr virus (EBV), herpes simplex virus (HSV), SARS-CoV-2 and human immunodeficiency virus (HIV) can also be excreted or transmitted through the saliva (Corstjens et al., 2000; Lee et al., 2009; Lima et al., 2010; Huang et al., 2021).

In this study, we utilized transcriptomics to characterize the mRNA expression in parotid tissues of mice infected with two RABV strains: the challenge virus standard (CVS)-11 strain and the street rabies virus PB4 strain, and we analyzed the biological functions of the differentially expressed mRNA. The accuracy of the RNA-seq data was confirmed by a quantitative real-time polymerase chain reaction (RT-qPCR). This research represents the first transcriptome analysis of salivary glands in RABV-infected mice and identifies potentially important salivary gland transcripts.

2 Materials and methods

2.1 Experimental animals

The Kunming mice used in this experiment were purchased from Liaoning Changsheng biotechnology co., Led. (License No. SCXK (Liao) 2015–0001). All mice were 3-week-old SPF Kunming mice, weighing approximately 15 g. Male and female mice were equally divided and provided with free access to food and water during the feeding period. A total of 36 mice were split into three groups: mock, CVS-11 and PB4-infected groups, with each group consisting of 12 mice. Within each group, four mice constituted a biological replicate.

2.2 Virus infection

The CVS-11 and PB4 strains used in this study were donated by Changchun Tu of the Diagnostic Laboratory of Rabies and Wildlife-associated Viral Zoonoses of Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences. The PB4 strain was isolated from a rabid pig that was bitten by a rabid dog in Hunan province in 2006 by the National Reference Laboratory for Rabies in China. All virus-related experiments were conducted in a Biosafety Level III (BSL-3) laboratory, following the biosafety regulations of the International Guiding Principles for Biomedical Research Involving Animals. Mice were divided into two experimental groups, and each group was inoculated intramuscularly (into the right thigh muscle) with 100 μL of 10^5.5 TCID50/mL of the CVS-11 strain or PB4 strain. A mock group was also included that was inoculated with an equal amount of Dulbecco’s Modified Eagle’s Medium (DMEM) (Sigma-Aldrich Co. LLC, St Louis, America). When the mice reached the late stage of the disease and died, their PGs were collected for detection and RNA sequencing.

2.3 RNA extraction and quantitative real-time PCR

Parotid tissues were lysed using TRIzol (Takara Biomedical Technology Co, Ltd., Kusatsu, Japan) for total RNA extraction. Subsequently, cDNA was synthesized using the All-In-One 5X RT MasterMix (Applied Biological Materials Inc., Richmond, Canada). RT-qPCR was conducted in a 20 μL reaction volume, consisting of 10 μL FastStart Universal SYBR Green Master Mix (ROX) (F. Hoffmann-La Roche Ltd., Basel, Switzerland), 1 μL PCR primers (containing 10 μM forward and reverse primers, respectively), 4 μL nuclease-free water, and 5 μL cDNA. The reference household gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was utilized. Samples were incubated at 95°C for 10 min to activate DNA polymerase, followed by 40 cycles at 95°C for 15 s and at 60°C for 1 min for amplification and analysis. Primer sequences for the detection of RABV nucleoprotein (RABV N), GAPDH and other corresponding mRNAs are listed in Table 1. Each target was normalized to that of GAPDH using the 2^−ΔΔCt method.

2.4 RNA extraction library construction and sequencing

The quantity of total RNA was analyzed using the Bioanalyzer 2,100 and RNA 6000 Nano LabChip Kit (Agilent, CA, United States, 5067–1,511). High-quality samples with a RIN number > 7.0 were selected to construct the sequencing library. The mRNA was purified from the total RNA and cleaved into short fragments, which were then reverse-transcribed to generate the cDNA. Second-stranded DNAs were synthesized, and dual-index adapters were ligated to the fragments, whose size selection was performed using AMPureXP beads. The ligated products were amplified by PCR. Finally, 2 × 150 bp paired-end sequencing (PE150) was conducted on an Illumina Novaseq™ 6,000 (LC-Bio Technology CO., Ltd., Hangzhou, China).

2.5 Sequencing and filtering of clean reads

The cDNA library was sequenced and run with the Illumina Novaseq™ 6,000 sequencing platform. We employed the Illumina paired-end RNA-seq approach to sequence the transcriptome and generate 2 × 150 bp paired-end reads. The reads obtained from the sequencing machines include raw reads with adapters or low quality bases. A quality analysis was performed before sequencing and clean reads were filtered to check the quality of the sequencing data and rule out possible low-quality reads. We then filtered the reads with Cutadapt¹ to obtain high quality clean reads. We utilized Q30, GC content, low-quality read content and other quality control indicators to evaluate data quality. During this process, we removed reads containing adapters, polyA and polyG, reads containing more than 5% of unknown nucleotides (N) and low quality reads containing more than 20% of low quality (Q-value ≤20) bases. We then verified the per- base sequence quality by using FastQC.² As a result, a total of 47.92GB of cleaned, paired-end reads were produced. We successfully uploaded and published the clean raw sequence data to the NCBI Gene Expression Omnibus (GEO) datasets, under GEO accession number GSE248365.

2.6 Alignment to the reference genome

We aligned the reads of all samples to the reference genomes of Mus musculus (Ensembl release 101) using HISAT2.³ This first removes portions of the reads based on quality information and then maps the reads to the reference genome.

2.7 Quantification of gene abundance

The mapped reads of each sample were assembled with StringTie.⁴ We then merged all transcriptomes of all samples to construct a comprehensive transcriptome using gffcompare software.⁵ After the final transcriptome was generated, we estimated the expression levels of all transcripts and calculated the FPKM (fragment per kilobase of transcript per million mapped reads) value by using StringTie and ballgown⁶ for mRNA expression abundance analysis.

2.8 Differential gene expression analysis

Differential gene expression analysis was conducted between different groups using the DESeq2 software.⁷ Gene count was used as the input chosen for DESeq2. The fold change was calculated as the ratio of the mean expression value of the gene between the experimental groups and the mock groups. The genes with a false discovery rate (FDR) parameter below 0.05 and an absolute fold change (FC) ≥ 2 were considered as differentially expressed genes (DEGs).

2.9 Go and KEGG pathway enrichment analysis

Gene Ontology (GO) is a standardized classification system of gene functions. The basic unit of GO is a term corresponding to a specific attribute. KEGG is a prominent public database focused on pathways. Pathway enrichment analysis takes the KEGG pathways as a unit and applies a hypergeometric test to identify those that are significantly enriched in genes. DEGs with GO and KEGG annotations were enriched and analyzed. The number of genes was then calculated for each term, and GO and KEGG terms that showed significant enrichment in DEGs compared to the genomic background were identified using the hypergeometric test. The p-value is a statistical significance value used to measure the significance of DEGs. It is calculated within a hypothesis testing framework to assess the probability of observing the observed difference or a more extreme difference under the null hypothesis. The Q-value is a statistic used to control for the probability of false positive DEGs. It is a corrected p-value that can be used to assess whether DEGs are truly biologically significant. In GO and KEGG analyses, significantly enriched terms in DEGs were defined as those with a p-value less than 0.05 for GO and KEGG terms.

3 Results

3.1 RNA identification and classification

To obtain the mRNA expression change profiles caused by RABV infection, a transcriptomic approach was performed using salivary gland tissue from mice infected with CVS-11 and PB4 strains, as well as a mock group. The RABV-infected mice exhibited clinical signs of infection, such as hunchback, collapse, tremors and hind limb paralysis as early as 3–4 days post-infection (dpi). At the advanced stage of the disease, we collected tissues from the SMG, SLG, and PG to measure the viral load using RT-qPCR. Our results showed that PG had the highest viral load compared to SMG and SLG (SI Figure 1). Therefore, our transcriptomic study focused on the parotid gland. We identified a total of 353,926,970 raw data points in all three groups, and the sequencing quality met the required standards (Table 2). The distribution of RNA on the chromosomes of the RABV- and mock-infected groups is shown in the pie chart (Figure 1). Valid data were derived from exons, introns, and intergenic regions.

Figure 1

3.2 Transcriptome regulation of parotid tissue after RABV infection

We identified a total of 3,764 DEGs (678 up-regulated and 3,086 down-regulated) in the CVS-11- infected group. Similarly, we found 4,557 DEGs in the PB4- infected group, with 874 genes up-regulated and 3,683 genes down-regulated. The top 20 up- and down-regulated DEGs with the largest fold change in the CVS-11 and PB4-infected parotid glands compared to the mock group are shown in Tables 3, 4. The intersection of DEGs between the mock group and the experimental group was visually shown using the Venn diagram (Figure 2A). Furthermore, we utilized the volcano map (Figures 2B,C) and the heat map (Figure 2D) to depict the distribution of DEGs. Our findings demonstrate that both CVS-11 and PB4 infection significantly alter the mRNA expression profile in the mouse parotid gland, leading to the dysregulation of multiple genes.

Figure 2

3.3 Functional enrichment analysis of mice PG transcriptome

To gain insight into the potential function of the aforementioned DEGs, we conducted GO and KEGG analyses. The GO enrichment analysis revealed that both infected groups showed similar enrichment patterns (Figures 3A,B). The ‘biological processes’ were mainly enriched in signal transduction, positive regulation of transcription by RNA polymerase II and cell differentiation. The ‘cell component’ was mainly enriched in the cell membrane, cytoplasm and integral component of the membrane. The ‘Molecular functions’ included protein binding, metal ion binding and identical protein binding. KEGG enrichment analysis showed similar pathways in both infected groups, with enrichment in pathways such as cancer, focal adhesion, the PI3K-Akt pathway, the salivary secretion, the cyclic guanosine monophosphate-protein kinase G (cGMP-PKG) signaling pathway, and the extracellular matrix (ECM)-receptor pathway (Figures 3C–F). We focused on pathways related to classical RABV infection signaling pathways or salivary gland functions for further analysis, which are described in detail below.

Figure 3

3.4 Associated crucial signaling pathways

CVS-11 or PB4 infection leads to thousands of DEGs in mouse PG tissue, and the DEGs are highly enriched in the salivary secretion pathway and the PI3K-AKT pathway. RABV was found to replicate in the salivary gland, and the physiology of this tissue in turn affects viral secretion. Figures 4A,B display some of the DEGs involved in the salivary secretion pathway. Our sequencing data revealed that 6 genes were up-regulated and 32 genes were down-regulated in the CVS-11-infected group, while 5 genes were up-regulated and 36 genes were down-regulated in the PB4-infected group. Among them, aquaporin 5 (AQP5), K⁺ channel protein-potassium calcium-activated channel subfamily N member 4 (KCNN4), and ADP-ribosyl cyclase (ADPRC)/cyclic ADP-ribose (cADPR) hydrolase 1 (CD38) were significantly down-regulated in both CVS-11 and PB4-infected samples, whereas the agglutinin deleted in malignant brain tumor 1 (DMBT1) showed obvious up-regulation in both RABV-infected groups. The majority of proteins encoded by DEGs in this pathway are Na⁺/K⁺-ATPase (NKA) subunits, ATPase, adrenergic receptor, adenylate cyclase, and inositol 1,4,5-trisphosphate receptor. Besides, DEGs in the parotid gland were significantly enriched in the PI3K-AKT pathway (Figures 4C,D), which is a frequently targeted intracellular signal transduction pathway by various viruses, including the mouse polyomavirus (MPyV), simian virus 40 (SV40) and papillomaviruses (PVs) (Buchkovich et al., 2008; Saeed et al., 2008; Feng et al., 2011). In the CVS-11-infected group, 111 genes were significantly enriched in this pathway (2 up-regulated and 109 down-regulated genes), while the PB4-infected group included 124 genes that were significantly enriched in this pathway (2 up-regulated and 122 down-regulated genes).

Figure 4

3.5 Validation of differentially expressed mRNAs by RT-qPCR

To validate the reliability of the RNA-seq data, we performed RT-qPCR to detect DEGs in RABV-infected and mock samples. We selected some of the Top 20 genes with the lowest Q-value that exhibited significant changes in CVS-11 and PB4-infected samples compared to the mock group for RT-qPCR validation (Figures 5A,B). The RT-qPCR results showed high concordance between the two datasets. The expression of lipocalin 11 (Lcn11), C-X-C motif chemokine ligand 13 (Cxcl13), alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase (Abo), and phytanoyl-CoA 2-hydroxylase (Phyh) were upregulated, while acetyl-CoA carboxylase alpha (Acaca), fatty acid synthase (Fasn), and ELOVL fatty acid elongase 6 (Elovl6) were significantly decreased in both strain-infected groups. Additionally, several genes were selected that showed significant trends of change in each group. Secretoglobin, family 2B, member 7 (Scgb2b7) and exocrine gland secreted peptide 31 (Esp31) were found to be increased in CVS-11 infected-parotids, while exocrine gland secreted peptide 23 (Esp23), anterior gradient 2 (Agr2) and lipase member O1 (Lipo1) were down-regulated in PB4-infected parotids. Although the multiples of change are not necessarily very similar, the RT-qPCR results exhibited similar variation trends with the transcriptomic data in the infected groups compared to the mock group, suggesting that the results of the RNA-seq datasets are reliable.

Figure 5

4 Discussion

Rabies is a highly fatal zoonotic disease that is primarily transmitted through the bites of infected animals, like dogs in many developing regions of Asia and Africa, and wild animals such as bats and foxes in Western Europe and the Americas (Wilde et al., 2013; Anderson and Shwiff, 2015). Rabies often progresses to a fatal form of encephalomyelitis (Hemachudha et al., 2013). It exhibits neurophilic properties and causes damage to the nervous system (Hemachudha et al., 2013). Salivary excretion of the virus is necessarily influenced by the basic physiology of the salivary gland (Boonsriroj et al., 2016). Transcriptomics has been widely used to screen for potentially functional genes. Previous RABV-related transcriptomics studies have mainly focused on the nervous system, investigating changes in lncRNA, circRNA, and mRNA expression in the mouse brain caused by RABV infection (Zhang et al., 2016b; Sui et al., 2021; Zhao et al., 2021). The spatial and cellular distribution of RABV infection in the mouse brain has also been deeply revealed (Huang and Sabatini, 2020; Zhang et al., 2022). However, studies on the pathogenic mechanism of RABV rarely focus on the salivary glands. To further examine the mechanism of salivary gland function, we infected mice with CVS-11 and PB4 strains and performed RNA sequencing of the PG for the first time. We identified 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.

The salivary glands are mainly composed of SMG, SLG, and PG, along with minor salivary glands (MSG) scattered throughout the mouth. Among these glands, 85–90% of the saliva is produced by the PG and SMG. The PG is located on each side of the face anterior to the ears, and is the main source of stimulated salivary secretion, secreting watery saliva that plays a role in lubricating the oral cavity and maintaining homeostasis. Actions such as biting and chewing are accompanied by contraction of the masseter muscle, which leads to secretion of saliva from the PG into the oral cavity (Brazen and Dyer, 2023). RABV has been found to invade and replicate in the PG of infected hosts (Howard, 1981; Charlton et al., 1984; Silva et al., 2009), and it is possible that a bite from a rabid animal may stimulate the secretory function of the PG. The process of fluid secretion in the salivary glands relies on various ion-transporting proteins expressed in the acinar cells. These proteins are primarily targeted at the apical or basolateral membrane of the acinar cells. In the salivary secretion pathway, several genes that regulate the transport of water and important ion-transporting proteins across membranes were found to be dysregulated in our RABV-infected samples. AQP5, a protein primarily expressed in exocrine glands like salivary and lacrimal glands, plays a crucial role in water transport and is decreased following RABV infection (Agre et al., 1993; Agre, 2004). AQP5 is critical for maintaining the normal physiology of the salivary glands, and its expression in salivary gland tissue could be affected by sialolithiasis and chronic salivary gland inflammation (D'Agostino et al., 2020). Research in AQP5 mutant rats (Murdiastuti et al., 2006), AQP5 knockout mice (Ma et al., 1999), and Sjögren’s syndrome patients (Steinfeld et al., 2001; Ichiyama et al., 2018; Hosoi et al., 2020) has demonstrated that abnormal expression of AQP5 leads to decreased salivary secretion and hypertonic and viscous saliva, which induces thirst or dryness in the oral cavity. Additionally, AQP5 has been found to be associated with viral infections. For instance, PRRSV and adenovirus infection could lead to the downregulation of AQP5, hindering water clearance and edema resolution, and causing abnormal fluid flux during pulmonary inflammation in infected animals (Generous et al., 2014; Zhang et al., 2018). The reduction of AQP5 expression in the infected parotid glands suggests that RABV may inhibit the expression and function of AQP5, thus affecting the intracellular and extracellular fluid balance and the function of water channels. This may be related to the invasion and spread of RABV and the development of the disease. CD38 is a transmembrane glycoprotein that acts as a Ca²⁺ releaser from intracellular Ca²⁺-stores (Lee, 1997; Lee, 2001). Furthermore, the Ca²⁺-activated K⁺ channel of the intermediate single channel conductance encoded by KCNN4 has been identified in exocrine salivary glands (Catalán et al., 2014). In our study, CD38 and KCNN4 were both decreased in mRNA levels, suggesting that RABV may affect the ion transporting system in salivary glands.

There are also dysregulated sequenced genes involved in the salivary secretion signaling pathway. DMBT1, which belongs to the scavenger receptor cysteine-rich (SRCR) superfamily (Generous et al., 2014; Reichhardt et al., 2017) and interacts with various pathogens in saliva (Ligtenberg et al., 2007), was significantly upregulated in both CVS-11 and PB4-infected groups. This may be due to the increased agglutination of virions in the salivary glands. NOS generates nitric oxide (NO) which is thought to contribute to protective and cytotoxic antiviral responses (Wink et al., 1996). Inducible nitric oxide synthase (iNOS) expression and activity have been shown to increase markedly in the brain of rats infected with RABV, indicating its importance in rabies neuropathogenesis (Hooper et al., 1995; Van Dam et al., 1995). Endothelial nitric oxide synthase (eNOS) expression was increased in cerebellar Purkinje cells and brainstem neurons in the brains of rabid cattle. Furthermore, eNOS and RABV were observed to co-localize in Negri bodies, suggesting the involvement of eNOS in the formation of RABV inclusion bodies (Akaike et al., 1995). NOS is widely distributed in various parts of the salivary glands, and is probably activated leading to the generation of cGMP, which opens ion channels to initiate the salivary secretory process (Bi and Reiss, 1995; Lomniczi et al., 1998). In our sequencing results, we observed down-regulation of NOS in the RABV-infected group, indicating that RABV may cause salivary gland damage similar to autoimmune inflammation, resulting in decreased NOS and NO expression levels to facilitate viral infection. Moreover, other signaling pathways are also involved in CVS-11 and PB4 infection, like the PI3K-Akt signaling pathway, the janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling pathway, and the retinoic acid-inducible gene I (RIG-I)-like receptor signaling pathway. Although it has been reported that RABV could stimulate these immune signaling pathways (Der et al., 1998; Goodbourn et al., 2000; Hornung et al., 2006; Faul et al., 2010), most immunity-related genes in our transcriptome data was restrained, like interferon gamma receptors (Ifngrs), interferon alpha receptors (Ifnars), Janus kinases (Jaks), signal transducers and activators of transcription (Stats), stimulator of interferon response cGAMP interactors (Stings), and TNF receptor associated factor 5 (Traf5). It has been speculated that the host immune system remains in a suppressed state as rabies progresses to late-stage disease, in order to persist the transmission and secretion of virions through saliva. However, the detailed mechanism by which the genes regulating salivary virion production and salivary gland function are dysregulated remains unclear. Further research is needed. Whether these genes could be potential targets for blocking virus transmission through saliva requires further validation. Whether the gene expression map of parotid tissue in the early stage of rabies differs from that in the later stage may merit intensive identification.

In summary, we demonstrated differentially expressed mRNA affected by RABV infection in the parotid tissues by employing the transcriptomic technique, and established GO and KEGG functional enrichment analyses. It is revealed that RABV invasion caused drastic regulation of the saliva secretion process in mice PG along with various host response processes. This study was the first to analyze the gene expression profiles and biological functions of the DEGs derived from RABV-infected mouse salivary tissues, which may have implications for the development of therapeutics, diagnostics and sanitation measures aimed at preventing the spread of rabies through saliva.

5 Conclusion

RNA-seq was utilized to establish mRNA expression profiles in PG mice infected with CVS-11 and PB4, and the biological functions of DEGs were analyzed. GO and KEGG enrichment analyses showed that most of the dysregulated genes were associated with the salivary secretion pathway, the PI3K-Akt signaling pathway, and the cGMP-PKG signaling pathway. Although some work has been undertaken, the detailed mechanism underlying RABV-induced salivary gland dysfunction remains elusive. These findings will contribute to a deeper understanding of the biology of viral spread.

References

• 1

  AgreP. (2004). Aquaporin water channels (Nobel lecture). Angew. Chem. Int. Ed. Engl.43, 4278–4290. doi: 10.1002/anie.200460804

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AgreP.PrestonG. M.SmithB. L.JungJ. S.RainaS.MoonC.et al. (1993). Aquaporin CHIP: the archetypal molecular water channel. Am. J. Phys.265, F463–F476. doi: 10.1152/ajprenal.1993.265.4.F463

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AkaikeT.WeiheE.SchaeferM.FuZ. F.ZhengY. M.VogelW.et al. (1995). Effect of neurotropic virus infection on neuronal and inducible nitric oxide synthase activity in rat brain. J. Neurovirol.1, 118–125. doi: 10.3109/13550289509111016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  AndersonA.ShwiffS. A. (2015). The cost of canine rabies on four continents. Transbound. Emerg. Dis.62, 446–452. doi: 10.1111/tbed.12168

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  Anfruns-EstradaE.SabriàA.FuentesC.SabatéS.RazquinE.CornejoT.et al. (2020). Detection of norovirus in saliva samples from acute gastroenteritis cases and asymptomatic subjects: association with age and higher shedding in stool. Viruses12:1369. doi: 10.3390/v12121369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  BeauregardM.CaseyG. A. (1969). Demonstration of rabies antigen in salivary glands of rabies suspected animals. Can. J. Comp. Med.33, 55–58. PMID:

  • Pubmed Abstract
  • Google Scholar

• 7

  BiZ.ReissC. S. (1995). Inhibition of vesicular stomatitis virus infection by nitric oxide. J. Virol.69, 2208–2213. doi: 10.1128/jvi.69.4.2208-2213.1995

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  BoonsrirojH.ManaloD. L.KimitsukiK.ShimatsuT.ShiwaN.ShinozakiH.et al. (2016). A pathological study of the salivary glands of rabid dogs in the Philippines. J. Vet. Med. Sci.78, 35–42. doi: 10.1292/jvms.15-0308

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  BrazenB.DyerJ. (2023). Histology, salivary glands. Treasure Island, FL: StatPearls Publishing LLC.

  • Google Scholar

• 10

  BuchkovichN. J.YuY.ZampieriC. A.AlwineJ. C. (2008). The TORrid affairs of viruses: effects of mammalian DNA viruses on the PI3K-Akt-mTOR signalling pathway. Nat. Rev. Microbiol.6, 266–275. doi: 10.1038/nrmicro1855

  • CrossRef
  • Google Scholar

• 11

  CatalánM. A.Peña-MunzenmayerG.MelvinJ. E. (2014). Ca²⁺-dependent K⁺ channels in exocrine salivary glands. Cell Calcium55, 362–368. doi: 10.1016/j.ceca.2014.01.005

  • CrossRef
  • Google Scholar

• 12

  CharltonK. M.CaseyG. A.CampbellJ. B. (1983). Experimental rabies in skunks: mechanisms of infection of the salivary glands. Can. J. Comp. Med.47, 363–369. PMID:

  • Pubmed Abstract
  • Google Scholar

• 13

  CharltonK. M.CaseyG. A.WebsterW. A. (1984). Rabies virus in the salivary glands and nasal mucosa of naturally infected skunks. Can. J. Comp. Med.48, 338–339. PMID:

  • Pubmed Abstract
  • Google Scholar

• 14

  CorstjensP. L.AbramsW. R.MalamudD. (2000). (2016) saliva and viral infections. Periodontol.70, 93–110. doi: 10.1111/prd.12112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  D'AgostinoC.ElkashtyO. A.ChivassoC.PerretJ.TranS. D.DelporteC. (2020). Insight into salivary gland Aquaporins. Cell9:1547. doi: 10.3390/cells9061547

  • CrossRef
  • Google Scholar

• 16

  DerS. D.ZhouA.WilliamsB. R.SilvermanR. H. (1998). Identification of genes differentially regulated by interferon alpha, beta, or gamma using oligonucleotide arrays. Proc. Natl. Acad. Sci. USA95, 15623–15628. doi: 10.1073/pnas.95.26.15623

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  DierksR. E.MurphyF. A.HarrisonA. K. (1969). Extraneural rabies virus infection. Virus development in fox salivary gland. Am. J. Pathol.54, 251–273. PMID:

  • Pubmed Abstract
  • Google Scholar

• 18

  DietzscholdB.LiJ.FaberM.SchnellM. (2008). Concepts in the pathogenesis of rabies. Future Virol.3, 481–490. doi: 10.2217/17460794.3.5.481

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  FaulE. J.WanjallaC. N.SutharM. S.GaleM.WirblichC.SchnellM. J. (2010). Rabies virus infection induces type I interferon production in an IPS-1 dependent manner while dendritic cell activation relies on IFNAR signaling. PLoS Pathog.6:e1001016. doi: 10.1371/journal.ppat.1001016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  FekaduM.ShaddockJ. H. (1984). Peripheral distribution of virus in dogs inoculated with two strains of rabies virus. Am. J. Vet. Res.45, 724–729. PMID:

  • Pubmed Abstract
  • Google Scholar

• 21

  FengS. Z.CaoW. S.LiaoM. (2011). The PI3K/Akt pathway is involved in early infection of some exogenous avian leukosis viruses. J. Gen. Virol.92, 1688–1697. doi: 10.1099/vir.0.030866-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  FooksA. R.BanyardA. C.HortonD. L.JohnsonN.McElhinneyL. M.JacksonA. C. (2014). Current status of rabies and prospects for elimination. Lancet384, 1389–1399. doi: 10.1016/s0140-6736(13)62707-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  FooksA. R.CliquetF.FinkeS.FreulingC.HemachudhaT.ManiR. S.et al. (2017). Rabies. Nat. Rev. Dis. Primers3:17091. doi: 10.1038/nrdp.2017.91

  • CrossRef
  • Google Scholar

• 24

  GenerousA.ThorsonM.BarcusJ.JacherJ.BuschM.SleisterH. (2014). Identification of putative interactions between swine and human influenza a virus nucleoprotein and human host proteins. Virol. J.11:228. doi: 10.1186/s12985-014-0228-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  GhoshS.KumarM.SantianaM.MishraA.ZhangM.LabayoH.et al. (2022). Enteric viruses replicate in salivary glands and infect through saliva. Nature607, 345–350. doi: 10.1038/s41586-022-04895-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  GoodbournS.DidcockL.RandallR. E. (2000). Interferons: cell signalling, immune modulation, antiviral response and virus countermeasures. J. Gen. Virol.81, 2341–2364. doi: 10.1099/0022-1317-81-10-2341

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  HemachudhaT.LaothamatasJ.RupprechtC. E. (2002). Human rabies: a disease of complex neuropathogenetic mechanisms and diagnostic challenges. Lancet Neurol.1, 101–109. doi: 10.1016/s1474-4422(02)00041-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  HemachudhaT.UgoliniG.WacharapluesadeeS.SungkaratW.ShuangshotiS.LaothamatasJ. (2013). Human rabies: neuropathogenesis, diagnosis, and management. Lancet Neurol.12, 498–513. doi: 10.1016/s1474-4422(13)70038-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  HooperD. C.OhnishiS. T.KeanR.NumagamiY.DietzscholdB.KoprowskiH. (1995). Local nitric oxide production in viral and autoimmune diseases of the central nervous system. Proc. Natl. Acad. Sci. USA92, 5312–5316. doi: 10.1073/pnas.92.12.5312

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  HornungV.EllegastJ.KimS.BrzózkaK.JungA.KatoH.et al. (2006). 5′-triphosphate RNA is the ligand for RIG-I. Science314, 994–997. doi: 10.1126/science.1132505

  • CrossRef
  • Google Scholar

• 31

  HosoiK.YaoC.HasegawaT.YoshimuraH.AkamatsuT. (2020). Dynamics of salivary gland AQP5 under Normal and pathologic conditions. Int. J. Mol. Sci.21:182. doi: 10.3390/ijms21041182

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HowardD. R. (1981). Rabies virus titer from tissues of naturally infected skunks (Mephitis mephitis). Am. J. Vet. Res.42, 1595–1597. PMID:

  • Pubmed Abstract
  • Google Scholar

• 33

  HuangN.PérezP.KatoT.MikamiY.OkudaK.GilmoreR. C.et al. (2021). SARS-CoV-2 infection of the oral cavity and saliva. Nat. Med.27, 892–903. doi: 10.1038/s41591-021-01296-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  HuangK. W.SabatiniB. L. (2020). Single-cell analysis of Neuroinflammatory responses following intracranial injection of G-deleted rabies viruses. Front. Cell. Neurosci.14:65. doi: 10.3389/fncel.2020.00065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  IchiyamaT.NakataniE.TatsumiK.HideshimaK.UranoT.NariaiY.et al. (2018). Expression of aquaporin 3 and 5 as a potential marker for distinguishing dry mouth from Sjögren's syndrome. J. Oral Sci.60, 212–220. doi: 10.2334/josnusd.17-0150

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  JacksonA. C. (2016). Diabolical effects of rabies encephalitis. J. Neurovirol.22, 8–13. doi: 10.1007/s13365-015-0351-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  JogaiS.RadotraB. D.BanerjeeA. K. (2002). Rabies viral antigen in extracranial organs: a post-mortem study. Neuropathol. Appl. Neurobiol.28, 334–338. doi: 10.1046/j.1365-2990.2002.00400.x

  • CrossRef
  • Google Scholar

• 38

  KirbyA.DoveW.AshtonL.HopkinsM.CunliffeN. A. (2010). Detection of norovirus in mouthwash samples from patients with acute gastroenteritis. J. Clin. Virol.48, 285–287. doi: 10.1016/j.jcv.2010.05.009

  • CrossRef
  • Google Scholar

• 39

  LeeH. C. (1997). Mechanisms of calcium signaling by cyclic ADP-ribose and NAADP. Physiol. Rev.77, 1133–1164. doi: 10.1152/physrev.1997.77.4.1133

  • CrossRef
  • Google Scholar

• 40

  LeeH. C. (2001). Physiological functions of cyclic ADP-ribose and NAADP as calcium messengers. Annu. Rev. Pharmacol. Toxicol.41, 317–345. doi: 10.1146/annurev.pharmtox.41.1.317

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  LeeJ. M.GaronE.WongD. T. (2009). Salivary diagnostics. Orthod. Craniofac. Res.12, 206–211. doi: 10.1111/j.1601-6343.2009.01454.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LewisP.FuY.LentzT. L. (2000). Rabies virus entry at the neuromuscular junction in nerve-muscle cocultures. Muscle Nerve23, 720–730. doi: 10.1002/(sici)1097-4598(200005)23:5<720::aid-mus9>3.0.co;2-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  LiZ.FengZ.YeH. (1995). Rabies viral antigen in human tongues and salivary glands. J. Trop. Med. Hyg.98, 330–332. PMID:

  • Pubmed Abstract
  • Google Scholar

• 44

  LigtenbergA. J.VeermanE. C.Nieuw AmerongenA. V.MollenhauerJ. (2007). Salivary agglutinin/glycoprotein-340/DMBT1: a single molecule with variable composition and with different functions in infection, inflammation and cancer. Biol. Chem.388, 1275–1289. doi: 10.1515/bc.2007.158

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  LimaD. P.DinizD. G.MoimazS. A.SumidaD. H.OkamotoA. C. (2010). Saliva: reflection of the body. Int. J. Infect. Dis.14, e184–e188. doi: 10.1016/j.ijid.2009.04.022

  • CrossRef
  • Google Scholar

• 46

  LomnicziA.SuburoA. M.ElverdinJ. C.MastronardiC. A.DiazS.RettoriV.et al. (1998). Role of nitric oxide in salivary secretion. Neuroimmunomodulation5, 226–233. doi: 10.1159/000026342

  • CrossRef
  • Google Scholar

• 47

  MaT.SongY.GillespieA.CarlsonE. J.EpsteinC. J.VerkmanA. S. (1999). Defective secretion of saliva in transgenic mice lacking aquaporin-5 water channels. J. Biol. Chem.274, 20071–20074. doi: 10.1074/jbc.274.29.20071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  MaierT.SchwartingA.MauerD.RossR. S.MartensA.KliemV.et al. (2010). Management and outcomes after multiple corneal and solid organ transplantations from a donor infected with rabies virus. Clin. Infect. Dis.50, 1112–1119. doi: 10.1086/651267

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  MitrabhakdiE.ShuangshotiS.WannakrairotP.LewisR. A.SusukiK.LaothamatasJ.et al. (2005). Difference in neuropathogenetic mechanisms in human furious and paralytic rabies. J. Neurol. Sci.238, 3–10. doi: 10.1016/j.jns.2005.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  MurdiastutiK.PurwantiN.KarabasilM. R.LiX.YaoC.AkamatsuT.et al. (2006). A naturally occurring point mutation in the rat aquaporin 5 gene, influencing its protein production by and secretion of water from salivary glands. Am. J. Physiol. Gastrointest. Liver Physiol.291, G1081–G1088. doi: 10.1152/ajpgi.00449.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  PisanicN.BallardS. B.ColquechaguaF. D.FrançoisR.ExumN.YoriP. P.et al. (2019). Minimally invasive saliva testing to monitor norovirus infection in community settings. J. Infect. Dis.219, 1234–1242. doi: 10.1093/infdis/jiy638

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  ReichhardtM. P.HolmskovU.MeriS. (2017). SALSA-A dance on a slippery floor with changing partners. Mol. Immunol.89, 100–110. doi: 10.1016/j.molimm.2017.05.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  SaeedM. F.KolokoltsovA. A.FreibergA. N.HolbrookM. R.DaveyR. A. (2008). Phosphoinositide-3 kinase-Akt pathway controls cellular entry of Ebola virus. PLoS Pathog.4:e1000141. doi: 10.1371/journal.ppat.1000141

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  SilvaM. L.Lima FdaS.de Barros GomesA. A.de AzevedoS. S.AlvesC. J.BernardiF.et al. (2009). Isolation of rabies virus from the parotid salivary glands of foxes (Pseudalopex vetulus) from Paraíba state, Northeastern Brazil. Braz. J. Microbiol.40, 446–449. doi: 10.1590/s1517-83822009000300004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  SinghR.SinghK. P.CherianS.SaminathanM.KapoorS.Manjunatha ReddyG. B.et al. (2017). Rabies - epidemiology, pathogenesis, public health concerns and advances in diagnosis and control: a comprehensive review. Vet. Q.37, 212–251. doi: 10.1080/01652176.2017.1343516

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  SrinivasanA.BurtonE. C.KuehnertM. J.RupprechtC.SutkerW. L.KsiazekT. G.et al. (2005). Transmission of rabies virus from an organ donor to four transplant recipients. N. Engl. J. Med.352, 1103–1111. doi: 10.1056/NEJMoa043018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  SteinfeldS.CoganE.KingL. S.AgreP.KissR.DelporteC. (2001). Abnormal distribution of aquaporin-5 water channel protein in salivary glands from Sjögren's syndrome patients. Lab. Investig.81, 143–148. doi: 10.1038/labinvest.3780221

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SuiB.ChenD.LiuW.TianB.LvL.PeiJ.et al. (2021). Comparison of lncRNA and mRNA expression in mouse brains infected by a wild-type and a lab-attenuated rabies lyssavirus. J. Gen. Virol.102:98. doi: 10.1099/jgv.0.001538

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  TobiumeM.SatoY.KatanoH.NakajimaN.TanakaK.NoguchiA.et al. (2009). Rabies virus dissemination in neural tissues of autopsy cases due to rabies imported into Japan from the Philippines: immunohistochemistry. Pathol. Int.59, 555–566. doi: 10.1111/j.1440-1827.2009.02406.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  Van DamA. M.BauerJ.ManA. H. W. K.MarquetteC.TildersF. J.BerkenboschF. (1995). Appearance of inducible nitric oxide synthase in the rat central nervous system after rabies virus infection and during experimental allergic encephalomyelitis but not after peripheral administration of endotoxin. J. Neurosci. Res.40, 251–260. doi: 10.1002/jnr.490400214

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  WildeH.HemachudhaT.WacharapluesadeeS.LumlertdachaB.TepsumethanonV. (2013). Rabies in Asia: the classical zoonosis. Curr. Top. Microbiol. Immunol.365, 185–203. doi: 10.1007/82_2012_228

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  WinkD. A.HanbauerI.GrishamM. B.LavalF.NimsR. W.LavalJ.et al. (1996). Chemical biology of nitric oxide: regulation and protective and toxic mechanisms. Curr. Top. Cell. Regul.34, 159–187. doi: 10.1016/s0070-2137(96)80006-9

  • CrossRef
  • Google Scholar

• 63

  ZhangC. Z.ChengX. Q.LiJ. Y.ZhangP.YiP.XuX.et al. (2016a). Saliva in the diagnosis of diseases. Int. J. Oral Sci.8, 133–137. doi: 10.1038/ijos.2016.38

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  ZhangD.HeF.BiS.GuoH.ZhangB.WuF.et al. (2016b). Genome-wide transcriptional profiling reveals two distinct outcomes in central nervous system infections of rabies virus. Front. Microbiol.7:751. doi: 10.3389/fmicb.2016.00751

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  ZhangY.XingX.LongB.CaoY.HuS.LiX.et al. (2022). A spatial and cellular distribution of rabies virus infection in the mouse brain revealed by fMOST and single-cell RNA sequencing. Clin. Transl. Med.12:e700. doi: 10.1002/ctm2.700

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  ZhangJ.YanM.GuW.ChenA.LiuJ.LiL.et al. (2018). Downregulation of Aquaporins (AQP1 and AQP5) and Na,K-ATPase in porcine reproductive and respiratory syndrome virus-infected pig lungs. Inflammation41, 1104–1114. doi: 10.1007/s10753-018-0762-2

  • CrossRef
  • Google Scholar

• 67

  ZhaoW.SuJ.WangN.ZhaoN.SuS. (2021). Expression profiling and bioinformatics analysis of CircRNA in mice brain infected with rabies virus. Int. J. Mol. Sci.22:537. doi: 10.3390/ijms22126537

  • Pubmed Abstract
  • CrossRef
  • Google Scholar