Limb development starts with the formation of limb buds (LBs), which consist of tissues from two different germ layers; the lateral plate mesoderm-derived mesenchyme and ectoderm-derived surface epithelium. Here, we report means for induction of an LB-like mesenchymal/epithelial complex tissues from murine pluripotent stem cells (PSCs) in vitro. The LB-like tissues selectively differentiate into forelimb- or hindlimb-type mesenchymes, depending on a concentration of retinoic acid. Comparative transcriptome analysis reveals that the LB-like tissues show similar gene expression pattern to that seen in LBs. We also show that manipulating BMP signaling enables us to induce a thickened epithelial structure similar to the apical ectodermal ridge. Finally, we demonstrate that the induced tissues can contribute to endogenous digit tissue after transplantation. This PSC technology offers a first step for creating an artificial limb bud in culture and might open the door to inducing other mesenchymal/epithelial complex tissues from PSCs.

Tissue-level characterization of deformation dynamics is crucial for understanding organ morphogenetic mechanisms, especially the interhierarchical links among molecular activities, cellular behaviors and tissue/organ morphogenetic processes. Limb development is a well-studied topic in vertebrate organogenesis. Nevertheless, there is still little understanding of tissue-level deformation relative to molecular and cellular dynamics. This is mainly because live recording of detailed cell behaviors in whole tissues is technically difficult. To overcome this limitation, by applying a recently developed Bayesian approach, we here constructed tissue deformation maps for chick limb development with high precision, based on snapshot lineage tracing using dye injection. The precision of the constructed maps was validated with a clear statistical criterion. From the geometrical analysis of the map, we identified three characteristic tissue growth modes in the limb and showed that they are consistent with local growth factor activity and cell cycle length. In particular, we report that SHH signaling activity changes dynamically with developmental stage and strongly correlates with the dynamic shift in the tissue growth mode. We also found anisotropic tissue deformation along the proximal-distal axis. Morphogenetic simulation and experimental studies suggested that this directional tissue elongation, and not local growth, has the greatest impact on limb shaping. This result was supported by the novel finding that anisotropic tissue elongation along the proximal-distal axis occurs independently of cell proliferation. Our study marks a pivotal point for multi-scale system understanding in vertebrate development.

To understand the morphogenetic mechanisms of organ development and regeneration, it is essential to clarify the inter-hierarchical relationship between microscopic, molecular/cellular activities and organ-level tissue deformation dynamics. While the former have been studied for several decades, the latter - macroscopic geometrical information about physical tissue deformation - is often missing, especially for many vertebrates. This is mainly because live recording of detailed cell behaviors in whole tissues during vertebrate organogenesis is technically difficult. In this study, we have developed a novel method that combines snapshot lineage tracing with Bayesian statistical estimation to construct whole-organ deformation maps from landmark data on limited numbers of space-time points. Following the validation of the method using artificially generated data sets, we applied it to the analysis of tissue deformation dynamics in chick limb development. A quantitative tissue deformation map for St.23-St.24 has been constructed, and its precision has been proven by evaluating its predictive performance. Geometrical analyses of the map have revealed a spatially heterogeneous volume growth pattern that is consistent with the expression pattern of a major morphogen and anisotropic tissue deformation along an axis. Thus, our method enables deformation dynamics analysis in organogenesis using practical lineage marking techniques.