Epigenetic regulators controlling osteogenic lineage commitment and bone formation.

Bone formation and homeostasis are controlled by environmental factors and endocrine regulatory cues that initiate intracellular signaling pathways capable of modulating gene expression in the nucleus. Bone-related gene expression is controlled by nucleosome-based chromatin architecture that limits the accessibility of lineage-specific gene regulatory DNA sequences and sequence-specific transcription factors. From a developmental perspective, bone-specific gene expression must be suppressed during the early stages of embryogenesis to prevent the premature mineralization of skeletal elements during fetal growth in utero. Hence, bone formation is initially inhibited by gene suppressive epigenetic regulators, while other epigenetic regulators actively support osteoblast differentiation. Prominent epigenetic regulators that stimulate or attenuate osteogenesis include lysine methyl transferases (e.g., EZH2, SMYD2, SUV420H2), lysine deacetylases (e.g., HDAC1, HDAC3, HDAC4, HDAC7, SIRT1, SIRT3), arginine methyl transferases (e.g., PRMT1, PRMT4/CARM1, PRMT5), dioxygenases (e.g., TET2), bromodomain proteins (e.g., BRD2, BRD4) and chromodomain proteins (e.g., CBX1, CBX2, CBX5). This narrative review provides a broad overview of the covalent modifications of DNA and histone proteins that involve hundreds of enzymes that add, read, or delete these epigenetic modifications that are relevant for self-renewal and differentiation of mesenchymal stem cells, skeletal stem cells and osteoblasts during osteogenesis.

Protein arginine methyltransferases PRMT1, PRMT4/CARM1 and PRMT5 have distinct functions in control of osteoblast differentiation.

Osteogenic differentiation of mesenchymal cells is controlled by epigenetic enzymes that regulate post-translational modifications of histones. Compared to acetyl or methyltransferases, the physiological functions of protein arginine methyltransferases (PRMTs) in osteoblast differentiation remain minimally understood. Therefore, we surveyed the expression and function of all nine mammalian PRMT members during osteoblast differentiation. RNA-seq gene expression profiling shows that Prmt1, Prmt4/Carm1 and Prmt5 represent the most prominently expressed PRMT subtypes in mouse calvarial bone and MC3T3 osteoblasts as well as human musculoskeletal tissues and mesenchymal stromal cells (MSCs). Based on effects of siRNA depletion, it appears that PRMT members have different functional effects: (i) loss of Prmt1 stimulates and (ii) loss of Prmt5 decreases calcium deposition of mouse MC3T3 osteoblasts, while (iii) loss of Carm1 is inconsequential for calcium deposition. Decreased Prmt5 suppresses expression of multiple genes involved in mineralization (e.g., Alpl, Ibsp, Phospho1) consistent with a positive role in osteogenesis. Depletion of Prmt1, Carm1 and Prmt5 has intricate but modest time-dependent effects on the expression of a panel of osteoblast differentiation and proliferation markers but does not change mRNA levels for select epigenetic regulators (e.g., Ezh1, Ezh2, Brd2 and Brd4). Treatment with the Class I PRMT inhibitor GSK715 enhances extracellular matrix mineralization of MC3T3 cells, while blocking formation of H3R17me2a but not H4R3me2a marks. In sum, Prmt1, Carm1 and Prmt5 have distinct biological roles during osteoblast differentiation, and different types histone H3 and H4 arginine methylation may contribute to the chromatin landscape during osteoblast differentiation.

The lysine methyltransferases SET and MYND domain containing 2 (Smyd2) and Enhancer of Zeste 2 (Ezh2) co-regulate osteoblast proliferation and mineralization.

