Endochondral ossification initiates the growth of the majority of the mammalian skeleton and is tightly controlled through gene regulatory networks. The forkhead box transcription factors Foxc1 and Foxc2 regulate aspects of osteoblast function in the formation of the skeleton, but their roles in chondrocytes to control endochondral ossification are less clear. Here, we demonstrate that Foxc1 expression is directly regulated by the activity of SRY (sex-determining region Y)-box 9, one of the earliest transcription factors to specify the chondrocyte lineage. Moreover, we demonstrate that elevated expression of Foxc1 promotes chondrocyte differentiation in mouse embryonic stem cells and loss of Foxc1 function inhibits chondrogenesis in vitro . Using chondrocyte-targeted deletion of Foxc1 and Foxc2 in mice, we reveal a role for these factors in chondrocyte differentiation in vivo . Loss of both Foxc1 and Foxc2 caused a general skeletal dysplasia predominantly affecting the vertebral column. The long bones of the limbs were smaller, mineralization was reduced, and organization of the growth plate was disrupted; in particular, the stacked columnar organization of the proliferative chondrocyte layer was reduced in size and cell proliferation was decreased. Differential gene expression analysis indicated disrupted expression patterns of chondrogenesis and ossification genes throughout the entire process of endochondral ossification in chondrocyte-specific Foxc1/Foxc2 KO embryos. Our results suggest that Foxc1 and Foxc2 are required for normal chondrocyte differentiation and function, as loss of both genes results in disorganization of the growth plate, reduced chondrocyte proliferation, and delays in chondrocyte hypertrophy that prevents ossification of the skeleton.

Despite the tremendous increase in omics data generated by modern sequencing technologies, their analysis can be tricky and often requires substantial expertise in bioinformatics. To address this concern, we have developed a user-friendly pipeline to analyze (cancer) genomic data that takes in raw sequencing data (FASTQ format) as input and outputs insightful statistics. Our iCOMIC toolkit pipeline featuring many independent workflows is embedded in the popular Snakemake workflow management system. It can analyze whole-genome and transcriptome data and is characterized by a user-friendly GUI that offers several advantages, including minimal execution steps and eliminating the need for complex command-line arguments. Notably, we have integrated algorithms developed in-house to predict pathogenicity among cancer-causing mutations and differentiate between tumor suppressor genes and oncogenes from somatic mutation data. We benchmarked our tool against Genome In A Bottle benchmark dataset (NA12878) and got the highest F1 score of 0.971 and 0.988 for indels and SNPs, respectively, using the BWA MEM—GATK HC DNA-Seq pipeline. Similarly, we achieved a correlation coefficient of r = 0.85 using the HISAT2-StringTie-ballgown and STAR-StringTie-ballgown RNA-Seq pipelines on the human monocyte dataset (SRP082682). Overall, our tool enables easy analyses of omics datasets, significantly ameliorating complex data analysis pipelines.

Endochondral ossification forms and grows the majority of the mammalian skeleton and is tightly controlled through gene regulatory networks. The forkhead box transcription factors Foxc1 and Foxc2 have been demonstrated to regulate aspects of osteoblast function in the formation of the skeleton but their roles in chondrocytes to control endochondral ossification are less clear. We demonstrate that Foxc1 expression is directly regulated by SOX9 activity, one of the earliest transcription factors to specify the chondrocyte lineages. Moreover we demonstrate that elevelated expression of Foxc1 promotes chondrocyte differentiation in mouse embryonic stem cells and loss of Foxc1 function inhibits chondrogenesis in vitro. Using chondrocyte-targeted deletion of Foxc1 and Foxc2 in mice, we reveal a role for these factors in chondrocyte differentiation in vivo. Loss of both Foxc1 and Foxc2 caused a general skeletal dysplasia predominantly affecting the vertebral column. The long bones of the limb were smaller and mineralization was reduced and organization of the growth plate was disrupted. In particular, the stacked columnar organization of the proliferative chondrocyte layer was reduced in size and cell proliferation in growth plate chondrocytes was reduced. Differential gene expression analysis indicated disrupted expression patterns in chondrogenesis and ossification genes throughout the entire process of endochondral ossification in Col2-cre;Foxc1Δ/Δ;Foxc2Δ/Δ embryos. Our results suggest that Foxc1 and Foxc2 are required for correct chondrocyte differentiation and function. Loss of both genes results in disorganization of the growth plate, reduced chondrocyte proliferation and delays in chondrocyte hypertrophy that prevents correct ossification of the endochondral skeleton.

No abstract