Piximi - An Images to Discovery web tool for bioimages and beyond.

Deep learning has greatly accelerated research in biological image analysis yet it often requires programming skills and specialized tool installation. Here we present Piximi, a modern, no-programming image analysis tool leveraging deep learning. Implemented as a web application at Piximi.app, Piximi requires no installation and can be accessed by any modern web browser. Its client-only architecture preserves the security of researcher data by running all computation locally. Piximi offers four core modules: a deep learning classifier, an image annotator, measurement modules, and pre-trained deep learning segmentation modules. Piximi is interoperable with existing tools and workflows by supporting import and export of common data and model formats. The intuitive researcher interface and easy access to Piximi allows biological researchers to obtain insights into images within just a few minutes. Piximi aims to bring deep learning-powered image analysis to a broader community by eliminating barriers to entry.

Drug repurposing platform for deciphering the druggable SARS-CoV-2 interactome.

The coronavirus disease 2019 (COVID-19) pandemic has heavily challenged the global healthcare system. Despite the vaccination programs, the new virus variants are circulating. Further research is required for understanding of the biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and for discovery of therapeutic agents against the virus. Here, we took advantage of drug repurposing to identify if existing drugs could inhibit SARS-CoV-2 infection. We established an open high throughput platform for in vitro screening of drugs against SARS-CoV-2 infection. We screened ∼1000 drugs for their ability to inhibit SARS-CoV-2-induced cell death in the African green monkey kidney cell line (Vero-E6), analyzed how the hit compounds affect the viral N (nucleocapsid) protein expression in human cell lines using high-content microscopic imaging and analysis, determined the hit drug targets in silico, and assessed their ability to cause phospholipidosis, which can interfere with the viral replication. Duvelisib was found by in silico interaction assay as a potential drug targeting virus-host protein interactions. The predicted interaction between PARP1 and S protein, affected by Duvelisib, was further validated by immunoprecipitation. Our results represent a rapidly applicable platform for drug repurposing and evaluation of the new emerging viruses' responses to the drugs. Further in silico studies help us to discover the druggable host pathways involved in the infectious cycle of SARS-CoV-2.

Image-based and machine learning-guided multiplexed serology test for SARS-CoV-2.

We present a miniaturized immunofluorescence assay (mini-IFA) for measuring antibody response in patient blood samples. The method utilizes machine learning-guided image analysis and enables simultaneous measurement of immunoglobulin M (IgM), IgA, and IgG responses against different viral antigens in an automated and high-throughput manner. The assay relies on antigens expressed through transfection, enabling use at a low biosafety level and fast adaptation to emerging pathogens. Using severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the model pathogen, we demonstrate that this method allows differentiation between vaccine-induced and infection-induced antibody responses. Additionally, we established a dedicated web page for quantitative visualization of sample-specific results and their distribution, comparing them with controls and other samples. Our results provide a proof of concept for the approach, demonstrating fast and accurate measurement of antibody responses in a research setup with prospects for clinical diagnostics.

Phosphorylation of PACSIN2 at S313 Regulates Podocyte Architecture in Coordination with N-WASP.

Changes in the dynamic architecture of podocytes, the glomerular epithelial cells, lead to kidney dysfunction. Previous studies on protein kinase C and casein kinase 2 substrates in neurons 2 (PACSIN2), a known regulator of endocytosis and cytoskeletal organization, reveal a connection between PACSIN2 and kidney pathogenesis. Here, we show that the phosphorylation of PACSIN2 at serine 313 (S313) is increased in the glomeruli of rats with diabetic kidney disease. We found that phosphorylation at S313 is associated with kidney dysfunction and increased free fatty acids rather than with high glucose and diabetes alone. Phosphorylation of PACSIN2 emerged as a dynamic process that fine-tunes cell morphology and cytoskeletal arrangement, in cooperation with the regulator of the actin cytoskeleton, Neural Wiskott-Aldrich syndrome protein (N-WASP). PACSIN2 phosphorylation decreased N-WASP degradation while N-WASP inhibition triggered PACSIN2 phosphorylation at S313. Functionally, pS313-PACSIN2 regulated actin cytoskeleton rearrangement depending on the type of cell injury and the signaling pathways involved. Collectively, this study indicates that N-WASP induces phosphorylation of PACSIN2 at S313, which serves as a mechanism whereby cells regulate active actin-related processes. The dynamic phosphorylation of S313 is needed to regulate cytoskeletal reorganization.

