Hoxa11 and Hoxd11 regulate chondrocyte differentiation upstream of Runx2 and Shox2 in mice

The role of the PTHrP/Ihh feedback loop in the unusual growth plate location in mammalian metatarsals and pisiforms

BACKGROUND

Longitudinal skeletal growth takes place in the cartilaginous growth plates. While growth plates are found at either end of conventional long bones, they occur at a variety of locations in the mammalian skeleton. For example, the metacarpals and metatarsals (MT) in the hands and feet form only a single growth plate at one end, and the pisiform in the wrist is the only carpal bone to contain a growth plate. We take advantage of this natural anatomical variation to test which components of the PTHrP/Ihh feedback loop, a fundamental regulator of chondrocyte differentiation, are specific to growth plate function.

RESULTS

Parathyroid hormone-like hormone (Pthlh), the gene that transcribes parathyroid hormone-related peptide (PTHrP), is expressed in the reserve zone of the growth plate-forming end of the MT. At the opposite end, the absence of a PTHrP+ reserve zone results in premature chondrocyte differentiation and Indian hedgehog (Ihh) expression. Pthlh is expressed in the reserve zone of the developing pisiform, confirming the existence of a true growth plate.

CONCLUSION

A pool of PTHrP+ reserve zone chondrocytes is a defining characteristic of growth plates, and its patterning may be key to evolved differences in growth plate location in the mammalian skeleton.

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

Single-cell RNA sequencing analysis of the temporomandibular joint condyle in 3 and 4-month-old human embryos

BACKGROUND

The temporomandibular joint (TMJ) is a complex joint consisting of the condyle, the temporal articular surface, and the articular disc. Functions such as mastication, swallowing and articulation are accomplished by the movements of the TMJ. To date, the TMJ has been studied more extensively, but the types of TMJ cells, their differentiation, and their interrelationship during growth and development are still unclear and the study of the TMJ is limited. The aim of this study was to establish a molecular cellular atlas of the human embryonic temporomandibular joint condyle (TMJC) by single-cell RNA sequencing, which will contribute to understanding and solving clinical problems.

RESULTS

Human embryos at 3 and 4 months of age are an important stage of TMJC development. We performed a comprehensive transcriptome analysis of TMJC tissue from human embryos at 3 and 4 months of age using single-cell RNA sequencing. A total of 16,624 cells were captured and the gene expression profiles of 15 cell clusters in human embryonic TMJC were determined, including 14 known cell types and one previously unknown cell type, "transition state cells (TSCs)". Immunofluorescence assays confirmed that TSCs are not the same cell cluster as mesenchymal stem cells (MSCs). Pseudotime trajectory and RNA velocity analysis revealed that MSCs transformed into TSCs, which further differentiated into osteoblasts, hypertrophic chondrocytes and tenocytes. In addition, chondrocytes (CYTL1^high + THBS1^high) from secondary cartilage were detected only in 4-month-old human embryonic TMJC.

CONCLUSIONS

Our study provides an atlas of differentiation stages of human embryonic TMJC tissue cells, which will contribute to an in-depth understanding of the pathophysiology of the TMJC tissue repair process and ultimately help to solve clinical problems.

HOX genes in stem cells: Maintaining cellular identity and regulation of differentiation

Stem cells display a unique cell type within the body that has the capacity to self-renew and differentiate into specialized cell types. Compared to pluripotent stem cells, adult stem cells (ASC) such as mesenchymal stem cells (MSCs) and hematopoietic stem cells (HSCs) exhibit restricted differentiation capabilities that are limited to cell types typically found in the tissue of origin, which implicates that there must be a certain code or priming determined by the tissue of origin. HOX genes, a subset of homeobox genes encoding transcription factors that are generally repressed in undifferentiated pluripotent stem cells, emerged here as master regulators of cell identity and cell fate during embryogenesis, and in maintaining this positional identity throughout life as well as specifying various regional properties of respective tissues. Concurrently, intricate molecular circuits regulated by diverse stem cell-typical signaling pathways, balance stem cell maintenance, proliferation and differentiation. However, it still needs to be unraveled how stem cell-related signaling pathways establish and regulate ASC-specific HOX expression pattern with different temporal-spatial topography, known as the HOX code. This comprehensive review therefore summarizes the current knowledge of specific ASC-related HOX expression patterns and how these were integrated into stem cell-related signaling pathways. Understanding the mechanism of HOX gene regulation in stem cells may provide new ways to manipulate stem cell fate and function leading to improved and new approaches in the field of regenerative medicine.

Limb Mesoderm and Head Ectomesenchyme Both Express a Core Transcriptional Program During Chondrocyte Differentiation

To explain how cartilage appeared in different parts of the vertebrate body at discrete times during evolution, we hypothesize that different embryonic populations co-opted expression of a core gene regulatory network (GRN) driving chondrocyte differentiation. To test this hypothesis, laser-capture microdissection coupled with RNA-seq was used to reveal chondrocyte transcriptomes in the developing chick humerus and ceratobranchial, which are mesoderm- and neural crest-derived, respectively. During endochondral ossification, two general types of chondrocytes differentiate. Immature chondrocytes (IMM) represent the early stages of cartilage differentiation, while mature chondrocytes (MAT) undergo additional stages of differentiation, including hypertrophy and stimulating matrix mineralization and degradation. Venn diagram analyses generally revealed a high degree of conservation between chondrocyte transcriptomes of the limb and head, including SOX9, COL2A1, and ACAN expression. Typical maturation genes, such as COL10A1, IBSP, and SPP1, were upregulated in MAT compared to IMM in both limb and head chondrocytes. Gene co-expression network (GCN) analyses of limb and head chondrocyte transcriptomes estimated the core GRN governing cartilage differentiation. Two discrete portions of the GCN contained genes that were differentially expressed in limb or head chondrocytes, but these genes were enriched for biological processes related to limb/forelimb morphogenesis or neural crest-dependent processes, respectively, perhaps simply reflecting the embryonic origin of the cells. A core GRN driving cartilage differentiation in limb and head was revealed that included typical chondrocyte differentiation and maturation markers, as well as putative novel "chondrocyte" genes. Conservation of a core transcriptional program during chondrocyte differentiation in both the limb and head suggest that the same core GRN was co-opted when cartilage appeared in different regions of the skeleton during vertebrate evolution.

