Data-analysis strategies for image-based cell profiling

A high-throughput platform for single-molecule tracking identifies drug interaction and cellular mechanisms

The regulation of cell physiology depends largely upon interactions of functionally distinct proteins and cellular components. These interactions may be transient or long-lived, but often affect protein motion. Measurement of protein dynamics within a cellular environment, particularly while perturbing protein function with small molecules, may enable dissection of key interactions and facilitate drug discovery; however, current approaches are limited by throughput with respect to data acquisition and analysis. As a result, studies using super-resolution imaging are typically drawing conclusions from tens of cells and a few experimental conditions tested. We addressed these limitations by developing a high-throughput single-molecule tracking (htSMT) platform for pharmacologic dissection of protein dynamics in living cells at an unprecedented scale (capable of imaging >106 cells/day and screening >104 compounds). We applied htSMT to measure the cellular dynamics of fluorescently tagged estrogen receptor (ER) and screened a diverse library to identify small molecules that perturbed ER function in real time. With this one experimental modality, we determined the potency, pathway selectivity, target engagement, and mechanism of action for identified hits. Kinetic htSMT experiments were capable of distinguishing between on-target and on-pathway modulators of ER signaling. Integrated pathway analysis recapitulated the network of known ER interaction partners and suggested potentially novel, kinase-mediated regulatory mechanisms. The sensitivity of htSMT revealed a new correlation between ER dynamics and the ability of ER antagonists to suppress cancer cell growth. Therefore, measuring protein motion at scale is a powerful method to investigate dynamic interactions among proteins and may facilitate the identification and characterization of novel therapeutics.

ImageDoubler: image-based doublet identification in single-cell sequencing

Single-cell sequencing provides detailed insights into individual cell behaviors within complex systems based on the assumption that each cell is uniquely isolated. However, doublets-where two or more cells are sequenced together-disrupt this assumption and can lead to potential data misinterpretations. Traditional doublet detection methods primarily rely on simulated genomic data, which may be less effective in homogeneous cell populations and can introduce biases from experimental processes. Therefore, we introduce ImageDoubler in this study, an innovative image-based model that identifies doublets and missing samples leveraging the Fluidigm single-cell sequencing image data. Our approach showcases a notable doublet detection efficacy, achieving a rate up to 93.87% and registering a minimum improvement of 33.1% in F1 scores compared to existing genomic-based methods. This advancement highlights the potential of using imaging to glean insight into developing doublet detection algorithms and exposes the limitations inherent in current genomic-based techniques.

Neuromorphic-enabled video-activated cell sorting

Imaging flow cytometry allows image-activated cell sorting (IACS) with enhanced feature dimensions in cellular morphology, structure, and composition. However, existing IACS frameworks suffer from the challenges of 3D information loss and processing latency dilemma in real-time sorting operation. Herein, we establish a neuromorphic-enabled video-activated cell sorter (NEVACS) framework, designed to achieve high-dimensional spatiotemporal characterization content alongside high-throughput sorting of particles in wide field of view. NEVACS adopts event camera, CPU, spiking neural networks deployed on a neuromorphic chip, and achieves sorting throughput of 1000 cells/s with relatively economic hybrid hardware solution (~$10 K for control) and simple-to-make-and-use microfluidic infrastructures. Particularly, the application of NEVACS in classifying regular red blood cells and blood-disease-relevant spherocytes highlights the accuracy of using video over a single frame (i.e., average error of 0.99% vs 19.93%), indicating NEVACS' potential in cell morphology screening and disease diagnosis.

A versatile information retrieval framework for evaluating profile strength and similarity

In profiling assays, thousands of biological properties are measured in a single test, yielding biological discoveries by capturing the state of a cell population, often at the single-cell level. However, for profiling datasets, it has been challenging to evaluate the phenotypic activity of a sample and the phenotypic consistency among samples, due to profiles' high dimensionality, heterogeneous nature, and non-linear properties. Existing methods leave researchers uncertain where to draw boundaries between meaningful biological response and technical noise. Here, we developed a statistical framework that uses the well-established mean average precision (mAP) as a single, data-driven metric to bridge this gap. We validated the mAP framework against established metrics through simulations and real-world data applications, revealing its ability to capture subtle and meaningful biological differences in cell state. Specifically, we used mAP to assess both phenotypic activity for a given perturbation (or a sample) as well as consistency within groups of perturbations (or samples) across diverse high-dimensional datasets. We evaluated the framework on different profile types (image, protein, and mRNA profiles), perturbation types (CRISPR gene editing, gene overexpression, and small molecules), and profile resolutions (single-cell and bulk). Our open-source software allows this framework to be applied to identify interesting biological phenomena and promising therapeutics from large-scale profiling data.

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.

Decoding phenotypic screening: A comparative analysis of image representations

Biomedical imaging techniques such as high content screening (HCS) are valuable for drug discovery, but high costs limit their use to pharmaceutical companies. To address this issue, The JUMP-CP consortium released a massive open image dataset of chemical and genetic perturbations, providing a valuable resource for deep learning research. In this work, we aim to utilize the JUMP-CP dataset to develop a universal representation model for HCS data, mainly data generated using U2OS cells and CellPainting protocol, using supervised and self-supervised learning approaches. We propose an evaluation protocol that assesses their performance on mode of action and property prediction tasks using a popular phenotypic screening dataset. Results show that the self-supervised approach that uses data from multiple consortium partners provides representation that is more robust to batch effects whilst simultaneously achieving performance on par with standard approaches. Together with other conclusions, it provides recommendations on the training strategy of a representation model for HCS images.

Applying cell painting in non-tumorigenic breast cells to understand impacts of common chemical exposures

The general population is exposed to many chemicals which have putative, but incompletely understood, links to breast cancer. Cell Painting is a high-content imaging-based in vitro assay that allows for unbiased measurements of concentration-dependent effects of chemical exposures on cellular morphology. We used Cell Painting to measure effects of 16 human exposure relevant chemicals, along with 21 small molecules with known mechanisms of action, in non-tumorigenic mammary epithelial cells, the MCF10A cell line. Using CellProfiler image analysis software, we quantified 3042 morphological features across approximately 1.2 million cells. We used benchmark concentration modeling to identify features both conserved and different across chemicals. Benchmark concentrations were compared to exposure biomarker concentration measurements from the National Health and Nutrition Examination Survey to assess which chemicals induce morphological alterations at human-relevant concentrations. We found significant feature overlaps between chemicals, including similarities between the organochlorine pesticide DDT metabolite p,p'-DDE and an activator of Wnt signaling CHIR99201. We validated these findings by assaying the activation of Wnt, as reflected by translocation of ꞵ-catenin, following p'-p' DDE exposure. Consistent with Wnt signaling activation, low concentration p',p'-DDE (25 nM) significantly enhanced the nuclear translocation of ꞵ-catenin. Overall, these findings highlight the ability of Cell Painting to enhance mode-of-action studies for toxicants which are common in our environment but incompletely characterized with respect to breast cancer risk.

Computational Methods for Image Analysis in Craniofacial Development and Disease

