RNA-seq of shEZH2 cells

Downregulation of EZH2 Leads to Cellular Senescence with Features of SASP Cells were infected with a lentivirus vector expressing shRNA against EZH2 and harvested at 4 and 8 days after infection. Total RNA was harvested from cells using Trizol reagent (Invitrogen) and further purified using the Purelink RNA Mini kit (Invitrogen) with DNase I digestion. RNA library preparation with polyA selection and Illumina HiSeq 2x150bp sequencing was performed by GeneWiz Inc. Paired-end reads were quality trimmed using Trim galore v0.4.0 and subsequently aligned to the human reference genome, hg19, using HISAT2 v2.1.0. Reads mapping to annotated genes were quantified using featureCounts (Liao et al., 2014). Differential gene expression was determined using DESeq2 v1.12.4 (Love et al., 2014) and significance was defined as FDR-corrected p-values of <0.05. The log2 fold change for each gene was used to rank the list of genes for GSEAPreranked analysis (Subramanian et al., 2005). FPKM values were calculated using DESeq2 and Z‐scores were generated from FPKMs