Supplemental Information. lncrna Epigenetic Landscape Analysis Identifies EPIC1 as an Oncogenic lncrna that Interacts with MYC

Transcription

1 Cancer Cell, Volume 33 Supplemental Information lncrna Epigenetic Landscape Analysis Identifies as an Oncogenic lncrna that Interacts with and Promotes Cell-Cycle Progression in Cancer Zehua Wang, Bo Yang, Min Zhang, Weiwei Guo, Zhiyuan Wu, Yue Wang, Lin Jia, Song Li, The Cancer Genome Atlas Research Network, Wen Xie, and Da Yang

2 A Data collection lncrna expression data (MiTranscriptome) Identification Correlation of methylation and expression DNA methylation level of 33 TCGA cancer types DNA methylation in normal tissue Nearest lncrna and PCG to probes DNA methylation in cancer sample B C Functional analysis Expression pattern in different tissue, development phase, cancer type or cancer subtype 1. Shared. PCG 3. lncrna. Intergenic Distance to lncrna TSS lncrna Protein-coding gene Transcription direction DNA methylation alteration TSS TSS Distance to PCG TSS To categorize patients into epigenetic activation or epigenetic silencing for each gene Prioritization lncrnas according to percentage of regulated patients Experimental analysis Methylated D Distance to lncrna TSS Coexpressed protein-coding gene and functional gene set or network Survival analysis TSS TSS DNA methylation alteration ESCA Distance to PCG TSS Figure S1. LncRNAs are both epigenetically activated and silenced by DNA methylation alteration in the promoter region, Related to Figure 1. (A) Flow chart of identification and functional analysis of EA and ES lncrnas in cancer. (B) Schematic of the annotation of DNA methylation probes to protein-coding genes and lncrna genes. (C) Differential DNA methylation between breast cancer and normal tissues. Density plot of average differential DNA methylation (indicated by beta values) within 1 windows in ± 1 kb from TSS sites are shown. The windows are arranged based on their distances to protein coding gene TSS (x-axis) and lncrna gene TSS (y-axis). The hypermethylation region in tumor is shown as red, whereas the hypomethylation region is shown as blue. (D) Differential DNA methylation between cancer and matched normal tissues in nine cancer types. TSS TSS TSS TSS

3 A LINC668 RP11-539E17.5 LINC91 RP11-68O1.1 AK333 CTD-31D.3 BC581 LINC6 RP11-556E13.1 AC LINC9 Nomal Tumor B LncRNA expression value DLBC ESCA LCLL MM SARC CERS3-AS1 AC66.13 HAND-AS1 FGF1-IT1 LINC1158 RP3-1A9.16 DIO3OS CTD-98J1. DPP1-AS1 RP11-7E.3 HHIP-AS1 LINC1197 AK15737 CTB-1I6. RP11-59O.1 AF AC MINCR SNHG1 MFI-AS1 LINC6 BC581 LINC91 LINC88 LINC668 LINC9