Observation is at the center of all biological sciences. Advances in imaging technologies are therefore essential to derive novel biological insights to better understand the complex workings of living systems. Recent high-throughput sequencing and imaging techniques are allowing researchers to simultaneously address complex molecular variations spatially and temporarily in tissues and organs. The availability of increasingly large dataset sizes has allowed for the evolution of robust deep learning models, designed to interrogate biomedical imaging data. These models are emerging as transformative tools in diagnostic medicine. Combined, these advances allow for dynamic, quantitative, and predictive observations of entire organisms and tissues. Here, we address 3 main tasks of bioimage analysis, image restoration, segmentation, and tracking and discuss new computational tools allowing for 3-dimensional spatial genomics maps. Finally, we demonstrate how these advances have been applied in studies of craniofacial development and oral disease pathogenesis.

SPACe: an open-source, single-cell analysis of Cell Painting data

Phenotypic profiling by high throughput microscopy, including Cell Painting, has become a leading tool for screening large sets of perturbations in cellular models. To efficiently analyze this big data, available open-source software requires computational resources usually not available to most laboratories. In addition, the cell-to-cell variation of responses within a population, while collected and analyzed, is usually averaged and unused. We introduce SPACe (Swift Phenotypic Analysis of Cells), an open-source platform for analysis of single-cell image-based morphological profiles produced by Cell Painting. We highlight several advantages of SPACe, including processing speed, accuracy in mechanism of action recognition, reproducibility across biological replicates, applicability to multiple models, sensitivity to variable cell-to-cell responses, and biological interpretability to explain image-based features. We illustrate SPACe in a defined screening campaign of cell metabolism small-molecule inhibitors tested in seven cell lines to highlight the importance of analyzing perturbations across models.

Phenomics Demonstrates Cytokines Additive Induction of Epithelial to Mesenchymal Transition

Epithelial to mesenchymal transition (EMT) is highly plastic with a programme where cells lose adhesion and become more motile. EMT heterogeneity is one of the factors for disease progression and chemoresistance in cancer. Omics characterisations are costly and challenging to use. We developed single cell phenomics with easy to use wide-field fluorescence microscopy. We analyse over 70,000 cells and combined 53 features. Our simplistic pipeline allows efficient tracking of EMT plasticity, with a single statistical metric. We discriminate four high EMT plasticity cancer cell lines along the EMT spectrum. We test two cytokines, inducing EMT in all cell lines, alone or in combination. The single cell EMT metrics demonstrate the additive effect of cytokines combination on EMT independently of cell line EMT spectrum. The effects of cytokines are also observed at the front of migration during wound healing assay. Single cell phenomics is uniquely suited to characterise the cellular heterogeneity in response to complex microenvironment and show potential for drug testing assays.

Pervasive mislocalization of pathogenic coding variants underlying human disorders

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease causing. This creates a new bottleneck: determining the functional impact of each variant-typically a painstaking, customized process undertaken one or a few genes and variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,448 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of uncertain significance. Our publicly available resource extends our understanding of coding variation in human diseases.

An automated high-content screening and assay platform for the analysis of spheroids at subcellular resolution

The endomembrane system is essential for healthy cell function, with the various compartments carrying out a large number of specific biochemical reactions. To date, almost all of our understanding of the endomembrane system has come from the study of cultured cells growing as monolayers. However, monolayer-grown cells only poorly represent the environment encountered by cells in the human body. As a first step to address this disparity, we have developed a platform that allows us to investigate and quantify changes to the endomembrane system in three-dimensional (3D) cell models, in an automated and highly systematic manner. HeLa Kyoto cells were grown on custom-designed micropatterned 96-well plates to facilitate spheroid assembly in the form of highly uniform arrays. Fully automated high-content confocal imaging and analysis were then carried out, allowing us to measure various spheroid-, cellular- and subcellular-level parameters relating to size and morphology. Using two drugs known to perturb endomembrane function, we demonstrate that cell-based assays can be carried out in these spheroids, and that changes to the Golgi apparatus and endosomes can be quantified from individual cells within the spheroids. We also show that image texture measurements are useful tools to discriminate cellular phenotypes. The automated platform that we show here has the potential to be scaled up, thereby allowing large-scale robust screening to be carried out in 3D cell models.

Capturing cell heterogeneity in representations of cell populations for image-based profiling using contrastive learning

Image-based cell profiling is a powerful tool that compares perturbed cell populations by measuring thousands of single-cell features and summarizing them into profiles. Typically a sample is represented by averaging across cells, but this fails to capture the heterogeneity within cell populations. We introduce CytoSummaryNet: a Deep Sets-based approach that improves mechanism of action prediction by 30-68% in mean average precision compared to average profiling on a public dataset. CytoSummaryNet uses self-supervised contrastive learning in a multiple-instance learning framework, providing an easier-to-apply method for aggregating single-cell feature data than previously published strategies. Interpretability analysis suggests that the model achieves this improvement by downweighting small mitotic cells or those with debris and prioritizing large uncrowded cells. The approach requires only perturbation labels for training, which are readily available in all cell profiling datasets. CytoSummaryNet offers a straightforward post-processing step for single-cell profiles that can significantly boost retrieval performance on image-based profiling datasets.

Exploring heterogeneous cell population dynamics in different microenvironments by novel analytical strategy based on images

Understanding the dynamic states and transitions of heterogeneous cell populations is crucial for addressing fundamental biological questions. High-content imaging provides rich datasets, but it remains increasingly difficult to integrate and annotate high-dimensional and time-resolved datasets to profile heterogeneous cell population dynamics in different microenvironments. Using hepatic stellate cells (HSCs) LX-2 as model, we proposed a novel analytical strategy for image-based integration and annotation to profile dynamics of heterogeneous cell populations in 2D/3D microenvironments. High-dimensional features were extracted from extensive image datasets, and cellular states were identified based on feature profiles. Time-series clustering revealed distinct temporal patterns of cell shape and actin cytoskeleton reorganization. We found LX-2 showed more complex membrane dynamics and contractile systems with an M-shaped actin compactness trend in 3D culture, while they displayed rapid spreading in early 2D culture. This image-based integration and annotation strategy enhances our understanding of HSCs heterogeneity and dynamics in complex extracellular microenvironments.

Dissecting the immune infiltrate of primary luminal B-like breast carcinomas in relation to age

