Transcriptional profiling of the effects of 25-hydroxycholesterol on human hepatocyte metabolism and the antiviral state it conveys against the hepatitis C virus

Background Hepatitis C virus (HCV) infection is a global health problem. A number of studies have implicated a direct role of cellular lipid metabolism in the HCV life cycle and inhibitors of the mevalonate pathway have been demonstrated to result in an antiviral state within the host cell. Transcriptome profiling was conducted on Huh-7 human hepatoma cells bearing subgenomic HCV replicons with and without treatment with 25-hydroxycholesterol (25-HC), an inhibitor of the mevalonate pathway that alters lipid metabolism, to assess metabolic determinants of pro- and antiviral states within the host cell. These data were compared with gene expression profiles from HCV-infected chimpanzees. Results Transcriptome profiling of Huh-7 cells treated with 25-HC gave 47 downregulated genes, 16 of which are clearly related to the mevalonate pathway. Fewer genes were observed to be upregulated (22) in the presence of 25-HC and 5 genes were uniquely upregulated in the HCV replicon bearing cells. Comparison of these gene expression profiles with data collected during the initial rise in viremia in 4 previously characterized HCV-infected chimpanzees yielded 54 overlapping genes, 4 of which showed interesting differential regulation at the mRNA level in both systems. These genes are PROX1, INSIG-1, NK4, and UBD. The expression of these genes was perturbed with siRNAs and with overexpression vectors in HCV replicon cells, and the effect on HCV replication and translation was assessed. Both PROX1 and NK4 regulated HCV replication in conjunction with an antiviral state induced by 25-hydroxycholesterol. Conclusion Treatment of Huh-7 cells bearing HCV replicons with 25-HC leads to the downregulation of many key genes involved in the mevalonate pathway leading to an antiviral state within the host cell. Furthermore, dysregulation of a larger subset of genes not directly related to the mevalonate pathway occurs both in 25-HC-treated HCV replicon harbouring cells as well as during the initial rise in viremia in infected chimpanzees. Functional studies of 3 of these genes demonstrates that they do not directly act as antiviral gene products but that they indirectly contribute to the antiviral state in the host cell. These genes may also represent novel biomarkers for HCV infection, since they demonstrate an outcome-specific expression profile.


