Comprehensive Proteomics Identification of IFN-λ3-regulated Antiviral Proteins in HBV-transfected Cells*

Interferon lambda (IFN-λ) is a relatively unexplored, yet promising antiviral agent. IFN-λ has recently been tested in clinical trials of chronic hepatitis B virus infection (CHB), with the advantage that side effects may be limited compared with IFN-α, as IFN-λ receptors are found only in epithelial cells. To date, IFN-λ's downstream signaling pathway remains largely unelucidated, particularly via proteomics methods. Here, we report that IFN-λ3 inhibits HBV replication in HepG2.2.15 cells, reducing levels of both HBV transcripts and intracellular HBV DNA. Quantitative proteomic analysis of HBV-transfected cells was performed following 24-hour IFN-λ3 treatment, with parallel IFN-α2a and PBS treatments for comparison using a dimethyl labeling method. The depth of the study allowed us to map the induction of antiviral proteins to multiple points of the viral life cycle, as well as facilitating the identification of antiviral proteins not previously known to be elicited upon HBV infection (e.g. IFITM3, XRN2, and NT5C3A). This study also shows up-regulation of many effectors involved in antigen processing/presentation indicating that this cytokine exerted immunomodulatory effects through several essential molecules for these processes. Interestingly, the 2 subunits of the immunoproteasome cap (PSME1 and PSME2) were up-regulated whereas cap components of the constitutive proteasome were down-regulated upon both IFN treatments, suggesting coordinated modulation toward the antigen processing/presentation mode. Furthermore, in addition to confirming canonical activation of interferon-stimulated gene (ISG) transcription through the JAK-STAT pathway, we reveal that IFN-λ3 restored levels of RIG-I and RIG-G, proteins known to be suppressed by HBV. Enrichment analysis demonstrated that several biological processes including RNA metabolism, translation, and ER-targeting were differentially regulated upon treatment with IFN-λ3 versus IFN-α2a. Our proteomic data suggests that IFN-λ3 regulates an array of cellular processes to control HBV replication.

Chronic hepatitis B (CHB) 1 is a major health problem worldwide, affecting 240 million people throughout the world with a prevalence in Africa and South-East Asia (1). Chronic HBVinfected individuals mostly acquire the virus at a young age through vertical transmission or contact with the blood or other body fluids of an infected person. CHB eventually progresses to severe and high-mortality liver diseases including liver cirrhosis and hepatocellular carcinoma (HCC) resulting in 600,000 deaths annually (2,3).
Current therapeutic agents for CHB are interferon (IFN)-␣ and nucleos(t)ide analogs (NAs) (4 -6). IFN-␣ is a cytokine that possesses antiviral and immunomodulatory effects, which promotes control and eradication of viral infection. The advantages of using IFN-␣ for CHB treatment are the decreased incidence of viral resistance following this treatment and the finite duration of therapy with a higher rate of seroconversion of viral antigens, i.e. HBeAg and HBsAg, compared with NAs. However, the drawbacks of IFN-␣ treatment are its inconvenient route of administration and its adverse effects such as influenza-like symptoms, nausea, vomiting, cytopenia, and psychiatric disorders. Regarding the latter agent, NAs against CHB, the lack of a 3Ј-hydroxyl group in these compounds allows competitive incorporation into the viral genome resulting in termination of viral replication (7). Although NAs can directly suppress HBV replication and have fewer unfavorable side effects compared with IFN-␣, CHB patients usually require a long-term treatment of these agents, thus increasing the chance of the emergence of viral resistance. In addition, the low rate of HBeAg and HBsAg seroconversion is another limitation of NAs. Therefore, new drugs that overcome restrictions of current anti-HBV treatments are still needed.
Several lines of reasoning suggest that IFN-could have a superior combination of efficacy and reduced side-effects. IFNor type III IFN is a cytokine in the class II cytokine family that has recently been used in clinical trials of CHB treatment (8,9). IFNhas 4 subtypes, namely IFN-1, -2, -3, and -4. IFN-3 was shown to have superior antiviral potency compared with other IFN-subtypes (10). IFN-receptor is composed of an IFNLR1 and IL10R2 dimer; hence upon ligand binding, the receptor activates both IFN and IL-10-like signaling pathways. The type I IFN-receptors also form a heterodimer (IFNAR1 and IFNAR2), however, the sequences of all these receptors do diverge, particularly at the C terminus (11,12). IFNLR1 is expressed only in epithelial cells including hepatocytes in contrast to IFN-␣ receptor, which is widely expressed in many cell types throughout the body; thus IFNtreatment results in fewer side-effects when compared with IFN-␣ treatment (9,13,14). IFN-has been shown to activate the JAK-STAT pathway, inducing formation of the ISGF3 transcription complex, leading to expression of interferon-stimulated genes (ISGs) in similar fashion to type I IFN, but with a different temporal profile compared with those induced by type I IFN (15)(16)(17). A limited number of reports regarding IL-10-like signaling pathways have been published, with evidence of STAT3/5 biological activities (18,19). In fact, no reports have comprehensively investigated the downstream molecular signaling effects of IFN-.
To better understand molecular mechanisms underlying direct antiviral and immunomodulatory effects of IFN-3, a comprehensive catalogue of protein effectors regulated by IFN-3 was compiled using quantitative proteomics analysis in the well-established HBV-transfected hepatoblastoma cell line model viz. HepG2.2.15 cells (20). HepG2.2.15 cells have been widely used as a model of chronic hepatitis B because they support HBV replication and virion secretion (21)(22)(23)(24)(25)(26)(27). These new findings could improve the treatment strategy of HBV infection and expand potential applications of IFN-3.