The impact of aging on the immune landscape of luminal breast cancer (Lum-BC) is poorly characterized. Understanding the age-related dynamics of immune editing in Lum-BC is anticipated to improve the therapeutic benefit of immunotherapy in older patients. To this end, here we applied the 'multiple iterative labeling by antibody neo-deposition' (MILAN) technique, a spatially resolved single-cell multiplex immunohistochemistry method. We created tissue microarrays by sampling both the tumor center and invasive front of luminal breast tumors collected from a cohort of treatment-naïve patients enrolled in the prospective monocentric IMAGE (IMmune system and AGEing) study. Patients were subdivided into three nonoverlapping age categories (35-45 = 'young', n = 12; 55-65 = 'middle', n = 15; ≥70 = 'old', n = 26). Additionally, depending on localization and amount of cytotoxic T lymphocytes, the tumor immune types 'desert' (n = 22), 'excluded' (n = 19), and 'inflamed' (n = 12) were identified. For the MILAN technique we used 58 markers comprising phenotypic and functional markers allowing in-depth characterization of T and B lymphocytes (T&B-lym). These were compared between age groups and tumor immune types using Wilcoxon's test and Pearson's correlation. Cytometric analysis revealed a decline of the immune cell compartment with aging. T&B-lym were numerically less abundant in tumors from middle-aged and old compared to young patients, regardless of the geographical tumor zone. Likewise, desert-type tumors showed the smallest immune-cell compartment and were not represented in the group of young patients. Analysis of immune checkpoint molecules revealed a heterogeneous geographical pattern of expression, indicating higher numbers of PD-L1 and OX40-positive T&B-lym in young compared to old patients. Despite the numerical decline of immune infiltration, old patients retained higher expression levels of OX40 in T helper cells located near cancer cells, compared to middle-aged and young patients. Aging is associated with important numerical and functional changes of the immune landscape in Lum-BC. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

How scanning probe microscopy can be supported by artificial intelligence and quantum computing?

The impact of Artificial Intelligence (AI) is rapidly expanding, revolutionizing both science and society. It is applied to practically all areas of life, science, and technology, including materials science, which continuously requires novel tools for effective materials characterization. One of the widely used techniques is scanning probe microscopy (SPM). SPM has fundamentally changed materials engineering, biology, and chemistry by providing tools for atomic-precision surface mapping. Despite its many advantages, it also has some drawbacks, such as long scanning times or the possibility of damaging soft-surface materials. In this paper, we focus on the potential for supporting SPM-based measurements, with an emphasis on the application of AI-based algorithms, especially Machine Learning-based algorithms, as well as quantum computing (QC). It has been found that AI can be helpful in automating experimental processes in routine operations, algorithmically searching for optimal sample regions, and elucidating structure-property relationships. Thus, it contributes to increasing the efficiency and accuracy of optical nanoscopy scanning probes. Moreover, the combination of AI-based algorithms and QC may have enormous potential to enhance the practical application of SPM. The limitations of the AI-QC-based approach were also discussed. Finally, we outline a research path for improving AI-QC-powered SPM. RESEARCH HIGHLIGHTS: Artificial intelligence and quantum computing as support for scanning probe microscopy. The analysis indicates a research gap in the field of scanning probe microscopy. The research aims to shed light into ai-qc-powered scanning probe microscopy.

Characterization of fission and fusion mitochondrial dynamics in HD fibroblasts according to patient's severity status

Huntington's Disease (HD) is an inheritable neurodegenerative condition caused by an expanded CAG trinucleotide repeat in the HTT gene with a direct correlation between CAG repeats expansion and disease severity with earlier onset-of- disease. Previously we have shown that primary skin fibroblasts from HD patients exhibit unique phenotype disease features, including distinct nuclear morphology and perturbed actin cap linked with cell motility, that are correlated with the HD patient disease severity. Here we provide further evidence that mitochondrial fission-fusion morphology balance dynamics, classified using a custom image-based high-content analysis (HCA) machine learning tool, that improved correlation with HD severity status. This mitochondrial phenotype is supported by appropriate changes in fission-fusion biomarkers (Drp1, MFN1, MFN2, VAT1) levels in the HD patients' fibroblasts. These findings collectively point towards a dysregulation in mitochondrial dynamics, where both fission and fusion processes may be disrupted in HD cells compared to healthy controls. This study shows for the first time a methodology that enables identification of HD phenotype before patient's disease onset (Premanifest). Therefore, we believe that this tool holds a potential for improving precision in HD patient's diagnostics bearing the potential to evaluate alterations in mitochondrial dynamics throughout the progression of HD, offering valuable insights into the molecular mechanisms and drug therapy evaluation underlying biological differences in any disease stage.

Machine learning in microscopy - insights, opportunities and challenges

Machine learning (ML) is transforming the field of image processing and analysis, from automation of laborious tasks to open-ended exploration of visual patterns. This has striking implications for image-driven life science research, particularly microscopy. In this Review, we focus on the opportunities and challenges associated with applying ML-based pipelines for microscopy datasets from a user point of view. We investigate the significance of different data characteristics - quantity, transferability and content - and how this determines which ML model(s) to use, as well as their output(s). Within the context of cell biological questions and applications, we further discuss ML utility range, namely data curation, exploration, prediction and explanation, and what they entail and translate to in the context of microscopy. Finally, we explore the challenges, common artefacts and risks associated with ML in microscopy. Building on insights from other fields, we propose how these pitfalls might be mitigated for in microscopy.

Extracellular Matrix Components and Mechanosensing Pathways in Health and Disease

Glycosaminoglycans (GAGs) and proteoglycans (PGs) are essential components of the extracellular matrix (ECM) with pivotal roles in cellular mechanosensing pathways. GAGs, such as heparan sulfate (HS) and chondroitin sulfate (CS), interact with various cell surface receptors, including integrins and receptor tyrosine kinases, to modulate cellular responses to mechanical stimuli. PGs, comprising a core protein with covalently attached GAG chains, serve as dynamic regulators of tissue mechanics and cell behavior, thereby playing a crucial role in maintaining tissue homeostasis. Dysregulation of GAG/PG-mediated mechanosensing pathways is implicated in numerous pathological conditions, including cancer and inflammation. Understanding the intricate mechanisms by which GAGs and PGs modulate cellular responses to mechanical forces holds promise for developing novel therapeutic strategies targeting mechanotransduction pathways in disease. This comprehensive overview underscores the importance of GAGs and PGs as key mediators of mechanosensing in maintaining tissue homeostasis and their potential as therapeutic targets for mitigating mechano-driven pathologies, focusing on cancer and inflammation.

da_Tracker: Automated workflow for high throughput single cell and single phagosome tracking in infected cells

Time-lapse microscopy has emerged as a crucial tool in cell biology, facilitating a deeper understanding of dynamic cellular processes. While existing tracking tools have proven effective in detecting and monitoring objects over time, the quantification of signals within these tracked objects often faces implementation constraints. In the context of infectious diseases, the quantification of signals at localized compartments within the cell and around intracellular pathogens can provide even deeper insight into the interactions between the pathogen and host cell organelles. Existing quantitative analysis at a single-phagosome level remains limited and dependent on manual tracking methods. We developed a near-fully automated workflow that performs with limited bias, high-throughput cell segmentation and quantitative tracking of both single cell and single bacterium/phagosome within multi-channel, z-stack, time-lapse confocal microscopy videos. We took advantage of the PyImageJ library to bring Fiji functionality into a Python environment and combined deep-learning-based segmentation from Cellpose with tracking algorithms from Trackmate. The 'da_tracker' workflow provides a versatile toolkit of functions for measuring relevant signal parameters at the single-cell level (such as velocity or bacterial burden) and at the single-phagosome level (i.e. assessment of phagosome maturation over time). Its capabilities in both single-cell and single-phagosome quantification, its flexibility and open-source nature should assist studies that aim to decipher for example the pathogenicity of bacteria and the mechanism of virulence factors that could pave the way for the development of innovative therapeutic approaches.