4 C MIR666A LINC9 LINC668 LINC88 RP11 539E17.5 LINC91 RP11 68O1.1 AK333 CTD 31D.3 BC581 LINC6 MFI AS1 RP11 556E13.1 LOC1797 SNHG1 BOLA MINCR RP11 78L15. AC AF ZNF667 AS1 RP11 59O.1 CTB 1I6. AK15737 LINC1197 HHIP AS1 RP11 7E.3 DPP1 AS1 CTD 98J1. DIO3OS RP3 1A9.16 LINC1158 BOLA3 AS1 FGF1 IT1 LOC19716 SNHG18 HAND AS1 AC66.13 CERS3 AS1 F HR CTD-13A5.3 LMO7-AS1 RP11-68O1.1 BOLA LOC LINC6 CTD-31D.3 LINC668 BC35 AK333 MIR5 ZFAS1 LINC673 RP11-3J1.1 LOC1576 AC ATPB3 ATP11A-AS1 LINC1197 ERBB MINCR RP11-1N16.3 LOC LOC BRAF RP11-69O1. CTB-1I6. CERS3-AS1 MIR138-1 EPHA5-AS1 ZNF793-AS1 TP53 ZNF667-AS1 KRAS RP11-1C8. p value <.1 <.1 <.1 <1 BC581 RP11-59O.1 LINC7 SETD AC AF LINC665 ARID1A PTEN LINC898 ATRX ZNF79-AS1 LINC1158 SLC6A-AS1 CTNNB1 NOTCH CCNE1 CTD-98J1. RP11-16E13. RNF3 RP11-31P5. DLX6-AS1 LINC118 PIK3CA GATA3 GNAS LINC958 CIC LINC61 PRDM1 LINC1116 LINC865 NF1 EMXOS SOX MIR666A FGF1-IT1 MEOX-AS1 BC361 RP11-573D15.8 RP11-3B7.3 D E HHIP-AS1 BCL6 Proportion Survival Proportion Survival Proportion Survival Proportion Survival RP3-1A9.16 p=. UC (5) EA (8) OS (month) p=. ES (3) UC (16) EA (3) OS (months) p=.6 OS (months) p=.1 ES (15) UC () ES (6) UC (33) LOC OS (months) RP11-6H.1 TPTEP1 BC7363 APC RP11-867G3.1 IDH1 FLJ355 LINC886 AC LOC RP11-17I1. LINC91 TP73-AS1 MEIS1-AS AC936. RP11-5H. Proportion Survival Proportion Survival Proportion Survival Proportion Survival PDZRN3-AS1 AC66.13 CDKNA RP11-195F19.9 MIR18A p=. OS (month) p=.6 UC (56) EA (85) UC (85) EA (181) p=. OS (months) ES (79) UC (193) OS (months) p=. ES (196) UC (98) RP11-89I19.1 DPP1-AS1 HPN-AS1 HAGLR EGFR OS (months) Tumor suppressor gene Oncogene ES lncrna EA lncrna Co-occurrence Mutual exclusivity MDM RP11-7E.3 SNORA69 TCONS_9157 MIR19A DM11953 SNHG18 AC863. AP51.3 GAS1RR AL WT1-AS RP11-55F3.9

5 G p value <.5 <5e 6 <5e 1 Mutual exclusivity Co-occurrence <5e <5e 8 SNHG1 RP11 68O1.1 RP11 556E13.1 RP11 539E17.5 RP11 78L15. MIR666A MINCR MFI AS1 LOC1797 LINC9 LINC91 LINC88 LINC668 LINC6 CTD 31D.3 BOLA BC581 AK333 AC UVM UCS THYM TGCT SARC PCPG MESO LAML ESCA DLBC CHOL ACC Figure S. EA lncrnas exhibit an on or off pattern with completely no expression in normal tissues, and are associated with tumor survival and tumor gene alterations, Related to Figure. (A) Representative expression pattern of EA lncrnas in multiple cancer types, compared to normal tissues. Blue dot denotes normal tissue, and red dot denotes tumor. (B) Expression pattern of EA and ES lncrnas in cancer cell lines. White denotes low expression, and dark blue denotes high expression. (C) Correlation of EA lncrnas (top panel) and ES lncrnas (bottom panel) with survival. The size of each cycle indicates the p value calculated by Cox regression survival analysis. The heatmap indicates the hazard ratio (HR). (D, E) Overall survival (OS) of representative EA lncrna LINC91 (D) and ES lncrna AF (E) in multiple cancer types. UC, unchanged; EA, epigenetic activation; ES, epigenetic silencing. (F) Mutual exclusivity and co-occurrence network for EA and ES lncrnas with tumor suppressor genes and oncogenes. (G) Mutual exclusivity and co-occurrence between EA lncrnas and TP53 mutation in 33 cancer types. The size of each circle indicates p value calculated by Fisher s exact test. The color indicates mutual exclusivity (blue) or co-occurrence (red).

