Differential gene expression analysis by RNA-seq reveals the importance of actin cytoskeletal proteins in erythroleukemia cells

Development of drug resistance limits the effectiveness of anticancer treatments. Understanding the molecular mechanisms triggering this event in tumor cells may lead to improved therapeutic strategies. Here we used RNA-seq to compare the transcriptomes of a murine erythroleukemia cell line (MEL) and a derived cell line with induced resistance to differentiation (MEL-R). RNA-seq analysis identified a total of 596 genes (Benjamini–Hochberg adjusted p-value < 0.05) that were differentially expressed by more than two-fold, of which 81.5% (486/596) of genes were up-regulated in MEL cells and 110 up-regulated in MEL-R cells. These observations revealed that for some genes the relative expression of mRNA amount in the MEL cell line has decreased as the cells acquired the resistant phenotype. Clustering analysis of a group of genes showing the highest differential expression allowed identification of a sub-group among genes up-regulated in MEL cells. These genes are related to the organization of the actin cytoskeleton network. Moreover, the majority of these genes are preferentially expressed in the hematopoietic lineage and at least three of them, Was (Wiskott Aldrich syndrome), Btk (Bruton’s tyrosine kinase) and Rac2, when mutated in humans, give rise to severe hematopoietic deficiencies. Among the group of genes that were up-regulated in MEL-R cells, 16% of genes code for histone proteins, both canonical and variants. A potential implication of these results on the blockade of differentiation in resistant cells is discussed.