Tumour immune characterisation of primary triple-negative breast cancer using automated image quantification of immunohistochemistry-stained immune cells

The tumour immune microenvironment (TIME) in breast cancer is acknowledged with an increasing role in treatment response and prognosis. With a growing number of immune markers analysed, digital image analysis may facilitate broader TIME understanding, even in single-plex IHC data. To facilitate analyses of the latter an open-source image analysis pipeline, Tissue microarray MArker Quantification (TMArQ), was developed and applied to single-plex stainings for p53, CD3, CD4, CD8, CD20, CD68, FOXP3, and PD-L1 (SP142 antibody) in a 218-patient triple negative breast cancer (TNBC) cohort with complementary pathology scorings, clinicopathological, whole genome sequencing, and RNA-sequencing data. TMArQ's cell counts for analysed immune markers were on par with results from alternative methods and consistent with both estimates from human pathology review, different quantifications and classifications derived from RNA-sequencing as well as known prognostic patterns of immune response in TNBC. The digital cell counts demonstrated how immune markers are coexpressed in the TIME when considering TNBC molecular subtypes and DNA repair deficiency, and how combination of immune status with DNA repair deficiency status can improve the prognostic stratification in chemotherapy treated patients. These results underscore the value and potential of integrating TIME and specific tumour intrinsic alterations/phenotypes for the molecular understanding of TNBC.

Hit screening with multivariate robust outlier detection

Hit screening, which involves the identification of compounds or targets capable of modulating disease-relevant processes, is an important step in drug discovery. Some assays, such as image-based high-content screenings, produce complex multivariate readouts. To fully exploit the richness of such data, advanced analytical methods that go beyond the conventional univariate approaches should be employed. In this work, we tackle the problem of hit identification in multivariate assays. As with univariate assays, a hit from a multivariate assay can be defined as a candidate that yields an assay value sufficiently far away in distance from the mean or central value of inactives. Viewed another way, a hit is an outlier from the distribution of inactives. A method was developed for identifying multivariate hit in high-dimensional data sets based on principal components and robust Mahalanobis distance (the multivariate analogue to the Z- or T-statistic). The proposed method, termed mROUT (multivariate robust outlier detection), demonstrates superior performance over other techniques in the literature in terms of maintaining Type I error, false discovery rate and true discovery rate in simulation studies. The performance of mROUT is also illustrated on a CRISPR knockout data set from in-house phenotypic screening programme.

Single-detector multiplex imaging flow cytometry for cancer cell classification with deep learning

Imaging flow cytometry, which combines the advantages of flow cytometry and microscopy, has emerged as a powerful tool for cell analysis in various biomedical fields such as cancer detection. In this study, we develop multiplex imaging flow cytometry (mIFC) by employing a spatial wavelength division multiplexing technique. Our mIFC can simultaneously obtain brightfield and multi-color fluorescence images of individual cells in flow, which are excited by a metal halide lamp and measured by a single detector. Statistical analysis results of multiplex imaging experiments with resolution test lens, magnification test lens, and fluorescent microspheres validate the operation of the mIFC with good imaging channel consistency and micron-scale differentiation capabilities. A deep learning method is designed for multiplex image processing that consists of three deep learning networks (U-net, very deep super resolution, and visual geometry group 19). It is demonstrated that the cluster of differentiation 24 (CD24) imaging channel is more sensitive than the brightfield, nucleus, or cancer antigen 125 (CA125) imaging channel in classifying the three types of ovarian cell lines (IOSE80 normal cell, A2780, and OVCAR3 cancer cells). An average accuracy rate of 97.1% is achieved for the classification of these three types of cells by deep learning analysis when all four imaging channels are considered. Our single-detector mIFC is promising for the development of future imaging flow cytometers and for the automatic single-cell analysis with deep learning in various biomedical fields.

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.

Information-Distilled Generative Label-Free Morphological Profiling Encodes Cellular Heterogeneity

Image-based cytometry faces challenges due to technical variations arising from different experimental batches and conditions, such as differences in instrument configurations or image acquisition protocols, impeding genuine biological interpretation of cell morphology. Existing solutions, often necessitating extensive pre-existing data knowledge or control samples across batches, have proved limited, especially with complex cell image data. To overcome this, "Cyto-Morphology Adversarial Distillation" (CytoMAD), a self-supervised multi-task learning strategy that distills biologically relevant cellular morphological information from batch variations, is introduced to enable integrated analysis across multiple data batches without complex data assumptions or extensive manual annotation. Unique to CytoMAD is its "morphology distillation", symbiotically paired with deep-learning image-contrast translation-offering additional interpretable insights into label-free cell morphology. The versatile efficacy of CytoMAD is demonstrated in augmenting the power of biophysical imaging cytometry. It allows integrated label-free classification of human lung cancer cell types and accurately recapitulates their progressive drug responses, even when trained without the drug concentration information. CytoMAD  also allows joint analysis of tumor biophysical cellular heterogeneity, linked to epithelial-mesenchymal plasticity, that standard fluorescence markers overlook. CytoMAD can substantiate the wide adoption of biophysical cytometry for cost-effective diagnosis and screening.

Artificial intelligence for high content imaging in drug discovery

Artificial intelligence (AI) and high-content imaging (HCI) are contributing to advancements in drug discovery, propelled by the recent progress in deep neural networks. This review highlights AI's role in analysis of HCI data from fixed and live-cell imaging, enabling novel label-free and multi-channel fluorescent screening methods, and improving compound profiling. HCI experiments are rapid and cost-effective, facilitating large data set accumulation for AI model training. However, the success of AI in drug discovery also depends on high-quality data, reproducible experiments, and robust validation to ensure model performance. Despite challenges like the need for annotated compounds and managing vast image data, AI's potential in phenotypic screening and drug profiling is significant. Future improvements in AI, including increased interpretability and integration of multiple modalities, are expected to solidify AI and HCI's role in drug discovery.

Deep learning large-scale drug discovery and repurposing

Large-scale drug discovery and repurposing is challenging. Identifying the mechanism of action (MOA) is crucial, yet current approaches are costly and low-throughput. Here we present an approach for MOA identification by profiling changes in mitochondrial phenotypes. By temporally imaging mitochondrial morphology and membrane potential, we established a pipeline for monitoring time-resolved mitochondrial images, resulting in a dataset comprising 570,096 single-cell images of cells exposed to 1,068 United States Food and Drug Administration-approved drugs. A deep learning model named MitoReID, using a re-identification (ReID) framework and an Inflated 3D ResNet backbone, was developed. It achieved 76.32% Rank-1 and 65.92% mean average precision on the testing set and successfully identified the MOAs for six untrained drugs on the basis of mitochondrial phenotype. Furthermore, MitoReID identified cyclooxygenase-2 inhibition as the MOA of the natural compound epicatechin in tea, which was successfully validated in vitro. Our approach thus provides an automated and cost-effective alternative for target identification that could accelerate large-scale drug discovery and repurposing.

