Systematic exploration of cell morphological phenotypes associated with a transcriptomic query

Deep-DPC: Deep learning-assisted label-free temporal imaging discovery of anti-fibrotic compounds by controlling cell morphology

INTRODUCTION

Fibrosis can damage the normal function of many organs, such as cardiac function, for which no effective clinical therapies exist. However, traditional approaches to anti-fibrosis drug discovery have primarily focused on the final biological indicators, often overlooking the dynamic morphological changes during fibrosis progression. Here, we present a novel approach, deep-DPC, which integrates label-free, time-series digital phase contrast (DPC) imaging with cell morphology analysis and unsupervised machine learning to dynamically control and monitor cell morphology.

OBJECTIVES

This method enables discrimination between resting and activated fibrocytes and facilitates the discovery of non-invasive labeled anti-fibrotic lead compounds.

METHODS

The deep-DPC comprises two major steps: (1) preliminary analysis by Harmony 4.9 software and (2) image classification via a neural network. For the experiment dataset, label-free time-series imaging was acquired from each well at 10 × magnification using the high-content imaging system, equipped with a high-speed charge-coupled device (CCD) camera. Dual-channel output images were generated through the imaging system, with one channel for bright-field and the other for DPC imaging, captured at 30-minute intervals. Firstly, applying the anti-fibrotic cell model as a case, a label-free time-series DPC imaging was developed by combining cell morphological analysis and deep learning, and its stability was verified by training with 12,000 images. Furthermore, the application of deep-DPC in the discovery of anti-fibrotic lead compounds.

RESULTS

Using the deep-DPC platform, over 100,000 images generated from 1,400 compounds were processed, identifying Neo-Przewaquinone A as a potent anti-fibrosis agent. Neo-Przewaquinone A exerts its effects by inhibiting TGF-β receptor I, thereby maintaining cells in a resting state and arresting the cell cycle.

CONCLUSION

The deep-DPC offers a promising strategy for fibrosis assessment by combining deep learning with dynamic cell morphology analysis based on time-series DPC images. Additionally, the platform holds potential as a novel therapeutic approach for anti-myocardial fibrosis by regulating cell morphology.

Cell Painting: a decade of discovery and innovation in cellular imaging

Modern quantitative image analysis techniques have enabled high-throughput, high-content imaging experiments. Image-based profiling leverages the rich information in images to identify similarities or differences among biological samples, rather than measuring a few features, as in high-content screening. Here, we review a decade of advancements and applications of Cell Painting, a microscopy-based cell-labeling assay aiming to capture a cell's state, introduced in 2013 to optimize and standardize image-based profiling. Cell Painting's ability to capture cellular responses to various perturbations has expanded owing to improvements in the protocol, adaptations for different perturbations, and enhanced methodologies for feature extraction, quality control, and batch-effect correction. Cell Painting is a versatile tool that has been used in various applications, alone or with other -omics data, to decipher the mechanism of action of a compound, its toxicity profile, and other biological effects. Future advances will likely involve computational and experimental techniques, new publicly available datasets, and integration with other high-content data types.

Systematic data analysis pipeline for quantitative morphological cell phenotyping

Quantitative morphological phenotyping (QMP) is an image-based method used to capture morphological features at both the cellular and population level. Its interdisciplinary nature, spanning from data collection to result analysis and interpretation, can lead to uncertainties, particularly among those new to this actively growing field. High analytical specificity for a typical QMP is achieved through sophisticated approaches that can leverage subtle cellular morphological changes. Here, we outline a systematic workflow to refine the QMP methodology. For a practical review, we describe the main steps of a typical QMP; in each step, we discuss the available methods, their applications, advantages, and disadvantages, along with the R functions and packages for easy implementation. This review does not cover theoretical backgrounds, but provides several references for interested researchers. It aims to broaden the horizons for future phenome studies and demonstrate how to exploit years of endeavors to achieve more with less.

High-Throughput Transcriptomics Screen of ToxCast Chemicals in U-2 OS Cells

New approach methodologies (NAMs) aim to accelerate the pace of chemical risk assessment while simultaneously reducing cost and dependency on animal studies. High Throughput Transcriptomics (HTTr) is an emerging NAM in the field of chemical hazard evaluation for establishing in vitro points-of-departure and providing mechanistic insight. In the current study, 1201 test chemicals were screened for bioactivity at eight concentrations using a 24-h exposure duration in the human- derived U-2 OS osteosarcoma cell line with HTTr. Assay reproducibility was assessed using three reference chemicals that were screened on every assay plate. The resulting transcriptomics data were analyzed by aggregating signal from genes into signature scores using gene set enrichment analysis, followed by concentration-response modeling of signatures scores. Signature scores were used to predict putative mechanisms of action, and to identify biological pathway altering concentrations (BPACs). BPACs were consistent across replicates for each reference chemical, with replicate BPAC standard deviations as low as 5.6 × 10-3 μM, demonstrating the internal reproducibility of HTTr-derived potency estimates. BPACs of test chemicals showed modest agreement (R2 = 0.55) with existing phenotype altering concentrations from high throughput phenotypic profiling using Cell Painting of the same chemicals in the same cell line. Altogether, this HTTr based chemical screen contributes to an accumulating pool of publicly available transcriptomic data relevant for chemical hazard evaluation and reinforces the utility of cell based molecular profiling methods in estimating chemical potency and predicting mechanism of action across a diverse set of chemicals.