Bone formation is controlled by histone modifying enzymes that regulate post-translational modifications on nucleosomal histone proteins and control accessibility of transcription factors to gene promoters required for osteogenesis. Enhancer of Zeste homolog 2 (EZH2/Ezh2), a histone H3 lysine 27 (H3K27) methyl transferase, is a suppressor of osteoblast differentiation. Ezh2 is regulated by SET and MYND domain-containing protein 2 (SMYD2/Smyd2), a lysine methyltransferase that modifies both histone and non-histone proteins. Here, we examined whether Smyd2 modulates Ezh2 suppression of osteoblast differentiation. Musculoskeletal RNA-seq data show that SMYD2/Smyd2 is the most highly expressed SMYD/Smyd member in human bone tissues and mouse osteoblasts. Smyd2 loss of function analysis in mouse MC3T3 osteoblasts using siRNA depletion enhances proliferation and calcium deposition. Loss of Smyd2 protein does not affect alkaline phosphatase activity nor does it result in a unified expression response for standard osteoblast-related mRNA markers (e.g., Bglap, Ibsp, Spp1, Sp7), indicating that Smyd2 does not directly control osteoblast differentiation. Smyd2 protein depletion enhances levels of the osteo-suppressive Ezh2 protein and H3K27 trimethylation (H3K27me3), as expected from increased cell proliferation, while elevating the osteo-inductive Runx2 protein. Combined siRNA depletion of both Smyd2 and Ezh2 protein is more effective in promoting calcium deposition when compared to loss of either protein. Collectively, our results indicate that Smyd2 inhibits proliferation and indirectly the subsequent mineral deposition by osteoblasts. Mechanistically, Smyd2 represents a functional epigenetic regulator that operates in parallel to the suppressive effects of Ezh2 and H3K27 trimethylation on osteoblast differentiation.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life.

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Profiling of human epigenetic regulators using a semi-automated real-time qPCR platform validated by next generation sequencing.

Epigenetic mechanisms control phenotypic commitment of mesenchymal stromal/stem cells (MSCs) into osteogenic, chondrogenic or adipogenic lineages. To investigate enzymes and chromatin binding proteins controlling the epigenome, we developed a hybrid expression screening strategy that combines semi-automated real-time qPCR (RT-qPCR), next generation RNA sequencing (RNA-seq), and a novel data management application (FileMerge). This strategy was used to interrogate expression of a large cohort (n>300) of human epigenetic regulators (EpiRegs) that generate, interpret and/or edit the histone code. We find that EpiRegs with similar enzymatic functions are variably expressed and specific isoforms dominate over others in human MSCs. This principle is exemplified by analysis of key histone acetyl transferases (HATs) and deacetylases (HDACs), H3 lysine methyltransferases (e.g., EHMTs) and demethylases (KDMs), as well as bromodomain (BRDs) and chromobox (CBX) proteins. Our results show gender-specific expression of H3 lysine 9 [H3K9] demethylases (e.g., KDM5D and UTY) as expected and upregulation of distinct EpiRegs (n>30) during osteogenic differentiation of MSCs (e.g., HDAC5 and HDAC7). The functional significance of HDACs in osteogenic lineage commitment of MSCs was functionally validated using panobinostat (LBH-589). This pan-deacetylase inhibitor suppresses osteoblastic differentiation as evidenced by reductions in bone-specific mRNA markers (e.g., ALPL), alkaline phosphatase activity and calcium deposition (i.e., Alizarin Red staining). Thus, our RT-qPCR platform identifies candidate EpiRegs by expression screening, predicts biological outcomes of their corresponding inhibitors, and enables manipulation of the human epigenome using molecular or pharmacological approaches to control stem cell differentiation.

Enhancer of Zeste Homolog 2 Inhibition Stimulates Bone Formation and Mitigates Bone Loss Caused by Ovariectomy in Skeletally Mature Mice.

