Pluripotent stem cell (PSC)-derived organoids have emerged as novel multicellular models of human tissue development but display immature phenotypes, aberrant tissue fates, and a limited subset of cells. Here, we demonstrate that integrated analysis and engineering of gene regulatory networks (GRNs) in PSC-derived multilineage human liver organoids direct maturation and vascular morphogenesis in vitro. Overexpression of PROX1 and ATF5, combined with targeted CRISPR-based transcriptional activation of endogenous CYP3A4, reprograms tissue GRNs and improves native liver functions, such as FXR signaling, CYP3A4 enzymatic activity, and stromal cell reactivity. The engineered tissues possess superior liver identity when compared with other PSC-derived liver organoids and show the presence of hepatocyte, biliary, endothelial, and stellate-like cell populations in single-cell RNA-seq analysis. Finally, they show hepatic functions when studied in vivo. Collectively, our approach provides an experimental framework to direct organogenesis in vitro by systematically probing molecular pathways and transcriptional networks that promote tissue development.

Mechanosensitive hair cells in the cochlea are responsible for hearing but are vulnerable to damage by genetic mutations and environmental insults. The paucity of human cochlear tissues makes it difficult to study cochlear hair cells. Organoids offer a compelling platform to study scarce tissues in vitro; however, derivation of cochlear cell types has proven non-trivial. Here, using 3D cultures of human pluripotent stem cells, we sought to replicate key differentiation cues of cochlear specification. We found that timed modulations of Sonic Hedgehog and WNT signaling promote ventral gene expression in otic progenitors. Ventralized otic progenitors subsequently give rise to elaborately patterned epithelia containing hair cells with morphology, marker expression, and functional properties consistent with both outer and inner hair cells in the cochlea. These results suggest that early morphogenic cues are sufficient to drive cochlear induction and establish an unprecedented system to model the human auditory organ.

Adult haematopoiesis is the outcome of distinct haematopoietic stem cell (HSC) subtypes with self-renewable repopulating ability, but with different haematopoietic cell lineage outputs. The molecular basis for this heterogeneity is largely unknown. BMP signalling regulates HSCs as they are first generated in the aorta-gonad-mesonephros region, but at later developmental stages, its role in HSCs is controversial. Here we show that HSCs in murine fetal liver and the bone marrow are of two types that can be prospectively isolated--BMP activated and non-BMP activated. Clonal transplantation demonstrates that they have distinct haematopoietic lineage outputs. Moreover, the two HSC types differ in intrinsic genetic programs, thus supporting a role for the BMP signalling axis in the regulation of HSC heterogeneity and lineage output. Our findings provide insight into the molecular control mechanisms that define HSC types and have important implications for reprogramming cells to HSC fate and treatments targeting distinct HSC types.

Remodelling the methylome is a hallmark of mammalian development and cell differentiation. However, current knowledge of DNA methylation dynamics in human tissue specification and organ development largely stems from the extrapolation of studies in vitro and animal models. Here, we report on the DNA methylation landscape using the 450k array of four human tissues (amnion, muscle, adrenal and pancreas) during the first and second trimester of gestation (9,18 and 22 weeks). We show that a tissue-specific signature, constituted by tissue-specific hypomethylated CpG sites, was already present at 9 weeks of gestation (W9). Furthermore, we report large-scale remodelling of DNA methylation from W9 to W22. Gain of DNA methylation preferentially occurred near genes involved in general developmental processes, whereas loss of DNA methylation mapped to genes with tissue-specific functions. Dynamic DNA methylation was associated with enhancers, but not promoters. Comparison of our data with external fetal adrenal, brain and liver revealed striking similarities in the trajectory of DNA methylation during fetal development. The analysis of gene expression data indicated that dynamic DNA methylation was associated with the progressive repression of developmental programs and the activation of genes involved in tissue-specific processes. The DNA methylation landscape of human fetal development provides insight into regulatory elements that guide tissue specification and lead to organ functionality.

The human fetal immune system must protect the infant against the sudden exposure to a large variety of pathogens upon birth. While it is known that the fetal immune system develops in sequential waves, relatively little is known about the composition of the innate and adaptive immune system in the tissues. Here, we applied high-dimensional mass cytometry to profile the immune system in human fetal liver, spleen, and intestine. With Hierarchical Stochastic Neighbor Embedding (HSNE) we distinguished 177 distinct immune cell clusters, including both previously identified and novel cell clusters. PCA analysis indicated substantial differences between the compositions of the immune system in the different organs. Through dual t-SNE we identified tissue-specific cell clusters, which were found both in the innate and adaptive compartment. To determine the spatial location of tissue-specific subsets we developed a 31-antibody panel to reveal both the immune compartment and surrounding stromal elements through analysis of snap-frozen tissue samples with imaging mass cytometry. Imaging mass cytometry reconstructed the tissue architecture and allowed both the characterization and determination of the location of the various immune cell clusters within the tissue context. Moreover, it further underpinned the distinctness of the immune system in the tissues. Thus, our results provide evidence for early compartmentalization of the adaptive and innate immune compartment in fetal spleen, liver, and intestine. Together, our data provide a unique and comprehensive overview of the composition and organization of the human fetal immune system in several tissues.