Cell identity revealed by precise cell cycle state mapping links data modalities

Several methods for cell cycle inference from sequencing data exist and are widely adopted. In contrast, methods for classification of cell cycle state from imaging data are scarce. We have for the first time integrated sequencing and imaging derived cell cycle pseudo-times for assigning 449 imaged cells to 693 sequenced cells at an average resolution of 3.4 and 2.4 cells for sequencing and imaging data respectively. Data integration revealed thousands of pathways and organelle features that are correlated with each other, including several previously known interactions and novel associations. The ability to assign the transcriptome state of a profiled cell to its closest living relative, which is still actively growing and expanding opens the door for genotype-phenotype mapping at single cell resolution forward in time.

Evaluating batch correction methods for image-based cell profiling

High-throughput image-based profiling platforms are powerful technologies capable of collecting data from billions of cells exposed to thousands of perturbations in a time- and cost-effective manner. Therefore, image-based profiling data has been increasingly used for diverse biological applications, such as predicting drug mechanism of action or gene function. However, batch effects severely limit community-wide efforts to integrate and interpret image-based profiling data collected across different laboratories and equipment. To address this problem, we benchmark ten high-performing single-cell RNA sequencing (scRNA-seq) batch correction techniques, representing diverse approaches, using a newly released Cell Painting dataset, JUMP. We focus on five scenarios with varying complexity, ranging from batches prepared in a single lab over time to batches imaged using different microscopes in multiple labs. We find that Harmony and Seurat RPCA are noteworthy, consistently ranking among the top three methods for all tested scenarios while maintaining computational efficiency. Our proposed framework, benchmark, and metrics can be used to assess new batch correction methods in the future. This work paves the way for improvements that enable the community to make the best use of public Cell Painting data for scientific discovery.

Recent Advances in Graphene Oxide-Based on Organoid Culture as Disease Model and Cell Behavior - A Systematic Literature Review

Due to their ability to replicate the in vivo microenvironment through cell interaction and induce cells to stimulate cell function, three-dimensional cell culture models can overcome the limitations of two-dimensional models. Organoids are 3D models that demonstrate the ability to replicate the natural structure of an organ. In most organoid tissue cultures, matrigel made of a mouse tumor extracellular matrix protein mixture is an essential ingredient. However, its tumor-derived origin, batch-to-batch variation, high cost, and safety concerns have limited the usefulness of organoid drug development and regenerative medicine. Its clinical application has also been hindered by the fact that organoid generation is dependent on the use of poorly defined matrices. Therefore, matrix optimization is a crucial step in developing organoid culture that introduces alternatives as different materials. Recently, a variety of substitute materials has reportedly replaced matrigel. The purpose of this study is to review the significance of the latest advances in materials for cell culture applications and how they enhance build network systems by generating proper cell behavior. Excellence in cell behavior is evaluated from their cell characteristics, cell proliferation, cell differentiation, and even gene expression. As a result, graphene oxide as a matrix optimization demonstrated high potency in developing organoid models. Graphene oxide can promote good cell behavior and is well known for having good biocompatibility. Hence, advances in matrix optimization of graphene oxide provide opportunities for the future development of advanced organoid models.

Demultiplexing of single-cell RNA-sequencing data using interindividual variation in gene expression

Motivation

Pooled designs for single-cell RNA sequencing, where many cells from distinct samples are processed jointly, offer increased throughput and reduced batch variation. This study describes expression-aware demultiplexing (EAD), a computational method that employs differential co-expression patterns between individuals to demultiplex pooled samples without any extra experimental steps.

Results

We use synthetic sample pools and show that the top interindividual differentially co-expressed genes provide a distinct cluster of cells per individual, significantly enriching the regulation of metabolism. Our application of EAD to samples of six isogenic inbred mice demonstrated that controlling genetic and environmental effects can solve interindividual variations related to metabolic pathways. We utilized 30 samples from both sepsis and healthy individuals in six batches to assess the performance of classification approaches. The results indicate that combining genetic and EAD results can enhance the accuracy of assignments (Min. 0.94, Mean 0.98, Max. 1). The results were enhanced by an average of 1.4% when EAD and barcoding techniques were combined (Min. 1.25%, Median 1.33%, Max. 1.74%). Furthermore, we demonstrate that interindividual differential co-expression analysis within the same cell type can be used to identify cells from the same donor in different activation states. By analysing single-nuclei transcriptome profiles from the brain, we demonstrate that our method can be applied to nonimmune cells.