6 A B D Proportion survival Proportion survival % 5% Proportion survival Proportion survival GSE711 logrank p =.53 Expression low high Time (months) GSE685 logrank p =. Expression low high Time (months) % EA p =.1 x 1-7 Normal -like Lum A p =. EA, n = 86 Non EA, n = Time (months) p =.739 EA, n = 3 Non EA, n = 8 Lum B Non EA HER Time (months) Proportion survival Proportion survival Basal GSE1653 logrank p =.5 Expression low high Time (months) E F GSE166 logrank p =.51 Expression low high Time (months) DNA mehtylatoin beta value (HM5) expression (RNA-seq) C expression Normal -like Lum A 1. p p.5 Proportion survival Proportion survival Lum B GSE1797 logrank p =.17 Expression low high Time (months) GSE19615 logrank p =.17 Expression low high 6 8 Time (months) HER Basal MB361 ZR-75-1 MB68 MB15 MB175 EFM-19 HCC38 HCC15 T-7D HCC1569 SK-BR-3 HCC1187 EFM-19A BT BT-7 MB36 HCC195 MCF-7 HCC119 MB53 MB31 ZR-75-3 HCC1937 CAMA-1

7 G H Relative level Relative expression RP11-78L15. HeLa 8 hr A DAC (mm) RP11-78L HeLa 7 hr Relative level DAC (mm) HeLa LINC6 SK DAC (mm) 8 hr 3 1 LINC6 A78 Relative level 8 hr HeLa T-7D DAC (mm) 8 hr hr 7 hr Hs578T 8 hr Relative expression CTD-31D.3 6 A78 8 hr CTD-31D A78cis 7 hr CTD-31D.3 6 MCF-7 8 hr 15 1 RP11-539E A78 8 hr RP11-539E MCF-7 8 hr RP11-539E MCF-7 7 hr Relative expression RP11-68O1.1 5 MCF-7 8 hr 15 1 RP11-68O1.1 5 MCF-7 7 hr RP11-556E13.1 MCF-7 8 hr RP11-556E MCF-7 7 hr RP11-556E Hs578T 8 hr Relative expression 15 1 AC MCF-7 8 hr AC MCF-7 7 hr AC A78cis 7 hr AC A78cis 8 hr AC A78 8 hr

8 Figure S3. is robustly correlated with poor survival in 189 breast cancer samples and can be epigenetically activated by decitabine, Related to Figure 3. (A) The association between expression and breast cancer survival in six independent breast cancer cohorts. (B) The association between s EA status and five breast cancer subtypes. The overall p value is calculated by chi-square test. (C) Comparison of s expression in five subtypes (B). The p value is calculated by pair-wise Wilcoxon Rank Sum test. p.5, p.1. (D) s correlation with survival in luminal B (top) and HER (bottom) subtypes. (E, F) Association of DNA methylation status (E) and expression (F) in breast cancer cell lines based on lncrna expression (CCLE) and genome-wide DNA methylation data (GSE837). (G) qrt-pcr analysis of expression in A78, SK--3, and T-7D cells treated with decitabine (DAC) as indicated time and dosage. (H) qrt-pcr analysis of eight EA lncrna expression in five cell lines treated with DAC. Error bars indicate mean ± SD, n = 3 for technical replicates. p.5, p.1.

