Abstract
Paclitaxel is widely used in the combination chemotherapy for many cancers including esophageal squamous cell carcinoma (ESCC). However, the paclitaxel resistance occurs frequently in treating ESCC and the mechanism is not fully understood yet. The heterogeneity of gene expression within the drug-resistant cancer cells may be one of the major factors contributing to its resistance. In the present study, we successfully induced paclitaxel resistance in ESCC cell line KYSE-30 through low dose and long-term treatment of paclitaxel. Gene expression profiles were measured utilizing population RNA-seq and single-cell RNA-seq (scRNA-seq). 37 single cells from KYSE-30 cells and 73 single cells from paclitaxel resistant KYSE-30 cells (Taxol-R) were subjected to scRNA-seq. Weighted gene co-expression network analysis (WGCNA) of scRNA-seq data revealed two major subpopulations in both KYSE-30 and Taxol-R cancer cells. Two subpopulations based on the KRT19 expression levels in KYSE-30 cells exhibited different paclitaxel sensitivity, suggesting the existence of an intrinsic paclitaxel resistance in KYSE-30 cells. In addition, the Taxol-R cells that acquired the resistance to paclitaxel through induction were characterized with higher expressions of proteasomes but a lower expression of HIF-1 signaling genes. Furthermore, we showed that carfilzomib (CFZ), a proteasome inhibitor, could attenuate the paclitaxel resistance in Taxol-R cancer cells through activating the HIF-1 signaling. Our new finding may pave a way leading to an improvement in the treatment on cancers including ESCC by combining CFZ with paclitaxel as a novel approach for cancer therapy.