Zebrafish model for spondylo-megaepiphyseal-metaphyseal dysplasia reveals post-embryonic roles of Nkx3.2 in the skeleton.

The regulated expansion of chondrocytes within growth plates and joints ensures proper skeletal development through adulthood. Mutations in the transcription factor NKX3.2 underlie spondylo-megaepiphyseal-metaphyseal dysplasia (SMMD), which is characterized by skeletal defects including scoliosis, large epiphyses, wide growth plates and supernumerary distal limb joints. Whereas nkx3.2 knockdown zebrafish and mouse Nkx3.2 mutants display embryonic lethal jaw joint fusions and skeletal reductions, respectively, they lack the skeletal overgrowth seen in SMMD patients. Here, we report adult viable nkx3.2 mutant zebrafish displaying cartilage overgrowth in place of a missing jaw joint, as well as severe dysmorphologies of the facial skeleton, skullcap and spine. In contrast, cartilage overgrowth and scoliosis are absent in rare viable nkx3.2 knockdown animals that lack jaw joints, supporting post-embryonic roles for Nkx3.2. Single-cell RNA-sequencing and in vivo validation reveal increased proliferation and upregulation of stress-induced pathways, including prostaglandin synthases, in mutant chondrocytes. By generating a zebrafish model for the skeletal overgrowth defects of SMMD, we reveal post-embryonic roles for Nkx3.2 in dampening proliferation and buffering the stress response in joint-associated chondrocytes.

Integrated Analysis of Transcriptome Expression Profiles Reveals miRNA-326-NKX3.2-Regulated Porcine Chondrocyte Differentiation.

The porcine body length trait is an essential factor affecting meat production and reproductive performance. It is evident that the development/lengthening of individual vertebrae is one of the main reasons for increases in body length; however, the underlying molecular mechanism remains unclear. In this study, RNA-seq analysis was used to profile the transcriptome (lncRNA, mRNA, and miRNA) of the thoracic intervertebral cartilage (TIC) at two time points (1 and 4 months) during vertebral column development in Yorkshire (Y) and Wuzhishan pigs (W). There were four groups: 1- (Y1) and 4-month-old (Y4) Yorkshire pigs and 1- (W1) and 4-month-old (W4) Wuzhishan pigs. In total, 161, 275, 86, and 126 differentially expressed (DE) lncRNAs, 1478, 2643, 404, and 750 DE genes (DEGs), and 74,51, 34, and 23 DE miRNAs (DE miRNAs) were identified in the Y4 vs. Y1, W4 vs. W1, Y4 vs. W4, and Y1 vs. W1 comparisons, respectively. Functional analysis of these DE transcripts (DETs) demonstrated that they had participated in various biological processes, such as cellular component organization or biogenesis, the developmental process, the metabolic process, bone development, and cartilage development. The crucial bone development-related candidate genes NK3 Homeobox 2 (NKX3.2), Wnt ligand secretion mediator (WLS), gremlin 1 (GREM1), fibroblast growth factor receptor 3 (FGFR3), hematopoietically expressed homeobox (HHEX), (collagen type XI alpha 1 chain (COL11A1), and Wnt Family Member 16 (WNT16)) were further identified by functional analysis. Moreover, lncRNA, miRNA, and gene interaction networks were constructed; a total of 55 lncRNAs, 6 miRNAs, and 7 genes formed lncRNA-gene, miRNA-gene, and lncRNA-miRNA-gene pairs, respectively. The aim was to demonstrate that coding and non-coding genes may co-regulate porcine spine development through interaction networks. NKX3.2 was identified as being specifically expressed in cartilage tissues, and it delayed chondrocyte differentiation. miRNA-326 regulated chondrocyte differentiation by targeting NKX3.2. The present study provides the first non-coding RNA and gene expression profiles in the porcine TIC, constructs the lncRNA-miRNA-gene interaction networks, and confirms the function of NKX3.2 in vertebral column development. These findings contribute to the understanding of the potential molecular mechanisms regulating pig vertebral column development. They expand our knowledge about the differences in body length between different pig species and provide a foundation for future studies.

Comprehensive Analysis of NKX3.2 in Liver Hepatocellular Carcinoma by Bigdata.

