RNAseq data for MDA-MB-231 cells selected for high and low ALDH expression

MDA-MB-231 cells were sorted via Florescence Associated Cell Sorting (FACS) us-ing the ALDEFLUOR™ Kit (Stem Cell Technologies, 01700, Canada). Cells were stained with the ALDEFLUOR™ Kit according to the manufactures recommendations and then sorted using the Sony LE-SH800S Cell Sorter for cells displaying the 5% lowest fluores-cence intensity (ALDH Low) and the 5% highest fluorescence intensity (ALDH High). Total RNA was purified using Monarch® Total RNA Miniprep Kit (New England Biolabs) following manufacturer’s instructions. Novogene’s (Novogene (UK) Company Limited, Cambridge, UK) RNA sequencing services were used for all mRNA sequencing. Briefly, an mRNA library was generated from total RNA, quantified, and sequenced using the Illumina NovaSeq 6000 platform (Illumina, San Diego, CA). Clean data was obtained from raw reads, and all downstream analysis were based on clean data. Ref-erence genome index was built using Hisat2 v2.0.5. Differential expression analysis was completed using DESeq2, and P-values were adjusted using the Benjamini and Hochberg's approach for controlling the false discovery rate. Genes with an adjusted P-value <=0.05 found by DESeq2 were assigned as differentially expressed. Enrichment analyses were completed using the clusterProfiler R package to test the statistical en-richment of differential expression genes in GO, KEGG, Reactome and DO databases.