Capturing cell heterogeneity in representations of cell populations for image-based profiling using contrastive learning

Image-based cell profiling is a powerful tool that compares perturbed cell populations by measuring thousands of single-cell features and summarizing them into profiles. Typically a sample is represented by averaging across cells, but this fails to capture the heterogeneity within cell populations. We introduce CytoSummaryNet: a Deep Sets-based approach that improves mechanism of action prediction by 30-68% in mean average precision compared to average profiling on a public dataset. CytoSummaryNet uses self-supervised contrastive learning in a multiple-instance learning framework, providing an easier-to-apply method for aggregating single-cell feature data than previously published strategies. Interpretability analysis suggests that the model achieves this improvement by downweighting small mitotic cells or those with debris and prioritizing large uncrowded cells. The approach requires only perturbation labels for training, which are readily available in all cell profiling datasets. CytoSummaryNet offers a straightforward post-processing step for single-cell profiles that can significantly boost retrieval performance on image-based profiling datasets.

The Autophagic Activator GHF-201 Can Alleviate Pathology in a Mouse Model and in Patient Fibroblasts of Type III Glycogenosis

Glycogen storage disease type III (GSDIII) is a hereditary glycogenosis caused by deficiency of the glycogen debranching enzyme (GDE), an enzyme, encoded by Agl, enabling glycogen degradation by catalyzing alpha-1,4-oligosaccharide side chain transfer and alpha-1,6-glucose cleavage. GDE deficiency causes accumulation of phosphorylase-limited dextrin, leading to liver disorder followed by fatal myopathy. Here, we tested the capacity of the new autophagosomal activator GHF-201 to alleviate disease burden by clearing pathogenic glycogen surcharge in the GSDIII mouse model Agl-/-. We used open field, grip strength, and rotarod tests for evaluating GHF-201's effects on locomotion, a biochemistry panel to quantify hematological biomarkers, indirect calorimetry to quantify in vivo metabolism, transmission electron microscopy to quantify glycogen in muscle, and fibroblast image analysis to determine cellular features affected by GHF-201. GHF-201 was able to improve all locomotion parameters and partially reversed hypoglycemia, hyperlipidemia and liver and muscle malfunction in Agl-/- mice. Treated mice burnt carbohydrates more efficiently and showed significant improvement of aberrant ultrastructural muscle features. In GSDIII patient fibroblasts, GHF-201 restored mitochondrial membrane polarization and corrected lysosomal swelling. In conclusion, GHF-201 is a viable candidate for treating GSDIII as it recovered a wide range of its pathologies in vivo, in vitro, and ex vivo.

A Saturation Artifacts Inpainting Method Based on Two-Stage GAN for Fluorescence Microscope Images

Fluorescence microscopic images of cells contain a large number of morphological features that are used as an unbiased source of quantitative information about cell status, through which researchers can extract quantitative information about cells and study the biological phenomena of cells through statistical and analytical analysis. As an important research object of phenotypic analysis, images have a great influence on the research results. Saturation artifacts present in the image result in a loss of grayscale information that does not reveal the true value of fluorescence intensity. From the perspective of data post-processing, we propose a two-stage cell image recovery model based on a generative adversarial network to solve the problem of phenotypic feature loss caused by saturation artifacts. The model is capable of restoring large areas of missing phenotypic features. In the experiment, we adopt the strategy of progressive restoration to improve the robustness of the training effect and add the contextual attention structure to enhance the stability of the restoration effect. We hope to use deep learning methods to mitigate the effects of saturation artifacts to reveal how chemical, genetic, and environmental factors affect cell state, providing an effective tool for studying the field of biological variability and improving image quality in analysis.

Pancreatic cancer organoid-screening captures personalized sensitivity and chemoresistance suppression upon cytochrome P450 3A5-targeted inhibition

Cytochrome P450 3A5 (CYP3A5) has been proposed as a predictor of therapy response in subtypes of pancreatic ductal adenocarcinoma cancer (PDAC). To validate CYP3A5 as a therapeutic target, we developed a high-content image organoid-based screen to quantify the phenotypic responses to the selective inhibition of CYP3A5 enzymatic activity by clobetasol propionate (CBZ), using a cohort of PDAC-derived organoids (PDACOs). The chemoresistance of PDACOs to a panel of standard-of-care drugs, alone or in combination with CBZ, was investigated. PDACO pharmaco-profiling revealed CBZ to have anti-cancer activity that was dependent on the CYP3A5 level. In addition, CBZ restored chemo-vulnerability to cisplatin in a subset of PDACOs. A correlative proteomic analysis established that CBZ caused the suppression of multiple cancer pathways sustained by or associated with a mutant form of p53. Limiting the active pool of CYP3A5 enables targeted and personalized therapy to suppress pro-oncogenic mechanisms that fuel chemoresistance in some PDAC tumors.

Unleashing the potential of cell painting assays for compound activities and hazards prediction

The cell painting (CP) assay has emerged as a potent imaging-based high-throughput phenotypic profiling (HTPP) tool that provides comprehensive input data for in silico prediction of compound activities and potential hazards in drug discovery and toxicology. CP enables the rapid, multiplexed investigation of various molecular mechanisms for thousands of compounds at the single-cell level. The resulting large volumes of image data provide great opportunities but also pose challenges to image and data analysis routines as well as property prediction models. This review addresses the integration of CP-based phenotypic data together with or in substitute of structural information from compounds into machine (ML) and deep learning (DL) models to predict compound activities for various human-relevant disease endpoints and to identify the underlying modes-of-action (MoA) while avoiding unnecessary animal testing. The successful application of CP in combination with powerful ML/DL models promises further advances in understanding compound responses of cells guiding therapeutic development and risk assessment. Therefore, this review highlights the importance of unlocking the potential of CP assays when combined with molecular fingerprints for compound evaluation and discusses the current challenges that are associated with this approach.

Distinct profile of antiviral drugs effects in aortic and pulmonary endothelial cells revealed by high-content microscopy and cell painting assays

Antiretroviral therapy have significantly improved the treatment of viral infections and reduced the associated mortality and morbidity rates. However, highly effective antiretroviral therapy (HAART) may lead to an increased risk of cardiovascular diseases, which could be related to endothelial toxicity. Here, seven antiviral drugs (remdesivir, PF-00835231, ritonavir, lopinavir, efavirenz, zidovudine and abacavir) were characterized against aortic (HAEC) and pulmonary (hLMVEC) endothelial cells, using high-content microscopy. The colourimetric study (MTS test) revealed similar toxicity profiles of all antiviral drugs tested in the concentration range of 1 nM-50 μM in aortic and pulmonary endothelial cells. Conversely, the drugs' effects on morphological parameters were more pronounced in HAECs as compared with hLMVECs. Based on the antiviral drugs' effects on the cytoplasmic and nuclei architecture (analyzed by multiple pre-defined parameters including SER texture and STAR morphology), the studied compounds were classified into five distinct morphological subgroups, each linked to a specific cellular response profile. In relation to morphological subgroup classification, antiviral drugs induced a loss of mitochondrial membrane potential, elevated ROS, changed lipid droplets/lysosomal content, decreased von Willebrand factor expression and micronuclei formation or dysregulated cellular autophagy. In conclusion, based on specific changes in endothelial cytoplasm, nuclei and subcellular morphology, the distinct endothelial response was identified for remdesivir, ritonavir, lopinavir, efavirenz, zidovudine and abacavir treatments. The effects detected in aortic endothelial cells were not detected in pulmonary endothelial cells. Taken together, high-content microscopy has proven to be a robust and informative method for endothelial drug profiling that may prove useful in predicting the organ-specific endothelial toxicity of various drugs.

