Application of 3D MAPs pipeline identifies the morphological sequence chondrocytes undergo and the regulatory role of GDF5 in this process

Single-cell morphometrics reveals T-box gene-dependent patterns of epithelial tension in the Second Heart field

The vertebrate heart tube extends by progressive addition of epithelial second heart field (SHF) progenitor cells from the dorsal pericardial wall. The interplay between epithelial mechanics and genetic mechanisms during SHF deployment is unknown. Here, we present a quantitative single-cell morphometric analysis of SHF cells during heart tube extension, including force inference analysis of epithelial stress. Joint spatial Principal Component Analysis reveals that cell orientation and stress direction are the main parameters defining apical cell morphology and distinguishes cells adjacent to the arterial and venous poles. Cell shape and mechanical forces display a dynamic relationship during heart tube formation. Moreover, while the T-box transcription factor Tbx1 is necessary for cell orientation towards the arterial pole, activation of Tbx5 in the posterior SHF correlates with the establishment of epithelial stress and SHF deletion of Tbx5 relaxes the progenitor epithelium. Integrating findings from cell-scale feature patterning and mechanical stress provides new insights into cardiac morphogenesis.

Limited column formation in the embryonic growth plate implies divergent growth mechanisms during pre- and postnatal bone development

Chondrocyte columns, which are a hallmark of growth plate architecture, play a central role in bone elongation. Columns are formed by clonal expansion following rotation of the division plane, resulting in a stack of cells oriented parallel to the growth direction. In this work, we analyzed hundreds of Confetti multicolor clones in growth plates of mouse embryos using a pipeline comprising 3D imaging and algorithms for morphometric analysis. Surprisingly, analysis of the elevation angles between neighboring pairs of cells revealed that most cells did not display the typical stacking pattern associated with column formation, implying incomplete rotation of the division plane. Morphological analysis revealed that although embryonic clones were elongated, they formed clusters oriented perpendicular to the growth direction. Analysis of growth plates of postnatal mice revealed both complex columns, composed of ordered and disordered cell stacks, and small, disorganized clusters located in the outer edges. Finally, correlation between the temporal dynamics of the ratios between clusters and columns and between bone elongation and expansion suggests that clusters may promote expansion, whereas columns support elongation. Overall, our findings support the idea that modulations of division plane rotation of proliferating chondrocytes determines the formation of either clusters or columns, a multifunctional design that regulates morphogenesis throughout pre- and postnatal bone growth. Broadly, this work provides a new understanding of the cellular mechanisms underlying growth plate activity and bone elongation during development.

Molecular mechanisms and genetic factors contributing to the developmental dysplasia of the hip

The most prevalent hip disease in neonates is developmental dysplasia of the hip (DDH). A timely and accurate diagnosis is required to provide the most effective treatment for pediatric patients with DDH. Heredity and gene variation have been the subject of increased attention and research worldwide as one of the factors contributing to the pathogenesis of DDH. Genome-wide association studies (GWAS), genome-wide linkage analyses (GWLA), and exome sequencing (ES) have identified variants in numerous genes and single-nucleotide polymorphisms (SNPs) as being associated with susceptibility to DDH in sporadic and DDH family patients. Furthermore, the DDH phenotype can be observed in animal models that exhibit susceptibility genes or loci, including variants in CX3CR1, KANSL1, and GDF5. The dentification of noncoding RNAs and de novo gene variants in patients with DDH-related syndrome has enhanced our understanding of the genes implicated in DDH. This article reviews the most recent molecular mechanisms and genetic factors that contribute to DDH.

Leveraging single cell multiomic analyses to identify factors that drive human chondrocyte cell fate

Cartilage plays a crucial role in skeletal development and function, and abnormal development contributes to genetic and age-related skeletal disease. To better understand how human cartilage develops in vivo , we jointly profiled the transcriptome and open chromatin regions in individual nuclei recovered from distal femurs at 2 fetal timepoints. We used these multiomic data to identify transcription factors expressed in distinct chondrocyte subtypes, link accessible regulatory elements with gene expression, and predict transcription factor-based regulatory networks that are important for growth plate or epiphyseal chondrocyte differentiation. We developed a human pluripotent stem cell platform for interrogating the function of predicted transcription factors during chondrocyte differentiation and used it to test NFATC2 . We expect new regulatory networks we uncovered using multiomic data to be important for promoting cartilage health and treating disease, and our platform to be a useful tool for studying cartilage development in vitro .