Availability and implementation

EAD workflow is available at https://isarnassiri.github.io/scDIV/ as an R package called scDIV (acronym for single-cell RNA-sequencing data demultiplexing using interindividual variations).

A Decade in a Systematic Review: The Evolution and Impact of Cell Painting

High-content image-based assays have fueled significant discoveries in the life sciences in the past decade (2013-2023), including novel insights into disease etiology, mechanism of action, new therapeutics, and toxicology predictions. Here, we systematically review the substantial methodological advancements and applications of Cell Painting. Advancements include improvements in the Cell Painting protocol, assay adaptations for different types of perturbations and applications, and improved methodologies for feature extraction, quality control, and batch effect correction. Moreover, machine learning methods recently surpassed classical approaches in their ability to extract biologically useful information from Cell Painting images. Cell Painting data have been used alone or in combination with other - omics data to decipher the mechanism of action of a compound, its toxicity profile, and many other biological effects. Overall, key methodological advances have expanded Cell Painting's ability to capture cellular responses to various perturbations. Future advances will likely lie in advancing computational and experimental techniques, developing new publicly available datasets, and integrating them with other high-content data types.

A Decade in a Systematic Review: The Evolution and Impact of Cell Painting

High-content image-based assays have fueled significant discoveries in the life sciences in the past decade (2013-2023), including novel insights into disease etiology, mechanism of action, new therapeutics, and toxicology predictions. Here, we systematically review the substantial methodological advancements and applications of Cell Painting. Advancements include improvements in the Cell Painting protocol, assay adaptations for different types of perturbations and applications, and improved methodologies for feature extraction, quality control, and batch effect correction. Moreover, machine learning methods recently surpassed classical approaches in their ability to extract biologically useful information from Cell Painting images. Cell Painting data have been used alone or in combination with other -omics data to decipher the mechanism of action of a compound, its toxicity profile, and many other biological effects. Overall, key methodological advances have expanded Cell Painting's ability to capture cellular responses to various perturbations. Future advances will likely lie in advancing computational and experimental techniques, developing new publicly available datasets, and integrating them with other high-content data types.

The Genomic Intersection of Oligodendrocyte Dynamics in Schizophrenia and Aging Unravels Novel Pathological Mechanisms and Therapeutic Potentials

Schizophrenia is a significant worldwide health concern, affecting over 20 million individuals and contributing to a potential reduction in life expectancy by up to 14.5 years. Despite its profound impact, the precise pathological mechanisms underlying schizophrenia continue to remain enigmatic, with previous research yielding diverse and occasionally conflicting findings. Nonetheless, one consistently observed phenomenon in brain imaging studies of schizophrenia patients is the disruption of white matter, the bundles of myelinated axons that provide connectivity and rapid signalling between brain regions. Myelin is produced by specialised glial cells known as oligodendrocytes, which have been shown to be disrupted in post-mortem analyses of schizophrenia patients. Oligodendrocytes are generated throughout life by a major population of oligodendrocyte progenitor cells (OPC), which are essential for white matter health and plasticity. Notably, a decline in a specific subpopulation of OPC has been identified as a principal factor in oligodendrocyte disruption and white matter loss in the aging brain, suggesting this may also be a factor in schizophrenia. In this review, we analysed genomic databases to pinpoint intersections between aging and schizophrenia and identify shared mechanisms of white matter disruption and cognitive dysfunction.

Evaluating batch correction methods for image-based cell profiling

High-throughput image-based profiling platforms are powerful technologies capable of collecting data from billions of cells exposed to thousands of perturbations in a time- and cost-effective manner. Therefore, image-based profiling data has been increasingly used for diverse biological applications, such as predicting drug mechanism of action or gene function. However, batch effects pose severe limitations to community-wide efforts to integrate and interpret image-based profiling data collected across different laboratories and equipment. To address this problem, we benchmarked seven high-performing scRNA-seq batch correction techniques, representing diverse approaches, using a newly released Cell Painting dataset, the largest publicly accessible image-based dataset. We focused on five different scenarios with varying complexity, and we found that Harmony, a mixture-model based method, consistently outperformed the other tested methods. Our proposed framework, benchmark, and metrics can additionally be used to assess new batch correction methods in the future. Overall, this work paves the way for improvements that allow the community to make best use of public Cell Painting data for scientific discovery.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

Robust mapping of spatiotemporal trajectories and cell-cell interactions in healthy and diseased tissues

Spatial transcriptomics (ST) technologies generate multiple data types from biological samples, namely gene expression, physical distance between data points, and/or tissue morphology. Here we developed three computational-statistical algorithms that integrate all three data types to advance understanding of cellular processes. First, we present a spatial graph-based method, pseudo-time-space (PSTS), to model and uncover relationships between transcriptional states of cells across tissues undergoing dynamic change (e.g. neurodevelopment, brain injury and/or microglia activation, and cancer progression). We further developed a spatially-constrained two-level permutation (SCTP) test to study cell-cell interaction, finding highly interactive tissue regions across thousands of ligand-receptor pairs with markedly reduced false discovery rates. Finally, we present a spatial graph-based imputation method with neural network (stSME), to correct for technical noise/dropout and increase ST data coverage. Together, the algorithms that we developed, implemented in the comprehensive and fast stLearn software, allow for robust interrogation of biological processes within healthy and diseased tissues.

