Yeast Proteome Dynamics from Single Cell Imaging and Automated Analysis.

Proteomics has proved invaluable in generating large-scale quantitative data; however, the development of systems approaches for examining the proteome in vivo has lagged behind. To evaluate protein abundance and localization on a proteome scale, we exploited the yeast GFP-fusion collection in a pipeline combining automated genetics, high-throughput microscopy, and computational feature analysis. We developed an ensemble of binary classifiers to generate localization data from single-cell measurements and constructed maps of ∼3,000 proteins connected to 16 localization classes. To survey proteome dynamics in response to different chemical and genetic stimuli, we measure proteome-wide abundance and localization and identified changes over time. We analyzed >20 million cells to identify dynamic proteins that redistribute among multiple localizations in hydroxyurea, rapamycin, and in an rpd3Δ background. Because our localization and abundance data are quantitative, they provide the opportunity for many types of comparative studies, single cell analyses, modeling, and prediction. VIDEO ABSTRACT.

Exploring the yeast acetylome using functional genomics.

Lysine acetylation is a dynamic posttranslational modification with a well-defined role in regulating histones. The impact of acetylation on other cellular functions remains relatively uncharacterized. We explored the budding yeast acetylome with a functional genomics approach, assessing the effects of gene overexpression in the absence of lysine deacetylases (KDACs). We generated a network of 463 synthetic dosage lethal (SDL) interactions involving class I and II KDACs, revealing many cellular pathways regulated by different KDACs. A biochemical survey of genes interacting with the KDAC RPD3 identified 72 proteins acetylated in vivo. In-depth analysis of one of these proteins, Swi4, revealed a role for acetylation in G1-specific gene expression. Acetylation of Swi4 regulates interaction with its partner Swi6, both components of the SBF transcription factor. This study expands our view of the yeast acetylome, demonstrates the utility of functional genomic screens for exploring enzymatic pathways, and provides functional information that can be mined for future studies.

Local statistics allow quantification of cell-to-cell variability from high-throughput microscope images.

MOTIVATION: Quantifying variability in protein expression is a major goal of systems biology and cell-to-cell variability in subcellular localization pattern has not been systematically quantified. RESULTS: We define a local measure to quantify cell-to-cell variability in high-throughput microscope images and show that it allows comparable measures of variability for proteins with diverse subcellular localizations. We systematically estimate cell-to-cell variability in the yeast GFP collection and identify examples of proteins that show cell-to-cell variability in their subcellular localization. CONCLUSIONS: Automated image analysis methods can be used to quantify cell-to-cell variability in microscope images.

Unsupervised clustering of subcellular protein expression patterns in high-throughput microscopy images reveals protein complexes and functional relationships between proteins.

Protein subcellular localization has been systematically characterized in budding yeast using fluorescently tagged proteins. Based on the fluorescence microscopy images, subcellular localization of many proteins can be classified automatically using supervised machine learning approaches that have been trained to recognize predefined image classes based on statistical features. Here, we present an unsupervised analysis of protein expression patterns in a set of high-resolution, high-throughput microscope images. Our analysis is based on 7 biologically interpretable features which are evaluated on automatically identified cells, and whose cell-stage dependency is captured by a continuous model for cell growth. We show that it is possible to identify most previously identified localization patterns in a cluster analysis based on these features and that similarities between the inferred expression patterns contain more information about protein function than can be explained by a previous manual categorization of subcellular localization. Furthermore, the inferred cell-stage associated to each fluorescence measurement allows us to visualize large groups of proteins entering the bud at specific stages of bud growth. These correspond to proteins localized to organelles, revealing that the organelles must be entering the bud in a stereotypical order. We also identify and organize a smaller group of proteins that show subtle differences in the way they move around the bud during growth. Our results suggest that biologically interpretable features based on explicit models of cell morphology will yield unprecedented power for pattern discovery in high-resolution, high-throughput microscopy images.

Cellular pathways regulating responses to compatible and self-incompatible pollen in Brassica and Arabidopsis stigmas intersect at Exo70A1, a putative component of the exocyst complex.