The fetus is thought to be protected from exposure to foreign antigens, yet CD45RO+ T cells reside in the fetal intestine. Here we combined functional assays with mass cytometry, single-cell RNA sequencing and high-throughput T cell antigen receptor (TCR) sequencing to characterize the CD4+ T cell compartment in the human fetal intestine. We identified 22 CD4+ T cell clusters, including naive-like, regulatory-like and memory-like subpopulations, which were confirmed and further characterized at the transcriptional level. Memory-like CD4+ T cells had high expression of Ki-67, indicative of cell division, and CD5, a surrogate marker of TCR avidity, and produced the cytokines IFN-γ and IL-2. Pathway analysis revealed a differentiation trajectory associated with cellular activation and proinflammatory effector functions, and TCR repertoire analysis indicated clonal expansions, distinct repertoire characteristics and interconnections between subpopulations of memory-like CD4+ T cells. Imaging mass cytometry indicated that memory-like CD4+ T cells colocalized with antigen-presenting cells. Collectively, these results provide evidence for the generation of memory-like CD4+ T cells in the human fetal intestine that is consistent with exposure to foreign antigens.

Current methods to generate human primordial germ cell-like cells (hPGCLCs) from human pluripotent stem cells (hPSCs) can be inefficient, and it is challenging to generate sufficient hPGCLCs to optimize in vitro gametogenesis. We present a differentiation method that uses diluted basement membrane extract (BMEx) and low BMP4 concentration to efficiently induce hPGCLC differentiation in scalable 2D cell culture. We show that BMEx overlay potentiated BMP/SMAD signaling, induced lumenogenesis, and increased expression of key hPGCLC-progenitor markers such as TFAP2A and EOMES. hPGCLCs that were generated using the BMEx overlay method were able to upregulate more mature germ cell markers, such as DAZL and DDX4, in human fetal ovary reconstitution culture. These findings highlight the importance of BMEx during hPGCLC differentiation and demonstrate the potential of the BMEx overlay method to interrogate the formation of PGCs and amnion in humans, as well as to investigate the next steps to achieve in vitro gametogenesis.

Accumulating evidence demonstrates important roles for metabolism in cell fate determination. However, it is a challenge to assess metabolism at a spatial resolution that acknowledges both heterogeneity and cellular dynamics in its tissue microenvironment. Using a multi-omics platform to study cell-type-specific dynamics in metabolism in complex tissues, we describe the metabolic trajectories during nephrogenesis in the developing human kidney. Exploiting in situ analysis of isotopic labeling, a shift from glycolysis toward fatty acid β-oxidation was observed during the differentiation from the renal vesicle toward the S-shaped body and the proximal tubules. In addition, we show that hiPSC-derived kidney organoids are characterized by a metabolic immature phenotype that fails to use mitochondrial long-chain fatty acids for energy metabolism. Furthermore, supplementation of butyrate enhances tubular epithelial differentiation and maturation in cultured kidney organoids. Our findings highlight the relevance of understanding metabolic trajectories to efficiently guide stem cell differentiation.

Human induced pluripotent stem cell (hiPSC)-derived hair-bearing skin organoids offer exciting new possibilities for modeling diseases like epidermolysis bullosa (EB). These inherited diseases affect 1 in 30,000 people worldwide and result from perturbed expression and/or structure of components of the epidermal-dermal junction (EDJ). To establish whether hiPSC-derived skin organoids might be able to capture salient features of EB, it is thus important to characterize their EDJ. Here, we report successful generation of hair-bearing skin organoids from two hiPSC lines that exhibited fully stratified interfollicular epidermis. Using immunofluorescence and electron microscopy, we showed that basal keratinocytes in organoids adhere to laminin-332 and type IV collagen-rich basement membrane via type I hemidesmosomes and integrin β1-based adhesion complexes. Importantly, we demonstrated that EDJs in organoids are almost devoid of type VII collagen, a fibril that mediates anchorage of the epidermis to dermis. This should be considered when using skin organoids for EB modeling.

During human fetal development, sex differentiation occurs not only in the gonads but also in the adjacent developing reproductive tract. However, while the cellular composition of male and female human fetal gonads is well described, that of the adjacent developing reproductive tract remains poorly characterized. Here, we performed single-cell transcriptomics on male and female human fetal gonads together with the adjacent developing reproductive tract from first and second trimesters, highlighting the morphological and molecular changes during sex differentiation. We validated different cell populations of the developing reproductive tract and gonads and compared the molecular signatures between the first and second trimesters, as well as between sexes, to identify conserved and sex-specific features. Together, our study provides insights into human fetal sex-specific gonadogenesis and development of the reproductive tract beyond the gonads.