Reproducible image-based profiling with Pycytominer

Advances in high-throughput microscopy have enabled the rapid acquisition of large numbers of high-content microscopy images. Whether by deep learning or classical algorithms, image analysis pipelines then produce single-cell features. To process these single-cells for downstream applications, we present Pycytominer, a user-friendly, open-source python package that implements the bioinformatics steps, known as "image-based profiling". We demonstrate Pycytominer's usefulness in a machine learning project to predict nuisance compounds that cause undesirable cell injuries.

Spatial Dissection of the Immune Landscape of Solid Tumors to Advance Precision Medicine

Chemotherapeutics, radiation, targeted therapeutics, and immunotherapeutics each demonstrate clinical benefits for a small subset of patients with solid malignancies. Immune cells infiltrating the tumor and the surrounding stroma play a critical role in shaping cancer progression and modulating therapy response. They do this by interacting with the other cellular and molecular components of the tumor microenvironment. Spatial multi-omics technologies are rapidly evolving. Currently, such technologies allow high-throughput RNA and protein profiling and retain geographical information about the tumor microenvironment cellular architecture and the functional phenotype of tumor, immune, and stromal cells. An in-depth spatial characterization of the heterogeneous tumor immune landscape can improve not only the prognosis but also the prediction of therapy response, directing cancer patients to more tailored and efficacious treatments. This review highlights recent advancements in spatial transcriptomics and proteomics profiling technologies and the ways these technologies are being applied for the dissection of the immune cell composition in solid malignancies in order to further both basic research in oncology and the implementation of precision treatments in the clinic.

Mesoscale standing wave imaging

Standing wave (SW) microscopy is a method that uses an interference pattern to excite fluorescence from labelled cellular structures and produces high-resolution images of three-dimensional objects in a two-dimensional dataset. SW microscopy is performed with high-magnification, high-numerical aperture objective lenses, and while this results in high-resolution images, the field of view is very small. Here we report upscaling of this interference imaging method from the microscale to the mesoscale using the Mesolens, which has the unusual combination of a low-magnification and high-numerical aperture. With this method, we produce SW images within a field of view of 4.4 mm × 3.0 mm that can readily accommodate over 16,000 cells in a single dataset. We demonstrate the method using both single-wavelength excitation and the multi-wavelength SW method TartanSW. We show application of the method for imaging of fixed and living cells specimens, with the first application of SW imaging to study cells under flow conditions.

Made to measure: An introduction to quantifying microscopy data in the life sciences

Images are at the core of most modern biological experiments and are used as a major source of quantitative information. Numerous algorithms are available to process images and make them more amenable to be measured. Yet the nature of the quantitative output that is useful for a given biological experiment is uniquely dependent upon the question being investigated. Here, we discuss the 3 main types of information that can be extracted from microscopy data: intensity, morphology, and object counts or categorical labels. For each, we describe where they come from, how they can be measured, and what may affect the relevance of these measurements in downstream data analysis. Acknowledging that what makes a measurement 'good' is ultimately down to the biological question being investigated, this review aims at providing readers with a toolkit to challenge how they quantify their own data and be critical of conclusions drawn from quantitative bioimage analysis experiments.

Connecting metabolome and phenotype: recent advances in functional metabolomics tools for the identification of bioactive natural products

Covering: 1995 to 2023Advances in bioanalytical methods, particularly mass spectrometry, have provided valuable molecular insights into the mechanisms of life. Non-targeted metabolomics aims to detect and (relatively) quantify all observable small molecules present in a biological system. By comparing small molecule abundances between different conditions or timepoints in a biological system, researchers can generate new hypotheses and begin to understand causes of observed phenotypes. Functional metabolomics aims to investigate the functional roles of metabolites at the scale of the metabolome. However, most functional metabolomics studies rely on indirect measurements and correlation analyses, which leads to ambiguity in the precise definition of functional metabolomics. In contrast, the field of natural products has a history of identifying the structures and bioactivities of primary and specialized metabolites. Here, we propose to expand and reframe functional metabolomics by integrating concepts from the fields of natural products and chemical biology. We highlight emerging functional metabolomics approaches that shift the focus from correlation to physical interactions, and we discuss how this allows researchers to uncover causal relationships between molecules and phenotypes.

Anomaly detection for high-content image-based phenotypic cell profiling

High-content image-based phenotypic profiling combines automated microscopy and analysis to identify phenotypic alterations in cell morphology and provide insight into the cell's physiological state. Classical representations of the phenotypic profile can not capture the full underlying complexity in cell organization, while recent weakly machine-learning based representation-learning methods are hard to biologically interpret. We used the abundance of control wells to learn the in-distribution of control experiments and use it to formulate a self-supervised reconstruction anomaly-based representation that encodes the intricate morphological inter-feature dependencies while preserving the representation interpretability. The performance of our anomaly-based representations was evaluated for downstream tasks with respect to two classical representations across four public Cell Painting datasets. Anomaly-based representations improved reproducibility, Mechanism of Action classification, and complemented classical representations. Unsupervised explainability of autoencoder-based anomalies identified specific inter-feature dependencies causing anomalies. The general concept of anomaly-based representations can be adapted to other applications in cell biology.

Three million images and morphological profiles of cells treated with matched chemical and genetic perturbations

The identification of genetic and chemical perturbations with similar impacts on cell morphology can elucidate compounds' mechanisms of action or novel regulators of genetic pathways. Research on methods for identifying such similarities has lagged due to a lack of carefully designed and well-annotated image sets of cells treated with chemical and genetic perturbations. Here we create such a Resource dataset, CPJUMP1, in which each perturbed gene's product is a known target of at least two chemical compounds in the dataset. We systematically explore the directionality of correlations among perturbations that target the same protein encoded by a given gene, and we find that identifying matches between chemical and genetic perturbations is a challenging task. Our dataset and baseline analyses provide a benchmark for evaluating methods that measure perturbation similarities and impact, and more generally, learn effective representations of cellular state from microscopy images. Such advancements would accelerate the applications of image-based profiling of cellular states, such as uncovering drug mode of action or probing functional genomics.

A spatial architecture-embedding HLA signature to predict clinical response to immunotherapy in renal cell carcinoma