Perturbations in skeletal development and bone degeneration may result in reduced bone mass and quality, leading to greater fracture risk. Bone loss is mitigated by bone protective therapies, but there is a clinical need for new bone-anabolic agents. Previous work has demonstrated that Ezh2 (enhancer of zeste homolog 2), a histone 3 lysine 27 (H3K27) methyltransferase, suppressed differentiation of osteogenic progenitors. Here, we investigated whether inhibition of Ezh2 can be leveraged for bone stimulatory applications. Pharmacologic inhibition and siRNA knockdown of Ezh2 enhanced osteogenic commitment of MC3T3 preosteoblasts. Next generation RNA sequencing of mRNAs and real time quantitative PCR profiling established that Ezh2 inactivation promotes expression of bone-related gene regulators and extracellular matrix proteins. Mechanistically, enhanced gene expression was linked to decreased H3K27 trimethylation (H3K27me3) near transcriptional start sites in genome-wide sequencing of chromatin immunoprecipitations assays. Administration of an Ezh2 inhibitor modestly increases bone density parameters of adult mice. Furthermore, Ezh2 inhibition also alleviated bone loss in an estrogen-deficient mammalian model for osteoporosis. Ezh2 inhibition enhanced expression of Wnt10b and Pth1r and increased the BMP-dependent phosphorylation of Smad1/5. Thus, these data suggest that inhibition of Ezh2 promotes paracrine signaling in osteoblasts and has bone-anabolic and osteoprotective potential in adults.

Genome-Wide Studies Reveal that H3K4me3 Modification in Bivalent Genes Is Dynamically Regulated during the Pluripotent Cell Cycle and Stabilized upon Differentiation.

Stem cell phenotypes are reflected by posttranslational histone modifications, and this chromatin-related memory must be mitotically inherited to maintain cell identity through proliferative expansion. In human embryonic stem cells (hESCs), bivalent genes with both activating (H3K4me3) and repressive (H3K27me3) histone modifications are essential to sustain pluripotency. Yet, the molecular mechanisms by which this epigenetic landscape is transferred to progeny cells remain to be established. By mapping genomic enrichment of H3K4me3/H3K27me3 in pure populations of hESCs in G2, mitotic, and G1 phases of the cell cycle, we found striking variations in the levels of H3K4me3 through the G2-M-G1 transition. Analysis of a representative set of bivalent genes revealed that chromatin modifiers involved in H3K4 methylation/demethylation are recruited to bivalent gene promoters in a cell cycle-dependent fashion. Interestingly, bivalent genes enriched with H3K4me3 exclusively during mitosis undergo the strongest upregulation after induction of differentiation. Furthermore, the histone modification signature of genes that remain bivalent in differentiated cells resolves into a cell cycle-independent pattern after lineage commitment. These results establish a new dimension of chromatin regulation important in the maintenance of pluripotency.

Epigenetic Control of Skeletal Development by the Histone Methyltransferase Ezh2.

Epigenetic control of gene expression is critical for normal fetal development. However, chromatin-related mechanisms that activate bone-specific programs during osteogenesis have remained underexplored. Therefore, we investigated the expression profiles of a large cohort of epigenetic regulators (>300) during osteogenic differentiation of human mesenchymal cells derived from the stromal vascular fraction of adipose tissue (AMSCs). Molecular analyses establish that the polycomb group protein EZH2 (enhancer of zeste homolog 2) is down-regulated during osteoblastic differentiation of AMSCs. Chemical inhibitor and siRNA knockdown studies show that EZH2, a histone methyltransferase that catalyzes trimethylation of histone 3 lysine 27 (H3K27me3), suppresses osteogenic differentiation. Blocking EZH2 activity promotes osteoblast differentiation and suppresses adipogenic differentiation of AMSCs. High throughput RNA sequence (mRNASeq) analysis reveals that EZH2 inhibition stimulates cell cycle inhibitory proteins and enhances the production of extracellular matrix proteins. Conditional genetic loss of Ezh2 in uncommitted mesenchymal cells (Prrx1-Cre) results in multiple defects in skeletal patterning and bone formation, including shortened forelimbs, craniosynostosis, and clinodactyly. Histological analysis and mRNASeq profiling suggest that these effects are attributable to growth plate abnormalities and premature cranial suture closure because of precocious maturation of osteoblasts. We conclude that the epigenetic activity of EZH2 is required for skeletal patterning and development, but EZH2 expression declines during terminal osteoblast differentiation and matrix production.

Epigenetic pathways regulating bone homeostasis: potential targeting for intervention of skeletal disorders.