In the Brassicaceae, compatible pollen-pistil interactions result in pollen adhesion to the stigma, while pollen grains from unrelated plant species are largely ignored. There can also be an additional layer of recognition to prevent self-fertilization, the self-incompatibility response, whereby self pollen grains are distinguished from nonself pollen grains and rejected. This pathway is activated in the stigma and involves the ARM repeat-containing 1 (ARC1) protein, an E3 ubiquitin ligase. In a screen for ARC1-interacting proteins, we have identified Brassica napus Exo70A1, a putative component of the exocyst complex that is known to regulate polarized secretion. We show through transgenic studies that loss of Exo70A1 in Brassica and Arabidopsis thaliana stigmas leads to the rejection of compatible pollen at the same stage as the self-incompatibility response. A red fluorescent protein:Exo70A1 fusion rescues this stigmatic defect in Arabidopsis and is found to be mobilized to the plasma membrane concomitant with flowers opening. By contrast, increased expression of Exo70A1 in self-incompatible Brassica partially overcomes the self pollen rejection response. Thus, our data show that the Exo70A1 protein functions at the intersection of two cellular pathways, where it is required in the stigma for the acceptance of compatible pollen in both Brassica and Arabidopsis and is negatively regulated by Brassica self-incompatibility.

Proteome-wide screens in Saccharomyces cerevisiae using the yeast GFP collection.

