Homozygous inactivating mutations in the NKX3-2 gene result in spondylo-megaepiphyseal-metaphyseal dysplasia

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

Hemifacial microsomia is linked to a rare homozygous variant V162I in FRK and validated in zebrafish

OBJECTIVES

Hemifacial microsomia (HFM) is a common birth defect involving the first and second branchial arch derivatives. Although several chromosomal abnormalities and causal gene variants have been identified, genetic etiologies in a majority of cases with HFM remain unknown. This study aimed to identify genetic mutations in affected individuals with HFM.

METHODS

Whole-exome sequencing and bioinformatics analysis were performed for 16 affected individuals and their family members. Sanger sequencing was applied for confirmation of selected mutations. Zebrafish embryos were used for in situ hybridization of candidate gene, microinjection with antisense morpholino, and cartilage staining.

RESULTS

A homozygous missense mutation (c.484G > A; p.V162I) in the FRK gene was identified in an 18-year-old girl with HFM and dental abnormalities. Heterozygous mutation of this mutation was identified in her parents, who are first cousins in a consanguineous family. FRK is highly expressed in the Meckel's cartilage during embryonic development in mouse and zebrafish. Knockdown of frk in zebrafish showed a lower length and width ratio of Meckel's cartilage, abnormal mandibular jaw joint, and disorganized ceratobranchial cartilage and bone.

CONCLUSIONS

We identified a recessive variant in the FRK gene as a novel candidate gene for a patient with HFM and mandibular hypoplasia and revealed its effects on craniofacial and embryonic development in zebrafish.

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.

Chondrocyte Hypertrophy in Osteoarthritis: Mechanistic Studies and Models for the Identification of New Therapeutic Strategies

Articular cartilage shows limited self-healing ability owing to its low cellularity and avascularity. Untreated cartilage defects display an increased propensity to degenerate, leading to osteoarthritis (OA). During OA progression, articular chondrocytes are subjected to significant alterations in gene expression and phenotype, including a shift towards a hypertrophic-like state (with the expression of collagen type X, matrix metalloproteinases-13, and alkaline phosphatase) analogous to what eventuates during endochondral ossification. Present OA management strategies focus, however, exclusively on cartilage inflammation and degradation. A better understanding of the hypertrophic chondrocyte phenotype in OA might give new insights into its pathogenesis, suggesting potential disease-modifying therapeutic approaches. Recent developments in the field of cellular/molecular biology and tissue engineering proceeded in the direction of contrasting the onset of this hypertrophic phenotype, but knowledge gaps in the cause-effect of these processes are still present. In this review we will highlight the possible advantages and drawbacks of using this approach as a therapeutic strategy while focusing on the experimental models necessary for a better understanding of the phenomenon. Specifically, we will discuss in brief the cellular signaling pathways associated with the onset of a hypertrophic phenotype in chondrocytes during the progression of OA and will analyze in depth the advantages and disadvantages of various models that have been used to mimic it. Afterwards, we will present the strategies developed and proposed to impede chondrocyte hypertrophy and cartilage matrix mineralization/calcification. Finally, we will examine the future perspectives of OA therapeutic strategies.

Torticollis in Connection with Spine Phenotype

