Supporting Information

Canaan et al. 10.1073/pnas.0911676106 SI Text pooled scores from several arrays were averaged for each Tissue Culture. The human B cell lymphoma line, BJAB, and the and accepted if Ͼ2 for up-regulated or Ͻ0.5 for down- ϩ BJAB/EBNA1 cell lines, were kindly provided by the lab of Elliott regulated genes, because the scores of EBNA1 samples were Kieff (Harvard University). EBNA1 expression phenotype was superimposed on their cell line counterpart. verified by Western blot analysis before these studies and was reconfirmed by QRT-PCR while confirming microarray data. Quantitative Real-Time PCR (QRT-PCR). To verify and evaluate the Excluding the expression of EBNA1 gene, the cell lines in this study cDNA microarray data, samples were run using two independent are EBV-negative cell lines. The human epithelial cell line 293 and RNA preparations for BJAB and one RNA preparation for 293. the 293/EBNA1 were originated at the laboratory of Lawrence S. These samples were then reverse transcribed in triplicate, while Young (Institute for Cancer Studies, University of Birmingham each cDNA was analyzed in a single QRT-PCR. QRT-PCR Medical School, U.K.). reactions were run in ABI PRISM 7000 Sequence Detection System. The following oligo sequence pairs were used in the cDNA-Microarray Slides. The Peter MacCallum Cancer Institute QRT-PCR assays, sequences shown from 5Ј to 3Ј ends: Hs.833: (PMCI) microarrays used in this report contain approximately CAGCTCCATGTCGGTGTCAG and CCAATCTTCTGGGT- 10,500 cDNA clones representing 9,389 unique clones (UniGene GATCTGC; Hs.82396: CTATGCTTGGGAGCGAGGG and build 144) (25). CCCTGGGCTGTGTTGAAATG; Hs.62354: TACGATGG- GAAACTCTCTAGTGCC and AGCCCGTGGATTG- Purification of RNA, Amplification, Labeling, and Hybridization. Cell lines TACGTGA; Hs.75256: AATGTGCCAGTATGGCTCCC and were grown in four 750-mL tissue culture flasks (Falcon) to a GACACCCACAAAAGGCCAAG; Hs.77100: TCACA- density of 106 cells per 1 mL medium for BJAB and Ͻ70% GAAAAAGCGACGCTTT and CACTCCAGCCAAGTGT- confluence for the 293 under standard humidifying conditions. TCGTT; Hs.46423: AGACGCCGTCACCTATACGG and AT- RNA was purified in a two-step process involving extraction with an GGCTGTGACAGTTTTGCG; Hs.75450: AGCGTGG- organic solvent (TRIzol, Life Technologies) followed by column TGGCCATAGACA and CTTCACCAGATCCATGGCCT; chromatography (RNeasy, Qiagen). Microarray analysis was per- Hs.75573: GCCTTTGTCCAGAGGTGCAA and AGGACCTG- formed using two-color competitive hybridization of fluorescently GCTGAGAATCCAC; Hs.75703: CTCTCAGCACCAAT- labeled cDNA to the glass slide array, comparing EBNA1-positive GGGCTC and GTAAGAAAAGCAGCAGGCGG; Two pairs of cell line to its EBNA1-negative counterpart. RNA samples were oligos were used for Hs.79070: CCGTCCAAGCAGAG- converted to cDNA, amplified, labeled, and hybridized as was GAGC and TCCGCAACAAGTCCTCTTCAG; CGAC- described in ref. 25. GAGAACAGTTGAAACACAA and CGCACAAGAGTTC- CGTAGCTG; EBNA1: CGACATCACACCATATACCGCA cDNA Microarray Scanning and Analysis. Slides were scanned using and ATTTCCAGGTCCTGTACCTGGC; the last pair was used to a dual UV-laser GSI Luminomics scanner (Hewlett Packard), confirm the presence or absence of EBNA1 in samples subjected to and fluorescence intensity was analyzed using the Quantarray QRT-PCR analysis. To overcome differences in initial mRNA program (Hewlett Packard). Images were formed by superim- amounts and to better quantify changes, we used posing the Cy3 and Cy5 images for each slide using the Scanalyze the mitogen-activated kinase 1 gene (MAPK1, Hs.324473) software (Eisen Lab). as a reference with the oligo pair: CCCTCGCATGACTGTTA- Scanning settings were used to equate intensity records of Cy3 CAGC and GGCACTTGGACAGTCATCTCTG. MAPK1 was and Cy5. Typically this means that Cy5 is scanned at 60% PMT not affected by EBNA1 expression by the microarray analysis. Oligo and Cy3 is scanned in 85% PMT. pairs with a Tm of 60 °C were chosen to generate 100 PCR The numerical analysis of the scanned results was performed products. To verify and ensure QRT-PCR specificity, we screened in GeneSpring 5 (Silicon Genetics) and data were normalized for the for each oligo pair using the ‘‘BLAST’’ (Basic Cy3/Cy5 ratio by applying LOWESS method to correct for Local Alignment Search Tool) protocol at the National Center for dye-specific intensity-dependant bias. Although this normaliza- Biotechnology Information (NCBI) web server. At the end of each tion produces more reliable results, it reduces the scores’ inten- QRT-PCR the synthesized products were verified for their speci- sity level. Data analysis of the normalized scores was carried out ficity by examining their actual Tm value. The sequence calculated by Access software (Microsoft) to pool results from several Tm for each PCR product was compared to the Tm measured experiments and then transferred to Excel software (Microsoft) through the decline of SYBR Green fluorescence from the melting for threshold cut off. To further reduce background signals, the strands in a gradual heat increase.