Statement of Significance

The identity and integrity of the articular cartilage lining our joints are crucial to pain-free activities of daily living. Here we identified a gene regulatory landscape of human chondrogenesis at single cell resolution, which is expected to open new avenues of research aimed at mitigating cartilage diseases that affect hundreds of millions of individuals world-wide.

Chondrocyte hypertrophy in the growth plate promotes stress anisotropy affecting long bone development through chondrocyte column formation

The length of long bones is determined by column formation of proliferative chondrocytes and subsequent chondrocyte hypertrophy in the growth plate during bone development. Despite the importance of mechanical loading in long bone development, the mechanical conditions of the cells within the growth plate, such as the stress field, remain unclear owing to the difficulty in investigating spatiotemporal changes within dynamically growing tissues. In this study, the mechanisms of longitudinal bone growth were investigated from a mechanical perspective through column formation of proliferative chondrocytes within the growth plate before secondary ossification center formation using continuum-based particle models (CbPMs). A one-factor model, which simply describes essential aspects of a biological signaling cascade regulating cell activities within the growth plate, was developed and incorporated into CbPM. Subsequently, the developmental process and maintenance of the growth plate structure and resulting bone morphogenesis were simulated. Thus, stress anisotropy in the proliferative zone that affects bone elongation through chondrocyte column formation was identified and found to be promoted by chondrocyte hypertrophy. These results provide further insights into the mechanical regulation of multicellular dynamics during bone development.

Compressive stress gradients direct mechanoregulation of anisotropic growth in the zebrafish jaw joint

Mechanical stimuli arising from fetal movements are critical factors underlying joint growth. Abnormal fetal movements negatively affect joint shape features with important implications for joint health, but the mechanisms by which mechanical forces from fetal movements influence joint growth are still unclear. In this research, we quantify zebrafish jaw joint growth in 3D in free-to-move and immobilised fish larvae between four and five days post fertilisation. We found that the main changes in size and shape in normally moving fish were in the ventrodorsal axis, while growth anisotropy was lost in the immobilised larvae. We next sought to determine the cell level activities underlying mechanoregulated growth anisotropy by tracking individual cells in the presence or absence of jaw movements, finding that the most dramatic changes in growth rates due to jaw immobility were in the ventrodorsal axis. Finally, we implemented mechanobiological simulations of joint growth with which we tested hypotheses relating specific mechanical stimuli to mechanoregulated growth anisotropy. Different types of mechanical stimulation were incorporated into the simulation to provide the mechanoregulated component of growth, in addition to the baseline (non-mechanoregulated) growth which occurs in the immobilised animals. We found that when average tissue stress over the opening and closing cycle of the joint was used as the stimulus for mechanoregulated growth, joint morphogenesis was not accurately predicted. Predictions were improved when using the stress gradients along the rudiment axes (i.e., the variation in magnitude of compression to magnitude of tension between local regions). However, the most accurate predictions were obtained when using the compressive stress gradients (i.e., the variation in compressive stress magnitude) along the rudiment axes. We conclude therefore that the dominant biophysical stimulus contributing to growth anisotropy during early joint development is the gradient of compressive stress experienced along the growth axes under cyclical loading.

Self-organized BMP signaling dynamics underlie the development and evolution of digit segmentation patterns in birds and mammals

