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.

A reductionist approach to determine the effect of cell-cell contact on human epidermal stem cell differentiation.

The balance between stem cell renewal and differentiation is determined by the interplay between intrinsic cellular controls and extrinsic factors presented by the microenvironment, or 'niche'. Previous studies on cultured human epidermis have utilised suspension culture and restricted cell spreading to investigate regulation of differentiation in single keratinocytes. However, keratinocytes are typically adherent to neighbouring cells in vivo. We therefore developed experimental models to investigate the combined effects of cell-ECM adhesion and cell-cell contact. We utilized lipid-modified oligonucleotides to form clusters of keratinocytes which were subsequently placed in suspension to induce terminal differentiation. In this experimental model cell-cell contact had no effect on suspension-induced differentiation of keratinocytes. We next developed a high-throughput platform for robust geometrical confinement of keratinocytes to hexagonal ECM-coated islands permitting direct cell-cell contact between single cells. As in the case of circular islands, differentiation was stimulated on the smallest single hexagonal islands. However, the percentage of involucrin-positive cells on small bowtie islands was significantly lower than on single islands, demonstrating that cell-cell contact reduced differentiation in response to decreased substrate adhesion. None of the small bowtie islands contained two involucrin-positive cells. Rather, if one cell was involucrin-positive the other was involucrin-negative. This suggests that there is intrinsic asymmetry in the effect of cell-cell contact in decreasing differentiation. Thus, our reductionist approaches provide new insights into the effect of the niche on keratinocyte differentiation. STATEMENT OF SIGNIFICANCE: Stem cell behaviour is regulated by a combination of external signals, including the nature of the adhesive substrate and cell-cell interactions. An understanding of how different signals are integrated creates the possibility of developing new biomaterials to promote tissue regeneration and broaden our understanding of skin diseases such as eczema and psoriasis, in which stem cell proliferation and differentiation are perturbed. In this study we have applied two methods to engineer intercellular adhesion of human epidermal stem cells, one involving lipid-modified DNA and the other involving hexagonal micropatterns. We show that the effect of cell-cell adhesion depends on cell-substrate adhesion and uncover evidence that two cells in equivalent environments can nevertheless behave differently.

Embigin is a fibronectin receptor that affects sebaceous gland differentiation and metabolism.

Stem cell renewal and differentiation are regulated by interactions with the niche. Although multiple cell populations have been identified in distinct anatomical compartments, little is known about niche-specific molecular factors. Using skin as a model system and combining single-cell RNA-seq data analysis, immunofluorescence, and transgenic mouse models, we show that the transmembrane protein embigin is specifically expressed in the sebaceous gland and that the number of embigin-expressing cells is negatively regulated by Wnt. The loss of embigin promotes exit from the progenitor compartment and progression toward differentiation, and also compromises lipid metabolism. Embigin modulates sebaceous niche architecture by affecting extracellular matrix organization and basolateral targeting of monocarboxylate transport. We discover through ligand screening that embigin is a direct fibronectin receptor, binding to the N-terminal fibronectin domain without impairing integrin function. Our results solve the long-standing question of how embigin regulates cell adhesion and demonstrate a mechanism that couples adhesion and metabolism.

Spatial profiling of early primate gastrulation in utero.

Gastrulation controls the emergence of cellular diversity and axis patterning in the early embryo. In mammals, this transformation is orchestrated by dynamic signalling centres at the interface of embryonic and extraembryonic tissues1-3. Elucidating the molecular framework of axis formation in vivo is fundamental for our understanding of human development4-6 and to advance stem-cell-based regenerative approaches7. Here we illuminate early gastrulation of marmoset embryos in utero using spatial transcriptomics and stem-cell-based embryo models. Gaussian process regression-based 3D transcriptomes delineate the emergence of the anterior visceral endoderm, which is hallmarked by conserved (HHEX, LEFTY2, LHX1) and primate-specific (POSTN, SDC4, FZD5) factors. WNT signalling spatially coordinates the formation of the primitive streak in the embryonic disc and is counteracted by SFRP1 and SFRP2 to sustain pluripotency in the anterior domain. Amnion specification occurs at the boundaries of the embryonic disc through ID1, ID2 and ID3 in response to BMP signalling, providing a developmental rationale for amnion differentiation of primate pluripotent stem cells (PSCs). Spatial identity mapping demonstrates that primed marmoset PSCs exhibit the highest similarity to the anterior embryonic disc, whereas naive PSCs resemble the preimplantation epiblast. Our 3D transcriptome models reveal the molecular code of lineage specification in the primate embryo and provide an in vivo reference to decipher human development.