Background and Objectives: The gene NKX3.2 plays a role in determining cell fate during development, and mutations of NKX3.2 have been studied in relation to human skeletal diseases. However, due to the lack of studies on the link between NKX3.2 and cancer, we aimed to provide insights into NKX3.2 as a new prognostic biomarker for liver hepatocellular carcinoma (LIHC). Materials and Methods: The clinical significance of LIHC was investigated using open gene expression databases. We comprehensively analyzed NKX3.2 expression in LIHC using Gene Expression Profiling Interactive Analysis 2, Tumor Immune Estimation Resource (TIMER), and Kaplan-Meier plotter databases. Then, we investigated the association between NKX3.2 expression and tumor-infiltrating immune cells (TIICs). Results: NKX3.2 expression was higher in the primary tumor group compared to the normal group, and expression was higher in fibrolamellar carcinoma (FLC) compared to other subtypes. When the prognostic value of NKX3.2 was evaluated, highly expressed NKX3.2 significantly improved the overall survival and had an unfavorable prognosis. In addition, NKX3.2 expression was associated with immune cell infiltration. Patients with low gene expression and high macrophage expression had a poorer survival rate than those with low NKX3.2 and low macrophage expression (p = 0.0309). Conclusions: High NKX3.2 expression may induce poorer prognosis in LIHC. In addition, these findings can be used as basic data due to the lack of available related research. However, further in vivo studies are essential to gain a deeper understanding of the biological role of NKX3.2 in LIHC and its potential implications for cancer development and progression.

Suppression of Osteoarthritis progression by post-natal Induction of Nkx3.2.

Osteoarthritis (OA) is an incurable joint disease affecting 240 million elderly population, and major unmet medical needs exist for better therapeutic options for OA. During skeletal development, Nkx3.2 has been shown to promote chondrocyte differentiation and survival, but to suppress cartilage hypertrophy and blood vessel invasion. Here we show that Nkx3.2 plays a key role in osteoarthritis (OA) pathogenesis. Marked reduction of Nkx3.2 expression was observed in three different murine OA models. Consistent with these findings, analyses of surgery-induced and age-driven OA models revealed that cartilage-specific post-natal induction of Nkx3.2 can suppress OA progression in mice. These results suggest that Nkx3.2 may serve as a promising target for OA drug development.

Identification and in silico modeling of enhancers reveals new features of the cardiac differentiation network.

Developmental patterning and tissue formation are regulated through complex gene regulatory networks (GRNs) driven through the action of transcription factors (TFs) converging on enhancer elements. Here, as a point of entry to dissect the poorly defined GRN underlying cardiomyocyte differentiation, we apply an integrated approach to identify active enhancers and TFs involved in Drosophila heart development. The Drosophila heart consists of 104 cardiomyocytes, representing less than 0.5% of all cells in the embryo. By modifying BiTS-ChIP for rare cells, we examined H3K4me3 and H3K27ac chromatin landscapes to identify active promoters and enhancers specifically in cardiomyocytes. These in vivo data were complemented by a machine learning approach and extensive in vivo validation in transgenic embryos, which identified many new heart enhancers and their associated TF motifs. Our results implicate many new TFs in late stages of heart development, including Bagpipe, an Nkx3.2 ortholog, which we show is essential for differentiated heart function.

BMP-mediated induction of GATA4/5/6 blocks somitic responsiveness to SHH.

The relative timing of SHH and BMP signals controls whether presomitic mesoderm (PSM) cells will adopt either a chondrogenic or lateral plate mesoderm fate. Here we document that SHH-mediated induction of Nkx3.2 maintains the competence of somitic cells to initiate chondrogenesis in response to subsequent BMP signals by repressing BMP-dependent induction of GATA genes. Conversely, administration of BMP signals to PSM or forced expression of GATA family members in chick PSM explants blocks induction of hedgehog-dependent gene expression. We demonstrate that GATA factors can interact with Gli factors and can recruit the transcriptional co-factor FOG1 (ZFPM1) to the regulatory region of the mouse Gli1 gene, repressing the induction of Gli1 by SHH by binding to both GATA and Gli binding sites. Knockdown of FOG1 reverses the ability of GATA factors to repress Gli1 expression. Our findings uncover a novel role for GATA transcription factors as repressors of hedgehog signaling, and document that NKX3.2 maintains the ability of sclerotomal cells to express SHH transcriptional targets in the presence of BMP signals by repressing the induction of Gata4/5/6.

Pax1 acts as a negative regulator of chondrocyte maturation.

