Signature Analysis of High-Throughput Transcriptomics Screening Data for Mechanistic Inference and Chemical Grouping

• Published in: Toxicological sciences 2024-11, Vol.202 (1), p.103-122
• Main Authors: Harrill, Joshua A, Everett, Logan J, Haggard, Derik E, Word, Laura J, Bundy, Joseph L, Chambers, Bryant, Harris, Felix, Willis, Clinton, Thomas, Russell S, Shah, Imran, Judson, Richard
• Format: Article
• Language: English
• ISSN: 1096-6080; 1096-0929; 1096-0929
• DOI: 10.1093/toxsci/kfae108

High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1751 unique chemicals from the EPA's ToxCast collection in MCF7 cells using eight concentrations and an exposure duration of 6 hours. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms-of-action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally-related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points-of-departure, predicting chemicals' molecular mechanism(s)-of-action (MeOA) and grouping chemicals with similar bioactivity profiles.