Abstract
Recent advances in single-cell RNA sequencing (scRNA-seq) have endowed researchers with the ability to detect and analyze the transcriptomes of individual cancer cells. In the present study, 16,128 tumor cells from EGFR wild-type and EGFRvIII mutant cells were profiled by scRNA-seq. Analyses of scRNA-seq data from both U87MG and U87MG-EGFRvIII libraries revealed inherent heterogeneity in gene expression and biological processes. The cells stably expressing EGFRvIII showed enhanced transcriptional activities and a relatively homogeneous pattern, which manifested as less diverse distributions, gene expression levels and functional annotations compared with those of cells expressing the nonmutated version. Moreover, the differentially expressed genes between the U87MG and U87MG-EGFRvIII groups were mainly enriched in DNA replication, DNA repair and angiogenesis. We compared scRNA-seq data with bulk RNA-seq and EGFRvIII xenograft RNA-seq data. RAD51AP1 was shown to be upregulated in all three databases. Further analysis of RAD51AP1 revealed that it is an independent prognostic factor of glioma. Knocking down RAD51AP1 significantly inhibited tumor volume in an intracranial EGFRvIII-positive GBM model and prolonged survival time. Collectively, our microfluidic-based scRNA-seq driven by a single genetic event revealed a previously unappreciated implication of EGFRvIII in the heterogeneity of GBM and identified RAD51AP1 as an oncogene in glioma.