A multiplex single-cell RNA-Seq pharmacotranscriptomics pipeline for drug discovery

The gene-regulatory dynamics governing drug responses in cancer are yet to be fully understood. Here, we report a pipeline capable of producing high-throughput pharmacotranscriptomic profiling through live-cell barcoding using antibody–oligonucleotide conjugates. This pipeline combines drug screening with 96-plex single-cell RNA sequencing. We show the potential of this approach by exploring the heterogeneous transcriptional landscape of primary high-grade serous ovarian cancer (HGSOC) cells after treatment with 45 drugs, with 13 distinct classes of mechanisms of action. A subset of phosphatidylinositol 3-OH kinase (PI3K), protein kinase B (AKT) and mammalian target of rapamycin (mTOR) inhibitors induced the activation of receptor tyrosine kinases, such as the epithelial growth factor receptor (EGFR), and this was mediated by the upregulation of caveolin 1 (CAV1). This drug resistance feedback loop could be mitigated by the synergistic action of agents targeting PI3K–AKT–mTOR and EGFR for HGSOC with CAV1 and EGFR expression. Using this workflow could enable the personalized testing of patient-derived tumor samples at single-cell resolution.