Increasing evidence suggests the presence of minor cell subpopulations in prostate cancer that are androgen independent and poised for selection as dominant clones after androgen deprivation therapy. In this study, we investigated this phenomenon by stratifying cell subpopulations based on transcriptome profiling of 144 single LNCaP prostate cancer cells treated or untreated with androgen after cell-cycle synchronization. Model-based clustering of 397 differentially expressed genes identified eight potential subpopulations of LNCaP cells, revealing a previously unappreciable level of cellular heterogeneity to androgen stimulation. One subpopulation displayed stem-like features with a slower cell doubling rate, increased sphere formation capability, and resistance to G-M arrest induced by a mitosis inhibitor. Advanced growth of this subpopulation was associated with enhanced expression of 10 cell-cycle-related genes (, and ) and decreased dependence upon androgen receptor signaling. analysis of RNA-seq data from The Cancer Genome Atlas further demonstrated that concordant upregulation of these genes was linked to recurrent prostate cancers. Analysis of receiver operating characteristic curves implicates aberrant expression of these genes and could be useful for early identification of tumors that subsequently develop biochemical recurrence. Moreover, this single-cell approach provides a better understanding of how prostate cancer cells respond heterogeneously to androgen deprivation therapies and reveals characteristics of subpopulations resistant to this treatment. Illustrating the challenge in treating cancers with targeted drugs, which by selecting for drug resistance can drive metastatic progression, this study characterized the plasticity and heterogeneity of prostate cancer cells with regard to androgen dependence, defining the character or minor subpopulations of androgen-independent cells that are poised for clonal selection after androgen-deprivation therapy. .

Advanced prostate cancer is a very heterogeneous disease reflecting in diverse regulations of oncogenic signaling pathways. Aberrant spatial dynamics of epidermal growth factor receptor (EGFR) promote their dimerization and clustering, leading to constitutive activation in oncogenesis. The EphB2 and Src signaling pathways are associated with the reorganization of the cytoskeleton leading to malignancy, but their roles in regulating EGFR dynamics and activation are scarcely reported. Using single-particle tracking techniques, we found that highly phosphorylated EGFR in the advanced prostate cancer cell line, PC3, was associated with higher EGFR diffusivity, as compared with LNCaP and less aggressive DU145. The increased EGFR activation and biophysical dynamics were consistent with high proliferation, migration, and invasion. After performing single-cell RNA-seq on prostate cancer cell lines and circulating tumor cells from patients, we identified that upregulated gene expression in the EphB2 and Src pathways are associated with advanced malignancy. After dasatinib treatment or siRNA knockdowns of EphB2 or Src, the PC3 cells exhibited significantly lower EGFR dynamics, cell motility, and invasion. Partial inhibitory effects were also found in DU145 cells. The upregulation of parts of the EphB2 and Src pathways also predicts poor prognosis in the prostate cancer patient cohort of The Cancer Genome Atlas. Our results provide evidence that overexpression of the EphB2 and Src signaling pathways regulate EGFR dynamics and cellular aggressiveness in some advanced prostate cancer cells.

No abstract