Comparison of two supporting matrices for patient-derived cancer cells in 3D drug sensitivity and resistance testing assay (3D-DSRT).

Central to the success of functional precision medicine of solid tumors is to perform drug testing of patient-derived cancer cells (PDCs) in tumor-mimicking ex vivo conditions. While high throughput (HT) drug screening methods have been well-established for cells cultured in two-dimensional (2D) format, this approach may have limited value in predicting clinical responses. Here, we describe the results of the optimization of drug sensitivity and resistance testing (DSRT) in three-dimensional (3D) growth supporting matrices in a HT mode (3D-DSRT) using the hepatocyte cell line (HepG2) as an example. Supporting matrices included widely used animal-derived Matrigel and cellulose-based hydrogel, GrowDex, which has earlier been shown to support 3D growth of cell lines and stem cells. Further, the sensitivity of ovarian cancer PDCs, from two patients included in the functional precision medicine study, was tested for 52 drugs in 5 different concentrations using 3D-DSRT. Shortly, in the optimized protocol, the PDCs are embedded with matrices and seeded to 384-well plates to allow the formation of the spheroids prior to the addition of drugs in nanoliter volumes with acoustic dispenser. The sensitivity of spheroids to drug treatments is measured with cell viability readout (here, 72 h after addition of drugs). The quality control and data analysis are performed with openly available Breeze software. We show the usability of both matrices in established 3D-DSRT, and report 2D vs 3D growth condition dependent differences in sensitivities of ovarian cancer PDCs to MEK-inhibitors and cytotoxic drugs. This study provides a proof-of-concept for robust and fast screening of drug sensitivities of PDCs in 3D-DSRT, which is important not only for drug discovery but also for personalized ex vivo drug testing in functional precision medicine studies. These findings suggest that comparing results of 2D- and 3D-DSRT is essential for understanding drug mechanisms and for selecting the most effective treatment for the patient.

Fibroblast subsets in non-small cell lung cancer: Associations with survival, mutations, and immune features.

BACKGROUND: Cancer-associated fibroblasts (CAFs) are molecularly heterogeneous mesenchymal cells that interact with malignant cells and immune cells and confer anti- and protumorigenic functions. Prior in situ profiling studies of human CAFs have largely relied on scoring single markers, thus presenting a limited view of their molecular complexity. Our objective was to study the complex spatial tumor microenvironment of non-small cell lung cancer (NSCLC) with multiple CAF biomarkers, identify novel CAF subsets, and explore their associations with patient outcome. METHODS: Multiplex fluorescence immunohistochemistry was employed to spatially profile the CAF landscape in 2 population-based NSCLC cohorts (n = 636) using antibodies against 4 fibroblast markers: platelet-derived growth factor receptor-alpha (PDGFRA) and -beta (PDGFRB), fibroblast activation protein (FAP), and alpha-smooth muscle actin (αSMA). The CAF subsets were analyzed for their correlations with mutations, immune characteristics, and clinical variables as well as overall survival. RESULTS: Two CAF subsets, CAF7 (PDGFRA-/PDGFRB+/FAP+/αSMA+) and CAF13 (PDGFRA+/PDGFRB+/FAP-/αSMA+), showed statistically significant but opposite associations with tumor histology, driver mutations (tumor protein p53 [TP53] and epidermal growth factor receptor [EGFR]), immune features (programmed death-ligand 1 and CD163), and prognosis. In patients with early stage tumors (pathological tumor-node-metastasis IA-IB), CAF7 and CAF13 acted as independent prognostic factors. CONCLUSIONS: Multimarker-defined CAF subsets were identified through high-content spatial profiling. The robust associations of CAFs with driver mutations, immune features, and outcome suggest CAFs as essential factors in NSCLC progression and warrant further studies to explore their potential as biomarkers or therapeutic targets. This study also highlights multiplex fluorescence immunohistochemistry-based CAF profiling as a powerful tool for the discovery of clinically relevant CAF subsets.