Epigenetic regulation utilizes different mechanisms to convey heritable traits to progeny cells that are independent of DNA sequence, including DNA silencing, post-translational modifications of histone proteins, and the post-transcriptional modulation of RNA transcript levels by non-coding RNAs. Although long non-coding RNAs have recently emerged as important regulators of gene imprinting, their functions during osteogenesis are as yet unexplored. In contrast, microRNAs (miRNAs) are well characterized for their control of osteogenic and osteoclastic pathways; thus, further defining how gene regulatory networks essential for skeleton functions are coordinated and finely tuned through the activities of miRNAs. Roles of miRNAs are constantly expanding as new studies uncover associations with skeletal disorders. The distinct functions of epigenetic regulators and evidence for integrating their activities to control normal bone gene expression and bone disease will be presented. In addition, potential for using "signature miRNAs" to identify, manage, and therapeutically treat osteosarcoma will be discussed in this review.

Modeling acclimatization by hybrid systems: condition changes alter biological system behavior models.

In order to describe the dynamic behavior of a complex biological system, it is useful to combine models integrating processes at different levels and with temporal dependencies. Such combinations are necessary for modeling acclimatization, a phenomenon where changes in environmental conditions can induce drastic changes in the behavior of a biological system. In this article we formalize the use of hybrid systems as a tool to model this kind of biological behavior. A modeling scheme called strong switches is proposed. It allows one to take into account both minor adjustments to the coefficients of a continuous model, and, more interestingly, large-scale changes to the structure of the model. We illustrate the proposed methodology with two applications: acclimatization in wine fermentation kinetics, and acclimatization of osteo-adipo differentiation system linking stimulus signals to bone mass.

Reusing and composing models of cell fate regulation of human bone precursor cells.

In order to treat osteoporosis and other bone mass disorders it is necessary to understand the regulatory processes that control the cell fate decisions responsible for going from bone precursor cells to bone tissue. Many processes interact to regulate cell division, differentiation and apoptosis. There are models for these basic processes, but not for their interactions. In this work we use the theory of switched systems, reuse and composition of validated models to describe the cell fate decisions leading to bone and fat formation. We describe the differentiation of osteo-adipo progenitor cells by composing its model with differentiation stimuli. We use the activation of the Wnt pathway as stimulus to osteoblast lineage, including regulation of cell division and apoptosis. This model is our first step to simulate physiological responses in silico to treatments for bone mass disorders.

Genomic promoter occupancy of runt-related transcription factor RUNX2 in Osteosarcoma cells identifies genes involved in cell adhesion and motility.

Runt-related transcription factors (RUNX1, RUNX2, and RUNX3) are key lineage-specific regulators of progenitor cell growth and differentiation but also function pathologically as cancer genes that contribute to tumorigenesis. RUNX2 attenuates growth and stimulates maturation of osteoblasts during bone formation but is also robustly expressed in a subset of osteosarcomas, as well as in metastatic breast and prostate tumors. To assess the biological function of RUNX2 in osteosarcoma cells, we examined human genomic promoter interactions for RUNX2 using chromatin immunoprecipitation (ChIP)-microarray analysis in SAOS-2 cells. Promoter binding of both RUNX2 and RNA polymerase II was compared with gene expression profiles of cells in which RUNX2 was depleted by RNA interference. Many RUNX2-bound loci (1550 of 2339 total) exhibit promoter occupancy by RNA polymerase II and contain the RUNX consensus motif 5'-((T/A/C)G(T/A/C)GG(T/G). Gene ontology analysis indicates that RUNX2 controls components of multiple signaling pathways (e.g. WNT, TGFß, TNFα, and interleukins), as well as genes linked to cell motility and adhesion (e.g. the focal adhesion-related genes FAK/PTK2 and TLN1). Our results reveal that siRNA depletion of RUNX2, PTK2, or TLN1 diminishes motility of U2OS osteosarcoma cells. Thus, RUNX2 binding to diverse gene loci may support the biological properties of osteosarcoma cells.