Paired box gene 1 (Pax1) indirectly promotes the early stages of chondrogenic differentiation through induction and transactivation of Nk3 homeobox 2 (Nkx3.2), a transcriptional repressor. Later in chondrogenic differentiation, Nkx3.2 blocks chondrocyte hypertrophy by repressing Runt-related transcription factor 2 (Runx2). Here we report the inhibitory action of Pax1 on chondrocyte maturation, independently of Nkx3.2. Upon cartilage formation, Pax1 expression in the ventral sclerotome was gradually decreased except for the perichondrial region of the vertebral bodies and the intervertebral region, both of which express SRY-box containing gene 9 (Sox9). Forced expression of Pax1 in the chick forelimb resulted in the formation of shortened skeletal elements with a significant reduction of proteoglycans (PGs) accumulation in cartilage as well as a lack of the cortical bone formation and vascular invasion into the primary ossification center. Pax1-misexpressing chondrocytes exhibited aberrant cell morphology with a marked downregulation of Aggrecan (Agc1). Pax1-misexpressing cultured chondrocytes failed to accumulate cartilaginous PGs and became fibroblastic, in association with downregulation of the expression of Sox9, Nkx3.2, Indian hedgehog (Ihh), type II collagen (Col2a1), Chondromodulin-1 (Chm1), and Agc1. Accumulation of cartilaginous PGs in chondrocytes was also reduced by forced expression of Pax1 and Sox9. Thus, chondrocyte maturation driven by Sox9 is antagonized by Pax1 that is downregulated during chondrogenic differentiation.

Pseudogenization of NK3 homeobox 2 (Nkx3.2) in monotremes provides insight into unique gastric anatomy and physiology.

The enzymatic breakdown and regulation of food passage through the vertebrate antral stomach and pyloric sphincter (antropyloric region) is a trait conserved over 450 million years. Development of the structures involved is underpinned by a highly conserved signalling pathway involving the hedgehog, bone morphogenetic protein and Wingless/Int-1 (Wnt) protein families. Monotremes are one of the few vertebrate lineages where acid-based digestion has been lost, and this is consistent with the lack of genes for hydrochloric acid secretion and gastric enzymes in the genomes of the platypus (Ornithorhynchus anatinus) and short-beaked echidna (Tachyglossus aculeatus) . Furthermore, these species feature unique gastric phenotypes, both with truncated and aglandular antral stomachs and the platypus with no pylorus. Here, we explore the genetic underpinning of monotreme gastric phenotypes, investigating genes important in antropyloric development using the newest monotreme genomes (mOrnAna1.pri.v4 and mTacAcu1) together with RNA-seq data. We found that the pathway constituents are generally conserved, but surprisingly, NK3 homeobox 2 (Nkx3.2) was pseudogenized in both platypus and echidna. We speculate that the unique sequence evolution of Grem1 and Bmp4 sequences in the echidna lineage may correlate with their pyloric-like restriction and that the convergent loss of gastric acid and stomach size genotypes and phenotypes in teleost and monotreme lineages may be a result of eco-evolutionary dynamics. These findings reflect the effects of gene loss on phenotypic evolution and further elucidate the genetic control of monotreme stomach anatomy and physiology.

Enhanced contrast synchrotron X-ray microtomography for describing skeleton-associated soft tissue defects in zebrafish mutants.

Detailed histological analyses are desirable for zebrafish mutants that are models for human skeletal diseases, but traditional histological techniques are limited to two-dimensional thin sections with orientations highly dependent on careful sample preparation. On the other hand, techniques that provide three-dimensional (3D) datasets including µCT scanning are typically limited to visualizing the bony skeleton and lack histological resolution. We combined diffusible iodine-based contrast enhancement (DICE) and propagation phase-contrast synchrotron radiation micro-computed tomography (PPC-SRµCT) to image late larval and juvenile zebrafish, obtaining high-quality 3D virtual histology datasets of the mineralized skeleton and surrounding soft tissues. To demonstrate this technique, we used virtual histological thin sections and 3D segmentation to qualitatively and quantitatively compare wild-type zebrafish and nkx3.2 -/- mutants to characterize novel soft-tissue phenotypes in the muscles and tendons of the jaw and ligaments of the Weberian apparatus, as well as the sinus perilymphaticus associated with the inner ear. We could observe disrupted fiber organization and tendons of the adductor mandibulae and protractor hyoideus muscles associated with the jaws, and show that despite this, the overall muscle volumes appeared unaffected. Ligaments associated with the malformed Weberian ossicles were mostly absent in nkx3.2 -/- mutants, and the sinus perilymphaticus was severely constricted or absent as a result of the fused exoccipital and basioccipital elements. These soft-tissue phenotypes have implications for the physiology of nkx3.2 -/- zebrafish, and demonstrate the promise of DICE-PPC-SRµCT for histopathological investigations of bone-associated soft tissues in small-fish skeletal disease models and developmental studies more broadly.