Role of lncRNAs into Mesenchymal Stromal Cell Differentiation

Currently, findings support that 75% of the human genome is actively transcribed, but only 2% is translated into a protein, according to databases such as ENCODE (Encyclopedia of DNA Elements) [1]. The development of high-throughput sequencing technologies, computational methods for genome assembly and biological models have led to the realization of the importance of the previously unconsidered non-coding fraction of the genome. Along with this, noncoding RNAs have been shown to be epigenetic, transcriptional and post-transcriptional regulators in a large number of cellular processes [2]. Within the group of non-coding RNAs, lncRNAs represent a fascinating field of study, given the functional versatility in their mode of action on their molecular targets. In recent years, there has been an interest in learning about lncRNAs in MSC differentiation. The aim of this review is to address the signaling mechanisms where lncRNAs are involved, emphasizing their role in either stimulating or inhibiting the transition to differentiated cell. Specifically, the main types of MSC differentiation are discussed: myogenesis, osteogenesis, adipogenesis and chondrogenesis. The description of increasingly new lncRNAs reinforces their role as players in the well-studied field of MSC differentiation, allowing a step towards a better understanding of their biology and their potential application in the clinic.

Stimulation of the Runx2 P1 promoter by collagen-derived dipeptide prolyl-hydroxyproline bound to Foxg1 and Foxo1 in osteoblasts

Collagen-derived dipeptide prolyl-hydroxyproline (Pro-Hyp) directly binds to the forkhead box g1 (Foxg1) protein and causes it to undergo structural alteration. Pro-Hyp also promotes the production of a regulator of osteoblast differentiation, Runt-related transcription factor 2 (Runx2), through Foxg1, inducing osteoblast differentiation. In addition, Pro-Hyp disrupts the interaction between Foxg1 and Runx2, and Foxg1 appears to interact with Runx2 in the absence of Pro-Hyp. To elucidate the mechanism of Pro-Hyp that promotes osteoblast differentiation, we investigated whether Pro-Hyp regulates the Runx2 P1 promoter together with Foxg1. The present study revealed that Pro-Hyp is taken up by osteoblastic cells via the solute carrier family 15 member (Slc15a) 4. In the presence of Pro-Hyp, Runx2 is translocated from the nucleus to the cytoplasm and Foxg1 is translocated from the cytoplasm to the nucleus. We also found that Pro-Hyp promoted the interaction between Forkhead box o1 (Foxo1) and Runx2 and the dissociation of Foxg1 from Runx2. Moreover, we identified the Pro-Hyp response element in the Runx2 distal P1 promoter at nt −375 to −316, including the Runx2 binding sites and Fox core sequence. In the presence of Pro-Hyp, Runx2 is dissociated from the Pro-Hyp response element in the Runx2 distal P1 promoter. Subsequently, Foxg1 and Foxo1 activated the Runx2 promoter by binding to the Pro-Hyp response element. In summary, we delineated the mechanism by which Pro-Hyp stimulates the bone-related Runx2 distal P1 promoter activity in osteoblastic cells through Foxg1, Foxo1, and Runx2.

The Embryological Landscape of Mayer-Rokitansky-Kuster-Hauser Syndrome: Genetics and Environmental Factors

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is a disorder caused by Müllerian ducts dysgenesis affecting 1 in 5000 women with a typical 46,XX karyotype. The etiology of MRKH syndrome is complex and largely unexplained. Familial clustering suggests a genetic component and the spectrum of clinical presentations seems consistent with an inheritance pattern characterized by incomplete penetrance and variable expressivity. Mutations of several candidate genes have been proposed as possible causes based on genetic analyses of human patients and animal models. In addition, studies of monozygotic twins with discordant phenotypes suggest a role for epigenetic changes following potential exposure to environmental compounds. The spectrum of clinical presentations is consistent with intricate disruptions of shared developmental pathways or signals during early organogenesis. However, the lack of functional validation and translational studies have limited our understanding of the molecular mechanisms involved in this condition. The clinical management of affected women, including early diagnosis, genetic testing of MRKH syndrome, and the implementation of counseling strategies, is significantly impeded by these knowledge gaps. Here, we illustrate the embryonic development of tissues and organs affected by MRKH syndrome, highlighting key pathways that could be involved in its pathogenesis. In addition, we will explore the genetics of this condition, as well as the potential role of environmental factors, and discuss their implications to clinical practice.

Molecular fingerprint of female bovine embryos produced in vitro with high competence to establish and maintain pregnancy†