Prognostic Role of Tumor Immune Microenvironment in Pleural Epithelioid Mesothelioma.

Background: Pleural mesothelioma (MPM) is an aggressive malignancy with an average patient survival of only 10 months. Interestingly, about 5%-10% of the patients survive remarkably longer. Prior studies have suggested that the tumor immune microenvironment (TIME) has potential prognostic value in MPM. We hypothesized that high-resolution single-cell spatial profiling of the TIME would make it possible to identify subpopulations of patients with long survival and identify immunophenotypes for the development of novel treatment strategies. Methods: We used multiplexed fluorescence immunohistochemistry (mfIHC) and cell-based image analysis to define spatial TIME immunophenotypes in 69 patients with epithelioid MPM (20 patients surviving ≥ 36 months). Five mfIHC panels (altogether 21 antibodies) were used to classify tumor-associated stromal cells and different immune cell populations. Prognostic associations were evaluated using univariate and multivariable Cox regression, as well as combination risk models with area under receiver operating characteristic curve (AUROC) analyses. Results: We observed that type M2 pro-tumorigenic macrophages (CD163+pSTAT1-HLA-DRA1-) were independently associated with shorter survival, whereas granzyme B+ cells and CD11c+ cells were independently associated with longer survival. CD11c+ cells were the only immunophenotype increasing the AUROC (from 0.67 to 0.84) when added to clinical factors (age, gender, clinical stage, and grade). Conclusion: High-resolution, deep profiling of TIME in MPM defined subgroups associated with both poor (M2 macrophages) and favorable (granzyme B/CD11c positivity) patient survival. CD11c positivity stood out as the most potential prognostic cell subtype adding prediction power to the clinical factors. These findings help to understand the critical determinants of TIME for risk and therapeutic stratification purposes in MPM.

Multiparametric platform for profiling lipid trafficking in human leukocytes.

Systematic insight into cellular dysfunction can improve understanding of disease etiology, risk assessment, and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validate this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrate lipid storage data into another readout parameter, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlates positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explains hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking and has broad possible applications in dissecting the cellular basis of metabolic disorders.

Nucleus segmentation: towards automated solutions.

Single nucleus segmentation is a frequent challenge of microscopy image processing, since it is the first step of many quantitative data analysis pipelines. The quality of tracking single cells, extracting features or classifying cellular phenotypes strongly depends on segmentation accuracy. Worldwide competitions have been held, aiming to improve segmentation, and recent years have definitely brought significant improvements: large annotated datasets are now freely available, several 2D segmentation strategies have been extended to 3D, and deep learning approaches have increased accuracy. However, even today, no generally accepted solution and benchmarking platform exist. We review the most recent single-cell segmentation tools, and provide an interactive method browser to select the most appropriate solution.

Spatial immunoprofiling of the intratumoral and peritumoral tissue of renal cell carcinoma patients.