Repeating patterns of synovial joints are a highly conserved feature of articulated digits, with variations in joint number and location resulting in diverse digit morphologies and limb functions across the tetrapod clade. During the development of the amniote limb, joints form iteratively within the growing digit ray, as a population of distal progenitors alternately specifies joint and phalanx cell fates to segment the digit into distinct elements. While numerous molecular pathways have been implicated in this fate choice, it remains unclear how they give rise to a repeating pattern. Here, using single-cell RNA sequencing and spatial gene expression profiling, we investigate the transcriptional dynamics of interphalangeal joint specification in vivo. Combined with mathematical modeling, we predict that interactions within the BMP signaling pathway-between the ligand GDF5, the inhibitor NOGGIN, and the intracellular effector pSMAD-result in a self-organizing Turing system that forms periodic joint patterns. Our model is able to recapitulate the spatiotemporal gene expression dynamics observed in vivo, as well as phenocopy digit malformations caused by BMP pathway perturbations. By contrasting in silico simulations with in vivo morphometrics of two morphologically distinct digits, we show how changes in signaling parameters and growth dynamics can result in variations in the size and number of phalanges. Together, our results reveal a self-organizing mechanism that underpins amniote digit segmentation and its evolvability and, more broadly, illustrate how Turing systems based on a single molecular pathway may generate complex repetitive patterns in a wide variety of organisms.

Roles of Local Soluble Factors in Maintaining the Growth Plate: An Update

The growth plate is a cartilaginous tissue found at the ends of growing long bones, which contributes to the lengthening of bones during development. This unique structure contains at least three distinctive layers, including resting, proliferative, and hypertrophic chondrocyte zones, maintained by a complex regulatory network. Due to its soft tissue nature, the growth plate is the most susceptible tissue of the growing skeleton to injury in childhood. Although most growth plate damage in fractures can heal, some damage can result in growth arrest or disorder, impairing leg length and resulting in deformity. In this review, we re-visit previously established knowledge about the regulatory network that maintains the growth plate and integrate current research displaying the most recent progress. Next, we highlight local secretary factors, such as Wnt, Indian hedgehog (Ihh), and parathyroid hormone-related peptide (PTHrP), and dissect their roles and interactions in maintaining cell function and phenotype in different zones. Lastly, we discuss future research topics that can further our understanding of this unique tissue. Given the unmet need to engineer the growth plate, we also discuss the potential of creating particular patterns of soluble factors and generating them in vitro.

Growth orientations, rather than heterogeneous growth rates, dominate jaw joint morphogenesis in the larval zebrafish

In early limb embryogenesis, synovial joints acquire specific shapes which determine joint motion and function. The process by which the opposing cartilaginous joint surfaces are moulded into reciprocal and interlocking shapes, called joint morphogenesis, is one of the least understood aspects of joint formation and the cell-level dynamics underlying it are yet to be unravelled. In this research, we quantified key cellular dynamics involved in growth and morphogenesis of the zebrafish jaw joint and synthesised them in a predictive computational simulation of joint development. Cells in larval zebrafish jaw joints labelled with cartilage markers were tracked over a 48-h time window using confocal imaging. Changes in distance and angle between adjacent cell centroids resulting from cell rearrangement, volume expansion and extracellular matrix (ECM) deposition were measured and used to calculate the rate and direction of local tissue deformations. We observed spatially and temporally heterogeneous growth patterns with marked anisotropy over the developmental period assessed. There was notably elevated growth at the level of the retroarticular process of the Meckel's cartilage, a feature known to undergo pronounced shape changes during zebrafish development. Analysis of cell dynamics indicated a dominant role for cell volume expansion in growth, with minor influences from ECM volume increases and cell intercalation. Cell proliferation in the joint was minimal over the timeframe of interest. Synthesising the dynamic cell data into a finite element model of jaw joint development resulted in accurate shape predictions. Our biofidelic computational simulation demonstrated that zebrafish jaw joint growth can be reasonably approximated based on cell positional information over time, where cell positional information derives mainly from cell orientation and cell volume expansion. By modifying the input parameters of the simulation, we were able to assess the relative contributions of heterogeneous growth rates and of growth orientation. The use of uniform rather than heterogeneous growth rates only minorly impacted the shape predictions, whereas isotropic growth fields resulted in altered shape predictions. The simulation results suggest that growth anisotropy is the dominant influence on joint growth and morphogenesis. This study addresses the gap of the cellular processes underlying joint morphogenesis, with implications for understanding the aetiology of developmental joint disorders such as developmental dysplasia of the hip and arthrogryposis.