The objective was to identify the transcriptomic profile of in vitro-derived embryos with high competence to establish and maintain gestation. Embryos produced with X-sorted sperm were cultured from day 5 to day 7 in serum-free medium containing 10 ng/ml recombinant bovine colony-stimulating factor 2 (CSF2) or vehicle. The CSF2 was administered because this molecule can increase blastocyst competence for survival after embryo transfer. Blastocysts were harvested on day 7 of culture and manually bisected. One demi-embryo from a single blastocyst was transferred into a synchronized recipient and the other half was used for RNA-seq analysis. Using P < 0.01 and a fold change >2-fold or <0.5 fold as cutoffs, there were 617 differentially expressed genes (DEG) between embryos that survived to day 30 of gestation vs those that did not, 470 DEG between embryos that survived to day 60 and those that did not, 432 DEG between embryos that maintained pregnancy from day 30 to day 60 vs those where pregnancy failed after day 30, and 635 DEG regulated by CSF2. Pathways and ontologies in which DEG were overrepresented included many related to cellular responses to stress and cell survival. It was concluded that gene expression in the blastocyst is different between embryos that are competent to establish and maintain pregnancy vs those that are not. The relationship between expression of genes related to cell stress and subsequent embryonic survival probably reflects cellular perturbations caused by embryonic development taking place in the artificial environment associated with cell culture.

Mesenchyme-specific loss of Dot1L histone methyltransferase leads to skeletal dysplasia phenotype in mice

Chromatin modifying enzymes play essential roles in skeletal development and bone maintenance, and deregulation of epigenetic mechanisms can lead to skeletal growth and malformation disorders. Here, we report a novel skeletal dysplasia phenotype in mice with conditional loss of Disruptor of telomeric silencing 1-like (Dot1L) histone methyltransferase in limb mesenchymal progenitors and downstream descendants. Phenotypic characterizations of mice with Dot1L inactivation by Prrx1-Cre (Dot1L-cKO^Prrx1) revealed limb shortening, abnormal bone morphologies, and forelimb dislocations. Our in vivo and in vitro data support a crucial role for Dot1L in regulating growth plate chondrocyte proliferation and differentiation, extracellular matrix production, and secondary ossification center formation. Micro-computed tomography analysis of femurs revealed that partial loss of Dot1L expression is sufficient to impair trabecular bone formation and microarchitecture in young mice. Moreover, RNAseq analysis of Dot1L deficient chondrocytes implicated Dot1L in the regulation of key genes and pathways necessary to promote cell cycle regulation and skeletal growth. Collectively, our data show that early expression of Dot1L in limb mesenchyme provides essential regulatory control of endochondral bone morphology, growth, and stability.

Hox genes maintain critical roles in the adult skeleton

Hox genes are indispensable for the proper patterning of the skeletal morphology of the axial and appendicular skeleton during embryonic development. Recently, it has been demonstrated that Hox expression continues from embryonic stages through postnatal and adult stages exclusively in a skeletal stem cell population. However, whether Hox genes continue to function after development has not been rigorously investigated. We generated a Hoxd11 conditional allele and induced genetic deletion at adult stages to show that Hox11 genes play critical roles in skeletal homeostasis of the forelimb zeugopod (radius and ulna). Conditional loss of Hox11 function at adult stages leads to replacement of normal lamellar bone with an abnormal woven bone-like matrix of highly disorganized collagen fibers. Examining the lineage from the Hox-expressing mutant cells demonstrates no loss of stem cell population. Differentiation in the osteoblast lineage initiates with Runx2 expression, which is observed similarly in mutants and controls. With loss of Hox11 function, however, osteoblasts fail to mature, with no progression to osteopontin or osteocalcin expression. Osteocyte-like cells become embedded within the abnormal bony matrix, but they completely lack dendrites, as well as the characteristic lacuno-canalicular network, and do not express SOST. Together, our studies show that Hox11 genes continuously function in the adult skeleton in a region-specific manner by regulating differentiation of Hox-expressing skeletal stem cells into the osteolineage.

lncRNAs: function and mechanism in cartilage development, degeneration, and regeneration

With the increasing incidence of cartilage-related diseases such as osteoarthritis (OA) and intervertebral disc degeneration (IDD), heavier financial and social burdens need to be faced. Unfortunately, there is no satisfactory clinical method to target the pathophysiology of cartilage-related diseases. Many gene expressions, signaling pathways, and biomechanical dysregulations were involved in cartilage development, degeneration, and regeneration. However, the underlying mechanism was not clearly understood. Recently, lots of long non-coding RNAs (lncRNAs) were identified in the biological processes, including cartilage development, degeneration, and regeneration. It is clear that lncRNAs were important in regulating gene expression and maintaining chondrocyte phenotypes and homeostasis. In this review, we summarize the recent researches studying lncRNAs’ expression and function in cartilage development, degeneration, and regeneration and illustrate the potential mechanism of how they act in the pathologic process. With continued efforts, regulating lncRNA expression in the cartilage regeneration may be a promising biological treatment approach.

Atoh8 acts as a regulator of chondrocyte proliferation and differentiation in endochondral bones

Atonal homolog 8 (Atoh8) is a transcription factor of the basic helix-loop-helix (bHLH) protein family, which is expressed in the cartilaginous elements of endochondral bones. To analyze its function during chondrogenesis we deleted Atoh8 in mice using a chondrocyte- (Atoh8flox/flox;Col2a1-Cre) and a germline- (Atoh8flox/flox;Prx1-Crefemale) specific Cre allele. In both strains, Atoh8 deletion leads to a reduced skeletal size of the axial and appendicular bones, but the stages of phenotypic manifestations differ. While we observed obviously shortened bones in Atoh8flox/flox;Col2a1-Cre mice only postnatally, the bones of Atoh8flox/flox;Prx1-Crefemale mice are characterized by a reduced bone length already at prenatal stages. Detailed histological and molecular investigations revealed reduced zones of proliferating and hypertrophic chondrocytes. In addition, Atoh8 deletion identified Atoh8 as a positive regulator of chondrocyte proliferation. As increased Atoh8 expression is found in the region of prehypertrophic chondrocytes where the expression of Ihh, a main regulator of chondrocyte proliferation and differentiation, is induced, we investigated a potential interaction of Atoh8 function and Ihh signaling. By activating Ihh signaling with Purmorphamine we demonstrate that Atoh8 regulates chondrocyte proliferation in parallel or downstream of Ihh signaling while it acts on the onset of hypertrophy upstream of Ihh likely by modulating Ihh expression levels.

