DRUG-seq for miniaturized high-throughput transcriptome profiling in drug discovery

Here we report Digital RNA with pertUrbation of Genes (DRUG-seq), a high-throughput platform for drug discovery. Pharmaceutical discovery relies on high-throughput screening, yet current platforms have limited readouts. RNA-seq is a powerful tool to investigate drug effects using transcriptome changes as a proxy, yet standard library construction is costly. DRUG-seq captures transcriptional changes detected in standard RNA-seq at 1/100th the cost. In proof-of-concept experiments profiling 433 compounds across 8 doses, transcription profiles generated from DRUG-seq successfully grouped compounds into functional clusters by mechanism of actions (MoAs) based on their intended targets. Perturbation differences reflected in transcriptome changes were detected for compounds engaging the same target, demonstrating the value of using DRUG-seq for understanding on and off-target activities. We demonstrate DRUG-seq captures common mechanisms, as well as differences between compound treatment and CRISPR on the same target. DRUG-seq provides a powerful tool for comprehensive transcriptome readout in a high-throughput screening environment.RNA-seq is a powerful tool to investigate how drugs affect the transcriptome but library construction can be costly. Here the authors introduce DRUG-seq, an automated platform for high-throughput transcriptome profiling.

[1]  Y. Furukawa,et al.  The FLT3 inhibitor PKC412 exerts differential cell cycle effects on leukemic cells depending on the presence of FLT3 mutations , 2008, Oncogene.

[2]  D. Haines,et al.  CDK9 Is Constitutively Expressed throughout the Cell Cycle, and Its Steady-State Expression Is Independent of SKP2 , 2003, Molecular and Cellular Biology.

[3]  Paul A Clemons,et al.  The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease , 2006, Science.

[4]  F. Tirode,et al.  Reconstitution of the transcription factor TFIIH: assignment of functions for the three enzymatic subunits, XPB, XPD, and cdk7. , 1999, Molecular cell.

[5]  I. Nakano,et al.  MELK—a conserved kinase: functions, signaling, cancer, and controversy , 2015, Clinical and Translational Medicine.

[6]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[7]  Jinsong Qiu,et al.  RASL‐seq for Massively Parallel and Quantitative Analysis of Gene Expression , 2012, Current protocols in molecular biology.

[8]  K. Eggan,et al.  Genetic Ablation of AXL Does Not Protect Human Neural Progenitor Cells and Cerebral Organoids from Zika Virus Infection. , 2016, Cell stem cell.

[9]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[10]  J. Strovel,et al.  BRD4 Structure-Activity Relationships of Dual PLK1 Kinase/BRD4 Bromodomain Inhibitor BI-2536. , 2015, ACS medicinal chemistry letters.

[11]  Baojie Li,et al.  LRRK2 interacts with ATM and regulates Mdm2–p53 cell proliferation axis in response to genotoxic stress , 2017, Human molecular genetics.

[12]  A Gordon,et al.  Classification, 2nd Edition , 1999 .

[13]  A. Minden,et al.  PAK4 is required for regulation of the cell‐cycle regulatory protein p21, and for control of cell‐cycle progression , 2011, Journal of cellular biochemistry.

[14]  E. Appella,et al.  Initiation of a G2/M checkpoint after ultraviolet radiation requires p38 kinase , 2001, Nature.

[15]  Tony Z. Jia,et al.  Digital RNA sequencing minimizes sequence-dependent bias and amplification noise with optimized single-molecule barcodes , 2012, Proceedings of the National Academy of Sciences.

[16]  M. Prunotto,et al.  Opportunities and challenges in phenotypic drug discovery: an industry perspective , 2017, Nature Reviews Drug Discovery.

[17]  Carlos L. Arteaga,et al.  Transforming Growth Factor (cid:2) Engages TACE and ErbB3 To Activate Phosphatidylinositol-3 Kinase/Akt in ErbB2-Overexpressing Breast Cancer and Desensitizes Cells to Trastuzumab (cid:1) † , 2008 .

[18]  Irina Neganova,et al.  An Important Role for CDK2 in G1 to S Checkpoint Activation and DNA Damage Response in Human Embryonic Stem Cells , 2011, Stem cells.

[19]  R. Blouin,et al.  Loss of DLK expression in WI-38 human diploid fibroblasts induces a senescent-like proliferation arrest. , 2011, Biochemical and biophysical research communications.

[20]  S. Soddu,et al.  CDKL5 localizes at the centrosome and midbody and is required for faithful cell division , 2017, Scientific Reports.

[21]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[22]  M. Eijken,et al.  Connectivity Map-based discovery of parbendazole reveals targetable human osteogenic pathway , 2015, Proceedings of the National Academy of Sciences.

[23]  A. Tan,et al.  Brusatol overcomes chemoresistance through inhibition of protein translation , 2017, Molecular carcinogenesis.

[24]  Andrea Califano,et al.  PLATE-Seq for genome-wide regulatory network analysis of high-throughput screens , 2017, Nature Communications.

[25]  L. Wodicka,et al.  Dual kinase-bromodomain inhibitors for rationally designed polypharmacology , 2014, Nature chemical biology.

[26]  D. Bojanic,et al.  Impact of high-throughput screening in biomedical research , 2011, Nature Reviews Drug Discovery.

[27]  Richard D Kolodner,et al.  An overview of Cdk1-controlled targets and processes , 2010, Cell Division.

[28]  M. Labow,et al.  Protein kinase D3 sensitizes RAF inhibitor RAF265 in melanoma cells by preventing reactivation of MAPK signaling. , 2011, Cancer research.

[29]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[30]  Angela N. Brooks,et al.  A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles , 2017, Cell.

[31]  P. Febbo,et al.  Application of a priori established gene sets to discover biologically important differential expression in microarray data , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[32]  J. Walter Evidence for Sequential Action of cdc7 and cdk2 Protein Kinases during Initiation of DNA Replication in Xenopus Egg Extracts* , 2000, The Journal of Biological Chemistry.

[33]  Qiang Yu,et al.  A p53-independent G1 Cell Cycle Checkpoint Induced by the Suppression of Protein Kinase C α and θ Isoforms* , 2003, Journal of Biological Chemistry.

[34]  N. Laird,et al.  Meta-analysis in clinical trials. , 1986, Controlled clinical trials.