Background
Hepatitis C virus (HCV) infection is a global health problem and a leading cause of liver disease in North America [1,2]. Previous studies have indicated that genes associated with cellular metabolic pathways are involved in the HCV life cycle. HCV infection causes the accumulation of lipid droplets (LDs) in patients' hepatocytes, a process called steatosis [3]. Both the HCV core protein and nonstructural protein NS5A are localized to LDs, replication complexes have been demonstrated to be located in LDassociated membranes, and that this recruitment is critical for producing infectious viruses [4][5][6]. Additionally, the expression of core protein in the liver of transgenic mice causes the development of steatosis [7].
Altering the host cell environment by using small molecules and other functional genomic approaches has been demonstrated to modulate HCV replication through diverse mechanisms [8][9][10][11][12][13][14][15]. Depletion of cholesterol using β-cyclodextrin disrupts the membranous web upon which HCV replication occurs [8]. Inhibitors of HMG-CoA reductase, the rate-controlling enzyme of the mevalonate pathway, like lovastatin and fluvastatin, also inhibit HCV replication [9,10]. In addition, small molecule inhibitors such as 25-hydroxycholesterol (25-HC) and cerulenin, both of which act at the mevalonate pathway, also inhibit HCV replication [11,12]. GGTI-286, an inhibitor of geranylgeranylation, inhibits HCV replication [13], and FBL2 has been implicated as a geranylgeranylated host protein required for HCV RNA replication [14]. The PPAR receptor antagonist 2-chloro-5-nitro-N-(pyridyl)benzamide also modulates HCV replication [15]. Collectively, these observations implicate a direct role of cellular lipid metabolism in the HCV life cycle and suggest that further insight into host cell metabolic determinants for the HCV life cycle may help to define new strategies for the development of antiviral therapeutics.
Gene expression profiling has been used to study HCVassociated liver pathologies, including cirrhosis [16], fibrosis [17], and hepatocellular carcinoma (HCC) [18]. Other researchers have used microarrays to elucidate the mechanism(s) underlying IFN resistance [19] and the molecular basis for treatment failure of pegylated IFNα plus ribavirin [20][21][22][23]. The progression of HCV infection has also been studied by transcriptional profiling of liver biopsies in acutely infected chimpanzees that developed persistent infection, transient viral clearance, or sustained clearance, and common as well as outcome-specific changes in gene expression were observed [11]. Transient and sustained viral clearance was associated with the induction of IFNγ-induced genes, antigen processing and presentation genes, and genes associated with the adaptive immune response [11]. Host genes involved in lipid metabolism were also differentially regulated during the early stages of infection and these genes correlated with the outcome of infection [11]. These observations are particularly important given that HCV-derived steatosis is correlated with more rapid disease progression and reduced response to antiviral therapy clinically [24]. Herein we have performed gene expression profiling of Huh-7 cells harbouring HCV replicons in the absence and presence of the mevalonate pathway inhibitor 25-HC ( Fig. 1) and examined which differentially regulated genes lead to an antiviral state in the replicon model. The changes in gene expression in the replicon model were compared with the gene expression profiles in HCV infected chimpanzees to evaluate their role in HCV pathogenesis.
Comparisons were conducted between each set of duplicate data to determine which genes were differentially expressed in 25-HC-treated Huh-7 cells and in the 25-HCtreated HCV replicon cells. The ratios of average expression for 25-HC treated Huh-7 cells/naïve Huh-7 cells were used to compute the fold change values for the 25-HC treated column Table 1. Similarly the 25-HC-treated HCV replicon column in Table 1 was computed from the ratio of average expression values from (replicon + 25-HC)/replicon samples. Transcriptome profiling of Huh-7 cells treated with 25-HC demonstrated that 69 genes were differentially regulated. Forty-seven genes were downregulated, 16 of which are clearly related to the mevalonate pathway (Table 1). A number of membrane proteins and oligonucleotide-binding proteins also appear to be downregulated (Table 1). By contrast, only 22 upregulated genes were observed upon treatment with 25-HC ( Table  2). None of these genes are directly involved in the mevalonate pathway and there appears to be no functional relationship between the entries within Table 2. These observations suggest that 25-HC acts at the transcriptional level as a negative regulator of fatty acid and cholesterol biosynthesis [8,11,23], and that other changes in gene expression observed occur as secondary effects of 25-HC treatment. It is interesting to note that Huh-7 cells treated with both HCV replicon RNA and 25-HC did not show identical transcriptional profiles to cells treated analogously without HCV replicon RNA or to those only treated with HCV replicon RNA or replication-defective replicon RNA (pFK-I 389 neo/NS3-3'/Δ5B replicon RNA). Specifically, 5 genes were uniquely upregulated and 22 genes were uniquely downregulated in the 25-HC-treated HCV replicon harbouring cells (Table 3). Virtually all the genes in Table 3 are hypothetical transcripts or genes of unknown function. Collectively the global gene expression changes represented by Tables 1, 2, 3, represent the transcriptional program that leads to an antiviral state within the host cell hepatocyte, at least for the subgenomic HCV replicon. One or more of these genes may serve as prognostic markers for HCV treatment. Therefore, we performed further studies to investigate the importance of the observed gene expression changes in Tables 1,   2, 3 in other models for HCV pathogenesis.
We examined the transcriptional profiles of the genes from this study in three previously studied chimpanzees that were infected with HCV resulting in different outcomes of infection [28] and for whom genome-wide transcriptional analyses were performed during the course of infection [11]. Briefly, one chimpanzee successfully cleared the virus (sustained clearance, SC), one became persistently infected but initially transiently cleared the virus (transient clearance, TC), and the third chimpanzee developed an unrestrained persistent HCV infection (persistence, PS) in the absence of both an intrahepatic antiviral immune response and IFN-γ induction [28]. Genomewide transcriptional profiling did suggest that the induction of both immunological and hepatocellular genes can influence the course and outcome of HCV infection [11]. Both the SC and TC chimpanzees displayed an altered lipid metabolism profile during the initial rise in viremia following inoculation with HCV as revealed by gene expression profiling, whereas the PS chimpanzee did not [11].
To determine which genes from our study were most relevant in the context of HCV infection in the chimpanzee model, we compared genes identified in Tables 1, 2, 3 The mevalonate pathway Figure 1 The mevalonate pathway. Schematic diagram of the mevalonate pathway with the targets of the inhibitor 25-hydroxycholesterol indicated.
with those of the chimpanzee data sets described previously [11]. To this end, we compared the replicon data sets (transient transfection) with the most representative chimpanzee data sets, those containing genes that changed significantly during the initial rise in viremia (described in the methods section), since both involve an initial host response to viral RNA and both involve the modulation of lipid metabolism genes. From the 69 dif-    indicate genes whose fold-change in gene expression did not exceed 1.5-fold and hence were omitted from the analysis. a Fold change was calculated from the average of two independent hybridizations on U133A Affymetrix GeneChip arrays compared with an average expression value from analogous hybridizations from control cells that were not treated with 25-HC. Fold change for 25-HC-treated column was determined from the ratio of expression 25-HC treated Huh-7 cells/naïve Huh-7 cells and for the 25-HC-treated HCV replicon column from the ratio of expression values from (replicon + 25-HC)/replicon. Dashes (---) indicate genes whose fold-change in gene expression did not exceed 1.5fold and hence were omitted from the analysis.

Table 1: GeneChip analysis of downregulated genes after 25-HC treatment of Huh-7 cells and Huh-7 cells bearing HCV subgenomic replicons. a (Continued)
ferentially regulated genes identified here and the 257 genes derived from transcriptome profiling of infected chimpanzees during the initial rise in viremia [11], we identified 54 overlapping genes. These genes were further analyzed based on outcome-specific gene expression profiles, fold-induction in chimpanzees during the initial rise in viremia, function as determined by bioinformatically and using the SymAtlas at http://symatlas.gnf.org/SymAt las, and using the Ingenuity interaction software. The four genes selected for further study were prospero-related homeobox 1 (PROX1), insulin induced gene 1 (INSIG1), natural killer cell transcript 4 (NK4), and diubiquitin (UBD), Table 4.
The PROX1 gene is a transcription factor expressed in the adult human liver that is involved in hepatocyte specification, proliferation, differentiation, and migration [29]. It is also involved in control of bile acid synthesis and cho-lesterol homeostasis in vivo [30]. Since HCV replication has been demonstrated to be linked with cellular lipid metabolism, and the PROX1 gene was found to have interesting expression patterns in previous gene expression studies carried out in chimpanzee [11,30], the PROX1 gene was chosen for further validation in the HCV replicon cells. In the replicon data set, PROX1 was downregulated in HCV replicon cells compared with naïve Huh-7 cells ( Table 4). The PROX1 gene is downregulated in replicon-harbouring cells compared to naïve Huh-7 cells and represents a class of genes that may be negatively regulated during HCV replication. Interestingly, PROX1 is upregulated in TC and SC during the initial rise in viremia, but negatively regulated in PS ( Table 4).
The INSIG-1 is highly expressed in the adult liver and plays a role in sterol homeostasis. INSIG-1 was downregulated in 25-HC-treated HCV replicon cells compared a Fold change were calculated from the average of two independent hybridizations on U133A Affymetrix GeneChip arrays compared with an average expression value from analogous hybridizations from control cells that were not treated with 25-HC. Fold change for 25-HC-treated column was determined from the ratio of expression 25-HC treated Huh-7 cells/naïve Huh-7 cells and for the 25-HC-treated HCV replicon column from the ratio of expression values from (replicon + 25-HC)/replicon. Dashes (---) indicate genes whose fold-change in gene expression did not exceed 1.5-fold and hence were omitted from the analysis.
with HCV replicon-harbouring cells (Table 4). INSIG-1 represents a class of genes that are negatively regulated during HCV replication under metabolic conditions that are induced by the inhibitor of the mevalonate pathway, 25-HC. This transcript is also upregulated in an outcomespecific manner in the chimpanzee who experienced a clearance episode (SC and TC) but not in that which resulted in a persistent (PS) infection ( Table 4).
The NK4 gene was chosen for further validation because its expression seems to be correlated with a clearance state in the chimpanzee and it was found to be upregulated in the 25-HC-treated HCV replicon cells compared with naïve Huh-7 cells treated with 25-HC (Table 4). Unlike PROX1 and INSIG1, NK4 represents a class of genes that are induced upon HCV replication and appears to be modulated by treatment with 25-HC.
Finally, the UBD gene was chosen for further validation because its expression was increased dramatically in the chimpanzees associated with sustained and transient clearance, 43-fold and 7-fold, respectively, and it was upregulated in 25-HC-treated HCV replicon cells compared with naïve Huh-7 cells treated with 25-HC (Table  4). UBD represents a class of genes that are induced during HCV replication when an antiviral state exists within the host cell, as demonstrated by the increase in expression in HCV replicon harbouring cells treated with 25-HC. This gene may serve as a prognostic marker for HCV infection.
Next we examined the overall gene expression patterns for these four genes (Table 4) in the TC, SC, and PS chimpanzees. The expression levels of the four genes in infected chimpanzees were plotted as a function of time alongside the HCV RNA levels measured in the chimpanzees [11] (Fig. 3). The comparative analysis between the chimpanzee and replicon data sets was done based on the first three chimpanzee time points, which corresponded to the initial rise in viremia [11]. However, the expression of the four genes of interest was determined and plotted for the entire measured time course of infection for the chimpanzees in Fig. 4. In general, the gene expression levels of the four genes paralleled the HCV levels in the chimpanzees that displayed sustained clearance (SC) and transient clearance (TC) of HCV. The overall gene expression levels were lowest in the persistently infected (PS) chimpanzee compared to the SC and TC chimpanzees. Also, the overall gene expression levels were higher for all four genes in the TC chimpanzee, which also had the highest viral titres (Fig. 4).
Expression of PROX1 in the HCV-infected chimpanzees early post-infection decreased slightly within the first six days for the SC chimpanzee, increased within the first eight days in the TC chimpanzee, and decreased within the first eight days in the PS chimpanzee (Fig. 4). Throughout the remaining course of infection, PROX1 levels roughly correlated with HCV levels in the SC and TC chimpanzees, while in the PS chimpanzee the PROX1 levels remained lower than the HCV levels. The expression level of INSIG-1 early post-infection increased in the SC   [11] for transcriptional profiles of 4 chimpanzees with different outcomes of viremia. See reference [11] for the details of the experiments for transcriptional profiling. The transcriptional profiling data from reference [11] was simply mined for the relevant data presented in this table; b Multiple INSIG1 Accessions correspond to the three INSIG1 probes; c Initial Rise in Viremia indicates gene expression levels at timepoint 3 after infection over timepoint 1 just prior to infection for each chimpanzee. Dashes (---) indicate genes whose fold-change in gene expression did not exceed 1.5-fold and hence were omitted from the analysis. N/A indicates transcripts that were not assessed in the chimpanzee gene expression profiles.
and TC chimpanzees while it varied only slightly within the first eighteen days post-infection in the PS chimpanzee (Fig. 4). Thereafter, the INSIG-1 levels closely paralleled the HCV RNA levels in the SC and TC chimpanzees and rose to meet the HCV level in the PS chimpanzee. The NK4 levels during the initial rise in viremia in the chimpanzees paralleled the rising HCV levels in the SC and TC chimpanzees (Fig. 4). In contrast, in the PS chimpanzee the NK4 levels decreased within the first eight days of infection and then increased along with the rising HCV RNA level. Finally, the UBD pattern of expression basically paralleled the levels of HCV RNA in the SC and TC chimpanzees with some variation in the absolute levels of expression relative to that of the HCV RNA (Fig. 4). In the PS chimpanzee, the UBD expression level remained close to baseline levels in the first eighteen days post-infection, and then the UBD expression steadily increased during the remaining time course of HCV infection.
To ascertain the importance of the four genes for HCV replication, siRNA knockdown and overexpression studies were carried out in Huh-7 cells harbouring HCV replicon RNA. RT-PCR and northern blot analyses were used to confirm siRNA and overexpression vector function in each case (data not shown). Huh-7 cells stably expressing neo/ luc/NS3-3'/5.1 HCV replicons [8,[25][26][27] were treated with siRNAs against the four genes in the absence or presence of 25-HC (Fig. 5A). As controls, the HCV replicon cells were also treated with GL2 and GL3 siRNAs, a positive control that targets the luciferase gene expressed from the HCV replicon RNA and a negative control, respectively, in addition to IFNγ treatment [15]. Interestingly, although the PROX1 siRNA did not affect HCV replication on its own, it enhanced the inhibitory effect of treatment with 5 μM 25-HC when compared to mock-transfected replicon cells treated with 5 μM 25-HC (p < 0.01) as indicated (Fig.  5A). Overexpression of PROX1 in replicon cells did not affect HCV replication (Fig. 5C). Knockdown of INSIG-1 with homologous siRNAs did not affect HCV replication in replicon cells either in the absence or presence of 25-HC (Fig. 5A). Similarly, overexpression of INSIG-1 in replicon cells did not affect HCV replication (Fig. 5C).
Although the transcriptional profiling showed that NK4 was upregulated in replicon cells treated with 25-HC, siRNA knockdown of NK4 did not affect HCV replication (Fig. 5A). However, overexpression of NK4 affected HCV replication in the presence of 25-HC treatment, causing a significant decrease in HCV replication compared to replicon cells that were transfected with the backbone β-gal control vector and treated with 25-HC (p < 0.05; Fig. 5C). From the profiling data, UBD was upregulated in 25-HCtreated replicon cells compared with naïve Huh-7 cells treated with 25-HC. However, neither knockdown of UBD using siRNA duplexes (Fig. 5A) nor overexpression of UBD (Fig. 5C) had an affect on HCV replication.
To determine if siRNA treatment or overexpression had any affect on the expression of the HCV nonstructural proteins NS3 and NS5A, Huh-7 cells stably expressing neo/ NS3-3'/5.1 replicons were transiently transfected with Treatment siRNA or overexpression constructs and western blot analyses were conducted. As expected, siRNA knockdown of PROX1 gene expression resulted in complete loss of PROX1 protein expression (Fig. 5B). However, knockdown of PROX1 expression did not have any affect on NS3 or NS5A protein levels (Fig. 5B). Corresponding to their negative effect on HCV replication, IFNγ and 25-HC treatment of replicon cells caused a decrease in the levels of NS3 and NS5A, with the effect being more pronounced for NS5A as we had previously reported [25]. Interestingly, overexpression of INSIG-1 caused an increase in the levels of both HCV non-structural proteins (Fig. 5C).

Inhibition of HCV RNA replication by 25-hydroxycholesterol in a dose-dependent manner
Overexpression of NK4 and UBD had no effect on the levels of HCV non-structural proteins (Fig. 5D).

Discussion
Several examples now exist whereby HCV replication can be modulated or even cured by the targeting of host cell pathways rather than through the direct targeting of the virus [8,11,[13][14][15][25][26][27]. Any perturbation of a host cell pathway that is harnessed or diverted by HCV during its life cycle has the potential to be an antiviral target, particularly if it is not an essential pathway to cell function. Small molecule probes have been used to block HCV replication by preventing the replication complexes from Gene expression patterns for PROX1, INSIG-1, NK4, and UBD in HCV-infected chimpanzees  forming [8,11,13,14], by the selective inhibition of host cell kinases [26], by targeting nuclear hormone receptors [15], and by blocking protein folding pathways [27]. 25-HC has been established to yield an antiviral state within Huh-7 cells by modulating the mevalonate pathway [8,11,13,14]. To better understand the nature of the changes induced to the host cell, we performed global transcriptional profiling of Huh-7 cells and Huh-7 cells bearing HCV replicons. We identified unique sets of genes that are both up-and downregulated by 25-HC and established the transcriptional changes that lead to the inhibition of HCV replication in Huh-7 cells. By comparing these transcriptional profiles (Tables 1, 2, 3) with the gene expression profiles of HCV-infected chimpanzees we identified a number of transcripts that may be important for HCV pathogenesis. These genes may also be prognostic markers for disease progression and outcome that are indicative of the antiviral state of the host cell. From the 54 overlapping genes, we identified four genes for further study (Table 4), in an effort to establish a marker for HCV outcome that was linked to the transcriptional state of the host cell.
PROX1 is a highly conserved transcription factor that is expressed in hepatocytes and is involved in proliferation and migration [29]. It is also involved in the control of bile acid synthesis and cholesterol homeostasis in vivo [30]. The down-regulation of PROX1 in replicon cells compared to naïve Huh-7 cells suggests that PROX1 may negatively influence HCV replication but that this effect may only occur under certain metabolic conditions. The down-regulation of PROX1 was observed in the HCV replicon cells treated with 25-HC. Contrary to this, PROX1 expression was observed to rise early post-infection in the SC and TC chimpanzees (Fig. 4). The decreased expression of PROX1 in the chimpanzee with a persistent HCV infection suggests that PROX1 may be a predictor of the outcome of viremia. Dysregulation of PROX1 was also independently observed in another HCV-infected chimpanzee [31]. However, the fold changes observed for PROX1 may in fact be too small to act as a reliable prognostic marker for clinical applications.
INSIG-1 is highly expressed in the liver and plays an essential role in sterol homeostasis. INSIG-1 inhibits the sterol-dependent retention of the SCAP/SREBP complex in the ER and accelerates the sterol-stimulated degradation of HMG-CoA reductase, thus regulating the expression of cholesterol, fatty acids, and triglycerides [32,33].
Early post-infection, the expression level of INSIG-1 increased in the SC and TC chimpanzees and slightly decreased in the PS chimpanzee (Fig. 4), suggesting that this gene acts as an early metabolic predictor of outcome. Like 25-HC, INSIG-1 is an HMG-CoA reductase inhibitor and, therefore, it would be expected to decrease HCV replication. Since we did not observe any change in HCV rep-lication with the INSIG-1 siRNA treatment, we can only conclude that the kinetics of siRNA knockdown of INSIG-1 were too slow to give rise to significant changes in HMG-CoA reductase during the course of the experiment, given that the INSIG-1 siRNA sequences had been previously validated [23]. Interestingly, overexpression of INSIG-1 caused an increase in the levels of HCV proteins NS3 and NS5A, which may be the result of a homeostatic response (Fig. 5D). The fold change values observed for INSIG1 may also limit its ability to serve as a prognostic marker.
NK4 was discovered by differential screening of a cDNA library derived from human NK cells that were activated by interleukin-2 [34]. Our observation that NK4 levels were increased early in viremia in the SC and TC chimpanzees (Fig. 4) is consistent with a previous report of increased NK4 expression in an HCV-infected chimpanzee [31]. NK4 was also upregulated in replicon cells treated with 25-HC ( Table 2), suggesting that it may contribute to the suppression of replication caused by 25-HC. Taken together, the data suggests that NK4 may be a cofactor in inhibiting HCV replication, but it is not sufficient on its own to clear the virus in a persistent infection.
Diubiquitin (UBD, FAT10) is an alternative and ubiquitinindependent targeting mechanism for degradation by the proteasome [35]. UBD expression is inducible with IFNγ and TNF-α [36], both of which are well known antiviral cytokines. Additionally, previous gene expression profiling in HCV-infected chimpanzees identified differential regulation of proteasome genes [11,31,37]. From the profiling data derived from the HCV replicon cells, UBD was upregulated in 25-HC-treated replicon cells compared with naïve Huh-7 cells treated with 25-HC. This suggests that upregulation of proteasome degradation of HCV proteins may be an ancillary effect of 25-HC treatment. Neither knockdown of UBD using siRNA complexes nor overexpression of UBD had an affect on HCV replication in Huh-7 cells (Fig. 5), indicating that UBD is not a target for the development of potential therapeutic strategies for HCV infection. Rather, the very large fold change values observed, particularly in the large animal model for infection, suggests that it is a good candidate prognostic marker for an antiviral state within the host cell.

Conclusion
Herein we have established the transcriptional changes that give rise to an antiviral state that inhibits HCV replication within Huh-7 cells after treatment with 25-HC. We identified four genes that can act as markers for HCV pathogenesis. The functional validation studies showed that none of these genes are likely to be useful as antiviral drug targets, they were identified as being linked to an antiviral state in the HCV replicon cells and in chimpanzee infections. Not surprisingly, INSIG1 and PROX1 transcriptional levels are related to regulation of the mevalonate pathway. The fact that these genes are differentially regulated in different models of HCV pathogenesis further supports their relevance. In the case of PROX1 and NK4, effects on HCV replication were observed in cell culture but only in certain metabolic states induced by the presence of 25-HC. The overall expression levels of all four genes were the lowest in the persistently infected chimpanzee which suggests that these genes may be predictors of the outcome of infection and might ultimately indicate a particular course of treatment on an individual patient basis. UBD appears to be the most promising since it was observed to have the largest changes in transcriptional levels in the chimpanzee model.

Cell culture
Cell monolayers of the human hepatoma cell line Huh-7 were grown in Dulbecco's modified minimal essential medium (DMEM) (Invitrogen, Burlington, Ontario) supplemented 10% with fetal bovine serum (FBS) (Cansera, Rexdale, Ontario) and with 100 nM non-essential amino acids (NEAA), 50 U/ml penicillin, and 50 μg/ml streptomycin. The medium was further supplemented with G418 to a concentration of 250 μg/ml for Huh-7 cells that stably expressed HCV replicons.

In vitro transcription and transient transfection of replicon RNAs for total RNA extraction
In vitro transcripts of HCV replicon RNAs were generated using the MEGAscript™ kit (Ambion, Austin, Texas) according to the manufacturer's protocol and as described previously [8]. Huh-7 cells were seeded at 5 × 10 5 cells in 60 mm dishes to obtain cells at 70 to 80% confluency the next day. The second day, the cells were washed twice with 1× PBS, pH 7.4 (Gibco/Invitrogen, Burlington, Ontario), and once with serum-and antibiotic-free DMEM-NEAA. Transfection complexes were prepared by mixing 5 μg of in vitro transcribed replicon RNA (either luc/NS3-3'/5.1 or neo/NS3-3'/Δ5B) and 30 μl of DMRIE-C (Gibco/Invitrogen) in 1.5 ml of serum-and antibiotic-free DMEM-NEAA for each dish. The lipid-RNA complexes were immediately added to the washed cells and incubated for 4 h before commencing treatments. Mock-transfected cells received only the DMRIE-C reagent without any RNA.

25-HC and IFN treatments
Four hours post-transfection, the medium containing the lipid:RNA complexes was removed and replaced by DMEM-NEAA supplemented 10% with FBS ± 2 μM 25-HC. The cells were treated for 24 h before the lysates were prepared. As a positive control, some cells were treated with IFNγ (PBL Biomedical Laboratories, Piscataway, NJ) at a concentration of 100 U/ml. Mock-transfected cells were exposed to the treatment solvent, typically ethanol.

Total RNA preparation for microarray experiments
Total RNA was prepared using the RNeasy ® Mini kit (QIA-GEN, Mississauga, Ontario) and following the manufacturer's protocol. The RNA was quantified using the NanoDrop ® ND-1000 Spectrophotometer, and the quality of the RNA was verified using the RNA 6000 PicoLabChip kit (Agilent Technologies, Mississauga, Ontario).

Gene expression analyses
Cellular cRNA samples were prepared and hybridized to high-density oligonucleotide arrays (Affymetrix U133A Human GeneChips) in the same manner as was done for the chimpanzeeanzee samples in our previous study [11]. Each sample was hybridized in duplicate. The neo/NS3-3'/ Δ5B profiles were used as a negative control for HCV replication. Primary analyses were performed as described previously [11]. Briefly, the 116 genes identified as being differentially regulated in the replicon data set represent the logarithm-transformed genes that were deemed "present" and whose average difference exceeded 150 in at least one of the samples described above. Comparative analysis between the chimpanzee and replicon data sets was done based on the first three chimpanzee time points, which corresponded to the initial rise in viremia [11]. The post-infection time points corresponded to weeks -1, 6, and 9 for the SC chimpanzee, weeks -1, 5, and 9 for the TC chimpanzee, and weeks 0, 8, and 14 for the PS chimpanzee. The microarray data was submitted to GEO http://ncbi.nih.gov/geo with accession numbers GSM348423-GSM348433.