Bone morphology is regulated modularly by global and regional genetic programs

Bone protrusions provide stable anchoring sites for ligaments and tendons and define the unique morphology of each long bone. Despite their importance, the mechanism by which superstructures are patterned is unknown. Here, we identify components of the genetic program that control the patterning of Sox9 ⁺/Scx ⁺ superstructure progenitors in mouse and show that this program includes both global and regional regulatory modules. Using light-sheet fluorescence microscopy combined with genetic lineage labeling, we mapped the broad contribution of the Sox9 ⁺/Scx ⁺ progenitors to the formation of bone superstructures. Then, by combining literature-based evidence, comparative transcriptomic analysis and genetic mouse models, we identified Gli3 as a global regulator of superstructure patterning, whereas Pbx1, Pbx2, Hoxa11 and Hoxd11 act as proximal and distal regulators, respectively. Moreover, by demonstrating a dose-dependent pattern regulation in Gli3 and Pbx1 compound mutations, we show that the global and regional regulatory modules work in a coordinated manner. Collectively, our results provide strong evidence for genetic regulation of superstructure patterning, which further supports the notion that long bone development is a modular process.This article has an associated 'The people behind the papers' interview.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Genome-wide DNA methylation profiling of articular cartilage reveals significant epigenetic alterations in Kashin-Beck disease and osteoarthritis

OBJECTIVE

To determine genome-wide DNA methylation profiles of knee cartilage from patients with Kashin-Beck disease (KBD) and osteoarthritis (OA).

METHOD

Knee cartilage was collected from 14 grade III KBD patients, 5 primary OA patients and 13 healthy subjects. The genome-wide methylation profiles of 5 KBD cartilage, 5 OA cartilage and 5 normal cartilage were determined by Illumina HumanMethylation450 array. Illumina Methylation Analyzer package was employed for identifying differentially methylated CpG sites. Functional annotation and enrichment analysis of differentially methylated genes (DMG) were conducted using GeneRIF database, Ingenuity Pathway Analysis (IPA) and The Database for Annotation, Visualization and Integrated Discovery (DAVID). Mass spectrometry (MS) and immunohistochemistry (IHC) were conducted to validate the functional relevance of identified KBD associated gene.

RESULTS

We identified a total of 1212 differentially methylated CpG sites in KBD vs Normal, annotated to 264 hypermethylated and 368 hypomethylated genes. Comparing the DNA methylation profiles of KBD vs Normal and OA vs Normal detected overlap of 367 differentially methylated CpG sites (annotated to 182 genes) as well as 845 KBD-specific differentially methylated CpG sites (annotated to 471 unique genes). MS and IHC confirmed the hypermethylation status and decreased protein expression of HAPLN1 gene in KBD cartilage.

CONCLUSION

Our data implicate epigenetic dysregulation of a host of genes in KBD and OA. Furthermore, we observed common causal epigenetic changes shared by KBD and OA.

Epigenetic inactivation of HOXA11, a novel functional tumor suppressor for renal cell carcinoma, is associated with RCC TNM classification

Epigenetic inactivation of HOXA11, a putative tumor suppressor, is frequently observed in a number of solid tumors, but has not been described in RCC (renal cell carcinoma). In this study, we investigated the expression, epigenetic changes and the function of HOXA11 in human renal cell carcinoma (RCC). HOXA11 was silenced or down-regulated in RCC cell lines and tissues. Methylation specific PCR (MSP) and bisulfite genomic sequencing (BGS) revealed that the HOXA11 promoter was hypermethylated in 5/6 RCC cell lines. Demethylation treatment resulted in demethylation of the promoter and increased HOXA11 expression in these cell lines. HOXA11 methylation was also detected in 68/95 (70.5%) primary RCC tumors, but only rare adjacent non-malignant renal tissues (13%, 3/23) showed hypermethylation of promoter. We also found that the methylation of HOXA11 was associated with higher TNM classification of RCC (p<0.05). Ectopic expression of HOXA11 led to significant inhibition of proliferation, colony formation, migration and invasion abilities and induced RCC cells apoptosis. Moreover, HOXA11 was found to inhibit Wnt signaling. Thus, our study demonstrated that HOXA11 function as a tumor suppressor in RCC, while it is frequently silenced by promoter methylation in RCC.

Hox11 Function Is Required for Region-Specific Fracture Repair

The processes that govern fracture repair rely on many mechanisms that recapitulate embryonic skeletal development. Hox genes are transcription factors that perform critical patterning functions in regional domains along the axial and limb skeleton during development. Much less is known about roles for these genes in the adult skeleton. We recently reported that Hox11 genes, which function in zeugopod development (radius/ulna and tibia/fibula), are also expressed in the adult zeugopod skeleton exclusively in PDGFRα+/CD51+/LepR+ mesenchymal stem/stromal cells (MSCs). In this study, we use a Hoxa11eGFP reporter allele and loss-of-function Hox11 alleles, and we show that Hox11 expression expands after zeugopod fracture injury, and that loss of Hox11 function results in defects in endochondral ossification and in the bone remodeling phase of repair. In Hox11 compound mutant fractures, early chondrocytes are specified but show defects in differentiation, leading to an overall deficit in the cartilage production. In the later stages of the repair process, the hard callus remains incompletely remodeled in mutants due, at least in part, to abnormal bone matrix organization. Overall, our data supports multiple roles for Hox11 genes following fracture injury in the adult skeleton. © 2017 American Society for Bone and Mineral Research.