The current understanding of mammalian kidney development is largely based on mouse models. Recent landmark studies revealed pervasive differences in renal embryogenesis between mouse and human. The scarcity of detailed gene expression data in humans therefore hampers a thorough understanding of human kidney development and the possible developmental origin of kidney diseases. In this paper, we present a single-cell transcriptomics study of the human fetal kidney. We identified 22 cell types and a host of marker genes. Comparison of samples from different developmental ages revealed continuous gene expression changes in podocytes. To demonstrate the usefulness of our data set, we explored the heterogeneity of the nephrogenic niche, localized podocyte precursors, and confirmed disease-associated marker genes. With close to 18,000 renal cells from five different developmental ages, this study provides a rich resource for the elucidation of human kidney development, easily accessible through an interactive web application.

Signaling by Bone Morphogenetic Proteins (BMP) has been implicated in early lung development, adult lung homeostasis and tissue-injury repair. However, the precise mechanism of action and the spatio-temporal pattern of BMP-signaling during these processes remains inadequately described. To address this, we have utilized a transgenic line harboring a BMP-responsive eGFP-reporter allele (BRE-eGFP) to construct the first detailed spatiotemporal map of canonical BMP-pathway activation during lung development, homeostasis and adult-lung injury repair. We demonstrate that during the pseudoglandular stage, when branching morphogenesis progresses in the developing lung, canonical BMP-pathway is active mainly in the vascular network and the sub-epithelial smooth muscle layer of the proximal airways. Activation of the BMP-pathway becomes evident in epithelial compartments only after embryonic day (E) 14.5 primarily in cells negative for epithelial-lineage markers, located in the proximal portion of the airway-tree, clusters adjacent to neuro-epithelial-bodies (NEBs) and in a substantial portion of alveolar epithelial cells. The pathway becomes activated in isolated E12.5 mesenchyme-free distal epithelial buds cultured in Matrigel suggesting that absence of reporter activity in these regions stems from a dynamic cross-talk between endoderm and mesenchyme. Epithelial cells with activated BMP-pathway are enriched in progenitors capable of forming colonies in three-dimensional Matrigel cultures.As lung morphogenesis approaches completion, eGFP-expression declines and in adult lung its expression is barely detectable. However, upon tissue-injury, either with naphthalene or bleomycin, the canonical BMP-pathways is re-activated, in bronchial or alveolar epithelial cells respectively, in a manner reminiscent to early lung development and in tissue areas where reparatory progenitor cells reside. Our studies illustrate the dynamic activation of canonical BMP-pathway during lung development and adult lung tissue-repair and highlight its involvement in two important processes, namely, the early development of the pulmonary vasculature and the management of epithelial progenitor pools both during lung development and repair of adult lung tissue-injury.

Hematopoietic stem cells (HSC), the self-renewing cells of the adult blood differentiation hierarchy, are generated during embryonic stages. The first HSCs are produced in the aorta-gonad-mesonephros (AGM) region of the embryo through endothelial to a hematopoietic transition. BMP4 and Hedgehog affect their production and expansion, but it is unknown whether they act to affect the same HSCs. In this study using the BRE GFP reporter mouse strain that identifies BMP/Smad-activated cells, we find that the AGM harbors two types of adult-repopulating HSCs upon explant culture: One type is BMP-activated and the other is a non-BMP-activated HSC type that is indirectly controlled by Hedgehog signaling through the VEGF pathway. Transcriptomic analyses demonstrate that the two HSC types express distinct but overlapping genetic programs. These results revealing the bifurcation in HSC types at early embryonic stages in the AGM explant model suggest that their development is dependent upon the signaling molecules in the microenvironment.

Pluripotent stem cells (PSCs) could provide a powerful system to model development of the human esophagus, whose distinct tissue organization compared to rodent esophagus suggests that developmental mechanisms may not be conserved between species. We therefore established an efficient protocol for generating esophageal progenitor cells (EPCs) from human PSCs. We found that inhibition of TGF-ß and BMP signaling is required for sequential specification of EPCs, which can be further purified using cell-surface markers. These EPCs resemble their human fetal counterparts and can recapitulate normal development of esophageal stratified squamous epithelium during in vitro 3D cultures and in vivo. Importantly, combining hPSC differentiation strategies with mouse genetics elucidated a critical role for Notch signaling in the formation of this epithelium. These studies therefore not only provide an efficient approach to generate EPCs, but also offer a model system to study the regulatory mechanisms underlying development of the human esophagus.