How vertebrates got their bite.

A newly discovered enhancer region may have allowed vertebrates to evolve the ability to open and close their jaws.

A novel cis-regulatory element drives early expression of Nkx3.2 in the gnathostome primary jaw joint.

The acquisition of movable jaws was a major event during vertebrate evolution. The role of NK3 homeobox 2 (Nkx3.2) transcription factor in patterning the primary jaw joint of gnathostomes (jawed vertebrates) is well known, however knowledge about its regulatory mechanism is lacking. In this study, we report a proximal enhancer element of Nkx3.2 that is deeply conserved in most gnathostomes but undetectable in the jawless hagfish and lamprey. This enhancer is active in the developing jaw joint region of the zebrafish Danio rerio, and was thus designated as jaw joint regulatory sequence 1 (JRS1). We further show that JRS1 enhancer sequences from a range of gnathostome species, including a chondrichthyan and mammals, have the same activity in the jaw joint as the native zebrafish enhancer, indicating a high degree of functional conservation despite the divergence of cartilaginous and bony fish lineages or the transition of the primary jaw joint into the middle ear of mammals. Finally, we show that deletion of JRS1 from the zebrafish genome using CRISPR/Cas9 results in a significant reduction of early gene expression of nkx3.2 and leads to a transient jaw joint deformation and partial fusion. Emergence of this Nkx3.2 enhancer in early gnathostomes may have contributed to the origin and shaping of the articulating surfaces of vertebrate jaws.

A novel NKX3-2 mutation associated with perinatal lethal phenotype of spondylo-megaepiphyseal-metaphyseal dysplasia in a neonate.

Spondylo-megaepiphyseal-metaphyseal dysplasia (SMMD) is an autosomal recessive skeletal dysplasia, characterized by disproportionate short stature with a short and stiff neck and trunk. SMMD is caused by inactivating mutations in NKX3-2, which encodes a homeobox-containing protein. Because of the rarity of the disorder, the diagnostic feature has not been fully established yet. We describe an affected newborn with dysmorphic facial features and severe short trunk. The patient required immediate intubation at the delivery room and duodenal atresia was detected during his course in neonatal intensive care unit. Skeletal survey revealed total absence of the ossification of the vertebral bodies, pubis, and ischia. Mainly the femora was short and broad with mild flaring of the metaphyses. The downward sloping or tented appearance of the ribs was distinctive. A diagnosis of SMMD was made on clinical and radiological grounds. Molecular analysis revealed homozygosity for a novel mutation, c.507-508delCA (p.Gly171Cysfs*55) in exon 2 of NKX3-2. The patient was operated on postnatal day 7 for duodenal atresia. In the post-operative period he developed sepsis and respiratory failure and he died on postnatal day 14. Although no neuroradiologic imaging could be performed, the findings of clubfoot, neuromuscular respiratory insufficiency requiring invasive mechanical ventilation and downward sloping or tented appearance of the ribs were suggestive of very early cervical cord compression leading to perinatal mortality. To our knowledge this patient yet represents one of the most severe postnatal phenotypes of SMMD.

Bapx1 upregulation is associated with ectopic mandibular cartilage development in amphibians.

BACKGROUND: The emergence of novel structures during evolution is crucial for creating variation among organisms, but the underlying processes which lead to the emergence of evolutionary novelties are poorly understood. The gnathostome jaw joint is such a novelty, and the incorporation of bapx1 expression into the intermediate first pharyngeal arch may have played a major role in the evolution of this joint. Knockdown experiments revealed that loss of bapx1 function leads to the loss of the jaw joint, because Meckel's cartilage and the palatoquadrate fuse during development. We used Xenopus laevis and Ambystoma mexicanum to further investigate the function of bapx1 in amphibians. Bapx1 expression levels were upregulated through the use of Ly-294,002 and we investigated the potential consequences of the enhanced bapx1 expression in amphibians to test the hypothesized joint inducing function of bapx1. RESULTS: We show that Ly-294,002 upregulates bapx1 expression in vivo. Additionally, ectopic mandibular arch derived cartilages develop after Ly-294,002 treatment. These ectopic cartilages are dorsoventrally oriented rods situated lateral to the palatoquadrate. The development of these additional cartilages did not change the muscular arrangement of mandibular arch-derived muscles. CONCLUSIONS: Development of additional mandibular cartilages is not unusual in larval anurans. Therefore, changes in the bapx1 expression during evolution may have been the reason for the development of several additional cartilages in the larval anuran jaw. Furthermore, our observations imply a joint-promoting function of bapx1, which further substantiates its hypothetical role in the evolution of the gnathostome jaw joint.