A Track Record on SHOX: From Basic Research to Complex Models and Therapy

SHOX deficiency is the most frequent genetic growth disorder associated with isolated and syndromic forms of short stature. Caused by mutations in the homeobox gene SHOX, its varied clinical manifestations include isolated short stature, Léri-Weill dyschondrosteosis, and Langer mesomelic dysplasia. In addition, SHOX deficiency contributes to the skeletal features in Turner syndrome. Causative SHOX mutations have allowed downstream pathology to be linked to defined molecular lesions. Expression levels of SHOX are tightly regulated, and almost half of the pathogenic mutations have affected enhancers. Clinical severity of SHOX deficiency varies between genders and ranges from normal stature to profound mesomelic skeletal dysplasia. Treatment options for children with SHOX deficiency are available. Two decades of research support the concept of SHOX as a transcription factor that integrates diverse aspects of bone development, growth plate biology, and apoptosis. Due to its absence in mouse, the animal models of choice have become chicken and zebrafish. These models, therefore, together with micromass cultures and primary cell lines, have been used to address SHOX function. Pathway and network analyses have identified interactors, target genes, and regulators. Here, we summarize recent data and give insight into the critical molecular and cellular functions of SHOX in the etiopathogenesis of short stature and limb development.

Unique functions of Gata4 in mouse liver induction and heart development

Gata4 and Gata6 are closely related transcription factors that are essential for the development of a number of embryonic tissues. While they have nearly identical DNA-binding domains and similar patterns of expression, Gata4 and Gata6 null embryos have strikingly different embryonic lethal phenotypes. To determine whether the lack of redundancy is due to differences in protein function or Gata4 and Gata6 expression domains, we generated mice that contained the Gata6 cDNA in place of the Gata4 genomic locus. Gata4(Gata6/Gata6) embryos survived through embryonic day (E)12.5 and successfully underwent ventral folding morphogenesis, demonstrating that Gata6 is able to replace Gata4 function in extraembryonic tissues. Surprisingly, Gata6 is unable to replace Gata4 function in the septum transversum mesenchyme or the epicardium, leading to liver agenesis and lethal heart defects in Gata4(Gata6/Gata6) embryos. These studies suggest that Gata4 has evolved distinct functions in the development of these tissues that cannot be performed by Gata6, even when it is provided in the identical expression domain. Our work has important implications for the respective mechanisms of Gata function during development, as well as the functional evolution of these essential transcription factors.

LncRNA-HIT Functions as an Epigenetic Regulator of Chondrogenesis through Its Recruitment of p100/CBP Complexes

Gene expression profiling in E 11 mouse embryos identified high expression of the long noncoding RNA (lncRNA), LNCRNA-HIT in the undifferentiated limb mesenchyme, gut, and developing genital tubercle. In the limb mesenchyme, LncRNA-HIT was found to be retained in the nucleus, forming a complex with p100 and CBP. Analysis of the genome-wide distribution of LncRNA-HIT-p100/CBP complexes by ChIRP-seq revealed LncRNA-HIT associated peaks at multiple loci in the murine genome. Ontological analysis of the genes contacted by LncRNA-HIT-p100/CBP complexes indicate a primary role for these loci in chondrogenic differentiation. Functional analysis using siRNA-mediated reductions in LncRNA-HIT or p100 transcripts revealed a significant decrease in expression of many of the LncRNA-HIT-associated loci. LncRNA-HIT siRNA treatments also impacted the ability of the limb mesenchyme to form cartilage, reducing mesenchymal cell condensation and the formation of cartilage nodules. Mechanistically the LncRNA-HIT siRNA treatments impacted pro-chondrogenic gene expression by reducing H3K27ac or p100 activity, confirming that LncRNA-HIT is essential for chondrogenic differentiation in the limb mesenchyme. Taken together, these findings reveal a fundamental epigenetic mechanism functioning during early limb development, using LncRNA-HIT and its associated proteins to promote the expression of multiple genes whose products are necessary for the formation of cartilage.

A fundamental problem studied by skeletal biologists is the development of regenerative therapies to replace cartilage tissues impacted by injury or disease, which for individuals affected by osteoarthritis represents nearly half of all of all adults over the age of sixty five. To date, no therapies exist to promote sustained cartilage regeneration, as we have not been able to recapitulate the programming events necessary to instruct cells to form articular cartilage without these cells continuing to differentiate into bone. Our analysis of the early programming events occurring during cartilage formation led to the identification of LncRNA-HIT a long noncoding RNA that is essential for the differentiation of the embryonic limb mesenchyme into cartilage. A genome wide analysis of LncRNA-HIT’s distribution in the mesenchyme revealed strong association between LncRNA-HIT and numerous genes whose products facilitate cartilage formation. In the absence of LncRNA-HIT, the expression of these chondrogenic genes is severely reduced, impacting the differentiation of these cells into cartilage. Mechanistically, LncRNA-HIT regulates these pro-chondrogenic genes by recruiting p100 and CBP to these loci, facilitating H3K27ac and transcriptional activation. LncRNA-HIT also appears to be present in most vertebrate species, suggesting that the epigenetic program regulated by this lncRNA may represent a fundamental mechanism used by many species to promote cartilage formation.

A potential role of homeobox transcription factors in osteoarthritis