Purpose: Torticollis is not of uncommon occurrence in orthopaedic departments. Various theories and studies concerning the pathogenesis of the deformity have been suggested. We aimed to highlight and discuss the underlying cervical and spine malformation complex in correlation with torticollis via radiographic and tomographic analysis and its connection with a specific syndromic entity. Methods: Torticollis has been recognised in six patients (2 boys and 4 girls with an age range of 14–18 years), in addition to a couple of parents manifested persistent backpain. A variable spine malformation complex was the main reason behind torticollis. In addition, some patients manifested plagiocephaly, facial asymmetry and scoliosis/kyphoscoliosis. In some patients, conventional radiographs were of limited value because of the overlapping anatomical structures. Three-dimensional reconstruction CT scanning was the modality of choice, which enlightens the path for the phenotypic characterisation. Results: A 16-year-old-boy presented with torticollis in correlation with pathologic aberration of the spine cartilaginous stage was analysed via 3DCT scan. Comprehensive clinical and radiological phenotypes were in favour of spondylomegepiphyseal dysplasia. The genotype showed a mutation of the NKX3-2 (BAPX1) gene compatible with the diagnosis of spondylo-meg-epiphyseal-metaphyseal dysplasia. His younger male sibling and parents were heterozygous carriers. In two patients with pseudoachondroplasia syndrome, in which odontoid hypoplasia associated with cervical spine synchondrosis causing life-threatening torticollis, Cartilage oligomeric matrix protein (COMP) gene mutation was identified. MURCS syndrome has been diagnosed in two unrelated girls. Torticollis associated with cervical kyphosis was the major presentation since early childhood. Interestingly, one girl showed omovertebral bones of the lower cervical and upper thoracic spine. Her karyotype manifested a balanced translocation of 46 XX, t (14q; 15q). Conclusion: To detect the underlying etiological diagnosis of torticollis, a skeletal survey was the primary diagnostic tool. Conventional radiographs of the craniocervical junction and spine resulted in confusing readings because of the overlapping anatomical structures. Cranio-cervical malformation complex could have serious neurological deficits, especially for children with indefinite diagnosis of torticollis. The widely used term of congenital muscular torticollis resulted in morbid or mortal consequences. Moreover, some patients received vigorous physical therapy on the bases of muscular torticollis. Sadly speaking, this resulted in grave complications. Understanding the imaging phenotype and the genotype in such patients is the baseline tool for precise and proper management. The value of this paper is to sensitise physicians and orthopaedic surgeons to the necessity of comprehensive clinical and radiological phenotypic characterisations in patients with long term skeletal pathology.

Single-cell RNA-seq reveals clonal diversity and prognostic genes of relapsed multiple myeloma

BACKGROUND

Multiple myeloma (MM) is a clinically and biologically heterogeneous plasma-cell malignancy. Despite extensive research, disease heterogeneity and relapse remain a big challenge in MM therapeutics. We tried to dissect this disease and identify novel biomarkers for patient stratification and treatment outcome prediction by applying single-cell technology.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) and variable-diversity-joining regions-targeted sequencing (scVDJ-seq) concurrently on bone marrow samples from a cohort of 18 patients with newly diagnosed MM (NDMM; n = 12) or refractory/relapsed MM (RRMM; n = 6). We analysed the malignant clonotypes using scVDJ-seq data and conducted data integration and cell-type annotation through the CCA algorithm based on gene expression profiling. Furthermore, we identified disease status-specific genes and modules by comparison of NDMM and RRMM datasets and explored the findings in a larger MM cohort from the MMRF CoMMpass study.

RESULTS

We found that all the myeloma cells in either diagnosed or relapsed samples were dominated by a major clone, with a few subclones in several samples (n = 5). Next, we investigated the universal transcriptional features of myeloma cells and identified eight meta-programs correlated with this disease, especially meta-programs 1 and 8 (M1 and M8), which were the most significant and related to cell cycle and stress response, respectively. Furthermore, we classified the malignant plasma cells into eight clusters and found that the cell numbers in clusters 2/6/7 were exclusively higher in relapsed samples. Besides, we identified several attractive candidates for biomarkers (e.g. SMAD1 and STMN1) associated with disease progression and relapse in our dataset and related to overall survival in the CoMMpass dataset.

CONCLUSIONS

Our data provide insights into the heterogeneity of MM as well as highlight the relevance of intra-tumour heterogeneity and discover novel biomarkers that might be a potent therapy.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

The broad role of Nkx3.2 in the development of the zebrafish axial skeleton

