|Home | About | Journals | Submit | Contact Us | Français|
High-throughput screening (HTS) of chemical libraries has become a critical tool in basic biology and drug discovery. However, its implementation and the adaptation of high content assays to human embryonic stem cells (hESCs) have been hampered by multiple technical challenges. Here we present a strategy to adapt hESCs to HTS conditions, resulting in an assay suitable for the discovery of small molecules that drive hESC self-renewal or differentiation. Use of this new assay has led to the identification of several marketed drugs and natural compounds promoting short-term hESC maintenance and compounds directing early lineage choice during differentiation. Global gene expression analysis upon drug treatment defines known and novel pathways correlated to hESC self-renewal and differentiation. Our results demonstrate feasibility of hESC-based HTS and enhance the repertoire of chemical compounds for manipulating hESC fate. The availability of high content assays should accelerate progress in basic and translational hESC biology.
The availability of libraries of chemical compounds is a critical tool in drug discovery. High-throughput chemical screens have also become an important strategy to address questions in basic biology including stem cell research (Ding and Schultz, 2004). Recent studies have identified novel small molecules that mediate cell fate acquisition and the differentiation behavior of a number of stem cell types. Examples include: self-renewal in mouse ES cells via dual inhibition of RasGAP and ERK1(Chen et al., 2006), neural (Ding et al., 2003) or cardiogenic (Wu et al., 2004a) specification of mouse ESCs, neurogenic response of neural stem/precursor cells by an orphan ligand, P-Ser (Saxe et al., 2007), or osteogenic differentiation in mesenchymal stem in response to purmorphamine (Wu et al., 2004b). Results from primary screens in mouse cells may often be applicable to similar cell populations of human origin. However, in the case of human ESCs this may not be the case given the distinct extrinsic signaling requirements to retain pluripotency in mouse (Ying et al., 2003) versus human (Vallier et al., 2005) ESCs.
Human ESCs (hESCs) self-renew extensively and give rise to various specialized cells relevant for modeling human development and for applications in regenerative medicine. While key transcription factors important for hESC renewal such as Oct4, Nanog, and Sox2 have been well characterized (Boyer et al., 2005;Ginis et al., 2004;Sato et al., 2003), extrinsic factors regulating hESC maintenance and early differentiation events are poorly understood (Hoffman and Carpenter, 2005). The development of high content assays in hESCs has been challenging due to the difficulties in establishing suitable growth and plating conditions. Here we report the first HTS assay in hESCs and identify molecules that promote short-term self-renewal or drive early differentiation. Our data pave the way for the routine application of HTS assays in hESC biology and expand the repertoire of compounds available for manipulating hESC fate.
A critical step in the development of a robust hESC based assay is the reproducible behavior of large scale hESC cultures in a multi-well format. Prior to seeding in 384 well plates, hESCs were maintained under feeder-free conditions (Xu et al., 2001) and assessed for the expression of pluripotency markers (Fig. 1A, Supplementary Fig. 1). To assure uniform cell distribution we developed a single cell dissociation protocol for hESCs by testing various enzymatic treatments and plating strategies (Supplementary Fig. 2). Accutase treatment resulted in fully dissociated cells suitable for plating on matrigel coated 384 well plates. Under these conditions, we observed a robust cell density-dependent growth response as measured by Alamar Blue at day 7 (Antczak et al., 2007;O’Brien et al., 2000) (Fig. 1B). Plating of 6000 cells/well led to a density of hESCs suitable for automated microscopy and image analysis at 7 days after plating. At the time of analysis nuclei were separated from each other and cells were distributed evenly over the area of each well without forming multilayered structures (Fig. 1C–E). The density was also suitable to withstand repeated wash cycles required for our immunostaining protocol at the time of analysis.
Increase and decrease in the expression of the pluripotency factor Oct4 was used as primary readout. FGF2 supplementation was withdrawn from hESC cultures at day 2 after plating. Continued exposure to FGF2 until day7 resulted in high Oct4 levels defined as high control (promoting self-renewal) while FGF2 withdrawal resulted in a sensitized state amenable to readout both self-renewal and differentiation compounds. Treatment with BMP4 induced a dramatic decrease in Oct4 expression, which was defined as low control (induction of differentiation) in our assay (Fig. 1C). HESCs were exposed to the chemicals from day 2 – 7 (Fig. 1D). Growth in conditioned medium for 48 hours prior to exposure to compounds provided sufficient time for dissociated hESCs to recover from the dissociation procedure and resume growth. The medium change at 48 hours after plating allowed removal of non-attached and dead cells. Compounds were diluted in KSR medium at a final concentration of 10μM and added to the cells. Cells were assayed at day 7 (5 days/120 hrs after exposure to the compounds) via high throughput laser scanning confocal microscopy (GE InCell Analyzer 3000) following automated immunocytochemistry. Oct4 levels were normalized to nuclear Hoechst 33342 staining in order to control for effects of the compounds on total cell number (Fig. 1D) and z’ values were calculated to establish statistical robustness of the assay (Supplementary Fig. 3). To probe the validity of the 384 well assay we tested a number of additional compounds with known activities in hESCs (Fig. 1E) including IGF2 (Bendall et al., 2007), BMP2 (Pera et al., 2004), and PD98059 (Li et al., 2007). While our primary screen was carried out in H9 cells (WA-09), the conditions for single cell dissociation, 384 well plating, and validations for high and low controls were fully reproduced in an independent hESC line (H1, WA-01; Supplementary Fig. 4).
Using this assay, we screened a library of 2,880 small molecules in multiple replicates. The library contained a broad set of biologically active and structurally diverse compounds including a set of 748 currently marketed drugs (for details see Experimental Procedures and (Antczak et al., 2007)). Compounds were screened at a final concentration of 10μM in a volume of 50μl/well containing 1% DMSO (v/v). Normalized Oct4 values for each of the replicate 384 well plates were ranked and analyzed for repeated low or high hits. Out of the 89 hits (~ 3% of total molecules screened), 65 compounds were selected for subsequent dose response assays based on their repeated high or low values in Oct4 expression intensity. A total of 22 compounds exhibited dose-dependent effects on normalized Oct4 values. Exposure to four of the compounds resulted in increased Oct4 intensity (activators) while 10 compounds caused a decrease in Oct4 intensity (inhibitors) (Fig. 2A). The remaining eight compounds showed dose-dependent toxicity in our assay without inducing differentiation (cytotoxic). A list of all 22 compounds is provided in Supplementary Table 1. Representative images for each of the three responses (self-renewal, differentiation, cell death) are shown in Fig. 2A, lower panel. Dose response curves demonstrated that several activators induced Oct4 levels comparable or higher than those achieved by FGF2 treatment. Similarly, several inhibitors led to a decrease in Oct4 levels comparable to BMP4 treatment (Fig. 2B). Several classes of drugs were represented multiple times within the list of 22 confirmed hits such as antiinflammatory drugs, antibiotics and cardiac glycosides (Fig. 2C). A significant overrepresentation of compounds in cardiac glycoside and antiinflammatory compound categories was observed compared with overall library composition: Cardiac glycosides: 18% of confirmed hits; < 0.1% of library, p < 0.000001; antiinflammatory compounds: 14% of confirmed hits; ~ 4% of library, p < 0.05; antibiotics: 18% of confirmed hits, < 10% of library, p not significant. Chemical identity was confirmed by mass spectrometry (data not shown).
Activator compounds were subjected to a number of biological validations. Effects on Oct4 were confirmed using Nanog expression as an alternative pluripotency marker and using H1 as an additional hESC line. All four activator compounds, Theanine (THEA), Sinomenine (SNM), Gatifloxacin (GTFX), and Flurbiprofen (FBP) induced expression of Nanog at levels comparable to Oct4 in both H9 and H1 (Fig. 3A). Maintenance of hESCs on matrigel in the presence of these activators was tested in a basic hESC medium (Perrier et al., 2004) that lacked FGF2 supplementation and conditioning in mouse embryonic fibroblasts (MEFs). The number of undifferentiated colonies was quantified after multiple passages in vitro. We observed that each of the four activator compounds yielded a strong increase in the percentage of undifferentiated colonies compared to cells maintained in hESC medium alone. The efficiency of these drugs in maintaining colonies undifferentiated was comparable to FGF2 exposure though it did not reach the levels achieved by combined MEF conditioning and FGF2 supplementation of the medium (Fig. 3B).
While key transcriptional regulators of pluripotency are shared between human and mouse ES cells, extrinsic conditions required to maintain an undifferentiated state are quite distinct. Here we tested whether the four activator compounds identified in our primary hESC based screen are capable of promoting self-renewal in mouse ES cells. To this end, an Oct4::eGFP line was generated by BAC transgenesis (Supplementary Fig. 5) following our previously published protocols (Tomishima et al., 2007). Using FACS analysis, we observed a dose-dependent effect of GTFX on Oct4 expression which was comparable to the effect of LIF treatment. Treatment with FBP or SNM resulted in a dose-dependent partial effect while supplementation with THEA did not increase Oct4 expression in mouse ESCs (Fig. 3C). We next assessed Oct4 expression in hESC cultures following initial exposure and subsequent passaging of cells under standard 6 well culture conditions. Oct4/DAPI ratio was determined following immunocytochemistry and quantitative imaging (ImageJ, for details see Experimental Procedures). These data confirmed a near complete loss of Oct4 in BMP4 treated cultures and maintenance of Oct4 expression in control, compound and FGF2 treated cultures (Fig. 3D,E). To further explore the impact of compound treatment on hESC growth, we performed cell cycle analysis via Flow Cytometry using cell loading with propidium iodide-citrate solution (Fig. 3F; for details see Experimental Procedures; Supplementary Fig. 6). The cell cycle data showed that FGF2 treatment results in a high percentage of cells in S- or G2/M phase indicating a high degree of cell division at both passages (1 and 2). In contrast, exposure to BMP4 led to a marked increase in the proportion of cells in G1/G0 compatible with an increase in cell differentiation at passage 2. Drug-treated cells showed a high proportion of cells in S- and G2/M phase, similar to IGF treatment. Cultures maintained under control conditions () yielded an increased proportion of cells at G1/G0 in passage 2.
A next set of studies addressed apoptotic cell death and maintenance of pluripotency in compound-treated hESC cultures. We observed low overall percentages of TUNEL+ cells at day 4 upon compound exposure without significant differences among the various treatment groups (Fig. 4A). Directed differentiation assays towards neurectoderm (Perrier et al., 2004), endoderm (Watanabe et al., 2007), and mesoderm (Laflamme et al., 2007) fate showed that differentiation potential towards all three germ-layers was retained (Fig. 4B). However, long-term maintenance revealed that the drugs were not capable of reproducing the full effect of conditioned medium. Spontaneous differentiation under feeder-free conditions upon long-term passaging was eventually observed in all conditions, even when combining multiple compounds (data not shown).
Global mRNA expression analysis identified a total of > 1000 transcripts (Supplementary Table 2; ) significantly regulated after drug treatment compared to basal hESC medium (without FGF2 or MEF conditioning). Gene expression analysis was performed immediately upon exposure to drugs (24 hours). At that stage of culture we did not observe any signs of differentiation or cell heterogeneity. A heat map of the changes showed co-regulation for many of the mRNAs induced by drug treatment with those induced by conditioned medium or conditioned medium supplemented with FGF2 (Fig. 4C). Both SNM and GTFX showed highly significant overlaps in gene expression with conditioned medium treatment (p < 10−9) and with combined conditioned medium + FGF2 treatment (p < 10−16) (Figs. 4D,E). Functional annotation of the most differentially regulated transcripts using DAVID http://david.abcc.ncifcrf.gov (Dennis, Jr. et al., 2003) confirmed strong regulation of TGFβ and Wnt signaling pathways in cultures treated with conditioned medium (CM). However, despite significant overlap in overall analysis (Fig. 4D, E), the top most differentially regulated categories in GTFX or SNM treated cultures were distinct (Fig. 4F, Supplementary Table 3). Regulation of heparin-binding in GTFX-treated cultures could potentiate endogenous FGF signaling in hESC cultures. Fibronectin-1 which was significantly regulated in SNM treated cultures has been a key component in a novel serum- and feeder-free culture system proposed for hESC maintenance (Amit and Itskovitz-Eldor, 2006). However, futures studies will have to directly address the contribution of these individual candidate pathways in mediating the self-renewal promoting response of each compound.
Interestingly, a chemical derivative of GTFX, Orbifloxacin, had comparable effects on Oct4 maintenance while another closely related compound, Ciprofloxacin, lacked Oct4 inducing activities and caused differentiation towards mesendodermal fates (Supplementary Fig. 7). Similarly, we identified a derivative of FBP, Felbinac, capable of promoting maintenance of Oct4 (Supplementary Fig. 7). The availability of related compounds with distinct action on hESC differentiation will provide an entry point towards structure-function studies and an understanding of mechanism of action. These data demonstrate that our HTS assay is suitable for the identification of compounds promoting hESC self-renewal, and that several of these compounds fall into the category of currently marketed drugs (GTFX, SNM, and FBP) or natural compounds such as THEA found in black tea. The fact that none of the activators was sufficient to induce long-term self-renewal indicates the need for screening larger libraries and further modifications in screening design. One promising strategy may be the use of larger primary screens followed by stringent secondary assays measuring long-term hESC self-renewal more directly.
We next tested the function of a selected set of 4 out of 10 inhibitor compounds including Tretinoin (all-trans retinoic acid, RA), Selegiline (SLG), Cymarin (CYM), and Sarmentogenin (SRM) that induced a dose-dependent reduction in Oct4 levels (Fig. 5A). We first validated the effect of these drugs on loss of pluripotency using Nanog as an alternative marker in WA-09 and confirmed loss of both Oct4 and Nanog in H1 (Fig. 5A). We hypothesized that loss of Oct4 expression in our high content assay can be used as a surrogate marker for the identification of compounds that induce early differentiation events in hESCs. Key lineage choices accessible to undifferentiated hESCs include extraembryonic fates such as trophoblast (Xu et al., 2002) or primitive endoderm differentiation (Pera et al., 2004), mesendodermal (Barberi et al., 2007;D’Amour et al., 2005), and ectodermal fates (Zhang et al., 2001). A selected repertoire of these fate choices were assessed at the time of Oct4 analysis by automated immunocytochemistry for Cytokeratin 18 (CK18, trophoblast marker), Sox17 (mesendoderm/endoderm) and Pax6 (neuroectoderm). We observed that exposure to Tretinoin induced CK18 at levels comparable to BMP4 treatment. Exposure to the cardiac glycosides (CYM or SARM) caused an increase in Sox17 expression (Fig. 5B). Quantification of the percentage of Sox17+ cells (Fig. 5C, D) was compatible with the Sox17 signal intensity data in Fig. 5B. Both CYM and SARM yielded an increase in the percentage of Sox17+ cells comparable to the effect observed upon ActivinA exposure. While SLG did not show any significant increase in the percentage of cells expressing differentiation markers at day 7 of differentiation, an increased number of Pax6+ clusters in longer term cultures is compatible with recent observations in mouse ES cells (Esmaeili et al., 2006).
Global gene expression analysis 24 hours after drug exposure (Fig. 6A) confirmed a large number of differentially expressed transcripts (Supplementary Table 4) that are significantly co-regulated in CYM and SRM treated cultures (Fig. 6B; p < 10−241). A large number of these transcripts were compatible with early differentiation events towards mesendodermal lineage choice. In contrast, exposure to Tretinoin (RA) treatment led to an upregulation of genes compatible with extraembryonic and neuroectodermal differentiation (Fig. 6C, Supplementary Table 4) while transcripts induced by SLG treatment were most compatible with neurectodermal and neural crest differentiation (Supplementary Fig. 8). Functional annotation analysis (DAVID: http://david.abcc.ncifcrf.gov, (Dennis, Jr. et al., 2003)) in RA treated cultures demonstrated enrichment in categories known to be associated with retinoid signaling including regulation of Hox genes, modulation of Id2 repressor and other bHLH related factors (Fig. 6D). SLG exposure showed significant enrichment in transcripts associated with signaling and cell differentiation. SRM and CYM exposure showed enrichment in the regulation of highly overlapping sets of functional category terms including lipoprotein, basic leucin zipper, and various categories associated with development and differentiation (Fig. 6D, Supplementary Table 5).
Our data demonstrate the feasibility of performing a primary HTS assay in hESCs and pave the way for establishing high content screening platforms ranging from assays that address questions in basic hESC biology to translational applications and drug discovery. The design of our primary HTS assay is suited for the screening of large compound libraries. The significant costs of such an effort may be offset by the identification of compounds promoting self-renewal as well as compounds driving differentiation. Putative differentiation compounds identified in the primary screen can be subsequently classified into compound classes promoting various lineage specific cell fates using secondary screens. However, it will be critical in the future to fully validate all secondary screens using z’ analysis up to the standards of a primary screen. Our current assay requires high initial plating densities to compensate for cell loss at the time of plating and to assure even cell distribution throughout the surface area of each well. At this cell density we did not observe stacking of multiple nuclei that would prevent reliable image analysis. One disadvantage of using such high cell densities is the relatively large number of total cells needed for a full scale HTS screen. For example screening a library of 1 million compounds without replicates but including internal high and low controls will require ~ 6 × 109 cells. Future studies should allow primary HTS assays in hESCs at lower cell densities by reducing initial cell death at plating. One interesting compound to be tested is the ROCK inhibitor Y- 27632 shown to increase efficiency of plating dissociated hESCs at low cell densities (Watanabe et al., 2007).
There is currently a very limited set of chemical compounds available that affect hESC behavior (Chen et al., 2006;Ding et al., 2003;Sato et al., 2004), and none of these compounds have been identified through primary HTS assays in hESCs. The newly discovered compounds enhance our repertoire of chemical modulators suitable for directing growth and differentiation of hESCs. Future studies will have to explore the mechanistic action of these compounds on the self-renewal and differentiation machinery of hESCs. Our gene expression data, functional annotation analysis, and data from the literature on drug action of several of these compounds will provide valuable entry points to this end. One issue to be addressed is whether any given compound works directly on undifferentiated hESCs or whether such an effect could be due to a bystander effect mediated by hESC derived differentiated progeny. Such a bystander effect could feedback on undifferentiated hESCs as proposed in the case of spontaneously emerging hESC derived fibroblast (Bendall et al., 2007). The lack of true long-term self-renewal factors among the activator compounds identified in the current study suggests that such molecules may be rare and require the screening of much larger libraries. Another critical question is whether the induction of Sox17 observed upon exposure to CYM and SRM can be translated into a protocol for generating endodermal derivatives from hESCs. While the percentage of Sox17+ cells was increased at levels comparable or higher than those achieved by ActivinA, CYM and SRM had a negative effect on overall cell growth. Future studies will have to address this issue in more detail by systematically mapping the expression of endodermal genes and cell fates in the presence of these compounds.
Several of the cytotoxic factors identified here (Fig. 2A, Supplementary Table 1) could be further explored for the selective elimination of undifferentiated hESCs intermixed in hESC progeny prior to transplantation and as a paradigm for cell toxicity studies. In addition, one direct application for such assays could be the identification of compounds with teratogenic potential. Differentiation of self-renewal promoting effects at the hESC stage may reflect the potential for similar effects during early human embryonic development in vivo. The identification of retinoic acid, a well known teratogenic compound in our study supports this hypothesis. The emergence of national HTS centers such as the 10 centers established through the NIH Roadmap initiative http://mli.nih.gov/mli/mlscn/screening-centers and additional centers sponsored by state initiatives in the US and multinational initiatives in Europe and Asia should provide access to HTS technology more readily. This will be particularly important for users in the stem cell field that do not have access to individual screening centers at their own institutions. Once an HTS assay has been optimized for use, which represents the longest step, the actual screening process is quite rapid. The cell numbers required for our assay (50 – 100 × 106 cells) can be generated starting from a single 6 cm plate of hESCs in approximately 4 weeks. From the time of plating to full data analysis of a 2880 compound library, the screen can be completed in less than two weeks. There was an overall failure rate of about 1 in 3 assays defined as not reaching acceptable standards in positive and negative controls that are included on every plate (value not comparable to original z’). This relatively large failure rate was cause mainly by the variability in the quality of undifferentiated hESCs maintained in bulk cultures. Future improvements in large scale hESC culturing should help reducing failure rates.
In conclusion, our data demonstrate the feasibility of performing primary HTS assays in hESCs. Our screen resulted in the identification of novel compounds for manipulating hESC fate. The use of HTS assays should dramatically enhance our understanding of hESC biology and facilitate the application of hESCs in drug development and cell therapy.
The hESC lines H9 (WA-09, passages 35 – 45) and H1 (WA-01, passages 30 – 45) were used for this study. HESCs were cultured on mitotically inactivated MEFs as described previously (Perrier et al., 2004) followed by feeder-free growth on matrigel in MEF-conditioned medium (Xu et al., 2001) prior to 384 well plating. For 384 well plating cells were harvested following Accutase dissociation for 20 minutes at 37°C. At that stage, single cell suspensions could be obtained without further mechanical dissociation, and dissociated cells displayed high levels of viability (> 95% based on trypan exclusion). hESCs were plated at 6,000 cells / well from a stirred single cell suspension (typically 50 – 100 × 106 cells / experiment) in MEF-conditioned hESC medium using a Multidrop dispenser (Thermo). Quality control (QC) steps were added with each set of experiments using Alamar blue assay. For the assay, the dye resazurin was added to living hESCs at 5 μM according to the specification of the manufacturer, and fluorescence was measured 24 hours later providing a quantitative estimate of the number of metabolically active cells. At the time of the assay (day 2 of assay), KSR medium (Perrier et al., 2004) was added to the plates of compounds diluted in 100% DMSO, and the KSR diluted compounds were transferred to the hESC plates at 10μM final concentration in 1% DMSO (v/v). The cells were maintained in compounds until day 7. DMSO is the preferred solvent for chemical libraries because of its stability and low temperature flash point making it amenable for use with robotic liquid handling systems. To address potential interference of DMSO with hESC self-renewal, we demonstrated robust long-term maintenance of undifferentiated hESCs in the presence of 1% DMSO when using conditioned medium + FGF2 (10ng/ml). These studies were based on measuring expression of pluripotency markers (Oct4, SSEA3), cell cycle analysis, and percentage of undifferentiated colonies over a period of 5–10 passages (data not shown)
Every 384 well plate in the screen included 16 wells each for high and for low internal controls. FGF2 (R&D Systems) was supplemented at 200ng/ml (Levenstein et al., 2006) in high control wells (promoting self-renewal). BMP4 (R&D Systems) was added at 200ng/ml in low control wells (induction of differentiation). These concentrations were selected based on dose titration assays for high and low Oct4 signal under HTS conditions taking into account stability and degradation pattern of the factors (data not shown). Additional control treatments were used to demonstrate the robustness of the assay: IGF (30ng/ml-R&D Systems), BMP2 (100ng/ml-R&D Systems) and PD98059 (50μM- Tocris Bioscience).
Fixation in 4% paraformaldehyde/0.15% picric acid was followed by standard fluorescent immunocytochemical techniques using the following primary antibodies: Monoclonal: BrdU (Sigma), Oct4 (Santa-Cruz), Polyclonal: Nanog (R&D Systems), Sox17 (R&D Systems), CK18 (Santa-Cruz), and Pax6 (Covance). All stainings were adapted for 384 well plates for automation. Liquid-handling systems were used to aspirate and dispense liquid from and in 384 well plates. Aspirations were performed using a Multidrop 384 (Thermo, Electron Corporation, Waltham, MA). Dispensions of fluids were performed using either the Multidrop 384, for rinsing steps, or a Flexdrop IV (PerkinElmer, Waltham, MA) for Matrigel (BD Biosciences, San Jose, CA), cell plating and antibody staining. Compounds were transferred on cells using a TPS-384 total pipetting platform (Apricot). Confocal image acquisition was carried out with an InCell Analyzer 3000 (GE Healthcare). An average of four images per well were collected; two images at 488-nm and two images for Hoechst 33342, using a 10X objective. Data analysis was performed using the InCell Investigator software. Images were first processed by a flat-field correction. The object intensity analysis module from the IN Cell Investigator software (GE Healthcare, USA) was used to measure the signal intensity of objects in the defined measurement region. The algorithm first identified the nucleus in the marker channel and then quantified the signal region within the area of the nucleus. Average signal intensity was calculated by measuring total intensity divided by the number of objects (i.e. cells.) and expressed as arbitrary units (AU). Screening data files obtained from the Amersham IN Cell Investigator software post processing were loaded into the HTS Core Screening Data Management System (Antczak et al., 2007). Hits were defined as the repeated top ranked compounds/plate that also exhibited at least a > 30% increase (inducing compounds) or > 30% decrease (inhibitory compounds) in normalized Oct4 values in ≥ 2 independent runs in duplicates. Data for some compounds data were confirmed in up to 6 independent runs.
The effect was tested over multiple passages on undifferentiated H9 hESCs passaged using a standard dispase (Worthington Biochemical Corp., Lakewood, NJ) protocol (Perrier et al., 2004) under feeder-free conditions on Matrigel in the presence of MEF-conditioned medium (CM) and FGF2 (10ng/ml). Two days after plating, CM was replaced with HES medium containing: THEA (50μM), SNM (50μM), FBP (25μM) or GTFX (50μM). Control conditions were matched with regard to DMSO concentration. At each passage, cells were dissociated with Accutase and subjected to cycle analysis by FACS after propidium iodide-Citrate solution treatment (50μg/ml propidium iodide, 0.1% sodium citrate dihydrate, 200μg/ml RNaseA and 0.1% Triton X-100). In brief, 2 × 105 dissociated hESCs were washed once with PBS, incubated with PI-Citrate solution for 10 min in the dark and analyzed. FACS analysis was performed on a FACSCalibur system (BD Biosciences). For Oct4/DAPI ratio quantification, cells were subjected to Oct4 immunocytochemistry and four 0.15 mm2 microscopic fields were imaged and analysed using ImageJ software (NIH, Bethesda, MD). In brief, Oct4 and DAPI channels were subjected to thresholding independently, the stain surface area was measured and finally, the ratio of the surface areas was calculated. Apoptotic cell death was quantified after plating of H9 hESCs on matrigel in CM + FGF2 (10ng/ml) for two days, followed by four days of drug treatment and subjection to fixation and TUNEL assay (Roche Diagnostics Corporation, Indianapolis, IN). The percentage of TUNEL+ cells was quantified in uniform randomly selected fields of 0.15 mm2 (n=5/group) together with DAPI for ratio measurement. Pluripotency was assessed by directed differentiation of H9 maintained over three passages in the presence of self-renewal compounds. For endoderm induction, hESCs were treated with Activin A (100 ng/ml) for 6 days in DMEM/F12 medium containing 1% FBS (Invitrogen, Carlsbad, CA) as previously described (Watanabe et al., 2007). For mesoderm induction, cells were cultured in RPMI medium supplemented with B27 (Invitrogen) as previously reported (Laflamme et al., 2007). In brief, Activin A (100 ng/ml; R&D Systems) was added for 1 day, followed by BMP4 (10ng/ml) for an additional next 4 days. After withdrawal of BMP4, cells were maintained in RPMI-B27 medium for 12 days before staining. Neuroectoderm induction was achieved using a co-culture system (Perrier et al., 2004) where hESCs are plated on stromal cells (MS5) in the presence of Noggin (250 ng/ml) added starting at day 3 of differentiation. Immunofluorescence analysis was performed at day 12 of differentiation. The following antibodies were used: Sox17 for endoderm induction, smooth muscle actin (SMA; Sigma) for mesodermal differentiation, and Pax6 (Covance) for neuroectoderm induction.
hESCs were plated in 384 well plates following the described HTS protocol. Each condition (untreated medium , BMP4, RA, SLG, CYM, SRM or ActivinA) had 6 replicates for each staining (CK18 or Sox17). Four images / well were captured by InCell 3000 scanner and quantified for intensity signal measurements. Percentages of Sox17+ cells were quantified by computer-assisted analysis using ImageJ with automated thresholding as described above. All data were normalized to Hoechst (intensity signal) or DAPI (Percentage of Sox17+ cells).
The library used for the current study combines chemicals obtained commercially from Prestwick and MicroSource. The MicroSource library (www.msdiscovery.com) contains 2000 biologically active and structurally diverse compounds from known drugs, experimental bioactives, and pure natural products. The library includes a collection of 720 natural products and their derivatives, a range of simple and complex oxygen-containing heterocycles, alkaloids, sesquiterpenes, diterpenes, pentacyclic triterpenes, sterols, and many other diverse representatives. The Prestwick Chemical Library (http://www.prestwickchemical.fr) is a collection of 880 high-purity chemical compounds selected for structural diversity, and representing drug classes with a broad spectrum of therapeutic uses. More than 85% of its compounds are marketed drugs.
To assess the performance of the assay, the Z′ factor value was determined. This factor represents a dimensionless parameter which value can range between <0 to 1. It is defined as 1−Z′ = (3σc+ + 3σc−)/|μ c+ − μcμ|, where σc+, σc−, μc+, and μc− are the standard deviations (σ) and averages (μ) of the high (c+) and low (c−) controls (Zhang et al., 1999). To obtain this value, three 384 plates of high control (1%DMSO+FGF2) and low controls (1%DMSO+BMP4) were analyzed (Supplementary Fig. 3). For dose-response analysis each 384 well plate of exposed to hits from the primary screen contained replicate wells for high and low controls. Each concentration was applied in duplicates over a range of 12 different doubling dilutions between 0 and 100 μM.
RNA was isolated (Qiagen), labeled and prepared for Illumina bead arrays following the instructions of the manufacturer (SKI Genomics Core Facility). All data were analyzed using the Bioconductor packages (www.bioconductor.com) for the R statistical system. The arrays were quantitated using the Illumina BeadStudio program and raw data (normalized and non-background subtracted signal levels) were exported. Normalization was carried out using the quantile method for the AFFT Bioconductor package. Probes that had detection p-values > 0.01 in 75% or more samples were removed as noise. Expression signal levels were then log 2 transformed. To determine genes that are differentially expressed between the various sample types the LIMMA package was used. This performs a variant of linear models (t-test or ANOVA) were the variance is corrected to deal with small sample numbers. To account for the multi-testing issue the False Discovery Rate (FDR) method was used with cutoffs of 0.20 for the HES samples comparison (self-renewal factors) and 0.05 for the KSR sample comparison (differentiation factors). To determine the significance of overlaps between different gene sets the standard Fisher test (hypergeometric distribution) was performed. All raw data are available online (GEO accession number: GSE11302).
Functional annotation analysis was performed using the DAVID platform (DAVID: http://david.abcc.ncifcrf.gov, (Dennis, Jr. et al., 2003) ). A list of differentially regulated transcripts for a given compound was established at the FDR levels described above. For analysis we selected only transcripts with a fold change of ≥ 2 compared with control. For those compounds with less than 200 differential transcripts at a fold change level of 2, the analysis was extended to those with a fold change of ≥ 1.5. Final selection of categories for presentation eliminated duplicate categories that contained identical gene lists. Furthermore all included terms exhibited significant p-values for enrichment and generally contained > 2 members / category.
We thank the members of the MSKCC HTS Core Facility for their help during the course of this study, especially C. Radu, D. Shum, and G.E. Sanchez. We also thank V. Tabar and H. López-Schier for critical comments on the manuscript and S. Gong and N. Heintz for providing the Oct4::eGFP BAC. Supported by NINDS grant (R21NS44231), The Starr Foundation, Mr. W.H. Goodwin and Mrs. A. Goodwin and the Commonwealth Foundation for Cancer Research, The William Randolph Hearst Foundation, and The Experimental Therapeutics Center of MSKCC.
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.