When a healthy joint progressively becomes osteoarthritic, the structures of the affected cartilage, bone and synovia undergo an initial phase of rearrangement. The exact molecular and cellular events occurring in this early osteoarthritic transition phase still remain elusive. Homeobox (Hox) genes encode for transcription factors that typically regulate limb morphogenesis and skeletal formation during development. More recently they were shown to be required for tissue remodelling and homeostasis in adults and to be modulated in a variety of pathologies. Here we present and discuss the hypothesis that dysregulation of specific Hox genes is associated with the onset and development of osteoarthritis (OA). Discovering mechanisms modulating Hox gene expression could not only provide important information in understanding OA pathology and its initiation, but also help to identify biomarkers reflecting the state of early OA. This knowledge would allow anticipating the time window for clinical treatment of the affected cartilage and assist in the development of innovative strategies to restore joint homeostasis, e.g., by cell or gene therapy.

Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development

Background

The 39 mammalian Hox genes show problematic patterns of functional overlap. In order to more fully define the developmental roles of Hox genes it is necessary to remove multiple combinations of paralogous and flanking genes. In addition, the downstream molecular pathways regulated by Hox genes during limb development remain incompletely delineated.

Results

In this report we examine limb development in mice with frameshift mutations in six Hox genes, Hoxa9,10,11 and Hoxd9,10,11. The mice were made with a novel recombineering method that allows the simultaneous targeting of frameshift mutations into multiple flanking genes. The Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant mice show a reduced ulna and radius that is more severe than seen in Hoxa11^−/−/Hoxd11^−/− mice, indicating a minor role for the flanking Hox9,10 genes in zeugopod development, as well as their primary function in stylopod development. The mutant mice also show severe reduction of Shh expression in the zone of polarizing activity, and decreased Fgf8 expression in the apical ectodermal ridge, thereby better defining the roles of these specific Hox genes in the regulation of critical signaling centers during limb development. Importantly, we also used laser capture microdissection coupled with RNA-Seq to characterize the gene expression programs in wild type and mutant limbs. Resting, proliferative and hypertrophic compartments of E15.5 forelimb zeugopods were examined. The results provide an RNA-Seq characterization of the progression of gene expression patterns during normal endochondral bone formation. In addition of the Hox mutants showed strongly altered expression of Pknox2, Zfp467, Gdf5, Bmpr1b, Dkk3, Igf1, Hand2, Shox2, Runx3, Bmp7 and Lef1, all of which have been previously shown to play important roles in bone formation.

Conclusions

The recombineering based frameshift mutation of the six flanking and paralogous Hoxa9,10,11 and Hoxd9,10,11 genes provides a resource for the analysis of their overlapping functions. Analysis of the Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant limbs confirms and extends the results of previous studies using mice with Hox mutations in single paralogous groups or with entire Hox cluster deletions. The RNA-Seq analysis of specific compartments of the normal and mutant limbs defines the multiple key perturbed pathways downstream of these Hox genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-015-0078-5) contains supplementary material, which is available to authorized users.

Modulation of chondrogenic differentiation of human mesenchymal stem cells in jellyfish collagen scaffolds by cell density and culture medium

Studies on tissue-engineering approaches for the regeneration of traumatized cartilage focus increasingly on multipotent human mesenchymal stem cells (hMSCs) as an alternative to autologous chondrocytes. The present study applied porous scaffolds made of collagen from the jellyfish Rhopilema esculentum for the in vitro chondrogenic differentiation of hMSCs. Culture conditions in those scaffolds differ from conditions in high-density pellet cultures, making a re-examination of these data necessary. We systematically investigated the influence of seeding density, basic culture media [Dulbecco's modified Eagle's medium (DMEM), α-minimum essential medium (α-MEM)] with varying glucose content and supplementation with fetal calf serum (FCS) or bovine serum albumin (BSA) on the chondrogenic differentiation of hMSCs. Gene expression analyses of selected markers for chondrogenic differentiation and hypertrophic development were conducted. Furthermore, the production of cartilage extracellular matrix (ECM) was analysed by quantification of sulphated glycosaminoglycan and collagen type II contents. The strongest upregulation of chondrogenic markers, along with the highest ECM deposition was observed in scaffolds seeded with 2.4 × 10⁶ cells/cm³ after cultivation in high-glucose DMEM and 0.125% BSA. Lower seeding densities compared to high-density pellet cultures were sufficient to induce in vitro chondrogenic differentiation of hMSCs in collagen scaffolds, which reduces the amount of cells required for the seeding of scaffolds and thus the monolayer expansion period. Furthermore, examination of the impact of FCS and α-MEM on chondrogenic MSC differentiation is an important prerequisite for the development of an osteochondral medium for simultaneous osteogenic and chondrogenic differentiation in biphasic scaffolds for osteochondral tissue regeneration. Copyright © 2015 John Wiley & Sons, Ltd.

Distinct Transcriptional Programs Underlie Sox9 Regulation of the Mammalian Chondrocyte

Sox9 encodes an essential transcriptional regulator of chondrocyte specification and differentiation. When Sox9 nuclear activity was compared with markers of chromatin organization and transcriptional activity in primary chondrocytes, we identified two distinct categories of target association. Class I sites cluster around the transcriptional start sites of highly expressed genes with no chondrocyte-specific signature. Here, Sox9 association reflects protein-protein association with basal transcriptional components. Class II sites highlight evolutionarily conserved active enhancers that direct chondrocyte-related gene activity through the direct binding of Sox9 dimer complexes to DNA. Sox9 binds through sites with sub-optimal binding affinity; the number and grouping of enhancers into super-enhancer clusters likely determines the levels of target gene expression. Interestingly, comparison of Sox9 action in distinct chondrocyte lineages points to similar regulatory strategies. In addition to providing insights into Sox family action, our comprehensive identification of the chondrocyte regulatory genome will facilitate the study of skeletal development and human disease.