9 A Copy number of 5 ( x 1 / μg RNA) C Cell line RNA copy # per cell MCF1A.5 ±.6 MCF ±.63 T-7D ± 3.1 ZR ±.3 MB ±.6 A78.56 ±.31 A78cis 19.6 ±.83 CAR-.87 ±.35 B - 39bp GAPDH - 15bp E. Marker MCF1A BT- BT-7 HCC1937 Hs578T MB31 MB361 MB68 MCF-7 T-7D ZR-75-1 a MCF1A BT- BT-7 HCC1937 Hs578T MB31 MB361 MB68 MCF-7 T-7D ZR-75-1 A78 A78cis IGR-1 CAR-3 CAR- CAR-8 SK--3 A78 A78cis IGR-1 CAR-3 CAR- CAR-8 SK--3 Marker HPDE AsPC-1 HPDE AsPC-1 BxPC-3 PANC-1 DU 15 BxPC-3 PANC-1 DU 15 PC-3 PC-3 A59 HeLa A-59 HeLa Hep G K56 Hep G K56 Water D kb chr: 7,, 7,5, 8,, 8,5, GENCODE v MCF-7 5 -Outer 5 -Inner T-7D 5 -Outer 5 -Inner MCF-7 T-7D 3 -Outer 3 -Inner 3 -Outer 3 -Inner ENST ENST65.1 b LOC893 uc3bik. uc6fgu.1 v3 RT-PCR v clones v1 ENST uc6fgv.1 c CpG: 7 Layered H3KMe3 Cons 1 Verts 5 RACE 3 RACE 13 - signal 1 _

10 F G > v1 AGTCCGCCATTGCAAACACGAAGCTCTTCCAGAAACGCCCTCACAGACACCCCGGAAGTCACGTACCCACTCTGTAGGTGCCCCGGGGCACAGGCAAGCG GACGAGCCAGTTATCCCTCAGAGCTCCTGCTGCCTCGCCCGCTTTCTCTCGGAAACGTGAAGTGTGGCCTCAGCTGAAAGTGAGGTGGGCCTCATTCAAT CAGTTGAATTCTTCAAGAGAGAAAAACTGAAGTCCCTTAGAAGGAAAGAGTTCTGCCTTCAGACTGTCTTTGAACTTAAGACTGTAGCGTCGACTCCTGC CGGAATTTCCAGCCTGCTGGCCAGCTCTGCAGATTCACACTTGCCAGCCTCCACAATCGTGTGAGCCAATTCCTTAACTTCTCTTTCTCCGTGTATCCCT TTGGTGCTGCCTCTCTGGGGAGCCCTGACTAATATGCATGCAGATGATACGGTGCCTGGCATTCTGAATACATGCACTAAATCCACCACTTTTCCCCATT TATAGATTTGGATTAACACACTAACTTACTCATATCTGCAAGTATAAATAAAAAAAATTGCTGGTGC > v AGTCCGCCATTGCAAACACGAAGCTCTTCCAGAAACGCCCTCACAGACACCCCGGAAGTCACGTACCCACTCTGTAGGTGCCCCGGGGCACAGGCAAGCG GACGAGCCAGTTATCCCTCAGAGCTCCTGCTGCCTCGCCCGCTTTCTCTCGGAAACGTGAAGTGTGGCCTCAGCTGAAAGTGAGGTGGGCCTCATTCAAT CAGTTGAATTCTTCAAGAGAGAAAAACTGAAGTCCCTTAGAAGGAAAGAGTTCTGCCTTCAGACTGTCTTTGAACTTAAGACTGTAGCGTCGACTCCTGC CGGAATTTCCAGCCTGCTGGCCAGCTCTGCAGATTCACACTTGCCAGCCTCCACAATCTTCCTGGATTTGAAACTGAAGAAGCAAGCAATCTGGAAATGT CAGTGGATGCACACAAAGAAACAACCGCAAAAGCCTGCTCGCTCTAGCCAAGGGACAAGAATAGGGGCAGTCCATCAAGACAGAATCCTTTTAAAAAATA ACCACTCCACTCCAGCAATACCACAGAAGAATCTGGCTGTACCCCAGGTACATCAGCAAAGATAACCTTTACCTAGCAGTAAAGAGGTCCCCCTTACACT GGGAGCCCTAGTGAAGAGCAGGGACTTTCACCCCCACTTAGCAGTGATGGGGCCCCACCCACCACAGTGCCAGCAGAGACCATGTGGGAGCCAGAATCCT CATCCCTACCCAGCAGTAACAAGGAGCCCTCCTCACTGCGGGCATCAAGGGTGAGTGAGTGCAAAACCTGGGTGTCACTCGGAAGGGAAGAATGGTGTCT CCTTCCTTCCCATCCCCTGCCAGAGTGATATCACTAGGAAAAAG > v3 AGTCCGCCATTGCAAACACGAAGCTCTTCCAGAAACGCCCTCACAGACACCCCGGAAGTCACGTACCCACTCTGTAGGTGCCCCGGGGCACAGGCAAGCG GACGAGCCAGTTATCCCTCAGAGCTCCTGCTGCCTCGCCCGCTTTCTCTCGGAAACGTGAAGTGTGGCCTCAGCTGAAAGTGAGGTGGGCCTCATTCAAT CAGTTGAATTCTTCAAGAGAGAAAAACTGAAGTCCCTTAGAAGGAAAGAGTTCTGCCTTCAGACTGTCTTTGAACTTAAGACTGTAGCGTCGACTCCTGC CGGAATTTCCAGCCTGCTGGCCAGCTCTGCAGATTCACACTTGCCAGCCTCCACAATCGCAGCTGAGGCGGAGGAACCCTAAGGGCTCATTGAGATCATG GATTTGCCCTTCTATGCATTGATGGAGCACCTGCTGCCCACAGCGTCTGTATTTGGTGCTGGGATGCTGAGCCTCCTTCTTTATGAATTTTTAAAAGGAC ACTGAGATCTTCAAACAGAGGCTGCCACTCTAAGCAAACAGATCCCGAGTCCTGGACTCTGAAGCTTGGGCCCAGTTCTCCTTTTCTCCGGGTTTCAGAT CCCACTGTAAAGTGAGGGGGCCCTTCTGATTCAGGACCCGGGGAAGCCAGGGGCATGAGCATCGGTGCCTCTTCTCTATTTCAAGGACCCTTCTGGGTGT AAAGTTCTCTGAGATGCCTTACATGGATTCCCACCACTGCAAGATAACCATCGTATGTAAAGTGTTATGACCAGCAGAGTGTAATTGAAGTGCATTCCAG AGGGAAAGACAGCGGCTCAGATTCTATTGAAAGAAACATGACATAATGATACCACAGCAAAAGCCAATCTTGCTCCTTTTTA Relative level sie3 sie sie5 sie6 Mock H Counts Counts MCF-7 PI ZR-75-1 I Relative level 1..5 shctrl she1 she J Relative growth rate 9 shctrl she1 6 she Day(s) K PI shctrl she1 she Relative # of colonies shctrl she1 she

