Chromatin Dynamics of Gene Activation and Repression in Response to Interferon α (IFNα) Reveal New Roles for Phosphorylated and Unphosphorylated Forms of the Transcription Factor STAT2*

Signal transducer and activator of transcription 2 (STAT2), the critical component of type I interferons signaling, is a prototype latent cytoplasmic signal-dependent transcription factor. Activated tyrosine-phosphorylated STAT2 associates with STAT1 and IRF9 to bind the ISRE elements in the promoters of a subset of IFN-inducible genes (ISGs). In addition to activate hundreds of ISGs, IFNα also represses numerous target genes but the mechanistic basis for this dual effect and transcriptional repression is largely unknown. We investigated by ChIP-chip the binding dynamics of STAT2 and “active” phospho(P)-STAT2 on 113 putative IFNα direct target promoters before and after IFNα induction in Huh7 cells and primary human hepatocytes (PHH). STAT2 is already bound to 62% of our target promoters, including most “classical” ISGs, before IFNα treatment. 31% of STAT2 basally bound promoters also show P-STAT2 positivity. By correlating in vivo promoter occupancy with gene expression and changes in histone methylation marks we found that: 1) STAT2 plays a role in regulating ISGs expression, independently from its phosphorylation; 2) P-STAT2 is involved in ISGs repression; 3) “activated” ISGs are marked by H3K4me1 and H3K4me3 before IFNα; 4) “repressed” genes are marked by H3K27me3 and histone methylation plays a dominant role in driving IFNα-mediated ISGs repression.