EXPERIMENTAL PROCEDURES
Experimental Design and Statistical Rationale-Mixing of dimethyllabeled samples corresponding to IFN-3, IFN-␣2a, and control (PBS) treatments were performed at n ϭ 5 versus the typical n ϭ 3 to emphasize depth. Only biological replicates were performed. The normality of all proteomics, western-blotting and qPCR data was evaluated with the Shapiro-Wilk test. In all cases, the unpaired Student's t test or one-way ANOVA was selected when the distribution of the data was normal; otherwise the Mann-Whitney U test was applied. For proteomic analysis, peak intensity log2 ratios (L/M, L/H, M/H) were compared against a value of 0 (no change, log2(1)). Significance was based on the following criterion: p value Ͻ 0.05 or cases where proteins were detected in only one condition making standard statistical analysis inapplicable, with the requirement that these proteins must be identified in at least 4 of 5 experiments. Proteins not fulfilling the above criteria as well as those appearing in fewer than 3 experiments, were excluded from downstream analysis. Regarding western-blotting and qPCR experiments, n ϭ 3 was set. One-way ANOVA was performed for these experiments with p value Ͻ 0.05 considered significant. In the case of DAVID enrichment analysis, all proteins identified by mass spectrometry were input as background. When Fisher-based enrichment analysis was performed against an in-house database of proteomic/transcriptomic studies, resulting log (p values) were conservatively adjusted by subtracting the log value corresponding to the total of all possible study/study combinations in the database.
Cell Culture and Cell Stimulation-HepG2.2.15 cells are stable HBV-transfected cells derived from a hepatoblastoma HepG2 cell line (20). These cells contain a plasmid which expresses the complete genome of HBV, which integrates into host cellular DNA (20,28). Because HBV within HepG2.2.15 replicates and secretes HBsAg, HBeAg and HBV DNA into the culture media, HepG2.2.15 has been widely used as a model of chronic hepatitis B (22,28). The HepG2.2.15 cell line was kindly provided from Professor Antonio Bertoletti (Singapore Institute for Clinical Sciences, A*Star). These cells were maintained in Dulbecco's Modified Eagle's Medium (DMEM; Gibco, MA) supplemented with 10% Fetal Bovine Serum (FBS; Gibco), 1% MEM Non-Essential Amino Acids (MEM-NEAA; Gibco), 1% Penicillin/Streptomycin (Gibco) and Geneticin (G418; Gibco) at a final concentration of 150 g/ml. The cultured cells were grown in a humidified incubator at 37°C with 5% CO 2 . One million HepG2.2.15 cells were seeded into 6-well plates with 1 ml media and maintained in complete DMEM for 24 h. For determining the effects of IFN-3 on HBV replication, these cells were left untreated or treated with 1, 10, 100 or 1000 ng/ml of IFN-3 (5259-IL-025, R&D Systems, MN) and incubated for another 24 h. For proteomics analysis, HepG2.2.15 cells were plated at 5 ϫ 10 6 cells in T-75 flasks and grown in complete DMEM for 24 h at 37°C. These cells were subsequently cultured in media with 100 ng/ml of IFN-3 or 100 ng/ml of IFN-␣2a (11100 -1, pbl assay science, NJ) or PBS (control) for another 24 h. For further qPCR experiments with an optimized concentration of IFN-3, HepG2.2.15 cells were stimulated with or without 100 ng/ml of IFN-3 for 0, 8, 16, and 24 h.
RNA Isolation, Reverse Transcription, and qPCR for Gene Expression-TRIzol Reagent (Thermo, MA) was used to extract total RNA as specified in the accompanying manual. Complementary DNA (cDNA) synthesis was carried out using the Taqman Reverse transcription kit (Applied Biosystems, MA). Conditions for reverse transcription were as specified in the manual. Relative gene expression was measured with the ABI Prism 7500 sequence detection system (Applied Biosystems). All primers and probes were designed with the "primer express 3" program (Applied Biosystems), and are shown in supplemental Table S1. To monitor gene amplification, the intensity of fluorescence from the 18S housekeeping gene was monitored via Taqman probe whereas all target gene levels were monitored via SYBR green dye (Applied Biosystems). To test the specificity of SYBR green dye, melting curve analysis was conducted for all target genes. The condition of amplification for both target and housekeeping genes was 1 cycle at 95°C for 5 min followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. Relative gene expression was calculated with the 2 ϪddCt method. The Student's t test and one-way ANOVA were used to compare the relative expressions of target genes in cells treated with various doses of IFN-3. A p value less than 0.05 was considered significant.
Mini Kit (Qiagen, MD) according to manufacturer's instructions. Quantification of HBV viral load was performed by absolute quantitative real-time PCR using the ABI Prism 7500 sequence detection system (Applied Biosystems). For ease of plasmid amplification, we used plasmids containing only the HBV gene preS1 as a surrogate marker for HBV DNA (preS1 was chosen because it is highly conserved across all HBV genotypes). First, preS1 plasmids were extracted from E. coli transformants with the GeneJET Plasmid Miniprep Kit (Fermentas, MA). The concentration of extracted plasmids was measured by spectrophotometer (Nanodrop, Thermo) and copy/l was determined. The plasmid concentration was adjusted and diluted in a range of 10 7 , 10 6 , 10 5 , 10 4 , 10 3 , and 10 2 copy/l. These concentrations were used to construct a standard curve. Both standard and sample preS1 were amplified at the same time using conditions described above. The fluorescent intensities were specified as C t values. For standards, C t values at the above concentrations were used to plot a standard curve. This curve and sample C t values were used to calculate the amount of HBV DNA in the samples. The Student's t test and one-way ANOVA were used to compare viral loads in cells treated with various levels of IFN-3, as well as untreated cells. A p value less than 0.05 was considered significant.
MTT Assay-HepG2.2.15 cells were seeded in 96-well plates at a density of 1 ϫ 10 4 cells per well and incubated for 24 h. The culture media was removed and replaced with fresh complete media in the absence or presence of IFN-3 (1, 10, 100 and 1000 ng/ml). The cells were further incubated for 24 h followed by addition of 10 l 5 mg/ml MTT (Sigma, MO) solution in each well with gentle shaking. After 4h incubation, the resulting purple formazan crystals were dissolved with dimethyl sulfoxide (DMSO, Riedel-deHaë n, NJ) and subsequently measured with ELISA plate reader (Thermo) at wavelength 570 nm. The absorbance values of cells treated with each concentration of drug were compared with that of control (untreated) cells and % cell viability was calculated. This experiment was performed in triplicate.
Protein Extraction and In-solution Digestion-After trypsinization, the cells were lysed with 5% sodium deoxycholate (SDC) and 1X protease inhibitor (Thermo) mixture, followed by sonication. All cell debris was removed by centrifugation and supernatant protein concentrations of each sample were measured via BCA Protein Assay (Thermo). Equal amounts of protein from treated and untreated HepG2.2.15 were reduced and alkylated by dithiothreitol (DTT) treatment for 30 min at 37°C and iodoacetamide (IA) treatment for 30 min at room temperature in the dark, respectively. These samples were further quenched with DTT at least 15 min at room temperature before incubating with trypsin at a ratio of 1:50 at 37°C overnight. These mixtures were incubated with 0.5% trifluoroacetic acid (TFA) for 30 min and then centrifuged to remove SDC precipitate. The amount of tryptic peptide of each sample was determined with the Pierce Quantitative Fluorometric Peptide Assay (Thermo).
Dimethyl Labeling and Fractionation-Peptides from untreated, IFN-␣2a treated, and IFN-3 treated HepG2.2.15 were labeled with light reagents (formaldehyde and cyanoborohydride), medium reagents (formaldehyde-d 2 and cyanoborohydride), and heavy reagents (deuterated and 13 C-labeled formaldehyde and cyanoborodeuteride), respectively, for an hour at room temperature. Ammonia solution and formic acid (FA) were sequentially used to stop the reaction. Labeling efficiency was tested, and we found that greater than 99% of peptides were labeled (data not shown). After combining these three samples, the mixed labeled-peptides were dried in a SpeedVac centrifuge at room temperature. Next, the pooled peptides were separated into 10 fractions to reduce complexity using the Pierce High pH Reversed-Phase Peptide Fractionation Kit (Thermo). Eluates of each fraction were dried in a SpeedVac centrifuge before LC-MS/MS analysis.

