Systematic Reconstruction of Molecular Cascades Regulating GP Development Using Single-Cell RNA-Seq

The growth plate (GP) comprising sequentially differentiated cell layers is a critical structure for bone elongation and regeneration. Although several key regulators in GP development have been identified using genetic perturbation, systematic understanding is still limited. Here, we used single-cell RNA-sequencing (RNA-seq) to determine the gene expression profiles of 217 single cells from GPs and developed a bioinformatics pipeline named Sinova to de novo reconstruct physiological GP development in both temporal and spatial high resolution. Our unsupervised model not only confirmed prior knowledge, but also enabled the systematic discovery of genes, potential signal pathways, and surface markers CD9/CD200 to precisely depict development. Sinova further identified the effective combination of transcriptional factors (TFs) that regulates GP maturation, and the result was validated using an in vitro EGFP-Col10a screening system. Our case systematically reconstructed molecular cascades in GP development through single-cell profiling, and the bioinformatics pipeline is applicable to other developmental processes.