High throughput microscopy and single cell phenotypic image-based analysis in toxicology and drug discovery

Measuring single cell responses to the universe of chemicals (drugs, natural products, environmental toxicants etc.) is of paramount importance to human health as phenotypic variability in sensing stimuli is a hallmark of biology that is considered during high throughput screening. One of the ways to approach this problem is via high throughput, microscopy-based assays coupled with multi-dimensional single cell analysis methods. Here, we will summarize some of the efforts in this vast and growing field, focusing on phenotypic screens (e.g., Cell Painting), single cell analytics and quality control, with particular attention to environmental toxicology and drug screening. We will discuss advantages and limitations of high throughput assays with various end points and levels of complexity.

Evolution and impact of high content imaging

The field of high content imaging has steadily evolved and expanded substantially across many industry and academic research institutions since it was first described in the early 1990's. High content imaging refers to the automated acquisition and analysis of microscopic images from a variety of biological sample types. Integration of high content imaging microscopes with multiwell plate handling robotics enables high content imaging to be performed at scale and support medium- to high-throughput screening of pharmacological, genetic and diverse environmental perturbations upon complex biological systems ranging from 2D cell cultures to 3D tissue organoids to small model organisms. In this perspective article the authors provide a collective view on the following key discussion points relevant to the evolution of high content imaging: • Evolution and impact of high content imaging: An academic perspective • Evolution and impact of high content imaging: An industry perspective • Evolution of high content image analysis • Evolution of high content data analysis pipelines towards multiparametric and phenotypic profiling applications • The role of data integration and multiomics • The role and evolution of image data repositories and sharing standards • Future perspective of high content imaging hardware and software.

Assessment of Drug-Induced Liver Injury through Cell Morphology and Gene Expression Analysis

Drug-induced liver injury (DILI) is a significant concern in drug development, often leading to drug withdrawal. Although many studies aim to identify biomarkers and gene/pathway signatures related to liver toxicity and aim to predict DILI compounds, this remains a challenge in drug discovery. With a strong development of high-content screening/imaging (HCS/HCI) for phenotypic screening, we explored the morphological cell perturbations induced by DILI compounds. In the first step, cell morphological signatures were associated with two datasets of DILI chemicals (DILIRank and eTox). The mechanisms of action were then analyzed for chemicals having transcriptomics data and sharing similar morphological perturbations. Signaling pathways associated with liver toxicity (cell cycle, cell growth, apoptosis, ...) were then captured, and a hypothetical relation between cell morphological perturbations and gene deregulation was illustrated within our analysis. Finally, using the cell morphological signatures, machine learning approaches were developed to predict chemicals with a potential risk of DILI. Some models showed relevant performance with validation set balanced accuracies between 0.645 and 0.739. Overall, our findings demonstrate the utility of combining HCI with transcriptomics data to identify the morphological and gene expression signatures related to DILI chemicals. Moreover, our protocol could be extended to other toxicity end points, offering a promising avenue for comprehensive toxicity assessment in drug discovery.

Using Transcriptomics and Cell Morphology Data in Drug Discovery: The Long Road to Practice

Gene expression and cell morphology data are high-dimensional biological readouts of much recent interest for drug discovery. They are able to describe biological systems in different states (e.g., healthy and diseased), as well as biological systems before and after compound treatment, and they are hence useful for matching both spaces (e.g., for drug repurposing) as well as for characterizing compounds with respect to efficacy and safety endpoints. This Microperspective describes recent advances in this direction with a focus on applied drug discovery and drug repurposing, as well as outlining what else is needed to advance further, with a particular focus on better understanding the applicability domain of readouts and their relevance for decision making, which is currently often still unclear.

Application of cell painting for chemical hazard evaluation in support of screening-level chemical assessments

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ~100 μM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Linking chemicals, genes and morphological perturbations to diseases

The progress in image-based high-content screening technology has facilitated high-throughput phenotypic profiling notably the quantification of cell morphology perturbation by chemicals. However, understanding the mechanism of action of a chemical and linking it to cell morphology and phenotypes remains a challenge in drug discovery. In this study, we intended to integrate molecules that induced transcriptomic perturbations and cellular morphological changes into a biological network in order to assess chemical-phenotypic relationships in humans. Such a network was enriched with existing disease information to suggest molecular and cellular profiles leading to phenotypes. Two datasets were used for this study. Firstly, we used the "Cell Painting morphological profiling assay" dataset, composed of 30,000 compounds tested on human osteosarcoma cells (named U2OS). Secondly, we used the "L1000 mRNA profiling assay" dataset, a collection of transcriptional expression data from cultured human cells treated with approximately 20,000 bioactive small molecules from the Library of Integrated Network-based Cellular Signatures (LINCS). Furthermore, pathways, gene ontology terms and disease enrichments were performed on the transcriptomics data. Overall, our study makes it possible to develop a biological network combining chemical-gene-pathway-morphological perturbation and disease relationships. It contains an ensemble of 9989 chemicals, 732 significant morphological features and 12,328 genes. Through diverse examples, we demonstrated that some drugs shared similar genes, pathways and morphological profiles that, taken together, could help in deciphering chemical-phenotype observations.