Canaan et al. www.pnas.org/cgi/content/short/0911676106 1of5 Fig. S1. QR-PCR analysis and agarose gel fractionation for CHUK and MUC13 promoters show the formation of specific amplicons in ChIP DNA samples. Samples of BJAB (Ϫ) and BJAB/EBNA1 (ϩ) from three biological replicates, which were amplified by QR-PCR have demonstrated similar melting curves characteristics (Left). These samples were also fractionated on an agarose gel to further demonstrate amplicon specificity and amplification differences. QR-PCR for the reference gene beta Actin is shown for these biological replicates for BJAB/EBNA1 (ϩ) and its control BJAB (Ϫ).

Canaan et al. www.pnas.org/cgi/content/short/0911676106 2of5 Fig. S2. Integrated genome browser view of ChIP-chip analysis for selected cellular promoters bound by EBNA1. The promoters for CHUK, MUC13, and OPN4 (Fig. 1, Fig. S1, and Table S2) are shown with the promoter for CCL18 (Table 1 and Fig. 1). The EBNA1 nonbound promoter for beta Actin (ACTB) is also shown as it was used as a reference for QR-PCR analysis and for the evaluation of enrichment folds. The EBNA1 nonbound PLAC9 and NRG3 promoter areas are also shown to further demonstrate selected promoter engagement by EBNA1.

Canaan et al. www.pnas.org/cgi/content/short/0911676106 3of5 Table S1. QR-PCR analysis of selected genes affected by EBNA1 in BJAB and 293 cell lines according to the microarray data Unigen Official gene symbol and name BJAB EBNA1/BJAB 293 EBNA1/293

Hs.202453 MYC, V-myc avian myelocytomatosis viral homolog 0.173 Ϯ 0.0132 0.451 Ϯ 0.062 Hs.46423 H4FG, H4 histone family, member G 3.272 Ϯ 0.22 8.938 Ϯ 0.331 Hs.75703 CCL4, Chemokine (C-C motif) ligand 4 (homologous to mouse Mip-1b) 49.58 Ϯ 0.721 1.989 Ϯ 0.66 Hs.524760 OAS1, 2Ј,5Ј-oligoadenylate synthetase 1 (40–46 kDa) 0.052 Ϯ 0.042 0.503 Ϯ 0.08 Hs.458485 ISG15, Interferon-stimulated protein, 15 kDa 0.1 Ϯ 0.052 0.615 Ϯ 0.05 Hs.480938 LRBA, LPS-responsive vesicle trafficking, beach and anchor 0.231 Ϯ 0.062 0.043 Ϯ 0.08 Hs.75256 RGS1, Regulator of G-protein signalling 1 5.897 Ϯ 0.7581 1.21 Ϯ 0.27 Hs.77100 GTF2E2, General IIE, polypeptide 2 (beta subunit, 34 kD) 4.777 Ϯ 1.01 0.558 Ϯ 0.05 Hs.75573 CENPE, Centromer protein E (312 kD) 0.575 Ϯ 0.112 2.049 Ϯ 0.01 Hs.522074 DSIPI, Delta sleep inducing peptide, immuno-reactor 0.817 Ϯ 0.02 2.412 Ϯ 0.221

Scores in table are given as a proportion of the curve slope values for EBNA1ϩ cell line over its cellular counterpart for the same PCR product. The slope was derived from the linear portion of the curve, which represents the fluorescent build up of the PCR products. Arrows indicate the score that initially was identified by microarray. Arrow direction represent up-regulation (1) and down-regulation (2) of the genes by microarray analysis. Dotted line separates the EBNA1-specific genes (above) from the tissue-specific genes (below). Scores are shown with their standard deviation that was generated from three QR-PCR samples. All samples were tested for EBNA1 expression by QR-PCR to verify EBNA1-positive and -negative samples.

