Abstract
Cell-to-cell expression heterogeneity within a single tumor is a common phenotype among various cancer types including squamous cell carcinoma. To further study the fundamentals and importance of heterogeneity of cell functions and its potential mechanisms, we performed single-cell RNA-seq (scRNA-seq) on human squamous cell carcinoma of the bladder (SCCB) and its corresponding physiologically normal epithelia. Extensive differentially expressed genes were uncovered by comparing cancer and normal single cells, which were preferentially enriched in cancer-correlated pathways, such as p53 signaling and bladder cancer pathway. Furthermore, the most diversely expressed genes were particularly enriched in MAPK signaling pathway, such as CACNG4, CACNA1E and CACNA1H, which involve in cancer evolution and heterogeneity formation. Co-expression network and hub-gene analyses revealed several remarkable "hub genes" of each regulatory module. Some of them are cancer related, such as POU2F3, NKD1 and CYP2C8, while LINC00189, GCC2 and OR9Q1 genes are rarely reported in human diseases. The genes within an interesting module are highly correlated with others, which could be treated as potential targets for SCCB patients. Our findings have fundamental implications for SCCB biology and therapeutic strategies.