While the abundance and phenotype of tumor-infiltrating lymphocytes are linked with clinical survival, their spatial coordination and its clinical significance remain unclear. Here, we investigated the immune profile of intratumoral and peritumoral tissue of clear cell renal cell carcinoma patients (n = 64). We trained a cell classifier to detect lymphocytes from hematoxylin and eosin stained tissue slides. Using unsupervised classification, patients were further classified into immune cold, hot and excluded topographies reflecting lymphocyte abundance and localization. The immune topography distribution was further validated with The Cancer Genome Atlas digital image dataset. We showed association between PBRM1 mutation and immune cold topography, STAG1 mutation and immune hot topography and BAP1 mutation and immune excluded topography. With quantitative multiplex immunohistochemistry we analyzed the expression of 23 lymphocyte markers in intratumoral and peritumoral tissue regions. To study spatial interactions, we developed an algorithm quantifying the proportion of adjacent immune cell pairs and their immunophenotypes. Immune excluded tumors were associated with superior overall survival (HR 0.19, p = 0.02) and less extensive metastasis. Intratumoral T cells were characterized with pronounced expression of immunological activation and exhaustion markers such as granzyme B, PD1, and LAG3. Immune cell interaction occurred most frequently in the intratumoral region and correlated with CD45RO expression. Moreover, high proportion of peritumoral CD45RO+ T cells predicted poor overall survival. In summary, intratumoral and peritumoral tissue regions represent distinct immunospatial profiles and are associated with clinicopathologic characteristics.

Regression plane concept for analysing continuous cellular processes with machine learning.

Biological processes are inherently continuous, and the chance of phenotypic discovery is significantly restricted by discretising them. Using multi-parametric active regression we introduce the Regression Plane (RP), a user-friendly discovery tool enabling class-free phenotypic supervised machine learning, to describe and explore biological data in a continuous manner. First, we compare traditional classification with regression in a simulated experimental setup. Second, we use our framework to identify genes involved in regulating triglyceride levels in human cells. Subsequently, we analyse a time-lapse dataset on mitosis to demonstrate that the proposed methodology is capable of modelling complex processes at infinite resolution. Finally, we show that hemocyte differentiation in Drosophila melanogaster has continuous characteristics.

BIAFLOWS: A Collaborative Framework to Reproducibly Deploy and Benchmark Bioimage Analysis Workflows.

Image analysis is key to extracting quantitative information from scientific microscopy images, but the methods involved are now often so refined that they can no longer be unambiguously described by written protocols. We introduce BIAFLOWS, an open-source web tool enabling to reproducibly deploy and benchmark bioimage analysis workflows coming from any software ecosystem. A curated instance of BIAFLOWS populated with 34 image analysis workflows and 15 microscopy image datasets recapitulating common bioimage analysis problems is available online. The workflows can be launched and assessed remotely by comparing their performance visually and according to standard benchmark metrics. We illustrated these features by comparing seven nuclei segmentation workflows, including deep-learning methods. BIAFLOWS enables to benchmark and share bioimage analysis workflows, hence safeguarding research results and promoting high-quality standards in image analysis. The platform is thoroughly documented and ready to gather annotated microscopy datasets and workflows contributed by the bioimaging community.

A functional genetic screen defines the AKT-induced senescence signaling network.

Exquisite regulation of PI3K/AKT/mTORC1 signaling is essential for homeostatic control of cell growth, proliferation, and survival. Aberrant activation of this signaling network is an early driver of many sporadic human cancers. Paradoxically, sustained hyperactivation of the PI3K/AKT/mTORC1 pathway in nontransformed cells results in cellular senescence, which is a tumor-suppressive mechanism that must be overcome to promote malignant transformation. While oncogene-induced senescence (OIS) driven by excessive RAS/ERK signaling has been well studied, little is known about the mechanisms underpinning the AKT-induced senescence (AIS) response. Here, we utilize a combination of transcriptome and metabolic profiling to identify key signatures required to maintain AIS. We also employ a whole protein-coding genome RNAi screen for AIS escape, validating a subset of novel mediators and demonstrating their preferential specificity for AIS as compared with OIS. As proof of concept of the potential to exploit the AIS network, we show that neurofibromin 1 (NF1) is upregulated during AIS and its ability to suppress RAS/ERK signaling facilitates AIS maintenance. Furthermore, depletion of NF1 enhances transformation of p53-mutant epithelial cells expressing activated AKT, while its overexpression blocks transformation by inducing a senescent-like phenotype. Together, our findings reveal novel mechanistic insights into the control of AIS and identify putative senescence regulators that can potentially be targeted, with implications for new therapeutic options to treat PI3K/AKT/mTORC1-driven cancers.