Canaan et al. www.pnas.org/cgi/content/short/0911676106 4of5 Table S2. Selected cellular promoters bound by EBNA1, ChIP-chip analysis # Chr Start End Signal Up-stream Annotation chr10 102295884 102296348 1.048 HIF1AN Hypoxia-inducible factor 1 alpha inhibitor chr11 88402235 88402812 1.021 OPN4* Opsin 4 chr10 102080666 102080996 1.014 PKD2L1 polycystic kidney disease 2-like 1 chr2 227897845 227897985 0.945 HRB Homo sapiens HIV-1 Rev binding protein (HRB) chr10 102035411 102035630 0.904 BLOC1S2 Biogenesis of lysosome-related organelles complex-1, subunit 2 chr10 102554161 102554629 0.888 PAX2 PAIRED BOX GENE 2 chr2 133896524 133896774 0.886 NAP5 Nck-associated protein 5 chr2 134753804 134753944 0.877 MGAT5 Mannosyl (alpha-1,6-)-glycoprotein beta-1,6-N-acetyl-glucosaminyltransferase chr2 88705057 88705197 0.861 RPIA Ribose 5-phosphate isomerase A chr8 128815762 128816430 0.858 MYC V-myc avian myelocytomatosis viral oncogene homolog chr2 121946079 121946219 0.856 CLASP1 Cytoplasmic linker associated protein 1 chr2 37308770 37308910 0.852 PRKD3 Protein kinase D3 chr9 5175954 5176476 0.851 INSL6 Insulin like 6 chr9 97252681 97252853 0.833 FAM120A Hypothetical protein chr6 42303600 42303884 0.833 MRPS10 Mitochondrial ribosomal protein S10 chr9 129928476 129928648 0.806 LCN2 No annotations chr10 102083738 102084067 0.792 PKD2L1 polycystic kidney disease 2-like 1 chr10 68355365 68355588 0.781 LRRTM3 Leucine rich repeat transmembrane neuronal 3 chr2 192769403 192769543 0.779 No annotations chr10 102255489 102255818 0.779 SEC31B Orthologue of Saccharomyces cerevisiae Sec31p chr10 101980039 101980808 0.768 CHUK* Conserved helix-loop-helix ubiquitous kinase chr3 55669351 55669665 0.765 No annotations chr3 158637297 158637721 0.751 PTX3 Pentraxin-related gene chr3 129252793 129253107 0.742 SEC61A1 Sec61 alpha 1 subunit chr5 158569846 158570745 0.739 FLJ31951 Hypothetical protein chr10 102270633 102270852 0.695 NDUFB8 NADH dehydrogenase (ubiquinone) 1 beta subcomplex 8 chr2 27651565 27652041 0.694 C2orf16 Hypothetical protein chr2 166638811 166638951 0.685 No annotations chr5 78405898 78406137 0.684 DMGDH Dimethylglycine dehydrogenase chr3 130599114 130599428 0.674 IFT122 Intraflagellar transport 122 homolog chr3 184464156 184464803 0.666 B3GNT5 Beta-1,3-N-acetylglucosaminyltransferase 5 chr5 180723411 180723536 0.661 OR4F3 Olfactory , family 4, subfamily F, member 3 chr3 184465180 184465605 0.616 B3GNT5 Beta-1,3-N-acetylglucosaminyltransferase 5 chr3 152546822 152547577 0.616 MED12L Mediator of RNA polymerase II transcription, subunit chr2 98175670 98175920 0.61 No annotations chr3 170316520 170316834 0.609 No annotations chr3 145053498 145053922 0.579 SLC9A9 Solute carrier family 9 (sodium/hydrogen exchanger), member 9 chr9 15472952 15473454 0.578 PSIP1 PC4 and SFRS1 interacting protein 1 chr3 158704359 158704899 0.574 VEPH1 Ventricular zone expressed PH domain homolog 1 chr6 151755678 151756201 0.574 ZBTB2 Zinc finger and BTB domain containing 2 chr3 142937866 142938510 0.571 RNF7 Ring finger protein 7 chr3 152638014 152638328 0.526 No annotations chr3 126136455 126137008 0.525 MUC13* Mucin 13, cell surface associated chr3 180762251 180762565 0.524 ACTL6A Actin-like 6A chr3 199170929 199171243 0.524 IQCG IQ motif containing G

Sequences that were successfully verified by QR-PCR are shown in bold and in italic sequences that failed. *Genes that were further analyzed.

Canaan et al. www.pnas.org/cgi/content/short/0911676106 5of5