An important challenge in the real-world management of patients with advanced clear-cell renal cell carcinoma (aRCC) is determining who might benefit from immune checkpoint blockade (ICB). Here we performed a comprehensive multiomics mapping of aRCC in the context of ICB treatment, involving discovery analyses in a real-world data cohort followed by validation in independent cohorts. We cross-connected bulk-tumor transcriptomes across >1,000 patients with validations at single-cell and spatial resolutions, revealing a patient-specific crosstalk between proinflammatory tumor-associated macrophages and (pre-)exhausted CD8+ T cells that was distinguished by a human leukocyte antigen repertoire with higher preference for tumoral neoantigens. A cross-omics machine learning pipeline helped derive a new tumor transcriptomic footprint of neoantigen-favoring human leukocyte antigen alleles. This machine learning signature correlated with positive outcome following ICB treatment in both real-world data and independent clinical cohorts. In experiments using the RENCA-tumor mouse model, CD40 agonism combined with PD1 blockade potentiated both proinflammatory tumor-associated macrophages and CD8+ T cells, thereby achieving maximal antitumor efficacy relative to other tested regimens. Thus, we present a new multiomics and spatial map of the immune-community architecture that drives ICB response in patients with aRCC.

Advancing materials science through next-generation machine learning

For over a decade, machine learning (ML) models have been making strides in computer vision and natural language processing (NLP), demonstrating high proficiency in specialized tasks. The emergence of large-scale language and generative image models, such as ChatGPT and Stable Diffusion, has significantly broadened the accessibility and application scope of these technologies. Traditional predictive models are typically constrained to mapping input data to numerical values or predefined categories, limiting their usefulness beyond their designated tasks. In contrast, contemporary models employ representation learning and generative modeling, enabling them to extract and encode key insights from a wide variety of data sources and decode them to create novel responses for desired goals. They can interpret queries phrased in natural language to deduce the intended output. In parallel, the application of ML techniques in materials science has advanced considerably, particularly in areas like inverse design, material prediction, and atomic modeling. Despite these advancements, the current models are overly specialized, hindering their potential to supplant established industrial processes. Materials science, therefore, necessitates the creation of a comprehensive, versatile model capable of interpreting human-readable inputs, intuiting a wide range of possible search directions, and delivering precise solutions. To realize such a model, the field must adopt cutting-edge representation, generative, and foundation model techniques tailored to materials science. A pivotal component in this endeavor is the establishment of an extensive, centralized dataset encompassing a broad spectrum of research topics. This dataset could be assembled by crowdsourcing global research contributions and developing models to extract data from existing literature and represent them in a homogenous format. A massive dataset can be used to train a central model that learns the underlying physics of the target areas, which can then be connected to a variety of specialized downstream tasks. Ultimately, the envisioned model would empower users to intuitively pose queries for a wide array of desired outcomes. It would facilitate the search for existing data that closely matches the sought-after solutions and leverage its understanding of physics and material-behavior relationships to innovate new solutions when pre-existing ones fall short.

Morphological profiling for drug discovery in the era of deep learning

Morphological profiling is a valuable tool in phenotypic drug discovery. The advent of high-throughput automated imaging has enabled the capturing of a wide range of morphological features of cells or organisms in response to perturbations at the single-cell resolution. Concurrently, significant advances in machine learning and deep learning, especially in computer vision, have led to substantial improvements in analyzing large-scale high-content images at high throughput. These efforts have facilitated understanding of compound mechanism of action, drug repurposing, characterization of cell morphodynamics under perturbation, and ultimately contributing to the development of novel therapeutics. In this review, we provide a comprehensive overview of the recent advances in the field of morphological profiling. We summarize the image profiling analysis workflow, survey a broad spectrum of analysis strategies encompassing feature engineering- and deep learning-based approaches, and introduce publicly available benchmark datasets. We place a particular emphasis on the application of deep learning in this pipeline, covering cell segmentation, image representation learning, and multimodal learning. Additionally, we illuminate the application of morphological profiling in phenotypic drug discovery and highlight potential challenges and opportunities in this field.

Nucleus-forming jumbophage PhiKZ therapeutically outcompetes non-nucleus-forming jumbophage Callisto

With the recent resurgence of phage therapy in modern medicine, jumbophages are currently under the spotlight due to their numerous advantages as anti-infective agents. However, most significant discoveries to date have primarily focused on nucleus-forming jumbophages, not their non-nucleus-forming counterparts. In this study, we compare the biological characteristics exhibited by two genetically diverse jumbophages: 1) the well-studied nucleus-forming jumbophage, PhiKZ; and 2) the newly discovered non-nucleus-forming jumbophage, Callisto. Single-cell infection studies further show that Callisto possesses different replication machinery, resulting in a delay in phage maturation compared to that of PhiKZ. The therapeutic potency of both phages was examined in vitro and in vivo, demonstrating that PhiKZ holds certain superior characteristics over Callisto. This research sheds light on the importance of the subcellular infection machinery and the organized progeny maturation process, which could potentially provide valuable insight in the future development of jumbophage-based therapeutics.

Genome-scale requirements for dynein-based transport revealed by a high-content arrayed CRISPR screen

The microtubule motor dynein plays a key role in cellular organization. However, little is known about how dynein's biosynthesis, assembly, and functional diversity are orchestrated. To address this issue, we have conducted an arrayed CRISPR loss-of-function screen in human cells using the distribution of dynein-tethered peroxisomes and early endosomes as readouts. From a genome-wide gRNA library, 195 validated hits were recovered and parsed into those impacting multiple dynein cargoes and those whose effects are restricted to a subset of cargoes. Clustering of high-dimensional phenotypic fingerprints revealed co-functional proteins involved in many cellular processes, including several candidate novel regulators of core dynein functions. Further analysis of one of these factors, the RNA-binding protein SUGP1, indicates that it promotes cargo trafficking by sustaining functional expression of the dynein activator LIS1. Our data represent a rich source of new hypotheses for investigating microtubule-based transport, as well as several other aspects of cellular organization captured by our high-content imaging.

Applying Cell Painting in Non-Tumorigenic Breast Cells to Understand Impacts of Common Chemical Exposures