Probabilistic embedding, clustering, and alignment for integrating spatial transcriptomics data with PRECAST

Spatially resolved transcriptomics involves a set of emerging technologies that enable the transcriptomic profiling of tissues with the physical location of expressions. Although a variety of methods have been developed for data integration, most of them are for single-cell RNA-seq datasets without consideration of spatial information. Thus, methods that can integrate spatial transcriptomics data from multiple tissue slides, possibly from multiple individuals, are needed. Here, we present PRECAST, a data integration method for multiple spatial transcriptomics datasets with complex batch effects and/or biological effects between slides. PRECAST unifies spatial factor analysis simultaneously with spatial clustering and embedding alignment, while requiring only partially shared cell/domain clusters across datasets. Using both simulated and four real datasets, we show improved cell/domain detection with outstanding visualization, and the estimated aligned embeddings and cell/domain labels facilitate many downstream analyses. We demonstrate that PRECAST is computationally scalable and applicable to spatial transcriptomics datasets from different platforms.

A Novel Nanosafety Approach Using Cell Painting, Metabolomics, and Lipidomics Captures the Cellular and Molecular Phenotypes Induced by the Unintentionally Formed Metal-Based (Nano)Particles

Additive manufacturing (AM) or industrial 3D printing uses cutting-edge technologies and materials to produce a variety of complex products. However, the effects of the unintentionally emitted AM (nano)particles (AMPs) on human cells following inhalation, require further investigations. The physicochemical characterization of the AMPs, extracted from the filter of a Laser Powder Bed Fusion (L-PBF) 3D printer of iron-based materials, disclosed their complexity, in terms of size, shape, and chemistry. Cell Painting, a high-content screening (HCS) assay, was used to detect the subtle morphological changes elicited by the AMPs at the single cell resolution. The profiling of the cell morphological phenotypes, disclosed prominent concentration-dependent effects on the cytoskeleton, mitochondria, and the membranous structures of the cell. Furthermore, lipidomics confirmed that the AMPs induced the extensive membrane remodeling in the lung epithelial and macrophage co-culture cell model. To further elucidate the biological mechanisms of action, the targeted metabolomics unveiled several inflammation-related metabolites regulating the cell response to the AMP exposure. Overall, the AMP exposure led to the internalization, oxidative stress, cytoskeleton disruption, mitochondrial activation, membrane remodeling, and metabolic reprogramming of the lung epithelial cells and macrophages. We propose the approach of integrating Cell Painting with metabolomics and lipidomics, as an advanced nanosafety methodology, increasing the ability to capture the cellular and molecular phenotypes and the relevant biological mechanisms to the (nano)particle exposure.

Using chemical and biological data to predict drug toxicity

Various sources of information can be used to better understand and predict compound activity and safety-related endpoints, including biological data such as gene expression and cell morphology. In this review, we first introduce types of chemical, in vitro and in vivo information that can be used to describe compounds and adverse effects. We then explore how compound descriptors based on chemical structure or biological perturbation response can be used to predict safety-related endpoints, and how especially biological data can help us to better understand adverse effects mechanistically. Overall, the described applications demonstrate how large-scale biological information presents new opportunities to anticipate and understand the biological effects of compounds, and how this can support predictive toxicology and drug discovery projects.

New roles of DNA-binding and forkhead-associated domains of Fkh1 and Fkh2 in cellular functions

Two yeast forkhead transcription factors Fkh1 and Fkh2 regulate the transcription of CLB2 cluster genes important for mitosis. Both proteins contain a DNA-binding domain (DBD) and a forkhead-associated domain (FHAD), which are essential for ternary complex formation with transcription factor Mcm1, the transcription of CLB2 cluster genes and the physical interaction with Ndd1 and Clb2. Fkh2 also contains an additional C' domain that contains six consensus Cdk phosphorylation sites, but the function of this domain is dispensable. Here, we found new roles of the DBD, the FHAD, and the C' domain of Fkh1 and Fkh2 in cellular functions. The Fkh2 DBD determines the genetic interaction with NDD1, while both the FHAD and DBD of Fkh1 or Fkh2 determine cell morphology and stability of their own transcripts. Both HFADs, but not DBDs, also mediate physical interaction between Fkh1 and Fkh2. DBD and HFAD of Fkh1 and DBD, but not HFAD, of Fkh2 are also fundamental for nuclear localization. However, the Fkh2-specific C' domain has no role in these aspects except in the stability of some fkh mutant transcripts, which is either increased or decreased in the presence of this domain. These findings reveal that Fkh1 and Fkh2 have multiple cellular functions and function mainly via their DBD and FHAD through a domain-controlled feedback regulation mechanism.