The transcription factor Nkx3.2 (Bapx1) is an important chondrocyte maturation inhibitor. Previous Nkx3.2 knockdown and overexpression studies in non-mammalian gnathostomes have focused on its role in primary jaw joint development, while the function of this gene in broader skeletal development is not fully described. We generated a mutant allele of nkx3.2 in zebrafish with CRISPR/Cas9 and applied a range of techniques to characterize skeletal phenotypes at developmental stages from larva to adult, revealing loss of the jaw joint, fusions in bones of the occiput, morphological changes in the Weberian apparatus, and the loss or deformation of bony elements derived from basiventral cartilages of the vertebrae. Axial phenotypes are reminiscent of Nkx3.2 knockout in mammals, suggesting that the function of this gene in axial skeletal development is ancestral to osteichthyans. Our results highlight the broad role of nkx3.2 in zebrafish skeletal development and its context-specific functions in different skeletal elements.

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.

Diagnostic Approach to Macrocephaly in Children

Macrocephaly affects up to 5% of the pediatric population and is defined as an abnormally large head with an occipitofrontal circumference (OFC) >2 standard deviations (SD) above the mean for a given age and sex. Taking into account that about 2-3% of the healthy population has an OFC between 2 and 3 SD, macrocephaly is considered as "clinically relevant" when OFC is above 3 SD. This implies the urgent need for a diagnostic workflow to use in the clinical setting to dissect the several causes of increased OFC, from the benign form of familial macrocephaly and the Benign enlargement of subarachnoid spaces (BESS) to many pathological conditions, including genetic disorders. Moreover, macrocephaly should be differentiated by megalencephaly (MEG), which refers exclusively to brain overgrowth, exceeding twice the SD (3SD-"clinically relevant" megalencephaly). While macrocephaly can be isolated and benign or may be the first indication of an underlying congenital, genetic, or acquired disorder, megalencephaly is most likely due to a genetic cause. Apart from the head size evaluation, a detailed family and personal history, neuroimaging, and a careful clinical evaluation are crucial to reach the correct diagnosis. In this review, we seek to underline the clinical aspects of macrocephaly and megalencephaly, emphasizing the main differential diagnosis with a major focus on common genetic disorders. We thus provide a clinico-radiological algorithm to guide pediatricians in the assessment of children with macrocephaly.

nkx3.2 mutant zebrafish accommodate jaw joint loss through a phenocopy of the head shapes of Paleozoic jawless fish

The vertebrate jaw is a versatile feeding apparatus. To function, it requires a joint between the upper and lower jaws, so jaw joint defects are often highly disruptive and difficult to study. To describe the consequences of jaw joint dysfunction, we engineered two independent null alleles of a single jaw joint marker gene, nkx3.2, in zebrafish. These mutations caused zebrafish to become functionally jawless via fusion of the upper and lower jaw cartilages (ankylosis). Despite lacking jaw joints, nkx3.2 mutants survived to adulthood and accommodated this defect by: (a) having a remodeled skull with a fixed open gape, reduced snout and enlarged branchial region; and (b) performing ram feeding in the absence of jaw-generated suction. The late onset and broad extent of phenotypic changes in the mutants suggest that modifications to the skull are induced by functional agnathia, secondarily to nkx3.2 loss of function. Interestingly, nkx3.2 mutants superficially resemble ancient jawless vertebrates (anaspids and furcacaudiid thelodonts) in overall head shape. Because no homology exists in individual skull elements between these taxa, the adult nkx3.2 phenotype is not a reversal but rather a convergence due to similar functional requirements of feeding without moveable jaws. This remarkable analogy strongly suggests that jaw movements themselves dramatically influence the development of jawed vertebrate skulls. Thus, these mutants provide a unique model with which to: (a) investigate adaptive responses to perturbation in skeletal development; (b) re-evaluate evolutionarily inspired interpretations of phenocopies generated by gene knockdowns and knockouts; and (c) gain insight into feeding mechanics of the extinct agnathans.

An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

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.

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.

Transcriptional control of chondrocyte specification and differentiation

