Table S2) Are Shown with the Promoter for CCL18 (Table 1 and Fig
Total Page:16
File Type:pdf, Size:1020Kb
Supporting Information Canaan et al. 10.1073/pnas.0911676106 SI Text pooled scores from several arrays were averaged for each gene Tissue Culture. The human B cell lymphoma line, BJAB, and the and accepted if Ͼ2 for up-regulated genes 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 MYC 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 gene expression changes, we used posing the Cy3 and Cy5 images for each slide using the Scanalyze the mitogen-activated protein 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 base pair 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 human genome 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 oncogene 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 transcription factor 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