There are a substantial number of chemicals to which individuals in the general population are exposed which have putative, but still poorly understood, links to breast cancer. Cell Painting is a high-content imaging-based in vitro assay that allows for rapid and unbiased measurements of the concentration-dependent effects of chemical exposures on cellular morphology. We optimized the Cell Painting assay and measured the effect of exposure to 16 human exposure relevant chemicals, along with 21 small molecules with known mechanisms of action, for 48 hours in non-tumorigenic mammary epithelial cells, the MCF10A cell line. Through unbiased imaging analyses using CellProfiler, we quantified 3042 morphological features across approximately 1.2 million cells. We used benchmark concentration modeling to quantify significance and dose-dependent directionality to identify morphological features conserved across chemicals and find features that differentiate the effects of toxicants from one another. Benchmark concentrations were compared to chemical exposure biomarker concentration measurements from the National Health and Nutrition Examination Survey to assess which chemicals induce morphological alterations at human-relevant concentrations. Morphometric fingerprint analysis revealed similar phenotypes between small molecules and prioritized NHANES-toxicants guiding further investigation. A comparison of feature fingerprints via hypergeometric analysis revealed significant feature overlaps between chemicals when stratified by compartment and stain. One such example was the similarities between a metabolite of the organochlorine pesticide DDT (p,p'-DDE) and an activator of canonical Wnt signaling CHIR99201. As CHIR99201 is a known Wnt pathway activator and its role in β-catenin translocation is well studied, we studied the translocation of β-catenin following p'-p' DDE exposure in an orthogonal high-content imaging assay. Consistent with activation of Wnt signaling, low dose p',p'-DDE (25nM) significantly enhances the nuclear translocation of β-catenin. Overall, these findings highlight the ability of Cell Painting to enhance mode-of-action studies for toxicants which are common exposures in our environment but have previously been incompletely characterized with respect to breast cancer risk.

Machine learning aided single cell image analysis improves understanding of morphometric heterogeneity of human mesenchymal stem cells

The multipotent stem cells of our body have been largely harnessed in biotherapeutics. However, as they are derived from multiple anatomical sources, from different tissues, human mesenchymal stem cells (hMSCs) are a heterogeneous population showing ambiguity in their in vitro behavior. Intra-clonal population heterogeneity has also been identified and pre-clinical mechanistic studies suggest that these cumulatively depreciate the therapeutic effects of hMSC transplantation. Although various biomarkers identify these specific stem cell populations, recent artificial intelligence-based methods have capitalized on the cellular morphologies of hMSCs, opening a new approach to understand their attributes. A robust and rapid platform is required to accommodate and eliminate the heterogeneity observed in the cell population, to standardize the quality of hMSC therapeutics globally. Here, we report our primary findings of morphological heterogeneity observed within and across two sources of hMSCs namely, stem cells from human exfoliated deciduous teeth (SHEDs) and human Wharton jelly mesenchymal stem cells (hWJ MSCs), using real-time single-cell images generated on immunophenotyping by imaging flow cytometry (IFC). We used the ImageJ software for identification and comparison between the two types of hMSCs using statistically significant morphometric descriptors that are biologically relevant. To expand on these insights, we have further applied deep learning methods and successfully report the development of a Convolutional Neural Network-based image classifier. In our research, we introduced a machine learning methodology to streamline the entire procedure, utilizing convolutional neural networks and transfer learning for binary classification, achieving an accuracy rate of 97.54%. We have also critically discussed the challenges, comparisons between solutions and future directions of machine learning in hMSC classification in biotherapeutics.

Image-based profiling and deep learning reveal morphological heterogeneity of colorectal cancer organoids

Patient-derived organoids have proven to be a highly relevant model for evaluating of disease mechanisms and drug efficacies, as they closely recapitulate in vivo physiology. Colorectal cancer organoids, specifically, exhibit a diverse range of morphologies, which have been analyzed with image-based profiling. However, the relationship between morphological subtypes and functional parameters of the organoids remains underexplored. Here, we identified two distinct morphological subtypes ("cystic" and "solid") across 31360 bright field images using image-based profiling, which correlated differently with viability and apoptosis level of colorectal cancer organoids. Leveraging object detection neural networks, we were able to categorize single organoids achieving higher viability scores as "cystic" than "solid" subtype. Furthermore, a deep generative model was proposed to predict apoptosis intensity based on a apoptosis-featured dataset encompassing over 17000 bright field and matched fluorescent images. Notably, a significant correlation of 0.91 between the predicted value and ground truth was achived, underscoring the feasibility of this generative model as a potential means for assessing organoid functional parameters. The underlying cellular heterogeneity of the organoids, i.e., conserved colonic cell types and rare immune components, was also verified with scRNA sequencing, implying a compromised tumor microenvironment. Additionally, the "cystic" subtype was identified as a relapse phenotype featuring intestinal stem cell signatures, suggesting that this visually discernible relapse phenotype shows potential as a novel biomarker for colorectal cancer diagnosis and prognosis. In summary, our findings demonstrate that the morphological heterogeneity of colorectal cancer organoids explicitly recapitulate the association of phenotypic features and exogenous perturbations through the image-based profiling, providing new insights into disease mechanisms.

Development of an Automated Morphometric Approach to Assess Vascular Outcomes following Exposure to Environmental Chemicals in Zebrafish

BACKGROUND

Disruptions in vascular formation attributable to chemical insults is a pivotal risk factor or potential etiology of developmental defects and various disease settings. Among the thousands of chemicals threatening human health, the highly concerning groups prevalent in the environment and detected in biological monitoring in the general population ought to be prioritized because of their high exposure risks. However, the impacts of a large number of environmental chemicals on vasculature are far from understood. The angioarchitecture complexity and technical limitations make it challenging to analyze the entire vasculature efficiently and identify subtle changes through a high-throughput in vivo assay.

OBJECTIVES

We aimed to develop an automated morphometric approach for the vascular profile and assess the vascular morphology of health-concerning environmental chemicals.

METHODS

High-resolution images of the entire vasculature in Tg(fli1a:eGFP) zebrafish were collected using a high-content imaging platform. We established a deep learning-based quantitative framework, ECA-ResXUnet, combined with MATLAB to segment the vascular networks and extract features. Vessel scores based on the rates of morphological changes were calculated to rank vascular toxicity. Potential biomarkers were identified by vessel-endothelium-gene-disease integrative analysis.

RESULTS

Whole-trunk blood vessels and the cerebral vasculature in larvae exposed to 150 representative chemicals were automatically segmented as comparable to human-level accuracy, with sensitivity and specificity of 95.56% and 95.81%, respectively. Chemical treatments led to heterogeneous vascular patterns manifested by 31 architecture indexes, and the common cardinal vein (CCV) was the most affected vessel. The antipsychotic medicine haloperidol, flame retardant 2,2-bis(chloromethyl)trimethylenebis[bis(2-chloroethyl) phosphate], and tert-butylphenyl diphenyl phosphate ranked as the top three in vessel scores. Pesticides accounted for the largest group, with a vessel score of ≥ 1 , characterized by a remarkable inhibition of subintestinal venous plexus and delayed development of CCV. Multiple-concentration evaluation of nine per- and polyfluoroalkyl substances (PFAS) indicated a low-concentration effect on vascular impairment and a positive association between carbon chain length and benchmark concentration. Target vessel-directed single-cell RNA sequencing of f l i 1 a + cells from larvae treated with λ -cyhalothrin , perfluorohexanesulfonic acid, or benzylbutyl phthalate, along with vessel-endothelium-gene-disease integrative analysis, uncovered potential associations with vascular disorders and identified biomarker candidates.

DISCUSSION

This study provides a novel paradigm for phenotype-driven screenings of vascular-disrupting chemicals by converging morphological and transcriptomic profiles at a high-resolution level, serving as a powerful tool for large-scale toxicity tests. Our approach and the high-quality morphometric data facilitate the precise evaluation of vascular effects caused by environmental chemicals. https://doi.org/10.1289/EHP13214.