High-throughput micropatterning platform reveals Nodal-dependent bisection of peri-gastrulation-associated versus preneurulation-associated fate patterning.

In vitro models of postimplantation human development are valuable to the fields of regenerative medicine and developmental biology. Here, we report characterization of a robust in vitro platform that enabled high-content screening of multiple human pluripotent stem cell (hPSC) lines for their ability to undergo peri-gastrulation-like fate patterning upon bone morphogenetic protein 4 (BMP4) treatment of geometrically confined colonies and observed significant heterogeneity in their differentiation propensities along a gastrulation associable and neuralization associable axis. This cell line-associated heterogeneity was found to be attributable to endogenous Nodal expression, with up-regulation of Nodal correlated with expression of a gastrulation-associated gene profile, and Nodal down-regulation correlated with a preneurulation-associated gene profile expression. We harness this knowledge to establish a platform of preneurulation-like fate patterning in geometrically confined hPSC colonies in which fates arise because of a BMPs signalling gradient conveying positional information. Our work identifies a Nodal signalling-dependent switch in peri-gastrulation versus preneurulation-associated fate patterning in hPSC cells, provides a technology to robustly assay hPSC differentiation outcomes, and suggests conserved mechanisms of organized fate specification in differentiating epiblast and ectodermal tissues.

A stepwise model of reaction-diffusion and positional information governs self-organized human peri-gastrulation-like patterning.

How position-dependent cell fate acquisition occurs during embryogenesis is a central question in developmental biology. To study this process, we developed a defined, high-throughput assay to induce peri-gastrulation-associated patterning in geometrically confined human pluripotent stem cell (hPSC) colonies. We observed that, upon BMP4 treatment, phosphorylated SMAD1 (pSMAD1) activity in the colonies organized into a radial gradient. We developed a reaction-diffusion (RD)-based computational model and observed that the self-organization of pSMAD1 signaling was consistent with the RD principle. Consequent fate acquisition occurred as a function of both pSMAD1 signaling strength and duration of induction, consistent with the positional-information (PI) paradigm. We propose that the self-organized peri-gastrulation-like fate patterning in BMP4-treated geometrically confined hPSC colonies arises via a stepwise model of RD followed by PI. This two-step model predicted experimental responses to perturbations of key parameters such as colony size and BMP4 dose. Furthermore, it also predicted experimental conditions that resulted in RD-like periodic patterning in large hPSC colonies, and rescued peri-gastrulation-like patterning in colony sizes previously thought to be reticent to this behavior.

Context-explorer: Analysis of spatially organized protein expression in high-throughput screens.

A growing body of evidence highlights the importance of the cellular microenvironment as a regulator of phenotypic and functional cellular responses to perturbations. We have previously developed cell patterning techniques to control population context parameters, and here we demonstrate context-explorer (CE), a software tool to improve investigation cell fate acquisitions through community level analyses. We demonstrate the capabilities of CE in the analysis of human and mouse pluripotent stem cells (hPSCs, mPSCs) patterned in colonies of defined geometries in multi-well plates. CE employs a density-based clustering algorithm to identify cell colonies. Using this automatic colony classification methodology, we reach accuracies comparable to manual colony counts in a fraction of the time, both in micropatterned and unpatterned wells. Classifying cells according to their relative position within a colony enables statistical analysis of spatial organization in protein expression within colonies. When applied to colonies of hPSCs, our analysis reveals a radial gradient in the expression of the transcription factors SOX2 and OCT4. We extend these analyses to colonies of different sizes and shapes and demonstrate how the metrics derived by CE can be used to asses the patterning fidelity of micropatterned plates. We have incorporated a number of features to enhance the usability and utility of CE. To appeal to a broad scientific community, all of the software's functionality is accessible from a graphical user interface, and convenience functions for several common data operations are included. CE is compatible with existing image analysis programs such as CellProfiler and extends the analytical capabilities already provided by these tools. Taken together, CE facilitates investigation of spatially heterogeneous cell populations for fundamental research and drug development validation programs.

Modeling signaling-dependent pluripotency with Boolean logic to predict cell fate transitions.

Pluripotent stem cells (PSCs) exist in multiple stable states, each with specific cellular properties and molecular signatures. The mechanisms that maintain pluripotency, or that cause its destabilization to initiate development, are complex and incompletely understood. We have developed a model to predict stabilized PSC gene regulatory network (GRN) states in response to input signals. Our strategy used random asynchronous Boolean simulations (R-ABS) to simulate single-cell fate transitions and strongly connected components (SCCs) strategy to represent population heterogeneity. This framework was applied to a reverse-engineered and curated core GRN for mouse embryonic stem cells (mESCs) and used to simulate cellular responses to combinations of five signaling pathways. Our simulations predicted experimentally verified cell population compositions and input signal combinations controlling specific cell fate transitions. Extending the model to PSC differentiation, we predicted a combination of signaling activators and inhibitors that efficiently and robustly generated a Cdx2+Oct4- cells from naïve mESCs. Overall, this platform provides new strategies to simulate cell fate transitions and the heterogeneity that typically occurs during development and differentiation.