Interferons are pleiotropic cytokines induced upon virus infection and other stimuli to modulate host immune response and are classified as type I interferon ␣ and ␤ (IFN␣ and IFN␤), 3 type II (IFN␥), and the recently discovered type III (IFN) (1)(2)(3)). Interferons exert their function by phosphorylating latent transcription factors belonging to the signal transducer and activator of transcription (STAT) family after a signaling cascade, which begins with the binding of interferon to their membrane receptors and involves Janus kinases (JAKs). Receptor dimerization or oligomerization leads to Jak apposition and the trans-phosphorylation on tyrosine residues that releases their intrinsic catalytic activity. Tyrosine-phosphorylated cytokine-receptor cytoplasmic domains then provides binding sites for the Src homology-2 (SH2) domain of the STAT proteins, which are recruited to the JAKs and phosphorylated on a single tyrosine residue (Tyr-689 in the case of STAT2). The interaction between phosphorylated-SH2 domains on STAT proteins leads to homo-or hetero-dimerization and nuclear translocation (4 -6). STAT2 Tyr-689 and STAT1 Tyr-701 can also be phosphorylated by non-receptor TKs, including SRC and ABL in the absence of ligand-induced receptor signaling (7). STAT dimers directly activate genes containing the IFN␥ activation site (GAS) element, while the association of STATs with the DNA-binding protein interferon regulatory factor (IRF) 9 expands the range of DNA response elements that can be targeted by the JAK-STAT pathway to interferon-stimulate response element (ISRE) and IRF response element (IRE) (8). Un-phosphorylated STAT2 binds IRF9 and constitutively shuttles in and out of the nucleus (9). STAT2 tyrosine phosphorylation promotes the association of STAT2/IRF9 with STAT1 to form the IFN-stimulate gene factor 3 (ISGF3). ISGF3 enters the nucleus, binds to the ISRE elements present in the promoter regions of a subset of IFNinducible genes and triggers transcription (10). In the ISGF3 complex, STAT1 and IRF9 are essential to mediate the DNA binding, whereas STAT2 provides a potent TAD (11). Although all STATs have been found to be activated by different type I IFNs in specific cell types and to participate in the formation of complexes with different intrinsic DNA binding specificities to GAS, ISRE, or IRE sites (11)(12)(13), STAT2 is absolutely required for a fully effective response to IFN␣. Indeed, STAT2 may form stable homodimers in response to IFN␣ that interact with IRF9 and activate transcription of ISRE-containing genes (14) and cells lacking STAT2 do not respond at all to IFN␣ activation of ISGs whose transcription is ISRE-dependent (15). In addition, a hybrid IRF9-STAT2 protein is able to recapitulate type I IFN-stimulated gene expression and antiviral response in the absence of STAT1 (16).
Significantly less progress has been made toward identifying the specific STAT-regulated genes responsible for the different biological effects of IFNs and determining how these target genes are selected and regulated at the chromatin level. The repertoire of ISGs whose expression is modulated by the different IFNs has been greatly expanded by DNA microarrays studies (17)(18)(19). However, very few of these genes were originally validated as ISGF3 direct target genes (20,21). Hartman et al. (20) have detected by ChIP-chip a large number of genomic sites in chromosome 22 to which STAT1 and/or STAT2 bind after IFN␥ and/or IFN␣ stimulation and identified a number of new candidate direct target genes for STAT1 homodimers and STAT1:STAT2 heterodimers. Similarly, STAT1 binding sites have been located on the whole genome using several approaches, including ChIPSeq (22)(23)(24). Expression profiling studies have also shown that, in addition to activate hundreds of genes, IFN␣ is also able to repress numerous target genes (21,(25)(26)(27). The basis for this dual effect, and in particular for transcriptional repression, is still unclear as the composition of STAT containing complexes in genes activated or repressed in response to IFNs is largely unknown. To generate further mechanistic knowledge on ISREs-bound protein complexes, we investigated the binding dynamics of STAT2 and its "active" phospho-STAT2 form to a large panel of putative IFN␣ direct target promoters by ChIP-chip on a custom oligonucleotide promoter array. We found that Stat2 is already bound to 62% of investigated target promoters, including most "classical" ISGs, before IFN␣ treatment. Only a proportion of STAT2 basally bound promoters also show P-STAT2 positivity before and after IFN␣. Combined ChIP and expression analysis indicates that STAT2 plays a role in regulating ISGs expression, independently from its phosphorylation. By correlating in vivo promoter occupancy and gene expression with RNA Pol II recruitment and changes in histone methylation, we could show that "activated" ISGs are marked by H3K4me and H3K4me3 before IFN␣ whereas IFN␣ "repressed" genes either have STAT2 bound before IFN␣ or recruit both STAT2 and P-STAT2, thus linking STAT2 with ISGs repression. Finally, "repressed" genes are consistently marked by H3K27me3, suggesting a dominant role of histone methylation in driving IFN␣-mediated ISGs repression.

EXPERIMENTAL PROCEDURES
IFN␣/STAT2 Array Design, Probe Generation, and Microarray Hybridization-We aimed at designing an oligonucleotidebased array containing a large repertoire of genes targeted by IFN-␣/␤, but not IFN-␥ (21). Selection criteria and detail about the oligos can be found in supplemental materials. Probe generation was done essentially as previously described in the Whole Genome Amplification Method (28) using 10 ng of the input and the whole ChIP samples. For a detailed description of the amplification and the array hybridization see supplemental materials.
Chromatin Immunoprecipitation (ChIP)-ChIP analysis was performed as described in (29). Briefly, cells (0.5/1 ϫ 10 8 ) were washed in phosphate-buffered saline (PBS) and incubated for 10 min with 1% formaldehyde; after quenching the reaction with glycine 0.125 M, cells were sonicated, and chromatin fragments of an average length of 1 kb recovered by centrifugation. Immunoprecipitations were performed with ProtG-Sepharose (KPL, 2235101) and 3-5 g of the indicated antibodies. After immunoprecipitation, washes, and reverse cross-linking, the samples were extracted twice with phenol/chloroform, once with chloroform and ethanol precipitated in the presence of 30 g of glycogen. In the Re-ChIp protocol, after the first immunoprecipitation, the samples were eluted in 10 mM DTT at 37°C for 30Ј, the supernatant was collected, diluted 1:20 in Re-ChIP buffer (1% Triton, 2 mM EDTA, 150 mM NaCl, 20 mM Tris-HCl, pH 8) and incubated with the second antibody. Washes, reverse crosslinking, and extraction were the same as previously described. Quantitative PCR was performed using the SYBR-Green ROX mix, and an ABI Prism 7900 Sequence Detection System according to manufacturer's instructions (Applied Biosystems). Primers used are listed in supplemental Table S1. Dissociation curves after amplification showed that all primer pairs generated single products. The amount of PCR product amplified was calculated as % of a standard curve of the input. A genomic region between the GAPDH gene and the chromosome condensation-related SMC-associated protein (CNAP1) gene was used as negative internal control (CTL) (ChIP-IT kit -Active Motif) and the % of input values on target sites were divided for % of input values on the CTL region to obtain the specific enrichment.
Data Analysis-Data analysis was conducted essentially as described in Ref. 30. For each array, background-subtracted signal intensities for ChIP Cy5 and the Input Cy3 channels were obtained. The STAT2 (or the STAT2ϩP-STAT2 in the re-ChIP experiments) ChIP-derived Cy5 values were analyzed both by considering the enrichment versus the Input Cy3 and the NoAb (or the Stat2ϩNoAb in the re-ChIP experiments) ChIP Cy5 channels. In the first analysis, as Cy5-ChIPed chromatin is hybridized together with Cy3-Input chromatin on the same array, values corrected for the corresponding background signals of the NoAb controls were used to calculate the ratio between Cy5 and Cy3 channels. A value Ͼ2 was considered significant. Moreover, to compare the Cy5 intensities of the same spot on different arrays (STAT2 versus NoAb, STAT2ϩP-STAT2 versus STAT2ϩNoAb), raw background-subtracted data were normalized to mean value of elements present in the second quartile of lowest intensities in the Cy5 channel. As previously reported, this normalization approach equalizes negative array elements without dampening the enrichment of signals in positive array elements (31). After this normalization, the Cy5/Cy3 intensity ratio of each spot in the STAT2 or STAT2ϩP-STAT2 slides were divided by the Cy5/Cy3 intensity ratio of the corresponding spot of the NoAb or STAT2ϩNoAb slide. An enrichment ratio of 1.4 was scored positive. The three replicate experiments for both STAT2 and STAT2ϩP-STAT2 were scaled to one another using the ExpressYourself data processing platform. More than 90% of oligonucleotides with a Cy5/Cy3 ratio considered "positive" (Ͼ2) were confirmed by the enrichment analysis over the corresponding NoAb. Entire list of values for the 339 oligos is available upon request.
Cell Cultures and IFN␣ Treatment-Huh7 hepatoma cells were maintained in DMEM supplemented with 10% fetal bovine serum (Invitrogen), 1% penicillin/streptomycin and 1% glutamine (Sigma). IFN␣ (Roferon-A, Roche, Basel) was used at a final concentration of 1000 IU/ml and added directly to the culture medium.
Liver Tissue Samples and Primary Culture of Human Hepatocytes-Primary human hepatocytes (PHH) were prepared from adult patients undergoing lobectomy or segmental liver resection for medically required purposes unrelated to this research program. The use of these human hepatic specimens for scientific purposes has been approved by the French National Ethic Committee. The PHH used in this study were obtained from two patients: FT304, 52-year-old female transplant donor and FT310, 60-year-old female, undergoing liver resection for a cystadenoma. Both patients were negative for HBV, HCV, and HIV serologic markers. PHH were prepared and cultured as described elsewhere (32). The cells were plated into collagen-coated dishes (BD Biosciences) at 1.7 ϫ 105 cells/ cm 2 in a hormonally and chemically defined medium (32). Forty-eight hours after plating, PHH were exposed to IFN␣ at 1000 units/ml for the indicated times.
RNA Extraction and Taqman Low Density Arrays-Subconfluent Huh7 cells were treated with 1000 IU/ml of IFN␣ and/or 300 nM TSA and total mRNAs were extracted at the indicated time points using the RNAeasy mini kit (Qiagen). Reverse transcription was performed using the random examers method (Superscript II kit, Invitrogen). RNA quality and quantity were monitored by ethidium bromide staining and by UV absorbance. Custom real-time PCR liquid arrays (Taqman Low Density Arrays or TLDAs -Applied Biosystems) were loaded with 200 ng of cDNA per lane and run in a 7900HT Fast Real-Time apparatus (Applied Biosystems), according to manufacturer's instructions. Data were analyzed using the SDS 2.2.2 program. 18 S RNA was used as internal control for normalizing equal loading of samples.

In Vivo Binding of STAT2 to Target Promoters in Huh7 Cells-
Aiming to get insights on IFN␣/STAT2 trascriptome regulation at the chromatin level, we designed a custom oligonucleotide array to probe by ChIP-chip 113 gene promoters known to be modulated by type I IFNs (and not IFN␥) in expression profiling experiments and/or to be specifically bound by STAT2 in vivo in ChIP assays (see supplemental Table S2 and "Experi-mental Procedures" for selection criteria). For each promoter, three 50-mers were designed, the first one located in the upstream regulatory region (Oligo 1; between Ϫ1500 and Ϫ500 from tss), the second in the proximal promoter (Oligo 2; between Ϫ500 and ϩ0) and the third downstream of the tss (oligo 3, between ϩ0 and ϩ1000). A complete description of the ISREs and 50-mers probes positions is given in supplemental Table S3. Cross-linked chromatin was prepared from untreated and IFN␣ treated Huh7 cells and immunoprecipitated with an ␣-STAT2-specific antibody. Amplified and Cy5labeled anti-STAT2 and control NoAb ChIPed-chromatins were hybridized on the array together with the Cy3-labeled Input chromatin. Both the enrichment of STAT2 Cy5 channel over the NoAb Cy5 channel hybridized in parallel and of Cy5-ChIPed chromatin over the corresponding Cy3-Input chromatin were considered for data analysis (see "Experimental Procedures"). As expected, positive spots corresponded to the ISRE-positive oligonucleotides for each promoter set, the flanking negative oligonucleotides serving as internal controls for hybridization specificity (see supplemental Table S4). As shown in supplemental Fig. S1a (S1A) when oligonucleotides are divided into classes according to their enrichment over the Input before (left panel) and after (right panel) IFN␣ treatment, the negative internal control oligos fall into the group with Cy5/ Cy3 ratios between 1 and 2 (therefore scored "negative") and the frequency of negative internal control oligos scored "negative" is not significantly affected by IFN treatment. Notably, the vast majority of oligonucleotides scored as "positive" [i.e. with a Cy5/Cy3 ratio Ͼ2.0 (blue colums)] and display enrichment ratios Ͼ2,5 and only a marginal number of them falls near the cut-off values (2.0 to 2,5).
We found STAT2 already bound to 70 promoters, corresponding to 62% of total, before IFN␣ stimulation, and to 90 promoters, corresponding to 80% of total, after treatment ( Fig.  1A). In particular, 60 (53%) promoters display STAT2 binding before and after IFN␣ stimulation, 10 (9%) lose a pre-existing binding after IFN␣ stimulation, while 33 (29%) acquire STAT2 binding only after IFN␣ treatment. 10 promoters (9%) show no STAT2 binding in both conditions (Fig. 1A). This latter observation may be related to the specific cell line used in the study or might reflect a binding kinetic different from the time points under investigation. The Gene Ontology categories most enriched in basal STAT2 binding are Immune modulation/ Host defense (24 out of 40 genes), Transcription Regulation (10/25), Protein/Cell metabolism (12/13), Ubiquitination/Protein degradation (5/6) cytoskeleton/cell adhesion (3/3). After IFN␣ treatment, while cytoskeleton/cell adhesion and protein/ cell metabolism genes tend to lose STAT2 binding, the immune modulation/host defense and transcription regulation categories dramatically increase the association with STAT2 (Fig. 1B). These data indicate that STAT2 binds in vivo to a number of target promoters, among which the large majority of "classical" ISGs before IFN␣ treatment and that the cytokine stimulation triggers a STAT2 re-distribution in addition to its well characterized recruitment on responsive genes. The notion of STAT2 occupancy on target promoters in unstimulated cells is indirectly supported by the observation that STAT2 shuttles in and out the nucleus without being phosphorylated (9), although the biological significance and the consequences on gene transcription were not established.
Phosphorylated STAT2 Binding to Target Promoters-Next, we investigated the binding of the putative "active" phosphorylated STAT2 (P-STAT2) to the target promoters before and after IFN␣ treatment. To this aim, we performed ChIP-chip experiments using chromatin sequentially chipped with ␣-STAT2 and then with ␣-P-STAT2 antibodies (see "Experimental Procedures"). Preliminary experiments showed that this approach permits to minimize false positive signals. supplemental Fig. S1B shows that, likewise for STAT2, also for P-STAT2 most oligonucleotides with a Cy5/Cy3 ratio evaluated as "positive" in our analysis (blue colums) display higher enrichment ratios (Ͼ2.5). The frequency of the negative internal control oligonucleotides is comparable to that found after STAT2 hybridization and is, again, not affected by IFN␣ stimulation.
The overall view of the results shows that the majority of the P-STAT2 positive promoters (49; 43%) acquire P-STAT2 binding after IFN␣ treatment, while only a small subset of targets loses (5; 5%) or displays P-STAT2 occupancy before IFN␣ stimulation (17; 15%), respectively ( Fig. 2A). 42 out of 113 promoters (37%) are not occupied by P-STAT2 at the investigated time points. All the promoters that are P-STAT2 positive also show STAT2 positivity. Conversely, all the promoters that are STAT2 negative do not show any P-STAT2 recruitment, thus confirming the reliability of the sequential ChIP protocol. It is worth to note that most well-known target genes belonging to the Immune modulation/Host defense Gene Ontology category acquire, as expected, P-STAT2 binding after IFN␣ stimulation. A similar behavior is shared by many target genes belonging to ubiquitination/protein degradation and transcription regulation GO categories (Fig. 2B).
Next, we crossed the STAT2 and P-STAT2 occupancy data sets obtained in all the six ChIP-chip experiments. The analysis of STAT2 and P-STAT2 dynamic recruitment on chromatin before and after IFN␣ revealed a rather complex scenario. On the whole, 22 of the 72 (30%) target promoters that are STAT2 positive before IFN␣ treatment are also P-STAT2 positive (i.e. STAT2 and P-STAT2 pos-pos ϩ pos-neg). This percentage raises to 74% after IFN treatment (i.e. STAT2 and P-STAT2 pos-pos ϩ neg-pos) ( Table 1). Among the genes that recruit STAT2 after IFN treatment (i.e. STAT2 neg-pos), about twothirds of them (23/33) display a concomitant recruitment of P-STAT2 (i.e. P-STAT2 neg-pos). Similarly, between the gene targets that are STAT2 positive before and after IFN stimulation (i.e. STAT2 pos-pos), STAT2 occupation is accompanied by a P-STAT2 binding in 19 out of 60 (i.e. P-STAT2 pos-pos ϩ pos-neg) before treatment, and in 46 out of 60 (i.e. P-STAT2 pos-pos ϩ neg-pos) after treatment. Finally, the 12 promoters that possess STAT2 binding before IFN stimulation and lose it after treatment show a concomitant P-STAT2 occupancy only FIGURE 1. ISRE occupancy by STAT2 before and after IFN␣ treatment. A, STAT2-ChIPped DNAs from untreated and IFN␣-treated Huh7 cells were labeled and hybridized to the STAT2/IFN␣ array as described under "Experimental Procedures." Genes were divided into four categories depending on whether they: 1) already display STAT2 binding to their ISRE site before IFN␣ treatment and maintain it after the treatment (pos-pos; 53%) 2) acquire STAT2 binding after treatment (neg-pos; 29%) 3) STAT2 is bound prior to IFN␣ treatment and detaches after treatment (pos-neg; 9%) 4) display no STAT2 binding before and after IFN␣ treatment (neg-neg; 9%). B, dynamics of ISRE/STAT2 occupancy according to Gene Ontology categories (within brackets the number of genes belonging to each category). The figures represent the percentage of genes for each GO category displaying ISRE STAT2 binding before (pos-pos and pos-neg groups, as defined in A) and after IFN␣ stimulation (pos-pos or neg-pos groups), respectively. The schemes representing the ISGF3 include all the known components of the complex (STAT1, STAT2, and IRF9). However, because our ChIP-chip experiments investigated only STAT2 direct binding to the ISREs, STAT1, and IRF9 are depicted as shaded with dotted line (STAT1) and as only dotted line (IRF9) (see the Introduction section for a detailed discussion of IRF9 requirement for STAT2 activity).
in three cases. From this data, it is evident that STAT2 and P-STAT2 binding is not always concomitant. In particular, P-STAT2 can be recruited to promoters that already show a pre-existing STAT2 binding (i.e. STAT2 pos-pos and P-STAT2 neg-pos); P-STAT2 could not be found on a number of STAT2 positive targets (i.e. P-STAT2 neg-neg) either before or after IFN, and, finally, it may be lost from regions that continue to display a STAT2 binding (i.e. STAT2 pos-pos and P-STAT2 pos-neg). Although some of these observations might be due to kinetic changes of P-STAT2 recruitment that we unable to reveal at the investigated time points, with mobile complexes that continuously bind and detach from DNA, they could also reflect a role of un-phosphorylated STAT2 in transcription regulation.
ChIP-chip Results Validation-We selected 18 gene promoters representative of the different combinations of STAT2/P-STAT2 occupation status before and after IFN␣ to independently test the array data generated in Huh7 cells by performing quantitative PCR reactions (Fig. 3A). USP18, CXCL10 and the zinc finger transcription factor ZBED4 are positive for both STAT2 and P-STAT2 before and after IFN␣ treatment. USP18 was originally identified as a protease cleaving ISG15 from its conjugated proteins, but recent works attribute to USP18 also a role in the long-term desensitization of IFN signal transduction pathway in the mouse liver (33,34). CXCL10 is a chemokine induced by both IFN␣ and IFN␥ responses. Among the gene promoters that are STAT2 positive before and after IFN treatment, but recruit P-STAT2 only after the stimulus, we chose the typical ISGs MX1 (35), IFI6 (36), and RSAD2 (VIPERIN; (37)), the ubiquitin ligase RBX1 (38) and "Suppressor of tumorigenity 13" (ST13; (39)), encoding for the HIP1␣ protein that interacts with and potentiates the chaperone functions of HSP70 in protein folding and repair. Whereas MX1, IFI6, and RSAD2 increase STAT2 enrichment at their promoters after IFN treatment, ST13 and RBX1 display a global reduction of STAT2 binding after treatment that is accompanied by the recruitment of P-STAT2. The tryptophanyl-tRNA synthetase  Fig. 1. Moreover, as a fraction of the promoters that are negative for P-STAT2 binding could be occupied by un-phosphorylated STAT2 in this categories the whole ISGF3 complex is drawn with dotted lines.  WARS (40) and the histone deacetylase inhibitor GHDC (LGP1; (41)) promoters are bound by STAT2 but do not display any P-STAT2 recruitment before or after IFN treatment. RFC2, HIF1A, and the interleukin-1␤-converting enzyme (ICE) family member LAP3 show STAT2 and P-STAT2 recruitment only after IFN stimulation. RFC2 encodes the third-largest subunit of the replication factor C complex, involved in clamp loading and DNA replication (42). HIF1a, in addition to its well known role in the cell response to hypoxia, has been recently involved in innate immune responses: IFN␣ induces HIF1A in endothelial cells under normoxic conditions (43) and the administration of double-stranded nucleic acids results in an increase of HIF1␣ protein levels and HIF1␣ target genes expression (44). HERC5 (also referred to as CEB1), a E3 ligase that mediates type I IFN induced ISGylation of protein targets (45), shows STAT2 recruitment after IFN treatment and no binding of P-STAT2 at any investigated time points. The splicing factor subunit SF3A1, the caveolae structural protein CAV1 (46) and the transcription factor involved in ␤-catenin/Wnt signaling PITX2 (47) display STAT2 binding before IFN stimulation and no STAT2 and P-STAT2 binding after the treatment. Finally, we performed additional ChIP experiments to confirm STAT2 and P-STAT2 occupancy on the same set of promoters in primary human hepatocytes (PHH) (Fig. 3B). Altogether the validation step fully confirmed the ChIP-chip results, not only in the hepatocellular carcinoma cell line Huh7, but also in the context of non transformed primary cells.
Correlation between STAT2/P-STAT2 Promoter Occupation and IFN␣-dependent Gene Expression Regulation-Next we investigated how STAT2/P-STAT2 promoter occupation and changes in ISRE-bound STAT2 complexes translated into gene expression in untreated and IFN␣-treated Huh7 cells and PHH. More in detail we wanted to assess: 1) the function of STAT2 and P-STAT2 bound to promoters under basal conditions 2) whether un-phosphorylated STAT2 could play a role in directing ISGs transcription independent from P-STAT2, and 3) whether transcriptionally activated and repressed genes cluster in different STAT2/P-STAT2 binding categories. To this aim,  Table S5. Standard deviation is always lower than 5% B, Venn representation of gene expression overlaps between Huh7 cells and PHH.
we correlated the expression profile of 76 genes belonging to the ChIP-chip array with the STAT2/P-STAT2 occupancy at their promoters. Total RNAs were prepared from Huh7 cells and PHH left untreated or treated with 1000 IU/ml of IFN␣ for 2, 4, and 8 h and the cDNA loaded onto custom TLDAs (see "Experimental Procedures" for details). The results of the real time transcripts quantification in both cell lines are shown in Fig. 4 and in supplemental Fig. S6. Six genes that are not expressed in Huh7 cells, but are modulated in PHH are shown separately in panel C. Four of them are genes involved in metabolism of cells belonging to the immunomodulatory system (BST1, CSF2RB, CYBB, and NCF1), while RASD2 is a Ras protein homolog probably involved in schizophrenia and SLC5A1 is a sodium/glucose co-transporter (see supplemental Table S2 for gene annotations). Genes displaying a fold enrichment over control (untreated cells) Ͼ2 were considered as "activated" and those which had more than 40% reduction respect to the control were labeled as "repressed." Supplemental Fig. S1C (S1C) shows the kinetics of IFN␣ effect on gene expression in Huh7 cells. Interestingly, activated genes (red columns) show a fold induction several folds higher than cut-off (Ͼ2), and their number reaches the maximum after 4 h of stimulation and remains the same also after 8 h post-IFN␣. Repressed genes show a different kinetics, in that their highest frequency is found 2 h after treatment and then diminishes progressively.
42 genes (55,3% of total) in Huh7 and 49 genes (64% of total) in PHH are activated in response to IFN␣, whereas repressed genes represent 30% of the total in Huh7 cells (23 genes) and 20% in PHH (15 genes) (Fig. 4A). The overlap between the two cell lines is about 80% for the activated genes and 70% for the repressed ones (Fig. 4B). In Fig. 4A, to allow a better visual comparison of the results in Huh7 cells (upper panel) and PHH (lower panel) genes were arranged according to their Gene Ontology category. Notably, immune modulation genes show little differences among the two cell lines, the more relevant being represented by IL8, which is, as expected, activated by IFN␣ in PHH but has an opposite behavior in Huh7 cells and MX2, a "typical" ISG that is unexpectedly only marginally activated by IFN␣ in Huh7. Transcription, cell metabolism, ubiquitination, and cell signaling categories show a very conserved expression profile in Huh7 and PHH. The few exceptions are the tryptophanyl tRNA synthetase WARS, the transporters SLC25A1 and SLC2A11, the E3 ligase RBX1 and the transcription factor HIF1a. Interestingly, the categories that display more differences in gene expression between Huh7 and PHH are apoptosis/cell proliferation and chromatin remodeling. The 3 chromatin modifiers HDAC10, GHDC, and SMARCB1 are all activated in PHH and repressed in Huh7 cells. Interestingly, HDAC10 CϾT polymorphism, associated to an augmented HDAC10 expression, has been linked to an accelerated HCC development (48). GHDC functions as an HDAC inhibitor and it has been shown to sensitize U937 cells to TRAIL-induced apoptosis (41). SMARCB1 encodes for the key subunit of the SWI/SNF-like complex BAF, which has been involved in the rapid induction of ISGs by IFN␣ (49). The apoptosis/growth suppression promoting genes BAX, DAD1, and SGSM3 are all repressed in Huh7 cells and activated in PHH. Conversely, BAG3, a member of the BAG co-chaperone protein family which sustains cell survival (50), shows an opposite behavior. These differences also account for most of the differences existing between repressed genes in Huh7 with respect to PHH. Because most of the genes displaying a differential behavior are involved in cell proliferation and survival, it is likely that these differences might be related to the transformed nature of Huh7 cells.
Next, we crossed the expression profile data with the ChIPchip results. As shown in Table 2, among the 42 activated genes, 25 are STAT2 positive before and after IFN␣ treatment and almost all of them display also a P-STAT2 binding, either before and after IFN␣ (P-STAT2 pos-pos, 10 genes) or only after IFN␣ (P-STAT2 neg-pos, 12 genes). Except a few exceptions (i.e. IFIT1, IFIT3, IFITM1), all the "typical" immunomodulated ISGs present in our custom TLDA enter this category. A smaller number of activated genes show the recruitment of STAT2 alone (STAT2 neg-pos and P-STAT2 neg-neg, 4 genes) or of both STAT2 and P-STAT2 only after IFN␣ stimulation (STAT2 and P-STAT2 neg-pos, 8 genes). On the other hand, half of the repressed genes display a STAT2 binding before IFN treatment (STAT2 pos-pos, 10 genes, and STAT2 pos-neg, 2 genes. Finally, (6) of the 23 repressed genes show the recruitment of both STAT2 and P-STAT2 after IFN stimulation.
Scatter plots in supplemental Fig. S1D (S1D) show how activated (red spots) and repressed (green spots) genes divide according to their level of positivity for STAT2 or P-STAT2 enrichment before and after IFN␣ treatment. It is worth noting that both transcriptional activation and repression correlate with high levels of STAT2/P-STAT2 occupation. In addition, the scatter plot of P-STAT2 Cy5/Cy3 ratio in untreated Huh7 cells clearly shows that, differently from what observed with STAT2, there are no repressed genes having a basal P-STAT2 binding.
Altogether, these data suggest that many typical activated ISGs possess a basal STAT2 binding before IFN␣ stimulation, a Genes modulated in Huh7 cells exposed to 1000 IU/ml IFN␣ for 2, 4, and 8 h were divided into three categories: activated, repressed, and not changed (see "Experimental Procedures" for more details). b All combinations of STAT2/P-STAT2 occupancy at their ISRE sites are listed for activated and repressed ISGs. The frequency of each "ISRE occupancy group" is expressed both as absolute number of genes and as percentage within each category.
that is not always accompanied by a P-STAT2 binding either before or after IFN␣ and that, unexpectedly, IFN␣ repressed genes seem to require the presence of STAT2 before IFN␣ stimulation or the recruitment of both STAT2 and P-STAT2 after the treatment. Chromatin Dynamics at STAT2/P-STAT2 Target Promoters-To investigate the chromatin context associated to IFN␣ stimulation and STAT2/P-STAT2 binding dynamics, we performed ChIP assays followed by real time PCR analysis on the same set of 18 promoters analyzed for STAT2/P-STAT2 occupation (see Fig. 3) with antibodies against, monomethyl-and trimethyl-Lys4 of histone H3 and trimethyl-Lys27 of histone H3 (Fig. 5). IFN␣-activated genes show no association at all to the repres-sive mark H3K27me3. Conversely, the presence of H3K4me1 is restricted to activated promoters. Before IFN␣ stimulation, H3K4me1 can be found only on promoters that display a concomitant STAT2 occupation (RSAD2, USP18, CXCL10, MX1, IFI6 (in the STAT2 pos-pos group) and CAV1 (in the STAT2 pos-neg group)). Conversely, IFN␣-repressed genes display a strong correlation with H3K27me3. In particular, H3K27me3 is associated to P-STAT2 presence/recruitment on target promoters (ZBED4, RBX1 (in the P-STAT2 pos-pos group), ST13, HIF1a and RFC2 (in the P-STAT2 neg-pos group)). SF3A1 and PITX2 represent an exception: both are STAT2 positive before IFN␣ and become STAT2 negative after the treatment; they display H3K27me3 but do not show any P-STAT2 binding either before or after IFN. The H3K4me3 modification is present both in activated and repressed genes. In repressed genes, it seems to mirror H3K27me3 and it is found also on those promoters that do not recruit P-STAT2 (GHDC and WARS).

DISCUSSION
In this work we investigated by ChIP-chip the binding dynamics of STAT2 and its phosphorylated form P-STAT2 to the ISRE sequences present on a large repertoire of IFN␣/␤ specific target genes before and after IFN␣ induction in hepatocellular carcinoma Huh7 cells and in primary human hepatocytes (PHH). The STAT2 and P-STAT2 occupancy and the chromatin context around the STAT2-bound ISREs were correlated with activation and repression of genes expression. We found that the 62% of the investigated promoters show STAT2 binding prior to IFN␣ treatment. The large majority of these genes belongs to the Immune Modulation/Host Defense pathway and include most of the "typical" ISGs. Although the presence of STAT2 on target promoters before its activation by the IFN␣ signaling cascade might be unexpected, it has to be underlined that recent work indicates that p53, known to be also "signal" activated, binds to its target sequences prior to "activation" by DNA damaging agents or other stimuli (30,(51)(52)(53)(54). TCF4, the nuclear effector of Wnt signaling, has been shown to bind to its DNA target sequences on Wnt-responsive genes in a stimulus-independent manner and TCF4 chromatin recruitment seems to correlate with a specific epigenetic structure (55). A note of caution in the interpretation of chromatin occupancy by supposedly "latent" and "signal activated" transcription factors came from the observation that this phenomenon, in the case of p53 (Shaked et al.,54), might be restricted to transformed cells and not observed in non transformed cells. Our results do not support this view, and, indeed, we found a good correlation between STAT2 basal occupancy of target ISREs in Huh7 cells and in primary hepatocytes. Moreover, STAT2 is able to shuttle constitutively in and out of the nucleus without being phosphorylated (9). We found that only 30% of STAT2 basally bound promoters also show a P-STAT2 positivity and the percentage raises to 74% after IFN␣ stimulation. Although we cannot formally exclude that false-negatives results for P-STAT2 occupancy may have occurred for one or few genes, the combined analysis of the STAT2/P-STAT2 binding and target genes expression, strongly support the notion that STAT2 plays a role in regulating ISGs expression in a phoshorylation-independent manner. Similarly, un-phosphory- FIGURE 5. Histone marks at ISRE sites and transcriptional modulation of selected array genes. Chromatin from untreated (K) and IFN␣-treated (1000 UI/ml for 1 h) Huh7 cells was immunoprecipitated with ␣-H3K4me1, ␣-H3K4me3 or an anti-H3K27me3 antibodies and analyzed by real-time PCR with primers amplifying the regions corresponding to the oligonucleotide specifically enriched according to the ChIP-chip arrays data analysis (see supplemental Tables S1 and S4). Results are expressed as specific enrichment (% of input values on selected regions divided for % of input values on the control (CTL) region) (see "Experimental Procedures" for more details). The figure shows mean values derived from three independent ChIPs experiments. Bars indicate standard deviation. lated STAT1 has been recently found to drive the constitutive expression of some target genes and to prolong the expression of specific ISGs without phosphorylation (56). Finally, Lou et al. (57) have shown that STAT2-IRF9 heterodimers can drive the activation of RIG-G, a "typical" ISG, without requiring STAT2 phosphorylation. Thus, different STATs or P-STATs complexes could alternate each other on specific targets and exert both distinct and overlapping functions.
Our data also show that P-STAT2 can be directly involved in repressing some interferon target genes. Indeed, we have recently reported that the recruitment of a P-STAT2 containing ISGF3 complex is required for transcriptional repression of the DNp73 promoter in response to IFN␣ (58). 26% of our repressed genes recruit both STAT2 and P-STAT2 after IFN␣ treatment, thus behaving as the DNp73 promoter. Moreover, 6 out of 12 genes repressed by IFN␣ that display STAT2 binding before interferon recruit P-STAT2 after treatment. These observations might be consistent with a role of basal STAT2 in the maintenance of a pre-existing "chromatin repressive mark," whereas the "de novo" recruitment of a repressive complex would require P-STAT2. Notably, by the analysis of multiple ISRE sequence alignment, we could not find any significant differences in nucleotide composition or conservation between activated and repressed genes that might explain the recruitment of transcriptionally repressive versus active ISGF3 complexes(supplemental Fig. S2). Interestingly, in activated target genes, ISREs are mostly located between 0 and 200 from the tss, while in repressed genes the ISREs are far upstream or downstream.
We also investigated the chromatin context associated with STAT2/P-STAT2 binding by analyzing selected histone marks (i.e. H3K4me1, H3K4me3, and H3K27me3) around the ISREs involved in STAT2-mediated modulation of ISGs. It is well established that all monomethylated lysines are associated with active transcription (59). H3K4me3 is also widely associated to active transcription, but the results are more controversial (59). On the other hand, H3K27me3 appears to be a dominant mark and invariably mediates transcriptional repression (60). We found H3K4me1 only on activated promoters that display a concomitant STAT2 occupation before IFN␣ stimulation. Moreover, H3K4me1 is also found together with H3K4me3, confirming the frequent association of these two modifications at transcription start sites (61). These observations could suggest that basally bound STAT2 and H3K4me1/H3K4me3 may contribute together to "prime" chromatin loci that will be actively transcribed after IFN␣ stimulation. This would be also consistent with the finding that for the large majority of STAT1 binding sites after IFN␥ stimulation, dominant H3K4me1/ H3K4me3 combinations were already established before activation (62). Conversely, repressed genes display a strong correlation with the presence of H3K27me3, that is also associated to P-STAT2 presence/recruitment on target promoters. In repressed genes, H3K4me3 mirrors H3K27me3 pattern. The co-occurrence of the "positive" H3K4me3 and the "negative" H3K27me3 marks on the same locus has been already reported but it is not proven that both modifications occur on the same nucleosome.
Acetylation has been reported to play a dual role in the activation of ISGs and IFN␣ signaling, as it is required for a correct ISGF3 assembly (63, 64) but, after transcriptional complexes are recruited on the chromatin, it assumes a repressive role and HDACs activity is needed (65). When Huh7 cells were treated with both IFN␣ and TSA we could fully reproduce these findings. On the other hand, TSA treatment has no impact on IFN␣-induced repression, 4 indicating that acetylation does not play a pivotal role in IFN-repressive activity and probably histone methylation assumes a dominant role in driving IFN␣mediated ISGs repression.
In conclusion, we provide extensive new mechanistic knowledge on the transcriptional regulation at/around the ISRE regulatory sites of a large subset of ISGs and we define the molecular basis for IFN␣-induced repression. Altogether, our observations further challenge the classical paradigm of STAT2 as a latent transcription factor activated only in response to cytokine stimulation. In this regard, we clearly show that STAT2 may regulate ISGs expression independently from its phosphorylation and that, on the other hand, P-Stat2 is involved in ISGs repression. We also found that histone methylation, rather than histones deacetylation, seems to have a dominant role in driving IFN␣-mediated ISGs repression.
Our results may also translate into new insights for the clinical use of class I interferons. IFN␣ is the backbone drug for chronic hepatitis C patient (CHC) and an important option for chronic hepatitis B patients (CHB) (66,67). Current IFN␣based therapies cure 50% of CHC patients but the molecular mechanisms that differentiate IFN responder and non-responder CHC patients are still unclear. Genome wide expression profiling studies have shown that the up-regulation of a specific set of IFN␣ responsive genes, including a number of "non classical" ISGs, prior to IFN␣ treatment in chronic HCV hepatitis patients predicts a subsequent non response to exogenous IFN␣ therapy (68,69). A number of these "predictive" ISGs (namely USP18, CXCL10, MX1, IFI6, RSAD2, GHDC, LAP3, and HERC5) are included in this study and share some regulatory features and, in particular, the ability to be repressed by HDACs inhibitors. It will be interesting to expand the characterization of these "predictive" genes and to test whether the detailed knowledge of their chromatin dynamics might help in establishing meaningful ways to modulate their expression in CHC patients.