High-dimensional gene expression and morphology profiles of cells across 28,000 genetic and chemical perturbations

Cells can be perturbed by various chemical and genetic treatments and the impact on gene expression and morphology can be measured via transcriptomic profiling and image-based assays, respectively. The patterns observed in these high-dimensional profile data can power a dozen applications in drug discovery and basic biology research, but both types of profiles are rarely available for large-scale experiments. Here, we provide a collection of four datasets with both gene expression and morphological profile data useful for developing and testing multimodal methodologies. Roughly a thousand features are measured for each of the two data types, across more than 28,000 chemical and genetic perturbations. We define biological problems that use the shared and complementary information in these two data modalities, provide baseline analysis and evaluation metrics for multi-omic applications, and make the data resource publicly available ( https://broad.io/rosetta/ ).

Morphology and gene expression profiling provide complementary information for mapping cell state

Morphological and gene expression profiling can cost-effectively capture thousands of features in thousands of samples across perturbations by disease, mutation, or drug treatments, but it is unclear to what extent the two modalities capture overlapping versus complementary information. Here, using both the L1000 and Cell Painting assays to profile gene expression and cell morphology, respectively, we perturb human A549 lung cancer cells with 1,327 small molecules from the Drug Repurposing Hub across six doses, providing a data resource including dose-response data from both assays. The two assays capture both shared and complementary information for mapping cell state. Cell Painting profiles from compound perturbations are more reproducible and show more diversity but measure fewer distinct groups of features. Applying unsupervised and supervised methods to predict compound mechanisms of action (MOAs) and gene targets, we find that the two assays not only provide a partially shared but also a complementary view of drug mechanisms. Given the numerous applications of profiling in biology, our analyses provide guidance for planning experiments that profile cells for detecting distinct cell types, disease phenotypes, and response to chemical or genetic perturbations.

Discrete protein metric (DPM): A new image similarity metric to calculate accuracy of deep learning-generated cell focal adhesion predictions

Understanding cell behaviors can provide new knowledge on the development of different pathologies. Focal adhesion (FA) sites are important sub-cellular structures that are involved in these processes. To better facilitate the study of FA sites, deep learning (DL) can be used to predict FA site morphology based on limited microscopic datasets (e.g., cell membrane images). However, calculating the accuracy score of these predictions can be challenging due to the discrete/point pattern like nature of FA sites. In the present work, a new image similarity metric, discrete protein metric (DPM), was developed to calculate FA prediction accuracy. This metric measures differences in distribution (d), shape/size (s), and angle (a) of FA sites between predicted and ground truth microscopy images. Performance of the DPM was evaluated by comparing it to three other commonly used image similarity metrics: Pearson correlation coefficient (PCC), feature similarity index (FSIM), and Intersection over Union (IoU). A sensitivity analysis was performed by comparing changes in each metric value due to quantifiable changes in FA site location, number, aspect ratio, area, or orientation. Furthermore, accuracy score of DL-generated predictions was calculated using all four metrics to compare their ability to capture variation across samples. Results showed better sensitivity and range of variation for DPM compared to the other metrics tested. Most importantly, DPM had the ability to determine which FA predictions were quantitatively more accurate and consistent with qualitative assessments. The proposed DPM hence provides a method to validate DL-generated FA predictions and has the potential to be used for investigation of other sub-cellular protein aggregates relevant to cell biology.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Accelerating drug repurposing for COVID-19 treatment by modeling mechanisms of action using cell image features and machine learning

The novel coronavirus disease, COVID-19, has rapidly spread worldwide. Developing methods to identify the therapeutic activity of drugs based on phenotypic data can improve the efficiency of drug development. Here, a state-of-the-art machine-learning method was used to identify drug mechanism of actions (MoAs) based on the cell image features of 1105 drugs in the  LINCS database. As the multi-dimensional features of cell images are affected by non-experimental factors, the characteristics of similar drugs vary considerably, and it is difficult to effectively identify the MoA of drugs as there is substantial noise. By applying the supervised information theoretic metric-learning (ITML) algorithm, a linear transformation made drugs with the same MoA aggregate. By clustering drugs to communities and performing enrichment analysis, we found that transferred image features were more conducive to the recognition of drug MoAs. Image features analysis showed that different features play important roles in identifying different drug functions. Drugs that significantly affect cell survival or proliferation, such as cyclin-dependent kinase inhibitors, were more likely to be enriched in communities, whereas other drugs might be decentralized. Chloroquine and clomiphene, which block the entry of virus, were clustered into the same community, indicating that similar MoA could be reflected by the cell image. Overall, the findings of the present study laid the foundation for the discovery of MoAs of new drugs, based on image data. In addition, it provided a new method of drug repurposing for COVID-19.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-021-09727-5.

Building Organs Using Tissue-Specific Microenvironments: Perspectives from a Bioprosthetic Ovary

Recent research in tissue engineering and regenerative medicine has elucidated the importance of the matrisome. The matrisome, effectively the skeleton of an organ, provides physical and biochemical cues that drive important processes such as differentiation, proliferation, migration, and cellular morphology. Leveraging the matrisome to control these and other tissue-specific processes will be key to developing transplantable bioprosthetics. In the ovary, the physical and biological properties of the matrisome have been implicated in controlling the important processes of follicle quiescence and folliculogenesis. This expanding body of knowledge is being applied in conjunction with new manufacturing processes to enable increasingly complex matrisome engineering, moving closer to emulating tissue structure, composition, and subsequent functions which can be applied to a variety of tissue engineering applications.

Image-based profiling for drug discovery: due for a machine-learning upgrade?

Image-based profiling is a maturing strategy by which the rich information present in biological images is reduced to a multidimensional profile, a collection of extracted image-based features. These profiles can be mined for relevant patterns, revealing unexpected biological activity that is useful for many steps in the drug discovery process. Such applications include identifying disease-associated screenable phenotypes, understanding disease mechanisms and predicting a drug's activity, toxicity or mechanism of action. Several of these applications have been recently validated and have moved into production mode within academia and the pharmaceutical industry. Some of these have yielded disappointing results in practice but are now of renewed interest due to improved machine-learning strategies that better leverage image-based information. Although challenges remain, novel computational technologies such as deep learning and single-cell methods that better capture the biological information in images hold promise for accelerating drug discovery.

Pheno-RNA, a method to associate genes with a specific phenotype, identifies genes linked to cellular transformation

Cellular transformation is associated with dramatic changes in gene expression, but it is difficult to determine which regulated genes are oncogenically relevant. Here we describe Pheno-RNA, a general approach to identifying candidate genes associated with a specific phenotype. Specifically, we generate a "phenotypic series" by treating a nontransformed breast cell line with a wide variety of molecules that induce cellular transformation to various extents. By performing transcriptional profiling across this phenotypic series, the expression profile of every gene can be correlated with the strength of the transformed phenotype. We identify ∼200 genes whose expression profiles are very highly correlated with the transformation phenotype, strongly suggesting their importance in transformation. Within biological categories linked to cancer, some genes show high correlations with the transformed phenotype, but others do not. Many genes whose expression profiles are highly correlated with transformation have never been associated with cancer, suggesting the involvement of heretofore unknown genes in cancer.

On the correlation between material-induced cell shape and phenotypical response of human mesenchymal stem cells

Learning rules by which cell shape impacts cell function would enable control of cell physiology and fate in medical applications, particularly, on the interface of cells and material of the implants. We defined the phenotypic response of human bone marrow-derived mesenchymal stem cells (hMSCs) to 2176 randomly generated surface topographies by probing basic functions such as migration, proliferation, protein synthesis, apoptosis, and differentiation using quantitative image analysis. Clustering the surfaces into 28 archetypical cell shapes, we found a very strict correlation between cell shape and physiological response and selected seven cell shapes to describe the molecular mechanism leading to phenotypic diversity. Transcriptomics analysis revealed a tight link between cell shape, molecular signatures, and phenotype. For instance, proliferation is strongly reduced in cells with limited spreading, resulting in down-regulation of genes involved in the G2/M cycle and subsequent quiescence, whereas cells with large filopodia are related to activation of early response genes and inhibition of the osteogenic process. In this paper we were aiming to identify a universal set of genes that regulate the material-induced phenotypical response of human mesenchymal stem cells. This will allow designing implants that can actively regulate cellular, molecular signalling through cell shape. Here we are proposing an approach to tackle this question.

A deep learning model to predict RNA-Seq expression of tumours from whole slide images

Deep learning methods for digital pathology analysis are an effective way to address multiple clinical questions, from diagnosis to prediction of treatment outcomes. These methods have also been used to predict gene mutations from pathology images, but no comprehensive evaluation of their potential for extracting molecular features from histology slides has yet been performed. We show that HE2RNA, a model based on the integration of multiple data modes, can be trained to systematically predict RNA-Seq profiles from whole-slide images alone, without expert annotation. Through its interpretable design, HE2RNA provides virtual spatialization of gene expression, as validated by CD3- and CD20-staining on an independent dataset. The transcriptomic representation learned by HE2RNA can also be transferred on other datasets, even of small size, to increase prediction performance for specific molecular phenotypes. We illustrate the use of this approach in clinical diagnosis purposes such as the identification of tumors with microsatellite instability.

RNA-sequencing of tumour tissue can provide important diagnostic and prognostic information but this is costly and not routinely performed in all clinical settings. Here, the authors show that whole slide histology slides—part of routine care—can be used to predict RNA-sequencing data and thus reduce the need for additional analyses.

Multivariate analysis reveals activation-primed fibroblast geometric states in engineered 3D tumor microenvironments

Fibroblasts are a heterogeneous group of cells comprising subpopulations that have been found to be activated in the stromal microenvironment that regulates tumor initiation and growth. The underlying mechanisms of such selective activation of fibroblasts are not understood. We propose that the intrinsic geometric heterogeneity of fibroblasts modulates the nuclear mechanotransduction of signals from the microenvironment, resulting in their selective activation. To test this, we developed an engineered 3D fibroblast tumor coculture system and used high resolution images to quantify multiple cell geometry sensitive nuclear morphological and chromatin organizational features. These features were then mapped to activation levels as measured by the nuclear abundance of transcription cofactor, megakaryoblastic leukemia, and protein levels of its target, αSMA. Importantly, our results indicate the presence of activation-“primed” cell geometries that present higher activation levels, which are further enhanced in the presence of stimuli from cancer cells. Further, we show that by enriching the population of activation-primed cell geometric states by either increasing matrix rigidity or micropatterning primed cell shapes, fibroblast activation levels can be increased. Collectively, our results reveal important cellular geometric states that select for fibroblast activation within the heterogenous tumor microenvironment.

Signatures of cell death and proliferation in perturbation transcriptomics data-from confounding factor to effective prediction

Transcriptional perturbation signatures are valuable data sources for functional genomics. Linking perturbation signatures to screenings opens the possibility to model cellular phenotypes from expression data and to identify efficacious drugs. We linked perturbation transcriptomics data from the LINCS-L1000 project with cell viability information upon genetic (Achilles project) and chemical (CTRP screen) perturbations yielding more than 90 000 signature–viability pairs. An integrated analysis showed that the cell viability signature is a major factor underlying perturbation signatures. The signature is linked to transcription factors regulating cell death, proliferation and division time. We used the cell viability–signature relationship to predict viability from transcriptomics signatures, and identified and validated compounds that induce cell death in tumor cell lines. We showed that cellular toxicity can lead to unexpected similarity of signatures, confounding mechanism of action discovery. Consensus compound signatures predicted cell-specific drug sensitivity, even if the signature is not measured in the same cell line, and outperformed conventional drug-specific features. Our results can help in understanding mechanisms behind cell death and removing confounding factors of transcriptomic perturbation screens. To interactively browse our results and predict cell viability in new gene expression samples, we developed CEVIChE (CEll VIability Calculator from gene Expression; https://saezlab.shinyapps.io/ceviche/).

Morphometric signatures of exposure to endocrine disrupting chemicals in zebrafish eleutheroembryos

Understanding the mode of action of the different pollutants in human and wildlife health is a key step in environmental risk assessment. The aim of this study was to determine signatures that could link morphological phenotypes to the toxicity mechanisms of four Endocrine Disrupting Chemicals (EDCs): bisphenol A (BPA), perfluorooctanesulfonate potassium salt (PFOS), tributyltin chloride (TBT), and 17-ß-estradiol (E2). Zebrafish (Danio rerio) eleutheroembryos were exposed from 2 to 5 dpf to a wide range of BPA, PFOS, TBT and E2 concentrations. At the end of the exposures several morphometric features were assessed. Common and non-specific effects on larvae pigmentation or swim bladder area were observed after exposures to all compounds. BPA specifically induced yolk sac malabsorption syndrome and altered craniofacial parameters, whereas PFOS had specific effects on the notochord formation presenting higher rates of scoliosis and kyphosis. The main effect of E2 was an increase in the body length of the exposed eleutheroembryos. In the case of TBT, main alterations on the morphological traits were related to developmental delays. When integrating all morphometrical parameters, BPA showed the highest rates of malformations in terms of equilethality, followed by PFOS and, distantly, by TBT and E2. In the case of BPA and PFOS, we were able to relate our results with effects on the transcriptome and metabolome, previously reported. We propose that methodized morphometric analyses in zebrafish embryo model can be used as an inexpensive and easy screening tool to predict modes of action of a wide-range number of contaminants.

Systems Pharmacogenomic Landscape of Drug Similarities from LINCS data: Drug Association Networks

Modern research in the biomedical sciences is data-driven utilizing high-throughput technologies to generate big genomic data. The Library of Integrated Network-based Cellular Signatures (LINCS) is an example for a large-scale genomic data repository providing hundred thousands of high-dimensional gene expression measurements for thousands of drugs and dozens of cell lines. However, the remaining challenge is how to use these data effectively for pharmacogenomics. In this paper, we use LINCS data to construct drug association networks (DANs) representing the relationships between drugs. By using the Anatomical Therapeutic Chemical (ATC) classification of drugs we demonstrate that the DANs represent a systems pharmacogenomic landscape of drugs summarizing the entire LINCS repository on a genomic scale meaningfully. Here we identify the modules of the DANs as therapeutic attractors of the ATC drug classes.