Doxorubicin resistance involves modulation of interferon signaling, transcriptional bursting, and gene co-expression patterns of U-ISGF3-related genes

Chemotherapy, although effective in treating cancer, can induce various cellular responses, including senescence and drug resistance. Here, we investigate the transcriptomic alterations induced by doxorubicin (DOX), a commonly used chemotherapeutic agent, in human colon cancer cells. Using single-cell RNA sequencing, we identified distinct cell populations and their transcriptional profiles following subtoxic DOX treatment, revealing cell clusters characterized by differential expression of genes involved in cell cycle regulation and interferon (IFN) signaling. DOX-persisting proliferating cells exhibited upregulation of genes reported to be linked to the unphosphorylated form of ISGF3 (U-ISGF3) transcription factor. Furthermore, we found that HSH2D, a poor prognostic marker, was highly upregulated in doxorubicin-surviving proliferative cells, and its expression was correlated with U-ISGF3-related genes. Analysis of transcription kinetics via mathematical modeling revealed that the number of mRNA molecules produced per transcriptional burst was increased for U-ISGF3-related genes. We also observed altered gene co-expression patterns of U-ISGF3-related genes and others upon DOX treatment, which potentially contributes to chemoresistance of DOX-surviving proliferative cells and may influence cancer cell fate after chemotherapy. Our findings highlight U-ISGF3-related genes and the JAK/STAT pathway as potential therapeutic targets for overcoming chemoresistance in colon cancer.