Genes uniquely expressed in human growth plate chondrocytes uncover a distinct regulatory network.

BACKGROUND: Chondrogenesis is the earliest stage of skeletal development and is a highly dynamic process, integrating the activities and functions of transcription factors, cell signaling molecules and extracellular matrix proteins. The molecular mechanisms underlying chondrogenesis have been extensively studied and multiple key regulators of this process have been identified. However, a genome-wide overview of the gene regulatory network in chondrogenesis has not been achieved. RESULTS: In this study, employing RNA sequencing, we identified 332 protein coding genes and 34 long non-coding RNA (lncRNA) genes that are highly selectively expressed in human fetal growth plate chondrocytes. Among the protein coding genes, 32 genes were associated with 62 distinct human skeletal disorders and 153 genes were associated with skeletal defects in knockout mice, confirming their essential roles in skeletal formation. These gene products formed a comprehensive physical interaction network and participated in multiple cellular processes regulating skeletal development. The data also revealed 34 transcription factors and 11,334 distal enhancers that were uniquely active in chondrocytes, functioning as transcriptional regulators for the cartilage-selective genes. CONCLUSIONS: Our findings revealed a complex gene regulatory network controlling skeletal development whereby transcription factors, enhancers and lncRNAs participate in chondrogenesis by transcriptional regulation of key genes. Additionally, the cartilage-selective genes represent candidate genes for unsolved human skeletal disorders.

MED resulting from recessively inherited mutations in the gene encoding calcium-activated nucleotidase CANT1.

Multiple Epiphyseal Dysplasia (MED) is a relatively mild skeletal dysplasia characterized by mild short stature, joint pain, and early-onset osteoarthropathy. Dominantly inherited mutations in COMP, MATN3, COL9A1, COL9A2, and COL9A3, and recessively inherited mutations in SLC26A2, account for the molecular basis of disease in about 80-85% of the cases. In two families with recurrent MED of an unknown molecular basis, we used exome sequencing and candidate gene analysis to identify homozygosity for recessively inherited missense mutations in CANT1, which encodes calcium-activated nucleotidase 1. The MED phenotype is thus allelic to the more severe Desbuquois dysplasia phenotype and the results identify CANT1 as a second locus for recessively inherited MED.

Heredity and cardiometabolic risk: naturally occurring polymorphisms in the human neuropeptide Y(2) receptor promoter disrupt multiple transcriptional response motifs.

OBJECTIVES: The neuropeptide Y(2) G-protein-coupled receptor (NPY2R) relays signals from PYY or neuropeptide Y toward satiety and control of body mass. Targeted ablation of the NPY2R locus in mice yields obesity, and studies of NPY2R promoter genetic variation in more than 10,000 human participants indicate its involvement in control of obesity and BMI. Here we searched for genetic variation across the human NPY2R locus and probed its functional effects, especially in the proximal promoter. METHODS AND RESULTS: Twin pair studies indicated substantial heritability for multiple cardiometabolic traits, including BMI, SBP, DBP, and PYY, an endogenous agonist at NPY2R. Systematic polymorphism discovery by resequencing across NPY2R uncovered 21 genetic variants, 10 of which were common [minor allele frequency (MAF) >5%], creating one to two linkage disequilibrium blocks in multiple biogeographic ancestries. In vivo, NPY2R haplotypes were associated with both BMI (P = 3.75E-04) and PYY (P = 4.01E-06). Computational approaches revealed that proximal promoter variants G-1606A, C-599T, and A-224G disrupt predicted IRF1 (A>G), FOXI1 (T>C), and SNAI1 (A>G) response elements. In neuroendocrine cells transfected with NPY2R promoter/luciferase reporter plasmids, all three variants and their resulting haplotypes influenced transcription (G-1606A, P < 2.97E-06; C-599T, P < 1.17E-06; A-224G, P < 2.04E-06), and transcription was differentially augmented or impaired by coexpression of either the cognate full-length transcription factors or their specific siRNAs at each site. Endogenous expression of transcripts for NPY2R, IRF1, and SNAI1 was documented in neuroendocrine cells, and the NPY2R mRNA was differentially expressed in two neuroendocrine tissues (adrenal gland, brainstem) of a rodent model of hypertension and the metabolic syndrome, the spontaneously hypertensive rat. CONCLUSION: We conclude that common genetic variation in the proximal NPY2R promoter influences transcription factor binding so as to alter gene expression in neuroendocrine cells, and consequently cardiometabolic traits in humans. These results unveil a novel control point, whereby cis-acting genetic variation contributes to control of complex cardiometabolic traits, and point to new transcriptional strategies for intervention into neuropeptide actions and their cardiometabolic consequences.