11 Figure S. functions as an oncogenic lncrna regulating cell cycle progression, Related to Figure. (A, B) copy numbers by qrt-pcr (A) and RT-PCR analysis of the products (B) in multiple cancer-type cell lines. Copy numbers were calculated according to a standard curve of a serial dilutions of cdna from in vitro-transcribed ranging from 1 to 1 9 molecules. GAPDH was used as the internal control. (C) Quantification of RNA copy number/cell according to a standard curve of in vitro-transcribed. (D) 5 -RACE and 3 - RACE cloning of in MCF-7 and T-7D cells. (E) Alignment of with UCSC browser. a, Genomic location of in GENCODE is highlighted in background color of light yellow. The nearest protein coding genes (upstream TBC1DA and downstream FAM19A5) are also shown in the two ends. b, 's gene structure, isoforms from GENCODE (black), RefSeq (green) and UCSC (blue) annotation are enlarged. isoforms, i.e. v1, v, and v3, are also listed in the red window. c, The CpG island, H3KMe3 signal from ENCODE project and conservation tracks are presented at top. Sequences derived from 5'RACE and 3'RACE are listed in red window. signal from CCLE breast cancer cell lines are shown at bottom. (F) Sequences of isoforms cloned are shown. (G) qrt-pcr analysis of knockdown efficiency of sirnas in MCF-7 cells. (H) Cell cycle profiles of MCF-7 cells and ZR-75-1 cells treated with sirnas. (I-K) qrt-pcr analysis of expression (I), MTT assay (J), and anchorage-independent colony formation assays and representative images (K) of MCF-7 cells stably expressing shctrl and sh RNA, respectively. Error bars indicate mean ± SD, n = 3 for technical replicates. p.5, p.1.