NANOG is a key player in pluripotency and its expression is restricted to pluripotent cells of the inner cell mass, the epiblast and to primordial germ cells. Spermatogenesis is closely associated with pluripotency, because through this process highly specialized sperm cells are produced that contribute to the formation of totipotent zygotes. Nevertheless, it is unknown if NANOG plays a role in this process.

During healthy pregnancy, a balanced microenvironment at the maternal-fetal interface with coordinated interaction between various immune cells is necessary to maintain immunological tolerance. While specific decidual immune cell subsets have been investigated, a system-wide unbiased approach is lacking. Here, mass cytometry was applied for data-driven, in-depth immune profiling of the total leukocyte population isolated from first, second, and third trimester decidua, as well as maternal peripheral blood at time of delivery. The maternal-fetal interface showed a unique composition of immune cells, different from peripheral blood, with significant differences between early and term pregnancy samples. Profiling revealed substantial heterogeneity in the decidual lymphoid and myeloid cell lineages that shape gestational-specific immune networks and putative differentiation trajectories over time during gestation. Uncovering the overall complexity at the maternal-fetal interface throughout pregnancy resulted in a human atlas that may serve as a foundation upon which comprehension of the immune microenvironment and alterations thereof in pregnancy complications can be built.

DNA methylation (5mC) is central to cellular identity. The global erasure of 5mC from the parental genomes during preimplantation mammalian development is critical to reset the methylome of gametes to the cells in the blastocyst. While active and passive modes of demethylation have both been suggested to play a role in this process, the relative contribution of these two mechanisms to 5mC erasure remains unclear. Here, we report a single-cell method (scMspJI-seq) that enables strand-specific quantification of 5mC, allowing us to systematically probe the dynamics of global demethylation. When applied to mouse embryonic stem cells, we identified substantial cell-to-cell strand-specific 5mC heterogeneity, with a small group of cells displaying asymmetric levels of 5mCpG between the two DNA strands of a chromosome suggesting loss of maintenance methylation. Next, in preimplantation mouse embryos, we discovered that methylation maintenance is active till the 16-cell stage followed by passive demethylation in a fraction of cells within the early blastocyst at the 32-cell stage of development. Finally, human preimplantation embryos qualitatively show temporally delayed yet similar demethylation dynamics as mouse embryos. Collectively, these results demonstrate that scMspJI-seq is a sensitive and cost-effective method to map the strand-specific genome-wide patterns of 5mC in single cells.

Human ovarian folliculogenesis is a highly regulated and complex process. Characterization of follicular cell signatures during this dynamic process is important to understand follicle fate (to grow, become dominant, or undergo atresia). The transcriptional signature of human oocytes and granulosa cells (GCs) in early-growing and ovulatory follicles have been previously described; however, that of oocytes with surrounding GCs in small antral follicles have not been studied yet. Here, we have generated a unique dataset of single-cell transcriptomics (SmartSeq2) consisting of the oocyte with surrounding GCs from several individual (non-dominant) small antral follicles isolated from adult human ovaries. We have identified two main types of (healthy) follicles, with a distinct oocyte and GC signature. Using the CellphoneDB algorithm, we then investigated the bi-directional ligand-receptor interactions regarding the transforming growth factor-β (TGFβ)/bone morphogenetic protein (BMP), wingless-type (MMTV)-integration site (WNT), NOTCH, and receptor tyrosine kinases (RTK) signaling pathways between oocyte and GCs within each antral follicle type. Our work not only revealed the diversity of small antral follicles, but also contributes to fill the gap in mapping the molecular landscape of human folliculogenesis and oogenesis.

Trophoblasts are essential in fetal-maternal interaction during pregnancy. The goal was to study HLA profiles of primary trophoblasts derived from placentas, and to investigate their usefulness in studying interaction with immune cells.

Optimization of CRISPR/Cas9-mediated genome engineering has resulted in base editors that hold promise for mutation repair and disease modeling. Here, we demonstrate the application of base editors for the generation of complex tumor models in human ASC-derived organoids. First we show efficacy of cytosine and adenine base editors in modeling CTNNB1 hot-spot mutations in hepatocyte organoids. Next, we use C > T base editors to insert nonsense mutations in PTEN in endometrial organoids and demonstrate tumorigenicity even in the heterozygous state. Moreover, drug sensitivity assays on organoids harboring either PTEN or PTEN and PIK3CA mutations reveal the mechanism underlying the initial stages of endometrial tumorigenesis. To further increase the scope of base editing we combine SpCas9 and SaCas9 for simultaneous C > T and A > G editing at individual target sites. Finally, we show that base editor multiplexing allow modeling of colorectal tumorigenesis in a single step by simultaneously transfecting sgRNAs targeting five cancer genes.