A milestone in the evolutionary emergence of vertebrates was the invention of cartilage, a tissue that has key roles in modeling, protecting and complementing the bony skeleton. Cartilage is elaborated and maintained by chondrocytes. These cells derive from multipotent skeletal progenitors and they perform highly specialized functions as they proceed through sequential lineage commitment and differentiation steps. They form cartilage primordia, the primary skeleton of the embryo. They then transform these primordia either into cartilage growth plates, temporary drivers of skeletal elongation and endochondral ossification, or into permanent tissues, namely articular cartilage. Chondrocyte fate decisions and differentiated activities are controlled by numerous extrinsic and intrinsic cues, and they are implemented at the gene expression level by transcription factors. The latter are the focus of this review. Meritorious efforts from many research groups have led over the last two decades to the identification of dozens of key chondrogenic transcription factors. These regulators belong to all types of transcription factor families. Some have master roles at one or several differentiation steps. They include SOX9 and RUNX2/3. Others decisively assist or antagonize the activities of these masters. They include TWIST1, SOX5/6, and MEF2C/D. Many more have tissue-patterning roles and regulate cell survival, proliferation and the pace of cell differentiation. They include, but are not limited to, homeodomain-containing proteins and growth factor signaling mediators. We here review current knowledge of all these factors, one superclass, class, and family at a time. We then compile all knowledge into transcriptional networks. We also identify remaining gaps in knowledge and directions for future research to fill these gaps and thereby provide novel insights into cartilage disease mechanisms and treatment options.

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^-5) and 88 SNPs had suggestive (P < 10^-4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.

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 chondrocytic journey in endochondral bone growth and skeletal dysplasia

The endochondral bones of the skeleton develop from a cartilage template and grow via a process involving a cascade of chondrocyte differentiation steps culminating in formation of a growth plate and the replacement of cartilage by bone. This process of endochondral ossification, driven by the generation of chondrocytes and their subsequent proliferation, differentiation, and production of extracellular matrix constitute a journey, deviation from which inevitably disrupts bone growth and development, and is the basis of human skeletal dysplasias with a wide range of phenotypic severity, from perinatal lethality to progressively deforming. This highly coordinated journey of chondrocyte specification and fate determination is controlled by a myriad of intrinsic and extrinsic factors. SOX9 is the master transcription factor that, in concert with varying partners along the way, directs the different phases of the journey from mesenchymal condensation, chondrogenesis, differentiation, proliferation, and maturation. Extracellular signals, including bone morphogenetic proteins, wingless-related MMTV integration site (WNT), fibroblast growth factor, Indian hedgehog, and parathyroid hormone-related peptide, are all indispensable for growth plate chondrocytes to align and organize into the appropriate columnar architecture and controls their maturation and transition to hypertrophy. Chondrocyte hypertrophy, marked by dramatic volume increase in phases, is controlled by transcription factors SOX9, Runt-related transcription factor, and FOXA2. Hypertrophic chondrocytes mediate the cartilage to bone transition and concomitantly face a live-or-die situation, a subject of much debate. We review recent insights into the coordination of the phases of the chondrocyte journey, and highlight the need for a systems level understanding of the regulatory networks that will facilitate the development of therapeutic approaches for skeletal dysplasia.

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.

Genome-wide linkage study suggests a susceptibility locus for isolated bilateral microtia on 4p15.32-4p16.2

Microtia is a congenital deformity where the external ear is underdeveloped. Genetic investigations have identified many susceptibility genes of microtia-related syndromes. However, no causal genes were reported for isolated microtia, the main form of microtia. We conducted a genome-wide linkage analysis on a 5-generation Chinese pedigree with isolated bilateral microtia. We identified a suggestive linkage locus on 4p15.32-4p16.2 with parametric LOD score of 2.70 and nonparametric linkage score (Zmean) of 12.28 (simulated occurrence per genome scan equal to 0.46 and 0.47, respectively). Haplotype reconstruction analysis of the 4p15.32-4p16.2 region further confined the linkage signal to a 10-Mb segment located between rs12505562 and rs12649803 (9.65-30.24 cM; 5.54-15.58 Mb). Various human organ developmental genes reside in this 10-Mb susceptibility region, such as EVC, EVC2, SLC2A9, NKX3-2, and HMX1. The coding regions of three genes, EVC known for cartilage development and NKX3-2, HMX1 involved in microtia, were selected for sequencing with 5 individuals from the pedigree. Of the 38 identified sequence variants, none segregates along with the disease phenotype. Other genes or DNA sequences of the 10-Mb region warrant for further investigation. In conclusion, we report a susceptibility locus of isolated microtia, and this finding will encourage future studies on the genetic basis of ear deformity.

Interplay of Nkx3.2, Sox9 and Pax3 regulates chondrogenic differentiation of muscle progenitor cells

Muscle satellite cells make up a stem cell population that is capable of differentiating into myocytes and contributing to muscle regeneration upon injury. In this work we investigate the mechanism by which these muscle progenitor cells adopt an alternative cell fate, the cartilage fate. We show that chick muscle satellite cells that normally would undergo myogenesis can be converted to express cartilage matrix proteins in vitro when cultured in chondrogenic medium containing TGFß3 or BMP2. In the meantime, the myogenic program is repressed, suggesting that muscle satellite cells have undergone chondrogenic differentiation. Furthermore, ectopic expression of the myogenic factor Pax3 prevents chondrogenesis in these cells, while chondrogenic factors Nkx3.2 and Sox9 act downstream of TGFß or BMP2 to promote this cell fate transition. We found that Nkx3.2 and Sox9 repress the activity of the Pax3 promoter and that Nkx3.2 acts as a transcriptional repressor in this process. Importantly, a reverse function mutant of Nkx3.2 blocks the ability of Sox9 to both inhibit myogenesis and induce chondrogenesis, suggesting that Nkx3.2 is required for Sox9 to promote chondrogenic differentiation in satellite cells. Finally, we found that in an in vivo mouse model of fracture healing where muscle progenitor cells were lineage-traced, Nkx3.2 and Sox9 are significantly upregulated while Pax3 is significantly downregulated in the muscle progenitor cells that give rise to chondrocytes during fracture repair. Thus our in vitro and in vivo analyses suggest that the balance of Pax3, Nkx3.2 and Sox9 may act as a molecular switch during the chondrogenic differentiation of muscle progenitor cells, which may be important for fracture healing.

Nkx3.2 promotes primary chondrogenic differentiation by upregulating Col2a1 transcription

BACKGROUND

The Nkx3.2 transcription factor promotes chondrogenesis by forming a positive regulatory loop with a crucial chondrogenic transcription factor, Sox9. Previous studies have indicated that factors other than Sox9 may promote chondrogenesis directly, but these factors have not been identified. Here, we test the hypothesis that Nkx3.2 promotes chondrogenesis directly by Sox9-independent mechanisms and indirectly by previously characterized Sox9-dependent mechanisms.

METHODOLOGY/PRINCIPAL FINDINGS

C3H10T1/2 pluripotent mesenchymal cells were cultured with bone morphogenetic protein 2 (BMP2) to induce endochondral ossification. Overexpression of wild-type Nkx3.2 (WT-Nkx3.2) upregulated glycosaminoglycan (GAG) production and expression of type II collagen α1 (Col2a1) mRNA, and these effects were evident before WT-Nkx3.2-mediated upregulation of Sox9. RNAi-mediated inhibition of Nkx3.2 abolished GAG production and expression of Col2a1 mRNA. Dual luciferase reporter assays revealed that WT-Nkx3.2 upregulated Col2a1 enhancer activity in a dose-dependent manner in C3H10T1/2 cells and also in N1511 chondrocytes. In addition, WT-Nkx3.2 partially restored downregulation of GAG production, Col2 protein expression, and Col2a1 mRNA expression induced by Sox9 RNAi. ChIP assays revealed that Nkx3.2 bound to the Col2a1 enhancer element.

CONCLUSIONS/SIGNIFICANCE

Nkx3.2 promoted primary chondrogenesis by two mechanisms: Direct and Sox9-independent upregulation of Col2a1 transcription and upregulation of Sox9 mRNA expression under positive feedback system.