Epigenetic changes and functional study of HOXA11 in human gastric cancer

AIM

To examine epigenetic changes and the function of HOXA11 in human gastric cancer (GC).

MATERIALS & METHODS

Seven GC cell lines, five cases of normal gastric mucosa and 112 cases primary GC samples were used in this study.

RESULTS

Expression of HOXA11 and lack of promoter region methylation were found in NCI-N87, MKN45, BGC823 and HGC27 cells. Loss of expression and complete methylation were found in AGS gastric cancer cells. Reduced expression and partial methylation were found in MGC803 and SGC7901 cells. Restoration of HOXA11 expression was induced by 5-aza-2'-deoxycytidine. HOXA11 was methylated in 81.25% (91/112) of primary GCs. The presence of methylation was associated with male gender, tumor size, tumor differentiation and lymph node metastasis (all p < 0.05). Restoration of HOXA11 expression reduced cell proliferation, invasion, migration and induced apoptosis and G2/M phase arrest. HOXA11 was found to inhibit Wnt signaling by upregulating NKD1 expression.

CONCLUSION

Epigenetic silencing of HOXA11 promotes GC proliferation, migration and invasion through activation of Wnt signaling.

H3K4me1 marks DNA regions hypomethylated during aging in human stem and differentiated cells

In differentiated cells, aging is associated with hypermethylation of DNA regions enriched in repressive histone post-translational modifications. However, the chromatin marks associated with changes in DNA methylation in adult stem cells during lifetime are still largely unknown. Here, DNA methylation profiling of mesenchymal stem cells (MSCs) obtained from individuals aged 2 to 92 yr identified 18,735 hypermethylated and 45,407 hypomethylated CpG sites associated with aging. As in differentiated cells, hypermethylated sequences were enriched in chromatin repressive marks. Most importantly, hypomethylated CpG sites were strongly enriched in the active chromatin mark H3K4me1 in stem and differentiated cells, suggesting this is a cell type–independent chromatin signature of DNA hypomethylation during aging. Analysis of scedasticity showed that interindividual variability of DNA methylation increased during aging in MSCs and differentiated cells, providing a new avenue for the identification of DNA methylation changes over time. DNA methylation profiling of genetically identical individuals showed that both the tendency of DNA methylation changes and scedasticity depended on nongenetic as well as genetic factors. Our results indicate that the dynamics of DNA methylation during aging depend on a complex mixture of factors that include the DNA sequence, cell type, and chromatin context involved and that, depending on the locus, the changes can be modulated by genetic and/or external factors.

Genetic interactions between Shox2 and Hox genes during the regional growth and development of the mouse limb

The growth and development of the vertebrate limb relies on homeobox genes of the Hox and Shox families, with their independent mutation often giving dose-dependent effects. Here we investigate whether Shox2 and Hox genes function together during mouse limb development by modulating their relative dosage and examining the limb for nonadditive effects on growth. Using double mRNA fluorescence in situ hybridization (FISH) in single embryos, we first show that Shox2 and Hox genes have associated spatial expression dynamics, with Shox2 expression restricted to the proximal limb along with Hoxd9 and Hoxa11 expression, juxtaposing the distal expression of Hoxa13 and Hoxd13. By generating mice with all possible dosage combinations of mutant Shox2 alleles and HoxA/D cluster deletions, we then show that their coordinated proximal limb expression is critical to generate normally proportioned limb segments. These epistatic interactions tune limb length, where Shox2 underexpression enhances, and Shox2 overexpression suppresses, Hox-mutant phenotypes. Disruption of either Shox2 or Hox genes leads to a similar reduction in Runx2 expression in the developing humerus, suggesting their concerted action drives cartilage maturation during normal development. While we furthermore provide evidence that Hox gene function influences Shox2 expression, this regulation is limited in extent and is unlikely on its own to be a major explanation for their genetic interaction. Given the similar effect of human SHOX mutations on regional limb growth, Shox and Hox genes may generally function as genetic interaction partners during the growth and development of the proximal vertebrate limb.

Long bone development requires a threshold of Hox function

The Hoxd(Del(11-13)) mutant is one of the animal models for human synpolydactyly, characterized by short and syndactylous digits. Here we have characterized in detail the cartilage and bone defects in these mutants. We report two distinct phenotypes: (i) a delay and change in pattern of chondrocyte maturation of metacarpals/metatarsals and (ii) formation of a poor and not centrally positioned primary ossification center in the proximal-intermediate phalanx. In the metacarpals of Hoxd(Del(11-13)) mutants, ossification occurs postnataly, in the absence of significant Ihh expression and without the establishment of growth plates, following patterns similar to those of short bones. The strong downregulation in Ihh expression is associated with a corresponding increase of the repressor form of Gli3. To evaluate the contribution of this alteration to the phenotype, we generated double Hoxd(Del(11-13));Gli3 homozygous mutants. Intriguingly, these double mutants showed a complete rescue of the phenotype in metatarsals but only partial phenotypic rescue in metacarpals. Our results support Hox genes being required in a dose-dependent manner for long bone cartilage maturation and suggest that and excess of Gli3R mediates a significant part of the Hoxd(Del(11-13)) chondrogenic phenotype.

Ezh2 is required for neural crest-derived cartilage and bone formation