12 A RNA distribution (relative fold) B E 1..5 sim1 sim 5 CDKN1A 3 1 sim1 sim 1..5 CCNA sim1 sim C RT-PCR WB 1..5 F G H Relative growth rate UCSC Gene Ensembl Gene Cyto1 Cyto Nuc1 Nuc Cyto1 Cyto Nuc1 Nuc Cyto1 Cyto Nuc1 Nuc GM1878 K56 MCF Hs578T GAPDH SNRP7 GAPDH sim1 sim U6 WCL Cyto Nuc Hs578T 1 3 Day(s) (chr:q13.31) Percentage of cells T-7D WCL Cyto Nuc GAPDH T-7D WCL Cyto Nuc U6 D Gene sim1 sim G/G1 S G/M Correlation coefficiency CDC WCL Cyto Nuc sim1 sim p value CDC.15 6.e-6 CDC6.1.e-1 CCNB. 1.e-9 FOXM1.1.1e-5 CDC5A.16.9e-7 CDC e-15 CCNA.18 1.e-8 EF1.8.e-16 EF e-7 MADL1. 1.e-11 RAD51.7.e Gene CDC5 sim1 sim sim1 sim Correlation coefficiency 1..5 Relative # of colonies U6 -GAPDH SNRP7 GAPDH p value RAD51AP e-5 FANCI.1.1e-11 CCNB1..1e-13 MAPK8IP..6e-1 MCM.1 1.3e-1 CDKN1A e-1 IL1R e-6 ITGA e-7 MAP e-3 NRP e-13 ENTPD -.1.e-5 sim1 sim sim1 sim

13 Figure S5. is a nuclear lncrna regulating targets expression, Related to Figure 5. (A) qrt-pcr analysis of expression and Western blot of subcellular fractionation in Hs578T (left panel) and T-7D (right panel) cells. (B) RNA-seq analysis of subcellular localization in MCF7, K56, and GM1878 cells. Subcellular localization RNA-seq data are downloaded from ENCODE. (C) RT-PCR analysis of expression levels in different subcellular fractionation of MCF-7 cells treated with sirnas. (D) Correlation between -regulated genes ( qrt-pcr validated genes in Figure 5F) and expression in TCGA tumors. (E) qrt-pcr analysis of targets expression in MCF-7 cells treated with sirnas. (F-H) MTT assay (F), cell cycle analysis (G), and anchorage-independent colony formation assays (H) of MCF-7 cells treated with sirnas (i.e., sim1 and sim). Error bars indicate mean ± SD, n = 3 for technical replicates. p.5, p.1.

14 A D E H % Input Vector Flag- v1 v v3 v1 v v3 IP: Anti-Flag IB: Anti-Flag Input IB: Anti-Flag Tumor volume (mm3) Enriched RNA (% Input) Percentage of cells(%)6.5 MCF GAPDH v1 v v3 v1 v v3 Vector Flag- G/G1 S G/M Vector Vector Week(s) DAC (μm) GAPDH β-actin B F Vector G Vector v1 v 1. 1 CDKN1A Enriched RNA(% Input) Enriched RNA (% Input) GAPDH CCNA Flag- Vector A78 C GAPDH CDC5 CDC Cyclin A p1 β-actin GAPDH Tumor weight (mg) v1 v vector IgG RIP + μm DAC RIP + 1 μm DAC RIP + 1 μm DAC RIP + 1 μm DAC DAC (μm) β-actin Vector IP Flag Input Flag IB: Flag 11 RT-PCR Western blot 1