An architectural genetic and epigenetic perspective.

The organization and intranuclear localization of nucleic acids and regulatory proteins contribute to both genetic and epigenetic parameters of biological control. Regulatory machinery in the cell nucleus is functionally compartmentalized in microenvironments (focally organized sites where regulatory factors reside) that provide threshold levels of factors required for transcription, replication, repair and cell survival. The common denominator for nuclear organization of regulatory machinery is that each component of control is architecturally configured and every component of control is embedded in architecturally organized networks that provide an infrastructure for integration and transduction of regulatory signals. It is realistic to anticipate emerging mechanisms that account for the organization and assembly of regulatory complexes within the cell nucleus can provide novel options for cancer diagnosis and therapy with maximal specificity, reduced toxicity and minimal off-target complications.

Mitotic bookmarking of genes: a novel dimension to epigenetic control.

Regulatory machinery is focally organized in the interphase nucleus. The information contained in these focal nuclear microenvironments must be inherited during cell division to sustain physiologically responsive gene expression in progeny cells. Recent results suggest that focal mitotic retention of phenotypic transcription factors at promoters together with histone modifications and DNA methylation--a mechanism collectively known as gene bookmarking--is a novel parameter of inherited epigenetic control that sustains cellular identity after mitosis. The epigenetic signatures imposed by bookmarking poise genes for activation or suppression following mitosis. We discuss the implications of phenotypic transcription factor retention on mitotic chromosomes in biological control and disease.

The human SWI/SNF complex associates with RUNX1 to control transcription of hematopoietic target genes.

The acute myeloid leukemia 1 (AML1, RUNX1) transcription factor is a key regulator of hematopoietic differentiation that forms multi-protein complexes with co-regulatory proteins. These complexes are assembled at target gene promoters in nuclear microenvironments to mediate phenotypic gene expression and chromatin-related epigenetic modifications. Here, immunofluorescence microscopy and biochemical assays are used to show that RUNX1 associates with the human ATP-dependent SWI/SNF chromatin remodeling complex. The SWI/SNF subunits BRG1 and INI1 bind in vivo to RUNX1 target gene promoters (e.g., GMCSF, IL3, MCSF-R, MIP, and p21). These interactions correlate with histone modifications characteristic of active chromatin, including acetylated H4 and dimethylated H3 lysine 4. Downregulation of RUNX1 by RNA interference diminishes the binding of BRG1 and INI1 at selected target genes. Taken together, our findings indicate that RUNX1 interacts with the human SWI/SNF complex to control hematopoietic-specific gene expression.

Subnuclear targeting of the Runx3 tumor suppressor and its epigenetic association with mitotic chromosomes.

Runx proteins are tissue-specific transcriptional scaffolds that organize and assemble regulatory complexes at strategic sites of target gene promoters and at intranuclear foci to govern activation or repression. During interphase, fidelity of intranuclear targeting supports the biological activity of Runx1 and Runx2 proteins. Both factors regulate genes involved in cell cycle control and cell growth (e.g., rRNA genes), as well as lineage commitment. Here, we have examined the subcellular regulatory properties of the third Runx member, the tumor suppressor protein Runx3, during interphase and mitosis. Using in situ cellular and biochemical approaches we delineated a subnuclear targeting signal that directs Runx3 to discrete transcriptional foci that are nuclear matrix associated. Chromatin immunoprecipitation results show that Runx3 occupies rRNA promoters during interphase. We also find that Runx3 remains associated with chromosomes during mitosis and localizes with nucleolar organizing regions (NORs), reflecting an interaction with epigenetic potential. Taken together, our study establishes that common mechanisms control the subnuclear distribution and activities of Runx1, Runx2, and Runx3 proteins to support RNA polymerase I and II mediated gene expression during interphase and mitosis.

Phenotypic transcription factors epigenetically mediate cell growth control.