nucleAIzer: A Parameter-free Deep Learning Framework for Nucleus Segmentation Using Image Style Transfer.

Single-cell segmentation is typically a crucial task of image-based cellular analysis. We present nucleAIzer, a deep-learning approach aiming toward a truly general method for localizing 2D cell nuclei across a diverse range of assays and light microscopy modalities. We outperform the 739 methods submitted to the 2018 Data Science Bowl on images representing a variety of realistic conditions, some of which were not represented in the training data. The key to our approach is that during training nucleAIzer automatically adapts its nucleus-style model to unseen and unlabeled data using image style transfer to automatically generate augmented training samples. This allows the model to recognize nuclei in new and different experiments efficiently without requiring expert annotations, making deep learning for nucleus segmentation fairly simple and labor free for most biological light microscopy experiments. It can also be used online, integrated into CellProfiler and freely downloaded at www.nucleaizer.org. A record of this paper's transparent peer review process is included in the Supplemental Information.

Association of tamoxifen resistance and lipid reprogramming in breast cancer.

BACKGROUND: Tamoxifen treatment of estrogen receptor (ER)-positive breast cancer reduces mortality by 31%. However, over half of advanced ER-positive breast cancers are intrinsically resistant to tamoxifen and about 40% will acquire the resistance during the treatment. METHODS: In order to explore mechanisms underlying endocrine therapy resistance in breast cancer and to identify new therapeutic opportunities, we created tamoxifen-resistant breast cancer cell lines that represent the luminal A or the luminal B. Gene expression patterns revealed by RNA-sequencing in seven tamoxifen-resistant variants were compared with their isogenic parental cells. We further examined those transcriptomic alterations in a publicly available patient cohort. RESULTS: We show that tamoxifen resistance cannot simply be explained by altered expression of individual genes, common mechanism across all resistant variants, or the appearance of new fusion genes. Instead, the resistant cell lines shared altered gene expression patterns associated with cell cycle, protein modification and metabolism, especially with the cholesterol pathway. In the tamoxifen-resistant T-47D cell variants we observed a striking increase of neutral lipids in lipid droplets as well as an accumulation of free cholesterol in the lysosomes. Tamoxifen-resistant cells were also less prone to lysosomal membrane permeabilization (LMP) and not vulnerable to compounds targeting the lipid metabolism. However, the cells were sensitive to disulfiram, LCS-1, and dasatinib. CONCLUSION: Altogether, our findings highlight a major role of LMP prevention in tamoxifen resistance, and suggest novel drug vulnerabilities associated with this phenotype.

Intelligent image-based in situ single-cell isolation.

Quantifying heterogeneities within cell populations is important for many fields including cancer research and neurobiology; however, techniques to isolate individual cells are limited. Here, we describe a high-throughput, non-disruptive, and cost-effective isolation method that is capable of capturing individually targeted cells using widely available techniques. Using high-resolution microscopy, laser microcapture microscopy, image analysis, and machine learning, our technology enables scalable molecular genetic analysis of single cells, targetable by morphology or location within the sample.

Systems pathology by multiplexed immunohistochemistry and whole-slide digital image analysis.