The budding yeast is a simple and genetically tractable eukaryotic organism. It remains a leading system for functional genomic work and has been the focus of many pioneering efforts, including the systematic construction and analysis of gene deletion mutants. Over the past decade, many large-scale studies have made use of the deletion and other mutant collections to assay genetic interactions, chemical sensitivities, and other phenotypes, contributing enormously to our understanding of gene function. The deletion mutant collection has also been used in cell biological surveys to identify genes that control cell and organelle morphology. One valuable approach for systematic definition of gene function and biological pathways involves global assessment of the localization patterns of the proteins they encode and how these patterns are altered in response to environmental or genetic perturbation. However, proteome-wide, cell biological screens are extremely challenging, from both a technical and computational perspective. The yeast GFP collection, an elegant and unique strain set, is ideal for studying both protein localization and abundance across the proteome ( http://yeastgfp.yeastgenome.org/ ). In this chapter, we outline how the yeast GFP collection has been used to date and discuss approaches for conducting future surveys of the proteome.

Discovery of JNJ-63576253, a Next-Generation Androgen Receptor Antagonist Active Against Wild-Type and Clinically Relevant Ligand Binding Domain Mutations in Metastatic Castration-Resistant Prostate Cancer.

Numerous mechanisms of resistance arise in response to treatment with second-generation androgen receptor (AR) pathway inhibitors in metastatic castration-resistant prostate cancer (mCRPC). Among these, point mutations in the ligand binding domain can transform antagonists into agonists, driving the disease through activation of AR signaling. To address this unmet need, we report the discovery of JNJ-63576253, a next-generation AR pathway inhibitor that potently abrogates AR signaling in models of human prostate adenocarcinoma. JNJ-63576253 is advancing as a clinical candidate with potential effectiveness in the subset of patients who do not respond to or are progressing while on second-generation AR-targeted therapeutics.

Characterization of the Arabidopsis thaliana exocyst complex gene families by phylogenetic, expression profiling, and subcellular localization studies.

*The exocyst is a complex of eight proteins (Sec3p, Sec5p, Sec6p, Sec8p, Sec10p, Sec15p, Exo70p and Exo84p) involved in tethering vesicles to the plasma membrane during regulated or polarized secretion. Here, the plant exocyst complex was explored in phylogenetic, expression, and subcellular localization studies. *Evolutionary relationships of predicted exocyst subunits were examined in the complete genomes of Arabidopsis thaliana, Oryza sativa, Populus trichocarpa and Physcomitrella patens. Furthermore, detailed expression profiling of the A. thaliana microarray databases was performed and subcellular localization patterns were studied. *Several plant exocyst subunit genes appear to have undergone gene expansion in a common ancestor and subsequent duplication events in independent plant lineages. Expression profiling revealed that the A. thaliana Exo70 gene family exhibits dynamic expression patterns, while the remaining exocyst subunit genes displayed more static profiles. Subcellular localization patterns for A. thaliana exocyst subunits ranged from cytosolic to endosomal compartments (with enrichment in the early endosomes and the trans-Golgi network). Interestingly, two endosomal-localized AtExo70 proteins also recruited other exocyst subunits to these compartments. *Overall subcellular localization patterns were observed that were also found in yeast and animal cells, and this, coupled with the evolutionary relationships, suggests that the exocyst may perform similar conserved functions in plants.

Tales of 1,008 small molecules: phenomic profiling through live-cell imaging in a panel of reporter cell lines.

Phenomic profiles are high-dimensional sets of readouts that can comprehensively capture the biological impact of chemical and genetic perturbations in cellular assay systems. Phenomic profiling of compound libraries can be used for compound target identification or mechanism of action (MoA) prediction and other applications in drug discovery. To devise an economical set of phenomic profiling assays, we assembled a library of 1,008 approved drugs and well-characterized tool compounds manually annotated to 218 unique MoAs, and we profiled each compound at four concentrations in live-cell, high-content imaging screens against a panel of 15 reporter cell lines, which expressed a diverse set of fluorescent organelle and pathway markers in three distinct cell lineages. For 41 of 83 testable MoAs, phenomic profiles accurately ranked the reference compounds (AUC-ROC ≥ 0.9). MoAs could be better resolved by screening compounds at multiple concentrations than by including replicates at a single concentration. Screening additional cell lineages and fluorescent markers increased the number of distinguishable MoAs but this effect quickly plateaued. There remains a substantial number of MoAs that were hard to distinguish from others under the current study's conditions. We discuss ways to close this gap, which will inform the design of future phenomic profiling efforts.

PREDECT Protocols for Complex 2D/3D Cultures.

PREDECT, a European IMI consortium, has assumed the task to generate robust 2D and 3D culture platforms. Protocols established for 2D and 3D monoculture and stromal coculture models of increasing complexity (spheroid, stirred-tank bioreactor, Matrigel- and collagen-embedded cultures) have been established between six laboratories within academia, biotech, and pharma. These models were tested using three tumor cell lines (MCF7, LNCaP, and NCI-H1437), covering three pathologies (breast, prostate, and lung), but should be readily transferable to other model systems. Fluorescent protein tagged cell lines were used for all platforms, allowing for online measurement of growth curves and drug responses to treatments. All methods, from culture setup to phenotypic characterization and gene expression profiling are described in this chapter.The adaptable methodologies and detailed protocols described here should help to include these models more readily to the drug discovery pipeline.

Integrating images from multiple microscopy screens reveals diverse patterns of change in the subcellular localization of proteins.

The evaluation of protein localization changes on a systematic level is a powerful tool for understanding how cells respond to environmental, chemical, or genetic perturbations. To date, work in understanding these proteomic responses through high-throughput imaging has catalogued localization changes independently for each perturbation. To distinguish changes that are targeted responses to the specific perturbation or more generalized programs, we developed a scalable approach to visualize the localization behavior of proteins across multiple experiments as a quantitative pattern. By applying this approach to 24 experimental screens consisting of nearly 400,000 images, we differentiated specific responses from more generalized ones, discovered nuance in the localization behavior of stress-responsive proteins, and formed hypotheses by clustering proteins that have similar patterns. Previous approaches aim to capture all localization changes for a single screen as accurately as possible, whereas our work aims to integrate large amounts of imaging data to find unexpected new cell biology.

Repurposing High-Throughput Image Assays Enables Biological Activity Prediction for Drug Discovery.

In both academia and the pharmaceutical industry, large-scale assays for drug discovery are expensive and often impractical, particularly for the increasingly important physiologically relevant model systems that require primary cells, organoids, whole organisms, or expensive or rare reagents. We hypothesized that data from a single high-throughput imaging assay can be repurposed to predict the biological activity of compounds in other assays, even those targeting alternate pathways or biological processes. Indeed, quantitative information extracted from a three-channel microscopy-based screen for glucocorticoid receptor translocation was able to predict assay-specific biological activity in two ongoing drug discovery projects. In these projects, repurposing increased hit rates by 50- to 250-fold over that of the initial project assays while increasing the chemical structure diversity of the hits. Our results suggest that data from high-content screens are a rich source of information that can be used to predict and replace customized biological assays.

Liquid Growth of Arrayed Fluorescently Tagged Saccharomyces cerevisiae Strains for Live-Cell High-Throughput Microscopy Screens.

This protocol describes culturing arrays of fluorescently tagged yeast strains to early log-phase in a 96-well format for imaging on a high-throughput (HTP) microscope. The method assumes the use of the synthetic genetic array (SGA) technique to create the array of marked strains. When this approach is coupled with automated image analysis, the subcellular distribution and abundance of tagged proteins can be systematically and quantitatively examined in different genetic backgrounds and/or under different growth regimes.

Ellipsoid Segmentation Model for Analyzing Light-Attenuated 3D Confocal Image Stacks of Fluorescent Multi-Cellular Spheroids.

In oncology, two-dimensional in-vitro culture models are the standard test beds for the discovery and development of cancer treatments, but in the last decades, evidence emerged that such models have low predictive value for clinical efficacy. Therefore they are increasingly complemented by more physiologically relevant 3D models, such as spheroid micro-tumor cultures. If suitable fluorescent labels are applied, confocal 3D image stacks can characterize the structure of such volumetric cultures and, for example, cell proliferation. However, several issues hamper accurate analysis. In particular, signal attenuation within the tissue of the spheroids prevents the acquisition of a complete image for spheroids over 100 micrometers in diameter. And quantitative analysis of large 3D image data sets is challenging, creating a need for methods which can be applied to large-scale experiments and account for impeding factors. We present a robust, computationally inexpensive 2.5D method for the segmentation of spheroid cultures and for counting proliferating cells within them. The spheroids are assumed to be approximately ellipsoid in shape. They are identified from information present in the Maximum Intensity Projection (MIP) and the corresponding height view, also known as Z-buffer. It alerts the user when potential bias-introducing factors cannot be compensated for and includes a compensation for signal attenuation.

CYCLoPs: A Comprehensive Database Constructed from Automated Analysis of Protein Abundance and Subcellular Localization Patterns in Saccharomyces cerevisiae.

Changes in protein subcellular localization and abundance are central to biological regulation in eukaryotic cells. Quantitative measures of protein dynamics in vivo are therefore highly useful for elucidating specific regulatory pathways. Using a combinatorial approach of yeast synthetic genetic array technology, high-content screening, and machine learning classifiers, we developed an automated platform to characterize protein localization and abundance patterns from images of log phase cells from the open-reading frame-green fluorescent protein collection in the budding yeast, Saccharomyces cerevisiae. For each protein, we produced quantitative profiles of localization scores for 16 subcellular compartments at single-cell resolution to trace proteome-wide relocalization in conditions over time. We generated a collection of ∼300,000 micrographs, comprising more than 20 million cells and ∼9 billion quantitative measurements. The images depict the localization and abundance dynamics of more than 4000 proteins under two chemical treatments and in a selected mutant background. Here, we describe CYCLoPs (Collection of Yeast Cells Localization Patterns), a web database resource that provides a central platform for housing and analyzing our yeast proteome dynamics datasets at the single cell level. CYCLoPs version 1.0 is available at http://cyclops.ccbr.utoronto.ca. CYCLoPs will provide a valuable resource for the yeast and eukaryotic cell biology communities and will be updated as new experiments become available.

A quantitative literature-curated gold standard for kinase-substrate pairs.

We describe the Yeast Kinase Interaction Database (KID, http://www.moseslab.csb.utoronto.ca/KID/), which contains high- and low-throughput data relevant to phosphorylation events. KID includes 6,225 low-throughput and 21,990 high-throughput interactions, from greater than 35,000 experiments. By quantitatively integrating these data, we identified 517 high-confidence kinase-substrate pairs that we consider a gold standard. We show that this gold standard can be used to assess published high-throughput datasets, suggesting that it will enable similar rigorous assessments in the future.

A novel mechanism for SUMO system control: regulated Ulp1 nucleolar sequestration.

The small ubiquitin-related modifiers (SUMOs) are evolutionarily conserved polypeptides that are covalently conjugated to protein targets to modulate their subcellular localization, half-life, or activity. Steady-state SUMO conjugation levels increase in response to many different types of environmental stresses, but how the SUMO system is regulated in response to these insults is not well understood. Here, we characterize a novel mode of SUMO system control: in response to elevated alcohol levels, the Saccharomyces cerevisiae SUMO protease Ulp1 is disengaged from its usual location at the nuclear pore complex (NPC) and sequestered in the nucleolus. We further show that the Ulp1 region previously demonstrated to interact with the karyopherins Kap95 and Kap60 (amino acids 150 to 340) is necessary and sufficient for nucleolar targeting and that enforced sequestration of Ulp1 in the nucleolus significantly increases steady-state SUMO conjugate levels, even in the absence of alcohol. We have thus characterized a novel mechanism of SUMO system control in which the balance between SUMO-conjugating and -deconjugating activities at the NPC is altered in response to stress via relocalization of a SUMO-deconjugating enzyme.