Plating human iPSC lines on micropatterned substrates reveals role for ITGB1 nsSNV in endoderm formation.

Quantitative analysis of human induced pluripotent stem cell (iPSC) lines from healthy donors is a powerful tool for uncovering the relationship between genetic variants and cellular behavior. We previously identified rare, deleterious non-synonymous single nucleotide variants (nsSNVs) in cell adhesion genes that are associated with outlier iPSC phenotypes in the pluripotent state. Here, we generated micropatterned colonies of iPSCs to test whether nsSNVs influence patterning of radially ordered germ layers. Using a custom-built image analysis pipeline, we quantified the differentiation phenotypes of 13 iPSC lines that harbor nsSNVs in genes related to cell adhesion or germ layer development. All iPSC lines differentiated into the three germ layers; however, there was donor-specific variation in germ layer patterning. We identified one line that presented an outlier phenotype of expanded endodermal differentiation, which was associated with a nsSNV in ITGB1. Our study establishes a platform for investigating the impact of nsSNVs on differentiation.

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.

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.

Identifying Extrinsic versus Intrinsic Drivers of Variation in Cell Behavior in Human iPSC Lines from Healthy Donors.

Large cohorts of human induced pluripotent stem cells (iPSCs) from healthy donors are a potentially powerful tool for investigating the relationship between genetic variants and cellular behavior. Here, we integrate high content imaging of cell shape, proliferation, and other phenotypes with gene expression and DNA sequence datasets from over 100 human iPSC lines. By applying a dimensionality reduction approach, Probabilistic Estimation of Expression Residuals (PEER), we extracted factors that captured the effects of intrinsic (genetic concordance between different cell lines from the same donor) and extrinsic (cell responses to different fibronectin concentrations) conditions. We identify genes that correlate in expression with intrinsic and extrinsic PEER factors and associate outlier cell behavior with genes containing rare deleterious non-synonymous SNVs. Our study, thus, establishes a strategy for examining the genetic basis of inter-individual variability in cell 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.

Stem cell bioengineering: building from stem cell biology.

New fundamental discoveries in stem cell biology have yielded potentially transformative regenerative therapeutics. However, widespread implementation of stem-cell-derived therapeutics remains sporadic. Barriers that impede the development of these therapeutics can be linked to our incomplete understanding of how the regulatory networks that encode stem cell fate govern the development of the complex tissues and organs that are ultimately required for restorative function. Bioengineering tools, strategies and design principles represent core components of the stem cell bioengineering toolbox. Applied to the different layers of complexity present in stem-cell-derived systems - from gene regulatory networks in single stem cells to the systemic interactions of stem-cell-derived organs and tissues - stem cell bioengineering can address existing challenges and advance regenerative medicine and cellular therapies.

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.

Engineering the haemogenic niche mitigates endogenous inhibitory signals and controls pluripotent stem cell-derived blood emergence.

Efforts to recapitulate haematopoiesis, a process guided by spatial and temporal inductive signals, to generate haematopoietic progenitors from human pluripotent stem cells (hPSCs) have focused primarily on exogenous signalling pathway activation or inhibition. Here we show haemogenic niches can be engineered using microfabrication strategies by micropatterning hPSC-derived haemogenic endothelial (HE) cells into spatially-organized, size-controlled colonies. CD34+VECAD+ HE cells were generated with multi-lineage potential in serum-free conditions and cultured as size-specific haemogenic niches that displayed enhanced blood cell induction over non-micropatterned cultures. Intra-colony analysis revealed radial organization of CD34 and VECAD expression levels, with CD45+ blood cells emerging primarily from the colony centroid area. We identify the induced interferon gamma protein (IP-10)/p-38 MAPK signalling pathway as the mechanism for haematopoietic inhibition in our culture system. Our results highlight the role of spatial organization in hPSC-derived blood generation, and provide a quantitative platform for interrogating molecular pathways that regulate human haematopoiesis.