High-throughput fingerprinting of human pluripotent stem cell fate responses and lineage bias.

Populations of cells create local environments that lead to emergent heterogeneity. This is particularly evident with human pluripotent stem cells (hPSCs): microenvironmental heterogeneity limits hPSC cell fate control. We developed a high-throughput platform to screen hPSCs in configurable microenvironments in which we optimized colony size, cell density and other parameters to achieve rapid and robust cell fate responses to exogenous cues. We used this platform to perform single-cell protein expression profiling, revealing that Oct4 and Sox2 costaining discriminates pluripotent, neuroectoderm, primitive streak and extraembryonic cell fates. We applied this Oct4-Sox2 code to analyze dose responses of 27 developmental factors to obtain lineage-specific concentration optima and to quantify cell line-specific endogenous signaling pathway activation and differentiation bias. We demonstrated that short-term responses predict definitive endoderm induction efficiency and can be used to rescue differentiation of cell lines reticent to cardiac induction. This platform will facilitate high-throughput hPSC-based screening and quantification of lineage-induction bias.

Blebbing as a physical force in cancer EMT - parallels with mitosis.

Epithelial to mesenchymal transition (EMT) during metastasis is initially a two-step process beginning with delamination of cells from the solid tumor followed by acquisition of a migratory phenotype. Several reports indicate that plasma membrane blebbing, associated with cell division, coincides with cell delamination during developmental EMT. This raises a speculative question if blebbing drives EMT in cancer cells in a similar way. Here, we review available data on factors and processes that may support such a connection.

Virtual cells in a virtual microenvironment recapitulate early development-like patterns in human pluripotent stem cell colonies.

The mechanism by which morphogenetic signals engage the regulatory networks responsible for early embryonic tissue patterning is incompletely understood. Here, we developed a minimal gene regulatory network (GRN) model of human pluripotent stem cell (hPSC) lineage commitment and embedded it into "cellular" agents that respond to a dynamic morphogenetic signaling microenvironment. Simulations demonstrated that GRN wiring had significant non-intuitive effects on tissue pattern order, composition, and dynamics. Experimental perturbation of GRN connectivities supported model predictions and demonstrated the role of OCT4 as a master regulator of peri-gastrulation fates. Our so-called GARMEN strategy provides a multiscale computational platform to understand how single-cell-based regulatory interactions scale to tissue domains. This foundation provides new opportunities to simulate the impact of network motifs on normal and aberrant tissue development.

Endogenous suppression of WNT signalling in human embryonic stem cells leads to low differentiation propensity towards definitive endoderm.

Low differentiation propensity towards a targeted lineage can significantly hamper the utility of individual human pluripotent stem cell (hPSC) lines in biomedical applications. Here, we use monolayer and micropatterned cell cultures, as well as transcriptomic profiling, to investigate how variability in signalling pathway activity between human embryonic stem cell lines affects their differentiation efficiency towards definitive endoderm (DE). We show that endogenous suppression of WNT signalling in hPSCs at the onset of differentiation prevents the switch from self-renewal to DE specification. Gene expression profiling reveals that this inefficient switch is reflected in NANOG expression dynamics. Importantly, we demonstrate that higher WNT stimulation or inhibition of the PI3K/AKT signalling can overcome the DE commitment blockage. Our findings highlight that redirection of the activity of Activin/NODAL pathway by WNT signalling towards mediating DE fate specification is a vulnerable spot, as disruption of this process can result in poor hPSC specification towards DE.

Functional arrays of human pluripotent stem cell-derived cardiac microtissues.

To accelerate the cardiac drug discovery pipeline, we set out to develop a platform that would be capable of quantifying tissue-level functions such as contractile force and be amenable to standard multiwell-plate manipulations. We report a 96-well-based array of 3D human pluripotent stem cell (hPSC)-derived cardiac microtissues - termed Cardiac MicroRings (CaMiRi) - in custom 3D-print-molded multiwell plates capable of contractile force measurement. Within each well, two elastomeric microcantilevers are situated above a circumferential ramp. The wells are seeded with cell-laden collagen, which, in response to the gradual slope of the circumferential ramp, self-organizes around tip-gated microcantilevers to form contracting CaMiRi. The contractile force exerted by the CaMiRi is measured and calculated using the deflection of the cantilevers. Platform responses were robust and comparable across wells, and we used it to determine an optimal tissue formulation. We validated the contractile force response of CaMiRi using selected cardiotropic compounds with known effects. Additionally, we developed automated protocols for CaMiRi seeding, image acquisition, and analysis to enable the measurement of contractile force with increased throughput. The unique tissue fabrication properties of the platform, and the consequent effects on tissue function, were demonstrated upon adding hPSC-derived epicardial cells to the system. This platform represents an open-source contractile force screening system useful for drug screening and tissue engineering applications.