Ribosomal RNA (rRNA) genes are down-regulated during osteogenesis, myogenesis, and adipogenesis, necessitating a mechanistic understanding of interrelationships between growth control and phenotype commitment. Here, we show that cell fate-determining factors [MyoD, myogenin (Mgn), Runx2, C/EBPbeta] occupy rDNA loci and suppress rRNA expression during lineage progression, concomitant with decreased rRNA expression and reciprocal loss of occupancy by c-Myc, a proliferation-specific activator of rRNA transcription. We find interaction of phenotypic factors with the polymerase I activator upstream binding factor UBF-1 at interphase nucleoli, and this interaction is epigenetically retained on mitotic chromosomes at nucleolar organizing regions. Ectopic expression and RNA interference establish that MyoD, Mgn, Runx2, and C/EBPbeta each functionally suppress rRNA genes and global protein synthesis. We conclude that epigenetic control of ribosomal biogenesis by lineage-specific differentiation factors is a general developmental mechanism for coordinate control of cell growth and phenotype.

Mitotic occupancy and lineage-specific transcriptional control of rRNA genes by Runx2.

Regulation of ribosomal RNA genes is a fundamental process that supports the growth of cells and is tightly coupled with cell differentiation. Although rRNA transcriptional control by RNA polymerase I (Pol I) and associated factors is well studied, the lineage-specific mechanisms governing rRNA expression remain elusive. Runt-related transcription factors Runx1, Runx2 and Runx3 establish and maintain cell identity, and convey phenotypic information through successive cell divisions for regulatory events that determine cell cycle progression or exit in progeny cells. Here we establish that mammalian Runx2 not only controls lineage commitment and cell proliferation by regulating genes transcribed by RNA Pol II, but also acts as a repressor of RNA Pol I mediated rRNA synthesis. Within the condensed mitotic chromosomes we find that Runx2 is retained in large discrete foci at nucleolar organizing regions where rRNA genes reside. These Runx2 chromosomal foci are associated with open chromatin, co-localize with the RNA Pol I transcription factor UBF1, and undergo transition into nucleoli at sites of rRNA synthesis during interphase. Ribosomal RNA transcription and protein synthesis are enhanced by Runx2 deficiency that results from gene ablation or RNA interference, whereas induction of Runx2 specifically and directly represses rDNA promoter activity. Runx2 forms complexes containing the RNA Pol I transcription factors UBF1 and SL1, co-occupies the rRNA gene promoter with these factors in vivo, and affects local chromatin histone modifications at rDNA regulatory regions. Thus Runx2 is a critical mechanistic link between cell fate, proliferation and growth control. Our results suggest that lineage-specific control of ribosomal biogenesis may be a fundamental function of transcription factors that govern cell fate.

Transcriptional induction of the osteocalcin gene during osteoblast differentiation involves acetylation of histones h3 and h4.

The remodeling of chromatin is required for tissue-specific gene activation to permit interactions of transcription factors and coregulators with their cognate elements. Here, we investigate the chromatin-mediated mechanisms by which the bone-specific osteocalcin (OC) gene is transcriptionally activated during cessation of cell growth in ROS 17/2.8 osteosarcoma cells and during normal osteoblast differentiation. Acetylation of histones H3 and H4 at the OC gene promoter was assayed during the proliferative and postproliferative stages of cell growth by using chromatin immunoprecipitation assays with antibodies that recognize different acetylated forms of histones H3 or H4. The results show that the promoter and coding regions of the OC gene contain very low levels of acetylated histones H3 and H4 during the proliferative period of osteoblast differentiation when the OC gene is inactive. Active expression of the OC gene in mature osteoblasts and confluent ROS 17/2.8 cells is functionally linked to preferential acetylation of histone H4 and, to a lesser extent, to acetylation of histone H3. Histone acetylation at the loci for RUNX2 (CBFA1), alkaline phosphatase, bone sialoprotein, osteopontin, and the cell growth regulator p21, which are expressed throughout osteoblast differentiation, is not altered postproliferatively. We conclude that acetylation of histones H3 and H4 is functionally coupled to the chromatin remodeling events that mediate the developmental induction of OC gene transcription in bone cells.