The emergence of craniofacial skeletal elements, and of the jaw in particular, was a crucial step in the evolution of higher vertebrates. Most facial bones and cartilage are generated during embryonic development by cranial neural crest cells, while an osteochondrogenic fate is suppressed in more posterior neural crest cells. Key players in this process are Hox genes, which suppress osteochondrogenesis in posterior neural crest derivatives. How this specific pattern of osteochondrogenic competence is achieved remains to be elucidated. Here we demonstrate that Hox gene expression and osteochondrogenesis are controlled by epigenetic mechanisms. Ezh2, which is a component of polycomb repressive complex 2 (PRC2), catalyzes trimethylation of lysine 27 in histone 3 (H3K27me3), thereby functioning as transcriptional repressor of target genes. Conditional inactivation of Ezh2 does not interfere with localization of neural crest cells to their target structures, neural development, cell cycle progression or cell survival. However, loss of Ezh2 results in massive derepression of Hox genes in neural crest cells that are usually devoid of Hox gene expression. Accordingly, craniofacial bone and cartilage formation is fully prevented in Ezh2 conditional knockout mice. Our data indicate that craniofacial skeleton formation in higher vertebrates is crucially dependent on epigenetic regulation that keeps in check inhibitors of an osteochondrogenic differentiation program.

Jellyfish collagen scaffolds for cartilage tissue engineering

Porous scaffolds were engineered from refibrillized collagen of the jellyfish Rhopilema esculentum for potential application in cartilage regeneration. The influence of collagen concentration, salinity and temperature on fibril formation was evaluated by turbidity measurements and quantification of fibrillized collagen. The formation of collagen fibrils with a typical banding pattern was confirmed by atomic force microscopy and transmission electron microscopy analysis. Porous scaffolds from jellyfish collagen, refibrillized under optimized conditions, were fabricated by freeze-drying and subsequent chemical cross-linking. Scaffolds possessed an open porosity of 98.2%. The samples were stable under cyclic compression and displayed an elastic behavior. Cytotoxicity tests with human mesenchymal stem cells (hMSCs) did not reveal any cytotoxic effects of the material. Chondrogenic markers SOX9, collagen II and aggrecan were upregulated in direct cultures of hMSCs upon chondrogenic stimulation. The formation of typical extracellular matrix components was further confirmed by quantification of sulfated glycosaminoglycans.

Tbx4 interacts with the short stature homeobox gene Shox2 in limb development

BACKGROUND

The short stature homeodomain transcription factors SHOX and SHOX2 play key roles in limb formation. To gain more insight into genes regulated by Shox2 during limb development, we analyzed expression profiles of WT and Shox2-/- mouse embryonic limbs and identified the T-Box transcription factor Tbx4 as a potential downstream target. Tbx4 is known to exert essential functions in skeletal and muscular hindlimb development. In humans, haploinsufficiency of TBX4 causes small patella syndrome, a skeletal dysplasia characterized by anomalies of the knee, pelvis, and foot.

RESULTS

Here, we demonstrate an inhibitory regulatory effect of Shox2 on Tbx4 specifically in the forelimbs. We also show that Tbx4 activates Shox2 expression in fore- and hindlimbs, suggesting Shox2 as a feedback modulator of Tbx4. Using EMSA studies, we find that Tbx4/TBX4 is able to bind to distinct T-box binding sites within the mouse and human Shox2/SHOX2 promoter.

CONCLUSIONS

Our data identifies Tbx4 as a novel transcriptional activator of Shox2 during murine fore- and hindlimb development. Tbx4 is also regulated by Shox2 specifically in the forelimb bud possibly via a feedback mechanism. These data extend our understanding of the role and regulation of Tbx4 and Shox2 in limb development and limb associated diseases.

Understanding developmental mechanisms in the context of osteoarthritis

Osteoarthritis (OA) is a joint disease that is highly related to aging. However, as OA development is the consequence of interplay between external stimuli, such as mechanical loading and the structure and physiology of the joint, it can be anticipated that variation in developmental processes early in life will affect OA development later in life. Genes involved in patterning processes, such as the Hox genes, but also genes that encode transcription factors, growth factors and cytokines and their respective receptors and those that encode molecules involved in formation of the extracellular matrix, will influence embryonic skeletal development and OA incidence and severity in the adult. The function of genes involved in patterning processes can be partly be understood by close analysis of inborn diseases that result in musculoskeletal syndromes, but a deeper understanding will be the result of specific gene knockouts or overexpression in transgenic mouse models.

Decoupling the function of Hox and Shh in developing limb reveals multiple inputs of Hox genes on limb growth

Limb development relies on an exquisite coordination between growth and patterning, but the underlying mechanisms remain elusive. Anterior-posterior and proximal-distal specification initiates in early limb bud concomitantly with the proliferative expansion of limb cells. Previous studies have shown that limb bud growth initially relies on fibroblast growth factors (FGFs) produced in the apical ectodermal ridge (AER-FGFs), the maintenance of which relies on a positive-feedback loop involving sonic hedgehog (Shh) and the BMP antagonist gremlin 1 (Grem1). The positive cross-regulation between Shh and the HoxA and HoxD clustered genes identified an indirect effect of Hox genes on the maintenance of AER-FGFs but the respective function of Shh and Hox genes in this process remains unknown. Here, by uncoupling Hox and Shh function, we show that HoxA and HoxD genes are required for proper AER-FGFs expression, independently of their function in controlling Shh expression. In addition, we provide evidence that the Hox-dependent control of AER-FGF expression is achieved through the regulation of key mesenchymal signals, namely Grem1 and Fgf10, ensuring proper epithelial-mesenchymal interactions. Notably, HoxA and HoxD genes contribute to both the initial activation of Grem1 and the subsequent anterior expansion of its expression domain. We propose that the intricate interactions between Hox genes and the FGF and Shh signaling pathways act as a molecular network that ensures proper limb bud growth and patterning, probably contributing to the coordination of these two processes.