Human Evolutionary Biology, Harvard University, Cambridge, United States
Steven Pregizer
Department of Orthopedic Research, Boston Children’s Hospital, Boston, United States; Department of Orthopedic Surgery, Harvard Medical School, Boston, United States
Divya Venkatasubramanian
Department of Orthopedic Research, Boston Children’s Hospital, Boston, United States; Department of Orthopedic Surgery, Harvard Medical School, Boston, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
Rosanne M Raftery
Department of Orthopedic Research, Boston Children’s Hospital, Boston, United States; Department of Orthopedic Surgery, Harvard Medical School, Boston, United States
Pushpanathan Muthuirulan
Human Evolutionary Biology, Harvard University, Cambridge, United States
Zun Liu
Human Evolutionary Biology, Harvard University, Cambridge, United States
Terence D Capellini
Human Evolutionary Biology, Harvard University, Cambridge, United States; Broad Institute of MIT and Harvard, Cambridge, United States
Department of Orthopedic Research, Boston Children’s Hospital, Boston, United States; Department of Orthopedic Surgery, Harvard Medical School, Boston, United States; Harvard Stem Cell Institute, Cambridge, United States
To address large gaps in our understanding of the molecular regulation of articular and growth plate cartilage development in humans, we used our directed differentiation approach to generate these distinct cartilage tissues from human embryonic stem cells. The resulting transcriptomic profiles of hESC-derived articular and growth plate chondrocytes were similar to fetal epiphyseal and growth plate chondrocytes, with respect to genes both known and previously unknown to cartilage biology. With the goal to characterize the regulatory landscapes accompanying these respective transcriptomes, we mapped chromatin accessibility in hESC-derived chondrocyte lineages, and mouse embryonic chondrocytes, using ATAC-sequencing. Integration of the expression dataset with the differentially accessible genomic regions revealed lineage-specific gene regulatory networks. We validated functional interactions of two transcription factors (TFs) (RUNX2 in growth plate chondrocytes and RELA in articular chondrocytes) with their predicted genomic targets. The maps we provide thus represent a framework for probing regulatory interactions governing chondrocyte differentiation. This work constitutes a substantial step towards comprehensive and comparative molecular characterizations of distinct chondrogenic lineages and sheds new light on human cartilage development and biology.