Nkx3.2 induces oxygen concentration-independent and lysosome-dependent degradation of HIF-1α to modulate hypoxic responses in chondrocytes.

Hypoxia-inducible factor 1-alpha (HIF-1α) is a DNA-binding transcription factor regulating hypoxic responses. It plays a key role in vascularization and angiogenesis as well as various metabolic pathways. Interestingly, during early phase endochondral ossification when HIF expression in chondrocytes is evident, developing cartilage primordia remains avascular until hypertrophic calcification commences. In this work, we uncovered a novel pathway causing oxygen concentration-independent and proteasome-independent degradation of HIF-1α protein. In this pathway, Nkx3.2, a chondrogenic factor, in conjunction with CHIP E3 ligase and p62/SQSTM1 adaptor, induces HIF-1α degradation via a macroautophagy pathway in a hypoxic environment. Consistent with these findings, Nkx3.2 was capable of suppressing HIF-dependent reporter gene activity as well as endogenous HIF target genes in in vitro cell culture. Furthermore, we observed that cartilage-specific Nkx3.2 overexpression in mice attenuates HIF-1α protein levels as well as vascularization in cartilage growth plates. Therefore, these results suggest that Nkx3.2-mediated HIF regulation may allow cartilage-specific avascularity under hypoxic conditions during endochondral skeleton development.

A post-translational modification cascade employing HDAC9-PIASy-RNF4 axis regulates chondrocyte hypertrophy by modulating Nkx3.2 protein stability.

While Nkx3.2/Bapx1 promotes chondrogenic differentiation and plays a role in maintaining chondrocyte viability and suppressing chondrocyte hypertrophy, the regulatory mechanisms of Nkx3.2 remain poorly understood. Here we show that p300- and HDAC9-induced Nkx3.2 acetylation and de-acetylation, respectively, play critical roles in controlling Nkx3.2 protein stability. In addition, we also found that HDAC9-dependent de-acetylation of Nkx3.2 triggers PIASy-mediated sumoylation and subsequent RNF4-mediated SUMO-targeted ubiquitination. Furthermore, we demonstrate that Nkx3.2 regulation by HDAC9 can be linked to the management of chondrocyte survival and hypertrophic maturation during cartilage development. Finally, our results together reveal a novel mechanism of protein stability control involving complex interplay between acetylation, de-acetylation, sumoylation, and ubiquitination, and suggest that this post-translational modification of Nkx3.2 employing HDAC9-PIASy-RNF4 axis plays a crucial role in controlling chondrocyte viability and hypertrophic maturation during skeletal development in vertebrates.

Suppression of Nkx3.2 by phosphatidylinositol-3-kinase signaling regulates cartilage development by modulating chondrocyte hypertrophy.

Phosphatidylinositol-3-kinase (PI3K) is a key regulator of diverse biological processes including cell proliferation, migration, survival, and differentiation. While a role of PI3K in chondrocyte differentiation has been suggested, its precise mechanisms of action are poorly understood. Here we show that PI3K signaling can down-regulate Nkx3.2 at both mRNA and protein levels in various chondrocyte cultures in vitro. In addition, we have intriguingly found that p85ß, not p85α, is specifically employed as a regulatory subunit for PI3K-mediated Nkx3.2 suppression. Furthermore, we found that regulation of Nkx3.2 by PI3K requires Rac1-PAK1, but not Akt, signaling downstream of PI3K. Finally, using embryonic limb bud cultures, ex vivo long bone cultures, and p85ß knockout mice, we demonstrated that PI3K-mediated suppression of Nkx3.2 in chondrocytes plays a role in the control of cartilage hypertrophy during skeletal development in vertebrates.

The role of Nkx3.2 in chondrogenesis.

Transcription factor, Nkx3.2, is a member of the NK family of developmental genes and is expressed during embryogenesis in a variety of mammalian model organisms, including chicken and mouse. It was first identified in Drosophila as the Bagpipe (bap) gene, where it has been demonstrated to be essential during formation of the midgut musculature. However, mammalian homolog Nkx3.2 has been shown to play a significant role in axial and limb skeletogenesis; in particular, the human skeletal disease, spondylo-megaepiphyseal-metaphyseal dysplasia (SMMD), is associated with mutations of the Nkx3.2 gene. In this review, we highlight the role of Nkx3.2 during musculoskeletal development, with an emphasis on the factor's role in determining chondrogenic cell fate and its subsequent role in endochondral ossification and chondrocyte survival.