15 Figure S6. directly binds with, Related to Figure 6. (A) Binding of isoforms with proteins in 93T cells. isoforms retrieved by exogenous Flag-tagged protein with Flag RIP were detected by qrt-pcr, and Western blot of Flag- is shown (right). (B, C) qrt-pcr analysis of targets (B), RT-PCR of products, and Western blot of targets in MCF-7 stably over-expressing isoforms (C). (D) Cell cycle analysis of MCF- 7 cells with stable overexpression of and an empty vector. (E, F) Tumor growth (E), tumor size and weight (F) in MCF-7 cells with stable overexpression of and an empty vector. Error bars indicate mean ± SD (n = 6). p.5, p.1. (G) Flag-tagged retrieves RNA (left). Western blot of Flag- is shown (right). (H) IP retrieves in MCF-7 and A78 cells after decitabine (DAC) treatment as indicated. Western blot of is also shown. Error bars indicate mean ± SD, n = 3 for technical replicates. p.5, p.1.

16 A CDC ChIP-seq Scale chr1: 1 kb hg38 3,359,5 3,36,5 3,361,5 3,36, rep1 _ rep _ si _ _ Scale CDC5 ChIP-seq rep1 _ rep _ si _ _ chr: 19,85, 1 kb hg38 19,95, 19,55, 19,515, ChIP-seq Scale chr: kb hg38 39,31, 39,316, 39,318, 39,3, RPL rep1 _ rep _ si _ _ Scale NR1D1 ChIP-seq rep1 _ rep _ si _ _ chr17: kb hg38,9,,96,,98,,1, ChIP-seq Scale 1 kb hg38 chr: 39,51, 39,51,5 39,5, 39,5,5 ATF rep1 _ rep _ si _ _ Scale APEX1 ChIP-seq 6.5- rep1 _ 6.5- rep _ si _ _ chr1:,55,5 1 kb hg38,56,,56,5,57,,57,5 RPL19 Scale chr17: 1 kb hg38 39,1, 39,, 39,3, 39,, RPL35 Scale chr9: ,859, 1 kb hg38 1,86,, 1,861, 1,86, ChIP-seq rep1 _ rep _ si _ _ ChIP-seq rep1 _ rep _ si _ _ B MCF APEX ATF 1..5 RPL RPL35 ZR APEX ATF 1..5 RPL RPL35

17 C D F % of Input %ofinput % of Input IP MAX.3 TRIP CD16 IgG Anti- Anti-MAX si si IgG band Vector si HA- Cyclin A p1 β-actin. MCM MXD AKAP EXO SCML PRKRIR. RGS si 1..7 G Relative level LNA.15 TROAP Input si Scr 1# # 6#.5 CKS MAFF # E. RPS6KA5.1 CACNB TFEB.6.. Relative Luc unit CCNA-Luc - - Δ- - - MAX -.5 RECQL5. RAD5L.15 APEX1.15 ATF RBBP RPL NS FOXM IgG. RPL35

18 Figure S7. is required for the regulatory role of in cancer, Related to Figure 7. (A) Alignment of two biological replicates of ChIP-seq in MCF-7 cells (green) and RNA-seq from (blue) and si (red) RNA treated MCF-7 cells. -regulated targets (CDC, CDC5, RPL3, and NR1D1) and non -regualted targets (ATF, APEX1, RPL19, and RPL35) are shown. (B) qrt-pcr validation of non -regualted targets expression in MCF-7 and ZR-75-1 cells with sirnas treatment. (C) ChIP-qPCR analysis of occupancy on the promoters of target genes in MCF-7 cells treated with sirnas. Non -regulated targets, i.e., APEX1, ATF, RPL19, and RPL35 are chosen as negative controls. (D) Effect of on and MAX complexes. Same quantity of lysates from MCF-7 cells with sirnas treatment was used for co-ip with anti- or anti-max antibodies. (E) Reporter assay of WT- and Δ- truncated mutant on CCNA-Luc. (F) Western blot of p1 and Cyclin A in MCF-7 cells transfected with sirnas followed by overexpression of indicated vectors. (G) qrt-pcr analysis of knockdown efficiency of LNA in MCF-7 cells. Error bars indicate mean ± SD, n = 3 for technical replicates. p.5, p.1. NS, not significant.