The paradigm of molecular histopathology is shifting from a single-marker immunohistochemistry towards multiplexed detection of markers to better understand the complex pathological processes. However, there are no systems allowing multiplexed IHC (mIHC) with high-resolution whole-slide tissue imaging and analysis, yet providing feasible throughput for routine use. We present an mIHC platform combining fluorescent and chromogenic staining with automated whole-slide imaging and integrated whole-slide image analysis, enabling simultaneous detection of six protein markers and nuclei, and automatic quantification and classification of hundreds of thousands of cells in situ in formalin-fixed paraffin-embedded tissues. In the first proof-of-concept, we detected immune cells at cell-level resolution (n = 128,894 cells) in human prostate cancer, and analysed T cell subpopulations in different tumour compartments (epithelium vs. stroma). In the second proof-of-concept, we demonstrated an automatic classification of epithelial cell populations (n = 83,558) and glands (benign vs. cancer) in prostate cancer with simultaneous analysis of androgen receptor (AR) and alpha-methylacyl-CoA (AMACR) expression at cell-level resolution. We conclude that the open-source combination of 8-plex mIHC detection, whole-slide image acquisition and analysis provides a robust tool allowing quantitative, spatially resolved whole-slide tissue cytometry directly in formalin-fixed human tumour tissues for improved characterization of histology and the tumour microenvironment.

UNC-45a promotes myosin folding and stress fiber assembly.

Contractile actomyosin bundles, stress fibers, are crucial for adhesion, morphogenesis, and mechanosensing in nonmuscle cells. However, the mechanisms by which nonmuscle myosin II (NM-II) is recruited to those structures and assembled into functional bipolar filaments have remained elusive. We report that UNC-45a is a dynamic component of actin stress fibers and functions as a myosin chaperone in vivo. UNC-45a knockout cells display severe defects in stress fiber assembly and consequent abnormalities in cell morphogenesis, polarity, and migration. Experiments combining structured-illumination microscopy, gradient centrifugation, and proteasome inhibition approaches revealed that a large fraction of NM-II and myosin-1c molecules fail to fold in the absence of UNC-45a. The remaining properly folded NM-II molecules display defects in forming functional bipolar filaments. The C-terminal UNC-45/Cro1/She4p domain of UNC-45a is critical for NM-II folding, whereas the N-terminal tetratricopeptide repeat domain contributes to the assembly of functional stress fibers. Thus, UNC-45a promotes generation of contractile actomyosin bundles through synchronized NM-II folding and filament-assembly activities.

Data-analysis strategies for image-based cell profiling.

Image-based cell profiling is a high-throughput strategy for the quantification of phenotypic differences among a variety of cell populations. It paves the way to studying biological systems on a large scale by using chemical and genetic perturbations. The general workflow for this technology involves image acquisition with high-throughput microscopy systems and subsequent image processing and analysis. Here, we introduce the steps required to create high-quality image-based (i.e., morphological) profiles from a collection of microscopy images. We recommend techniques that have proven useful in each stage of the data analysis process, on the basis of the experience of 20 laboratories worldwide that are refining their image-based cell-profiling methodologies in pursuit of biological discovery. The recommended techniques cover alternatives that may suit various biological goals, experimental designs, and laboratories' preferences.

Advanced Cell Classifier: User-Friendly Machine-Learning-Based Software for Discovering Phenotypes in High-Content Imaging Data.

High-content, imaging-based screens now routinely generate data on a scale that precludes manual verification and interrogation. Software applying machine learning has become an essential tool to automate analysis, but these methods require annotated examples to learn from. Efficiently exploring large datasets to find relevant examples remains a challenging bottleneck. Here, we present Advanced Cell Classifier (ACC), a graphical software package for phenotypic analysis that addresses these difficulties. ACC applies machine-learning and image-analysis methods to high-content data generated by large-scale, cell-based experiments. It features methods to mine microscopic image data, discover new phenotypes, and improve recognition performance. We demonstrate that these features substantially expedite the training process, successfully uncover rare phenotypes, and improve the accuracy of the analysis. ACC is extensively documented, designed to be user-friendly for researchers without machine-learning expertise, and distributed as a free open-source tool at www.cellclassifier.org.