The functions of Runx family transcription factors and Cbfb in skeletal development

Network model integrated with multi-omic data predicts MBNL1 signals that drive myofibroblast activation

RNA-binding protein muscleblind-like1 (MBNL1) was recently identified as a central regulator of cardiac wound healing and myofibroblast activation. To identify putative MBNL1 targets, we integrated multiple genome-wide screens with a fibroblast network model. We expanded the model to include putative MBNL1-target interactions and recapitulated published experimental results to validate new signaling modules. We prioritized 14 MBNL1 targets and developed novel fibroblast signaling modules for p38 MAPK, Hippo, Runx1, and Sox9 pathways. We experimentally validated MBNL1 regulation of p38 expression in mouse cardiac fibroblasts. Using the expanded fibroblast model, we predicted a hierarchy of MBNL1 regulated pathways with strong influence on αSMA expression. This study lays a foundation to explore the network mechanisms of MBNL1 signaling central to fibrosis.

Lineage-specific differences and regulatory networks governing human chondrocyte development

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.

Runx1 is a key regulator of articular cartilage homeostasis by orchestrating YAP, TGFβ, and Wnt signaling in articular cartilage formation and osteoarthritis

Runt-related transcription factor 1 (Runx1) plays a key role in cartilage formation, but its function in articular cartilage formation is unclear. We generated non-inducible and inducible Runx1-deficient mice (Runx1f/fCol2α1-Cre and Runx1f/fCol2α1-CreER mice) and found that chondrocyte-specific Runx1-deficient mice developed a spontaneous osteoarthritis (OA)-like phenotype and showed exacerbated articular cartilage destruction under OA, characterized by articular cartilage degradation and cartilage ossification, with decreased Col2α1 expression and increased Mmp13 and Adamts5 expression. RNA-sequencing analysis of hip articular cartilage from the Runx1f/fCol2α1-Cre mice compared to that from wild-type mice and subsequent validation analyses demonstrated that Runx1 is a central regulator in multiple signaling pathways, converging signals of the Hippo/Yap, TGFβ/Smad, and Wnt/β-catenin pathways into a complex network to regulate the expression of downstream genes, thereby controlling a series of osteoarthritic pathological processes. RNA-sequencing analysis of mutant knee joints showed that Runx1's role in signaling pathways in articular cartilage is different from that in whole knee joints, indicating that Runx1 regulation is tissue-specific. Histopathologic analysis confirmed that Runx1 deficiency decreased the levels of YAP and p-Smad2/3 and increased the levels of active β-catenin. Overexpression of Runx1 dramatically increased YAP expression in chondrocytes. Adeno-associated virus-mediated Runx1 overexpression in the knee joints of osteoarthritic mice showed the protective effect of Runx1 on articular cartilage damaged in OA. Our results notably showed that Runx1 is a central regulator of articular cartilage homeostasis by orchestrating the YAP, TGFβ, and Wnt signaling pathways in the formation of articular cartilage and OA, and targeting Runx1 and its downstream genes may facilitate the design of novel therapeutic approaches for OA.

Interactions and Feedbacks in E-Cadherin Transcriptional Regulation

Epithelial tissues rely on the adhesion between participating cells to retain their integrity. The transmembrane protein E-cadherin is the major protein that mediates homophilic adhesion between neighbouring cells and is, therefore, one of the critical components for epithelial integrity. E-cadherin downregulation has been described extensively as a prerequisite for epithelial-to-mesenchymal transition and is a hallmark in many types of cancer. Due to this clinical importance, research has been mostly focused on understanding the mechanisms leading to transcriptional repression of this adhesion molecule. However, in recent years it has become apparent that re-expression of E-cadherin is a major step in the progression of many cancers during metastasis. Here, we review the currently known molecular mechanisms of E-cadherin transcriptional activation and inhibition and highlight complex interactions between individual mechanisms. We then propose an additional mechanism, whereby the competition between adhesion complexes and heterochromatin protein-1 for binding to STAT92E fine-tunes the levels of E-cadherin expression in Drosophila but also regulates other genes promoting epithelial robustness. We base our hypothesis on both existing literature and our experimental evidence and suggest that such feedback between the cell surface and the nucleus presents a powerful paradigm for epithelial resilience.

Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

Significance of MEF2C and RUNX3 Regulation for Endochondral Differentiation of Human Mesenchymal Progenitor Cells

Guiding progenitor cell development between chondral versus endochondral pathways is still an unachieved task of cartilage neogenesis, and human mesenchymal progenitor cell (MPC) chondrogenesis is considered as a valuable model to better understand hypertrophic development of chondrocytes. Transcription factors Runx2, Runx3, and Mef2c play prominent roles for chondrocyte hypertrophy during mouse development, but little is known on the importance of these key fate-determining factors for endochondral development of human MPCs. The aim of this study was to unravel the regulation of RUNX2, RUNX3, and MEF2C during MPC chondrogenesis, the pathways driving their expression, and the downstream hypertrophic targets affected by their regulation. RUNX2, RUNX3, and MEF2C gene expression was differentially regulated during chondrogenesis of MPCs, but remained low and unregulated when non-hypertrophic articular chondrocytes were differentiated under the same conditions. RUNX3 and MEF2C mRNA and protein levels rose in parallel to hypertrophic marker upregulation, but surprisingly, RUNX2 gene expression changed only by trend and RUNX2 protein remained undetectable. While RUNX3 expression was driven by TGF-β and BMP signaling, MEF2C responded to WNT-, BMP-, and Hedgehog-pathway inhibition. MEF2C but not RUNX3 levels correlated significantly with COL10A1, IHH, and IBSP gene expression when hypertrophy was attenuated. IBSP was a downstream target of RUNX3 and MEF2C but not RUNX2 in SAOS-2 cells, underlining the capacity of RUNX3 and MEF2C to stimulate osteogenic marker expression in human cells. Conclusively, RUNX3 and MEF2C appeared more important than RUNX2 for human endochondral MPC chondrogenesis. Pathways altering the speed of chondrogenesis (FGF, TGF-β, BMP) affected RUNX2 or RUNX3, while pathways changing hypertrophy (WNT, PTHrP/HH) regulated mainly MEF2C. Taken together, reduction of MEF2C levels is a new goal to shift human cartilage neogenesis toward the chondral pathway.