Infection of BALB/c mice with the JX H5N6 virus
Previously, we have demonstrated that HPAIV H5N6 (A/duck/Hubei/WH18/2015, abbreviated as JX) is neurotoxic and highly pathogenic to mice. To establish animal models of brain damage caused by IAV, BALB/c mice were inoculated intranasally with 105 TCID50/50 µL inoculum/mouse of IAV JX, and the mock-treated mice were used as control. All mice infected with JX showed clearly detectable symptoms of the flu at 3 dpi, and showed drastic weight loss (Fig. 1A, F (1, 98) = 293.2, P < 0.0001 ). Strikingly, these symptoms had progressed to neurological symptoms such as balancing problems, hind limb weakness, and paralysis by 5 and 6 dpi and mice began to die by 7 dpi. At the end of the experimental period, only two out of five mice survived (Fig. 1B, Mantel-Cox log-rank test, P = 0.0511 ). Next, IHC was performed to investigate whether the brain was successfully infected by the JX virus. Results suggested that the brains of mice were successfully infected (Fig. 1C, 1D). Meanwhile,IHC staining for IAV viral antigen revealed virus-positive neurons in the brains of mice sacrificed on 3 dpi. Virus antigen was detected in multiple areas of the brain, including cerebral cortex, brainstem, occipital lobe. For representative regions, an image of the cerebral cortex with a single staining of IAV antigen is shown (Fig. 1C). Brain tissue injury caused by JX virus involves sub-meningeal hemorrhage, degeneration, and necrosis of neurons and local glial cell nodules(Fig. 1D). Collectively, these results demonstrated that the JX virus can infect, replicate, and spread throughout the CNS, thereby successfully establishing a mouse model for further functional studies.
Brain transcription profile after JX virus infection
To determine the differential response of the host brain tissue to infection with JX and mock infection, RNA samples were prepared from whole brains and the complete gene expression profile was compared after whole transcriptome sequencing (Fig. 2). Histopathological and immunohistochemical examinations showed neuron necrosis, proliferation of glial cells, and viral antigen positivity on 3 dpi. Therefore, brain samples were obtained on this day. A cut off of > 2.0 in fold change and P < 0.05 was used initially to examine the significantly differentially expressed genes that were up- or down-regulated following JX infection (Fig. 2). In total, 32,716 genes were detected in the brain, 240 of which were identified as differentially expressed mRNAs between the virus-infected and uninfected mice. Among these, 179 genes were up-regulated with log2 fold change values ranging from 1.01 to 5.08, while 61 genes were down-regulated with log2 fold change values ranging from 1 to 4.59. The top 10 up-regulated genes listed were Cxx1b/Cxx1a, Scrt1, Bmyc, Pcp2, Car8,Arl4c, Fat2, Mybpc3, Rn7sk, and Irf2bp1 (Table 1). The top 10 most down-regulated genes were Nudc-ps1, Glrx5, Pgrmc1, Rnf113a2, Irgq, Rgs4, Tmem97, Zik1, Zcchc3, and Gng4 (Table 2). To further assess the connection between gene expression pattern and IAV infection-induced biological processes in the brain, functional classification of mRNA transcripts and pathway analysis were performed. These analyses revealed that these genes were involved in regulating neuronal differentiation, ion transmembrane transport, and protein binding. Some of the significant pathways involving the differentially expressed genes were associated with herpes simplex infection, glutamatergic synapse, and chemokine signaling pathway. Among the significantly down-regulated genes, we focused on PGRMC1 to evaluate their role in IAV infection in the brain.
Table 1
The top 10 most up-regulated genes of virus-infected mice compared to control mice
Symbol
|
Gene name
|
Log2FC
|
p-value
|
Cxx1b;Cxx1a
|
CAAX box 1A ;CAAX box 1B
|
5.68
|
0.014
|
Scrt1
|
scratch family zinc finger 1
|
4.47
|
0.021
|
Bmyc
|
brain expressed myelocytomatosis oncogene
|
4.35
|
0.046
|
Pcp2
|
Purkinje cell protein 2 (L7)
|
4.31
|
0.037
|
Car8
|
carbonic anhydrase 8
|
4.25
|
0.044
|
Arl4c
|
ADP-ribosylation factor-like 4C
|
3.94
|
0.046
|
Fat2
|
FAT atypical cadherin 2
|
3.56
|
0.045
|
Mybpc3
|
myosin binding protein C, cardiac
|
3.51
|
0.038
|
Rn7sk
|
RNA, 7SK, nuclear
|
3.32
|
0.05
|
Irf2bp1
|
interferon regulatory factor 2 binding protein 1
|
3.27
|
0.046
|
Table 2
The top 10 most down-regulated genes of virus-infected mice compared to that in the control mice
Symbol
|
Gene name
|
Log2FC
|
p-value
|
Nudc-ps1
|
nuclear distribution gene C homolog (Aspergillus), pseudogene 1
|
4.59
|
0.001
|
Glrx5
|
glutaredoxin 5
|
3.87
|
0.043
|
Pgrmc1
|
progesterone receptor membrane component 1
|
3.35
|
0
|
Rnf113a2
|
ring finger protein 113A2
|
3.1
|
0.038
|
Irgq
|
immunity-related GTPase family, Q
|
2.92
|
0.014
|
Rgs4
|
regulator of G-protein signaling 4
|
2.6
|
0.007
|
Tmem97
|
transmembrane protein 97
|
2.38
|
0.031
|
Zik1
|
zinc finger protein interacting with K protein 1
|
2.16
|
0.034
|
Zcchc3
|
zin(Kimura et al., 2006)c finger, CCHC domain containing 3
|
2.11
|
0.019
|
Gng4
|
guanine nucleotide binding protein (G protein), gamma 4
|
2.1
|
0.014
|
H5N6 down-regulated the expression of PGRMC1 in vitro
It is well-known that the nervous tissue is composed of two primary cell types: neurons and glial cells. Neurons transmit nerve messages, while the surrounding glial cells are in direct contact with neurons (Kimura et al., 2006). To verify whether the H5N6 virus can down-regulate PGRMC1 in vitro, both human neuroblastoma cell line SK-N-SH and human glioma cell line U251 were used. SK-N-SH or U251 cells were infected with JX at multiplicity of infection (MOI) of 0.01. At 24 and 48 hpi, PGRMC1 mRNA and protein levels were assessed using qRT-PCR and western blotting, respectively (Fig. 3). Compared to the mock, JX significantly inhibited the production of PGRMC1 mRNA in SK-N-SH after 24 h, which was in agreement with the protein level (Fig. 3A, F (1, 23) = 201.3, P < 0.0001, 3C). At 48 hpi, the extent of PGRMC1 protein inhibition by the JX virus was lower than that at 24 hpi. Similar result was observed in U251 cells (Fig. 3B, F (1, 20) = 87.57, P < 0.0001, 3D). These results indicated that IAV JX may significantly inhibit the expression of PGRMC1 in agreement with those identified using RNA-seq.
PGRMC1 inhibited H5N6 replication
As IAV JX inhibited PGRMC1 expression, we speculated that PGRMC1 might participate in IAV JX replication. To better understand the biological role of PGRMC1 in H5N6 replication, AG205, a putative PGRMC1-specific inhibitor (Teakel et al., 2020), was used. SK-N-SH or U251 cells were seeded in 12-well plates and grown to 80% confluence. Then, the cells were treated with 15 µM AG-205 or vehicle (cell culture grade DMSO) for 24 h. Cells were infected with JX virus at MOI of 0.01. Viral NP mRNA levels and viral titers were evaluated at different time points (Fig. 4). NP mRNA levels in the two cell lines significantly increased in the AG205 group compared to that in the control (Fig. 4A, F (1, 21) = 116.0, P < 0.0001; 4B, F (1, 19) = 62.26, P < 0.0001). Next, we observed that AG205 significantly increased the viral titers in U251 12 h, 24 h and 36 h post-infection (Fig. 4C, F (1, 12) = 44.47,P < 0.0001; 4D, F (1, 12) = 37.78, P < 0.0001). However, viral titers significantly increased after AG205 treatment of SK-N-SH cells at 36 h, although no effect on virus replication was evident 24 h after AG205 treatment. Subsequently, we compared the effect of AG205 on viral proliferation in the two cell types and observed that AG205 consistently enhanced the viral titers of IAV JX in the U251 cells compared to that in the SK-N-SH cells. The effect of PGRMC1 overexpression on viral replication was also determined by transfecting U251 cells with Flag-PGRMC1. The efficiency of PGRMC1 overexpression was determined using western blotting (Fig. 4F). As observed in the TCID50 assays, virus titer at 36 h.p.i were lower in the PGRMC1-overexpressing groups than in the control group (Fig. 4E, F (1, 8) = 26.36, P = 0.0009). Taken together, these findings indicated that PGRMC1 was an IAV JX restriction factor.
AG205 inhibited the IAV-induced RIG-I-dependent antiviral response
AG205 affects virus multiplication more in U251 cells than in SK-N-SH cells. Consistently, the innate immune surveillance is mainly coordinated by glial cells in the CNS [27]. These CNS resident cells are assumed to orchestrate the immune response, which assists in combating infections in the brain. Hence, we hypothesized that AG205 may inhibit IAV H5N6-mediated innate immune response. To validate this hypothesis, we evaluated a series of immunological molecules in IAV JX-infected U251 cells, which were incubated with AG205 or DMSO (Fig. 5, Fig. 6). As expected, AG205 significantly down-regulated IFN-β induced by IAV JX (Fig. 5A, F (1, 8) = 126.1, P < 0.0001). In addition, the expression of the downstream effectors of IFN-β, namely, MX1 and ISG15, was suppressed (Fig. 6). We further observed that AG205 markedly inhibited poly (I:C)-mediated IFN-β induction (Fig. 5B-G). The molecules upstream of IFN-β, such as RIG-I, TBK1, and IRF3, were inhibited by AG205. MX1 and ISG15 were also suppressed. Overall, these observations suggested that AG205 inhibited the type I interferon response.
AG205 inhibited RIG-I-mediated IFN-β production
The RIG-I signaling pathway is essential for IAV recognition. Once activated, RIG-I/MAVS induces the activation of transcription factor IRF 3 and NF-κB, ultimately resulting in establishment of host IFN-mediated antiviral response (Iwasaki and Pillai, 2014). To investigate how the PGRMC1 inhibitor AG205 regulates IFN-β, the IFN-β signaling pathway was analyzed using IFN-β luciferase activity assay (Fig. 7). AG205 or DMSO-pretreated 293 cells were co-transfected with the signal molecule expression plasmid, IFN-β promoter luciferase reporter plasmid, and internal control pRL-TK. AG205 markedly inhibited RIG-I-mediated IFN-β production (Fig. 7A, F (1, 8) = 156.1, P < 0.0001). Interestingly, it did not affect other molecules (MDA5, F (1, 8) = 3.410, P = 0.1020; TBK1, F (1, 8) = 0.08015, P = 0.7843; IKK-ξ, F (1, 8) = 1.774, P = 0.2196, IRF3, F (1, 7) = 3.168, P = 0.1183; IRF3-5D, F (1, 8) = 3.705, P = 0.0905) that stimulate IFN-β promoter luciferase activity (Fig. 6B-F). Notably, as shown in Fig. 7H (one-way ANOVA, F = 46.24, P = 0.005), we observed a dose-dependent relationship between AG205 and IFN-β promoter luciferase activity. On the contrary, overexpression of PGRMC1 increased RIG-I-stimulated IFN-β luciferase activity in a dose-dependent manner (Fig. 7G, one-way ANOVA, F = 78.66, P < 0.0001). Consequently, we concluded that down-regulation or functional inhibition of PGRMC1 reduces the ability of RIG-I to induce IFN-β production in cells infected with IAV H5N6.
PGRMC1 knockout promoted IAV replication by inhibiting the RIG-I-mediated IFN-β signaling pathway
As the PGRMC1 inhibitor AG205 has been shown to significantly promote IAV replication, we speculated whether PGRMC1 knockdown similarly affected the proliferation of IAV H5N6 in U251 cells. Toward this, PGRMC1 was knocked out using the lentiviral CRISPR-Cas9 system containing a pair of suitable guide RNAs driven by two independent U6 promoters on the same plasmid ( Fig. 8). After CRISPR/Cas9 gene editing and isolation of individual clones, PGRMC1 expression was analyzed using western blot analysis. Our results suggested that CRISPR/Cas9 targeting resulted in complete loss of the PGRMC1 protein (Fig. 8A). We first evaluated the effect of PGRMC1 deficiency on IAV JX propagation using a TCID50 assay. Indeed, cells with PGRMC1 knockout produced more infectious viral particles than control U251 cells (Fig. 8B, F (1, 12) = 70.53, P < 0.0001), in agreement with the results obtained in cells treated with AG205. Furthermore, the expression of RIG-I (protein level, Fig. 6) and IFN-β (mRNA, Fig. 8C, F (3, 61) = 185.1, P < 0.0001) in JX -stimulated cells clearly decreased. Alternatively, the IFN-β upstream factors (RIG-I, and TBK1) and IFN-β downstream effectors (MX1 and ISG15) were also significantly down-regulated (Fig. 6). These results further confirmed that PGRMC1 knockout promoted virus replication by efficiently suppressing RIG-I-mediated IFN-β production.
PGRMC1 knockout antagonized RIG-I ubiquitination
RIG-I exists in an inactive closed conformation under normal conditions. RIG activation requires ubiquitination induced by the tripartite motif 25 (TRIM25) ubiquitin E3 ligase. Hence, we next addressed whether PGRMC1 knockout specifically suppresses RIG-I ubiquitination [28]. Toward this, we assessed RIG-I ubiquitination in PGRMC1-KO U251 cells or wild-type U251 cells transfected with pCAGGS-HA-Ub vector or the empty vector along with or without p3⋅ Flag-RIG-I (Fig. 9), followed by infection with JX. The results of IP experiments clearly indicated that PGRMC1 knockout markedly inhibited RIG-I ubiquitination. These results collectively revealed that the PGRMC1 knockout inhibited RIG-I ubiquitination in U251cells, thereby suppressing RIG-I signal transduction.