LC-MS/MS and Analysis-
The fractionated samples were resuspended in 0.1% FA (Sigma) to a final volume of 15 l prior to MS injection. The peptides were then analyzed via an EASY-nLC1000 system (Thermo) coupled to a Q-Exactive Orbitrap Plus mass spectrometer (Thermo) equipped with a nano-electrospray ion source (Thermo). The peptides were eluted in 5-40% acetonitrile in 0.1% FA for 70 min followed by 40 -95% acetonitrile in 0.1% FA for 20 min at a flow rate of 300 nl/min. The MS methods included a full MS scan at a resolution of 70,000 followed by 10 data-dependent MS2 scans at a resolution of 17,500. The normalized collision energy of HCD fragmentation was set at 32%. An MS scan range of 350 to 1400 m/z was selected and precursor ions with unassigned charge states, a charge state of ϩ1, or a charge state of greater than ϩ8 were excluded. A dynamic exclusion of 30 s was used. The peaklistgenerating software used in this study was Proteome Discoverer™ Software 2.1 (Thermo). The SEQUEST-HT search engine was employed in data processing. MS raw data files were searched against the Human Swiss-Prot Database (20,219 proteins, June 2017) and the Hepatitis B Virus Swiss-Prot Database (225 proteins, June 2017), as well as a list of common protein contaminants (www.thegpm.org/ crap/). The following parameters were set for the search: (1) digestion enzyme: trypsin; (2) maximum allowance for missed cleavages: 2; (3) maximum of modifications: 4; (4) fixed modifications: carbamidomethylation of cysteine (ϩ57.02146 Da), as well as light, medium, and heavy dimethylation of N termini and lysine (ϩ28.031300, ϩ32.056407, and ϩ36.075670 Da); (5) variable modifications: oxidation of methionine (ϩ15.99491 Da). The mass tolerances for precursor and fragment ions were set to 10 ppm and 0.02 Da, respectively. Known contaminant ions were excluded. The Proteome Discoverer decoy database together with the Percolator algorithm were used to calculate the false positive discovery rate of the identified peptides based on Q-values which were set to 1%. The Precursor Ions Quantifier node in Proteome Discoverer™ Software was employed to quantify the relative MS signal intensities of dimethyl labeled-peptides. The control channels were used as denominators to generate abundance ratios of IFN-3/control and IFN-␣2a/control. Log2 of the normalized ratio was used to calculate the mean and standard deviation of fold change across all five biological replicates. When these ratios were found in less than three experiments, the relevant proteins were excluded. Significantly differentially regulated proteins were determined by Mann-Whitney U test and unpaired t-tests with p value Ͻ 0.05 considered significant.
Bioinformatics-We compiled a list of defense response to virus using a variety of resources as follows. The online resource Database for Annotation, Visualization and Integrated Discovery (DAVID, v6.8, https://david.ncifcrf.gov/) and Reactome (https://reactome.org/) were employed to classify the proteins regulated by IFN-3 into functional categories using all proteins identified by MS as background (for DAVID). We used terms such as "antiviral," "antigen processing/ presentation" to help extract a custom list of broad antiviral proteins. Additionally, the list contains proteins involved in the HBV life-cycle that were derived from manual literature curation. Further analysis of up-and down-regulated proteins was performed against an in-house database currently under assembly. The goal is simply to expand on DAVID's method by emphasizing data sets from individual MS and RNA-seq studies. Probabilities presented here are generated using Fisher's exact test with a background proteome size of 10,000 and are unadjusted. For consistency, lists of proteins with altered expression in our own HBV work are generated by filtering out all cases where differential expression is not accompanied by a p value Ͻ 0.05 (i.e. fold-change is not a factor); the lists are then subjected to Fisher analysis.

Validation of Responses to IFN-3 Treatment in HBV-transfected Hepatoblastoma
Cell Line Model-HepG2.2.15 is a hepatoblastoma cell line that contains a stable HBV expression plasmid that has been validated as a chronic HBV infection model in previous reports (20,28). The response to type III IFN treatment in HepG2.2.15 cells has not been reported. To determine whether HepG2.2.15 cells respond to IFN-3, we performed qPCR to investigate the expression of the classical ISGs, namely OAS1, MX1 and ISG15 in HepG2.2.15 cells treated with various amounts of IFN-3 for 24h. Fig. 1 shows that IFN-3 could significantly increase the expression of these 3 ISGs in a dose-dependent manner. These results indicated that HepG2.2.15 cells responded to IFN-3 stimulation.
For thoroughness, we investigated the anti-HBV effects of IFN-3 by determining the differential changes at 3 points in the viral life-cycle including the levels of preS1 (typically used as a representative gene for HBV transcripts, given its high conservation across all genotypes), replicative intermediate pre-genomic RNA (pgRNA), and intracellular HBV DNA (both rcDNA and cccDNA). qPCR was performed on HepG2.2.15 RNA following treatment with various amounts of IFN-3 for 24 h. As shown in Fig. 2A and 2B, IFN-3 reduced both preS1 and pgRNA expression compared with control in a dose-dependent manner. Measurement of intracellular HBV DNA showed that copy numbers of virus in IFN3-treated HepG2. 2.15 cells were diminished in a dose-dependent manner relative to control (Fig. 2C). The reduction reached significant levels when the doses of IFN-3 were 100 ng/ml (p value ϭ 0.04) and 1,000 ng/ml (p value ϭ 0.0134). Collectively, these results indicated that IFN-3 inhibits HBV replication in HepG2.2.15 at the given time point.
Before we investigated the cellular response to IFN-3, the toxicity of this drug was considered. The MTT cytotoxicity assay was performed to determine HepG2.2.15 viability under distinct concentrations of IFN-3. The percentage of viable cells is illustrated in Fig. 3. The increasing doses of IFN-3 significantly promoted cell death only at the highest two doses in this experiment, where the viable HepG2.2.15 cells reduced to 94% and 84% in 100 ng/ml and 1000 ng/ml conditions, respectively. Thus, we settled on 100 ng/ml of IFN-3 for further experiments because this dose showed the ability to significantly inhibit HBV replication with minimal cytotoxicity on HepG2.2.15 cells.
Quantitative Proteomics Analysis of IFN-3 Responses in HepG2.2.15- Fig. 4 shows the schematic workflow of this study. Briefly, untreated, IFN-␣2a-treated, and IFN-3-treated HepG2.2.15 cells were lysed and digested with trypsin. The tryptic peptides of these groups were labeled with light, medium, and heavy dimethyl reagents and then fractionated for subsequent LC-MS/MS analysis. In total, 4670 proteins were identified at a false discovery rate (FDR) of less than 1%, with 1471 proteins identified based on a single peptide (all information regarding peptide sequences assigned and protein identifications are supplied in supplemental Table S2 and S3). For the IFN-3 treatment condition, 2904 proteins were identified in at least 3 of 5 replicates allowing evaluation of significance, shown in the corresponding volcano plot (Fig. 5). Seven hundred thirty-seven proteins showed significant changes in abundance, with a slight bias toward down-regulated proteins in response to IFN-3 stimulation (see Fig. 5). Table I displays a list of significantly regulated proteins with ͉log 2 (IFN-3/Ctrl)͉ ratios of greater than 1 . We should point out that our primary intention in this work is to explore IFN-3 effects on HepG2.2.15 cells; the effects of IFN-␣2a would be secondary, for the sake of comparison with IFN-3 effects, and will be presented at the last part of the Results section. All data can be accessed from supplemental Table S4.
All HBV proteins were identified in this study. However, calculation of significance could not be performed because of absence in multiple MS runs, except for putative X-Core fused protein (which showed no significant change), see supplemental Table S4. Based on spectral counts, all the HBV peptides were apparently expressed at very low levels relative to the entire HepG2.2.15 proteome (below 1 ppm, see in PRIDE partner repository as mentioned in Data Availabilty). This finding is in agreement with a previous study that specifically examined intracellular HBV proteins in HepG2.2.15 cells but found them to be undetectable (29). Note that all HBV proteins, except HBV pol and HBx, are secreted (30 -32), likely causing intracellular levels of these proteins to be scarce.
To further explore the possible roles of transcriptional regulation for significantly altered proteins, we selected several proteins involved in antiviral processes for qPCR analysis. Fig.  7 show that up-regulation was seen at both protein and RNA expression levels without exception. However, no clear pattern emerged when down-regulated proteins were examined for transcript levels.
Bioinformatics Analysis and Antiviral/Immunomodulatory Process Mapping-The DAVID bioinformatics tool was used to classify and cluster significant functions of all up-and down-regulated proteins. Upon IFN-3 treatment, we found several biological processes expected to be canonically regulated on general IFN stimulation, such as immune response and viral infection defense. Additionally, however, several biological processes not previously emphasized upon interferon stimulation emerged from this work, including metabolism of RNA, major pathway of rRNA processing in the nucleolus and cytosol, selenoamino acid metabolism, peptide chain elongation, translation, prefoldin mediated transfer of substrate to CCT/TriC, transcription-coupled nucleotide excision repair (TC-NER), unfolded protein response (UPR), interleukin-23 signaling, and hypusine synthesis from eIF5A-lysine.
Following enrichment analysis, we conducted a protein-byprotein search for relevance to the viral processes. We manually searched the viral literature, as well as using DAVID's antiviral defense groups and antigen processing/presentation groups, to construct a map of points at which these proteins may partake in viral processes (Fig. 8) (as described under Experimental Procedures). We also built a map displaying steps at which IFN-3 treatment could promote antigen processing and presentation of viral proteins. Finally, we constructed a map that highlights points at which IFN-3 stimulation may heighten expression of proteins in the RIG-I pathway that have been shown to be inhibited via HBV infec- tion. Unless otherwise stated, we found no proteins that contradict the patterns that we report below.
Regarding the HBV life cycle, we identified many proteins which may be involved at numerous points ( Fig. 8). At the early steps of viral attachment, entry, and cytosolic release from endosomes (Fig. 8, step 1 and 2), we detected an increased expression of IFITM3, known to alter intracellular cholesterol homeostasis, preventing viral fusion with endosomes in at least 4 out of 7 of the Baltimore viral groups (39). RAB5C and RAB7A, proteins previously implicated in viral trafficking (Fig. 8, step 2), were down-regulated upon IFN-3 treatment (40), though RAB7A did not reach statistical significance. Two other trafficking proteins not mentioned in the viral literature, AP2B1 and EEA1, were significantly downregulated. Regarding cytosolic-to-nuclear transport of the HBV capsid along microtubules (Fig. 8, step 4), we found that several molecules involved in microtubule assembly and function evinced decreased expression after IFN-3 treatment, consistent with the literature. These proteins were DCTN1, DCTN2, KIF5B, MAP4, and MACF1 (41,42). For nuclear import (Fig. 8, step 5), all identified proteins known to play roles herein (42,43) were found to decrease because of IFN-3 Within the nucleus (Fig. 8, step 7), all identified ISG products shown previously to limit or degrade viral mRNA (44 -46) were observed to be up-regulated: ADAR1, OAS3, ISG20, PML, XRN2, ZC3HAV1, and TRIM25. At the level of RNA export (Fig.  8, step 8), ENY2 and CALR were down-regulated. Regarding translation (Fig. 8, step 9), antiviral proteins (EIF2AK2, IFIT1, and IFIT2) involved in protein synthesis increased in treated cells, as shown in previous work with type I IFN treatment (44,46,47). In addition, almost all identified translation initiation factors (EIF1AY, EIF2S3, EIF2B1, EIF3B, EIF4G1, EIF4B, EIF4G2, and EIF5A) decreased in this study. Encapsidation follows translation of viral proteins (Fig. 8, step 10). Interestingly, RIG-I was up-regulated and has been shown to interfere with the interaction between pgRNA and HBV polymerase, required for encapsidation (48). Other effects of RIG-I upregulation are discussed below. Reverse transcription, which occurs within the capsid (Fig. 8, step 11), has been demonstrated to be inhibited by SAMHD1 (HBV and HIV) (35,49) and NT5C3A (IAV and CSFV) (50,51); these proteins were upregulated in our study upon treatment. Proteins that might interfere with viral assembly (Fig. 8, step 12.1) were not identified in this study. Finally, at the egress step (Fig. 8, step 13), two up-regulated proteins thought to play a role in inhibiting HBV export, BST2 and C4A, were detected after IFN-3 stimulation.
Regarding antigen processing and presentation, we identified 22 proteins that are involved in this process (Fig. 9). Among proteins in the degradation-related ubiquitination cascade, we found only up-regulation, specifically UBA7 (E1), UBE2L6 (E2), and SYVN1 (E3) after IFN-3 stimulation. Two ubiquitination enzymes were down-regulated but have not been reported to instigate degradation. All deubiquitinating enzymes (DUBs) identified in this study, i.e. USP14, UCHL3, UCHL5, and OTUB1, decreased because of IFN-3 treatment. Downstream of these processes, core proteasome components, PSMA7, PSMB7, and PSMG2, were up-regulated with no exceptions. Importantly, both cap components (PSME1 and PSME2) associated with the immunoproteasome were up-regulated, whereas all significantly altered cap components of the constitutive proteasome were down-regulated. Evidence of possible post-proteasome processing was seen in the up-regulation of a cytosol aminopeptidase, LAP3. Next, we saw up-regulation of proteins involved in transport of peptides from the cytosol to the ER as well as assembly of peptide-HLA class I complexes; these proteins included TAP2, TAPBP, HLA-A, HLA-B, and HLA-C. Other proteins involved in the above step did not reach statistical significance but tended to increase expression following IFN-3 treatment; TAP1, PDIA3 (ERP57), and B2M. CALR, a player in this process, was down-regulated, the reasons for which will be discussed later.
One of the important functions of the innate immune response against HBV is the sensing of viral RNA, leading to IFN activation, which includes components such as toll-like receptors, RIG-I, and MDA5. Here, we saw significant up-regulation of RIG-I upon treatment, which could counteract the known suppression of RIG-I by HBV. Enhancers of RIG-I activity, such as OASL, TRIM25, and IFIT3 (RIG-G) were also significantly up-regulated. Increased levels of type I and III IFNs, which are activated by RIG-I, would be expected, but were not observed via MS. However, qPCR work confirmed significant up-regulation of these IFNs at the transcript level (Fig. 10). It is likely that the increase in IFN I and III gene expression because of RIG-I up-regulation could lead to more production of IFNs at the protein level and in turn activate IFN receptors as positive feedback. All three canonical ISGF3 components (STAT1, STAT2, and IRF9) downstream from IFN receptor activation were up-regulated, probably as a result of both initial IFN-3 treatment and further generation of type I and III IFNs (Fig. 11).
We found several likely antiviral proteins that do not conveniently fit into the above viral processes. Firstly, three 14-3-3 proteins (YWHAZ, YWHAH, and SFN) were found to be down-regulated following IFN-3 treatment.
Another down-regulated protein of interest was cyclophilin D (PPIF). Finally, CHID1 was the single most up-regulated protein in our study on IFN-3 treatment. Though not prominent in the viral literature, CHID1 is known to bind LPS (52) and is seen to be down-regulated in several viral studies (53)(54)(55)(56).
After removal of all antiviral proteins, enrichment analysis of the remainder produced an interesting result. RNA-binding proteins, particularly spliceosome components, were both up-and down-regulated on IFN-3 treatment versus control. For example, the splicing factor U2AF1 was significantly upregulated on treatment (p ϭ 0.00165). Searching through individual data sets from an in-house database containing a  6. Validation of altered proteins by immunoblotting assay. SDS-PAGE was performed on treated and untreated lysates, followed by membrane transfer and incubation with anti-OAS3, anti-SAMHD1, anti-STAT1, and anti-GAPDH overnight. The proteins of interest were visualized using the LI-COR Odyssey system. Consistent with proteomic results, the expression of OAS3, SAMHD1 and STAT1 increased after IFN-3 treatment. These experiments were performed in triplicate. broad modes of action between the two IFNs, significant overlaps should be seen in sets of up-and down-regulated proteins upon treatment, offering a simple validation of IFN-3's activity. At the same time, and most intriguingly, the comparison would allow us to suggest possible differences in effect. To take advantage of the depth of our data, we constructed a tree depicting all 27 possible differential expression outcomes resulting from 3 pairs of treatment comparisons (i.e. IFN-3 versus Ctrl, IFN-␣2a versus Ctrl, IFN-3 versus IFN-␣2a, see Fig. 12). A simple comparison between IFN-3 versus IFN-␣2a could lead to errors in interpretation. For example, the observation that IFN-3 treatment results in upregulation or down-regulation of a protein versus IFN-␣2a renders an uninteresting result if IFN-3 does not have any significant effect versus Ctrl. This scenario can be seen in outcomes #12 and #16 in Fig. 12. Another problem in interpretation arises when both IFNs have the same directional effect on cells but one effect is significantly more potent than the other; this would not be noticed upon mere IFN-3 versus IFN-␣2 comparison. This scenario can be seen in outcomes #3 and #25 in Fig. 12. The tree shows the number of MSderived proteins at their appropriate outcome branch points and shows proteins associated with various biological processes (as described under "Experimental procedures") for each outcome. As one simple example of the power of this approach, note that IFITM3 was significantly up-regulated upon treatment with both IFNs, however this effect was significantly stronger upon IFN-3 treatment (outcome #1). Another example was in outcome #4 where IFN-3 treatment caused up-regulation of proteins involved in antigen processing and presentation (SYVN1, LAP3, PSMA7, and PSMB7), whereas IFN-␣2a did not produce any significant effects. Interestingly, the mirror image of this branch, outcome #24 where IFN-3 caused down-regulation whereas IFN-␣2a did not, was also populated with a subset of proteins involved in antigen processing and presentation, specifically those associated with the constitutive proteasome cap (PSMC4, PSMD4, PSMD7, PSMD13, and PSMD14). A priori, one would expect the most interesting outcomes to be found in cases in which IFN-␣2a and IFN-3 have clearly opposite effects (outcomes #7 and #21). We did not identify proteins clearly related to antiviral processes in these two outcome groups. This reason pointed to broader unbiased analysis. Hence, we per- Ubiquitinating enzymes that promote protein degradation were found to be up-regulated whereas deubiquitinating enzymes known to remove ubiquitin were found to decrease after IFN-3 treatment. Although all identified subunits of the constitutive proteasome cap were down-regulated, the 2 subunits of the immunoproteasome cap as well as some subunits of the proteasome core were elevated in expression. Several effector molecules involved in peptide loading on class I HLA as shown in red text were up-regulated with the exception of CALR, which was down-regulated as shown in green text. The italic red text represents proteins with increased expression levels that failed to reach statistical significance.

Up-regulation
formed "Reactome" enrichment analysis against all outcome groups for complete enrichment analysis beyond the aforementioned antiviral processes. The resulting Reactome enrichment groups were subjected to 2-D cluster analysis ( Fig.  13 and underlying data in supplemental Table S4). This analysis offers new insights. Firstly, we note the aforementioned cases in which IFN-3 and IFN-␣2a showed opposing effects. In outcome #7, where IFN-3 exclusively caused up-regulation, the "peptide chain elongation," "metabolism of RNA," "major pathway of rRNA processing in the nucleolus and cytosol," "selenoamino acid metabolism," and "translation" terms are highlighted, whereas in outcome #21, where IFN-3 exclusively caused down-regulation, the "interleukin-23 signaling" term is emphasized. Next, we observe that in most cases where proteins were up-regulated upon IFN-␣2a treatment (outcomes#1, #2, #3, #11, and #21), "interferon signaling" and/or "class I MHC mediated antigen processing & presentation" terms were apparent, reinforcing expectations regarding IFN-␣'s potency on these pathways. Interestingly, specific terms related to "translation" and "metabolism of RNA" were significantly associated with 13 out of 17 different outcome groups. To our best knowledge, these biological processes have been not mentioned as consequences of general IFN activation. Finally, enrichment terms that are largely unexplored with regard to IFN treatment include "tran-scription-coupled nucleotide excision repair (TC-NER)" (1 significant occurrence), "interleukin-23 signaling" (1), "unfolded protein response (UPR) (2)," "prefoldin mediated transfer of substrate to CCT/TriC (4)," and "hypusine synthesis from eIF5A-lysine (1)." Interestingly, the final enrichment term was only associated with a single case, that of proteins downregulated in all treatment pairs (outcome #27). These results could help point toward differential effects of IFN-3 versus IFN-␣2a. DISCUSSION Type III IFN has been shown to exert antiviral effects on several viruses including EMCV, IAV, HSV, VSV, HIV, HCV, and HBV (17, 46, 59 -64). In clinical trials, peg IFNtreatment showed reduced HBV DNA, HBsAg and HBeAg in CHB patients at levels greater than or comparable to those treated with peg IFN-␣ during on-treatment and at end-of-treatment (9). Although the virological, serological, and biochemical responses of peg IFNat week 24 after treatment were not superior to peg-IFN-␣, adverse events associated with drug treatment were mostly seen in patients treated with peg-IFN-␣. Numerous HCV studies also illustrate this tendency toward diminished side-effects (13,14). Another clinical trial demonstrated that patients treated with peg IFNhad improved anti-HBV immunity through an increase in poly-func-  11. IFN-3 rescued the RIG-I signaling pathway. IFN-3 up-regulated the expression of RIG-I and IFIT3, which have been reported to be suppressed by HBV. IFN-3 also elevated the expression of OASL and TRIM25. These proteins promote type I and type III IFN production. These IFNs, in turn, activate the JAK-STAT pathway and induce the expression of ISGs. It is likely, then, that IFN-3 provides positive feedback to amplify ISG expression to control HBV replication. Proteins in red text represent up-regulated proteins, whereas the blue text represents the up-regulated transcripts.  tional NK cells and maintenance of both HBV-specific CD4 ϩ and CD8 ϩ T cells, important arms of immunity for viral elimination (8). The mechanisms involved in these immunomodulatory effects remain to be elucidated. At the level of cell culture, only two studies (59,65) have focused on the effects of IFN-treatment of HBV-transfected cells, one of which offers a comparison against a type I IFN, IFN-␤. Here, a comparison against IFN-␤ showed that IFN-␤ exhibited greater potency in viral elimination at the 6 h time point, but potency was equivalent at 24 h. Surprisingly, the effects of IFN-treatment have not been studied with high-resolution mass spectrometry to date, possibly because it is assumed that type I and type III IFNs share the same signaling pathways. In fact, modern, deep proteomic studies of the effects of IFN treatment in general are lacking, perhaps because of an underlying assumption that RNA-seq results should parallel those generated via MS. Here, we used a modern high-resolution LC-MS/MS system with upstream high-pH reversed phase fractionation to yield a total of 127,989 peptides (25,181 unique peptides) representing 4670 proteins. The peptides were labeled with dimethyl isotopes, allowing accurate quantification under different treatments. This study provides the most comprehensive examination of changes in protein abundance under IFN treatments to date.
To provide objective measures for study-study comparisons, in the discussion below we perform numerous comparisons between our studies and others with Fisher's exact test (see "Experimental procedures"). As expected, canonical IFNstimulated biological processes were prominent in IFN-3 treatment. Proteins significantly up-regulated upon IFN-3 treatment versus Ctrl (hereafter 3/c) overlapped with the GO "antiviral defense" term at p value ϭ 10 Ϫ14 (with the ␣/c condition giving p value ϭ 10 Ϫ12 ). Similarly, the 3/c list significantly overlapped (p value ϭ 10 Ϫ13 ) transcripts whose up/down-regulation patterns (i.e. "co-expression") mirror those of RIG-I transcripts reported in the ARCHS4 resource (66); the ␣/c list gave p value ϭ 10 Ϫ18 , perhaps illustrating the expected difference in potency. Comparing proteins significantly up-regulated in 3/c against those up-regulated in ␣/c, we derive p ϭ 10 Ϫ138 , a clear demonstration of similarities in effect versus control. Stepping outside the bounds of our own work, Bolen's transcriptomic comparison of the effects of 5 interferons, including IFN-3, on Huh7 cells can be examined against our proteomic results (16). Here, a list of highly upregulated transcripts derived from Bolen's 6, 12, and 24 h 3/c time-points significantly intersected with our own 3/c FIG. 12. Outcome tree. To compare the effects of IFN-3, IFN-␣2a, and PBS control, we generated a tree representing the 27 possible outcomes. The boxes in red and green represent proteins that are significantly up-and down-regulated, respectively. Groups of proteins with insignificantly changed expression are represented by yellow boxes. As an example, the topmost outcome represents proteins in which IFN-3 treatment causes significant up-regulation versus control, IFN-␣2a treatment causes significant up-regulation versus control, and IFN-3 treatment causes significant up-regulation versus IFN-␣2a. The number in parentheses refers to the number of identified proteins in the group. The proteins in each group that might be involved in suppressing HBV replication are shown in colored boxes specifying various processes, particularly antiviral.  13. Enrichment heatmap. Reactome enrichment analysis was performed on the tree outcomes that contained at least 9 members, excluding outcome #14 in which no significant alterations were seen on any treatment. Probabilities (Reactome's hypergeometric test, -log(p value)) were generated for every intersection of these tree outcomes with the specified Reactome groups. This matrix of probabilities was then submitted for two-dimensional cluster analysis using Cluster 3.0 and Java Treeview. The result is displayed as a heatmap. Red color intensity indicates the level of significance. results (p value ϭ 10 Ϫ25 ). Of interest, the degree of overlap (p value ϭ 10 Ϫ5 ) was less significant when performing the same exercise with down-regulated transcripts versus proteins, paralleling the tendency of lower correlation seen in our own qPCR results for several down-regulated proteins (Fig. 7). Other comparisons illustrating the expected canonical effects would be an ␣/c result in Huh7 (67) (a 3/c experiment was not performed) where up-regulated transcripts align with our own 3/c results (p value ϭ 10 Ϫ30 ), and a list of genes likely to be activated by STAT2 based on curated STAT2 chip-seq results that significantly overlaps with our 3/c results (p value ϭ 10 Ϫ25 ) (68).
In addition to acting on viral replication directly, we propose that the proteolytic pathway may be induced by IFN-3 to inhibit HBV replication. Robek et al. demonstrated that small molecule-based inhibition of proteasome activity could restrain the anti-HBV activity of type I IFN in vitro, suggesting that antiviral effects of IFN might relate to the proteasome (69). Yao et al. demonstrated that suppression of HBV replication was accompanied by an increase of 5 proteasome subunits in HepG2.2.15 cells in response to interleukin-4 (IL-4) treatment (70). In our study, we also note up-regulation of immunoproteasome members (cap components PSME1 and PSME2) on IFN-3 and IFN-␣2a treatment. The immunoproteasome is distinct from the constitutive proteasome in the composition of its catalytic subunits, its cap, and its peptide processing capacity (71). Interestingly, all significantly altered cap components of the constitutive proteasome were downregulated only upon IFN-3 treatment, suggesting coordinated modulation toward the antigen processing form of the proteasome on IFN-3 treatment. To our knowledge, this is the first report suggesting a shift from proteasome constitutive degradation mode to antigen processing/presentation mode in response to type III IFN treatment. HBV has been termed a "stealth" virus because innate immune responses in chimpanzees are not observed shortly after infection (72). The synthesis of type I IFNs is induced by the activation of RIG-I-like receptor (RLR) signaling pathways in response to viral infection. Previous studies have shown that HBV polymerase and HBx proteins interfere with the interaction of RIG-I and downstream signaling molecules, resulting in diminished type I IFN production (73)(74)(75).
HBV has also been shown to alter RIG-I expression by inducing up-regulation of miR146a, which may directly downregulate RIG-I and RIG-G transcripts (76). IFN treatment has not previously been shown to up-regulate RIG-I pathway components in HBV infected cells. We found, for the first time, that treating HepG2.2.15 with IFN-3 increased the expression of RIG-I. In addition, upon IFN-3 treatment, we observed elevated expression of RIG-G, known to enhance the activity of RIG-I, as well as TRIM25 and OASL, which promote the interaction between RIG-I and IPS1 on mitochondria or peroxisomes, leading to increased IFN production (77,78). We propose that IFN-3 restores RIG-I and associated down-stream signaling to produce type I IFN in HBV infected cells and likely other viruses. Although we did not find any IFNs in our proteomic work, possibly because of low abundances, qPCR verified increases in type I IFN (IFN-␤) and type III IFN  (IFN-1 and IFN-3) upon IFN-3 stimulation. Several biological processes not mentioned above emerged from this work, including the involvement of 14-3-3 proteins, cyclophilin D, calreticulin, and cell motility proteins. These subjects are explored in detail below. In our study, the 14-3-3 proteins YWHAZ, YWHAH, and SFN were found to be down-regulated following IFN-3 treatment. Several studies have shown that 14-3-3 proteins may facilitate the replication of viruses (79,80). The 14-3-3 proteins are highly conserved regulatory molecules having multiple functions, including involvement in signal transduction via kinases and phosphatases, the cell cycle, and apoptosis. These proteins have been shown to play important roles in viral infection. For example, complex formation of cdc25 with the HIV accessory protein Vpr, an event that alters the host cell life cycle, is facilitated by 14-3-3 proteins (79). Also, HCV core protein interacts with 14-3-3 proteins resulting in enhancement of Raf-1 kinase activity together with control of hepatocyte growth (80). Interestingly, sequence analysis revealed a 14-3-3 binding domain in HBx protein (81), suggesting that 14-3-3 proteins could directly interact with HBV. Notably, we could not observe significant alterations in 14-3-3 proteins upon IFN-␣2a treatment; the differential effects of IFN-versus IFN-␣ on 14-3-3 levels would be worthy of further investigation.
The down-regulation of cyclophilin D upon IFN-3 treatment observed in our work suggests that this molecule may play a role in inhibition of HBV replication. One of the pathways affected by HBx is calcium signaling, as Ca 2ϩ is necessary for viral replication and core assembly (82,83). In addition to interference with RIG-I signaling, the HBx protein can interact with and modulate the mitochondrial permeability transition pore (MPTP) leading to the release of mitochondrial calcium into the cytoplasm (82). The opening of the MPTP anion channel is regulated by cyclophilin D (84), thus downregulation of cyclophilin D would serve to counteract Ca 2ϩassisted viral replication and core assembly. Increasing levels of cytosolic calcium causes the stimulation of PYK2 kinase, in turn activating Src kinase signal transduction to promote HBV reverse transcription, replication, and core assembly (83,85).
In several studies, cyclosporine A (CsA) was used as an MPTP specific blocker, thus reducing core assembly and inhibiting HBV replication (30,84,85). Consistent with CsA treatment, the calcium ion chelating agent BAPTA-AM also suppressed HBV propagation (85).
Calreticulin (CALR), down-regulated upon IFN-3 treatment in our study, is an ER calcium-binding chaperone involved in the regulation of calcium homoeostasis, the folding of newly synthesized proteins, nuclear import, and peptide loading on MHC class I (86 -88). In HBV-transfected cells, Yue et al. demonstrated that the inhibition of IRF-7 translocation into the nucleus was induced by CALR, resulting in suppression of type I IFN production (89). Moreover, this group showed that CALR inhibited JAK-STAT pathway induction by IFN-␣ by inhibiting STAT1 phosphorylation, diminishing expression of PKR and OAS. Thus, the observation of down-regulation of CALR would be expected to have an antiviral effect. However, the down-regulation of CALR would also seem to interfere with antigen processing/presentation. We note that most antigen processing/presentation molecules were up-regulated upon treatment in our work, possibly compensating for this effect. Overall, the observed IFN-3-induced down-regulation of CALR in our study may have an inhibitory effect on HBV replication. Enrichment analysis (DAVID) showed that the most significantly down-regulated proteins following IFN-3 treatment were involved in cell motility and cell adhesion. We believe these processes could be relevant to the effects of IFN-3 treatment on liver cancer. Tan et al. demonstrated that HBVtransfected cells showed morphological changes because of formation of filopodia and lamellipodia, leading to cell migration. Of interest, CHB patients have 5-15-fold increased risk of HCC, a characteristic of which is increased cell motility (90,91). CDC42, down-regulated on IFN-3 treatment, has been reported to control cell proliferation, adhesion, and metastases. High expression of this protein was observed in several types of cancers including HBV-related HCC. Inhibition of CDC42 by CRISPR/Cas9 knockout and treatment with a specific CDC42 inhibitor in HBx-Huh7 cells (92) reduced cell proliferation and promoted apoptosis in these cells. In addition to CDC42 down-regulation, we also found that its downstream effectors, actin and IQGAP1, showed decreased expression following IFN-3 treatment. Also, proteins involved in cell movement and focal adhesion such as ACTN1, ACTN4, RDX, CFL1, PIP4K2C, and ROCK2 were down-regulated.
We have shown that IFN-3's effects were largely in accord with those expected of canonical IFN activity. Nevertheless, our tree and cluster analyses ( Fig. 12 and 13) highlight several differential effects between IFN-3 and IFN-␣2a worthy of further investigation. Most obviously, despite the division of treatment outcomes into 27 non-intersecting protein groups, unbiased Reactome enrichment analysis repeatedly elicited an RNA-metabolism/translation theme. This broad theme can be divided into Reactome sub-groups (e.g. metabolism of RNA, major pathway of rRNA processing in the nucleolus and cytosol, translation, transcription-coupled nucleotide excision repair (TC-NER), and peptide chain elongation). However, we could not specifically find IFN-3-or IFN-␣2a dominant patterns in these sub-groups, suggesting that investigation of this broad theme might best be undertaken at the level of individual proteins. Perhaps the strongest general observation would be a tendency for involvement of ER-associated proteins (GO "establishment of protein localization to endoplasmic reticulum," which is primarily composed of ribosomal proteins) in cases where IFN-3 treatment results in significantly greater protein abundance versus IFN-␣2a; in the 5 outcomes where such an effect would be discernable, -log(p) ϭ 4.7, 3.2, 4.8, 8.2, and 5.3. On the other hand, such significance (-log(p) ϭ 4.1) is seen only once out of six outcomes where IFN-␣2a treatment results in significantly more abundant protein expression. This result could imply differential regulation of these ER targeting proteins on IFN-3 versus IFN-␣2a treatment. Altered regulation of these ER proteins could interfere with viral replication (93,94). It is interesting to note that some viruses utilize the ER as a replication compartment, whereas others do not. Thus, our work suggests that in addition to IFNЈs role in regulating ER-targeting, one IFN treatment might be more appropriate than another for infections. Of other Reactome groups that emerged during clustering, we find the "hypusine synthesis from eIF5A-lysine" case especially interesting, as it is only associated with a single outcome group, that in which 3Ͻc, ␣2a Ͻc, 3Ͻ ␣2a (i.e. the two IFNs cause down-regulation, but IFN-3's effect is most potent). Hypusine, a rare amino acid, is apparently found only once in the human proteome, as a modification of eukaryotic initiation factor 5 (eIF5A), with proviral implications for HIV and Ebola (95). Thus, IFN-3's antiviral effect may be especially potent in down-regulating hypusination. The case of hypusination also illustrates how IFN-3-induced alterations to the RNA-metabolism/translation machinery could target viruses.
At the level of individual studies, we note that, in addition to the aforementioned 3/␣2a comparison, Bolen's data (16) also allows a 3/2 comparison; here, up-regulated transcripts aligned strongly with our own 3/c up-regulated proteins (p value ϭ 10 Ϫ27 ), suggesting possible differential effects even within the IFN-subgroup. As a side note, IFITM3, strongly up-regulated in our own 3/␣2a comparison (outcome #1), was the single most significantly up-regulated transcript in the 3/2 comparison (p value ϭ 10 Ϫ10 ). Having noted possible IFN-3 versus IFN-␣2a differences, we should reiterate the points that: 1) our experimental design did not allow for discrimination of kinetic effects (e.g. those caused over time by differing receptor/ligand affinities and half-lives) versus substantial alterations at the level of pathways and 2) the effects of IFN-3 and IFN-␣2a were indeed substantially similar, especially with regard to canonical IFN effects such as regulation of antiviral proteins, antigen processing/presentation, and the RIG-I signaling pathway. A recent work does make a strong case for a differential 3/␣2a effect on murine intestinal epithelial cells (96), but the effect is pronounced only on polarization of cells, and the underlying microarray data does not correlate with our proteomic data. Thorough studies focused on possible differential effects on PTMs may help resolve these questions; the immediate result of IFN-receptor/ ligand interaction is, after all, a phosphorylation event. Comprehensive identification of phosphorylation events and other kinds of PTMs will help us to elucidate possible differential effects in different IFN treatment.
In conclusion, our study found that IFN-3 exhibited anti-HBV activities by significant inhibition of HBV replication and expression of HBV RNA. We used high-throughput quantitative proteomics to obtain a comprehensive understanding of molecular events upon treatment of HepG2.2.15 with IFN-3. To our knowledge, in fact, this study is the most comprehensive proteomics-based analysis of IFN treatment to date. For the first time, we reported significant up-regulation of immunoproteasome components, restoration of HBV-inhibited RIG-I pathway proteins, as well as several proteins not previously associated with IFN-3 treatment. Further study of these proteins is required. Clearly, IFN-3 exhibited both antiviral and immunomodulatory effects to inhibit HBV replication; therefore, IFN-3 is an attractive novel candidate for CHB treatment and the altered proteins might be new therapeutic targets in CHB infection.

DATA AVAILABILITY
The mass spectrometry proteomics data, including annotated spectra for all modified peptides and proteins identified on the basis of a single peptide, have been deposited to: (1) the ProteomeXchange Consortium via the PRoteomics IDEntifications (PRIDE) partner repository with the data set identifier PXD007896 and (2) the MS-Viewer (http://msviewer. ucsf.edu/prospector/cgi-bin/msform.cgi?formmsviewer) with the following keys: dhjzinh2g0, tjnx2fzkzu, 3xzzxfanwm, ac4wmxx0tv, and ao9nga6qqm.