The bone-specific Runx2-P1 promoter displays conserved three-dimensional chromatin structure with the syntenic Supt3h promoter

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.

Regulation and Function of FOXC1 in Osteoblasts

Estrogens, which bind to estrogen receptor alpha (ERα), are important for proper bone mineral density. When women go through menopause, estrogen levels decrease, and there is a decrease in bone quality, along with an increased risk for fractures. We previously identified an enhancer near FOXC1 as the most significantly enriched binding site for estrogen receptor alpha (ERα) in osteoblasts. FOXC1 is a transcription factor belonging to a large group of proteins known as forkhead box genes and is an important regulator of bone formation. Here, we demonstrate that 17β-estradiol (E2) increases the mRNA and protein levels of FOXC1 in primary mouse and human osteoblasts. GATA4 is a pioneer factor for ERα and it is also recruited to enhancers near Foxc1. Knockdown of Gata4 in mouse osteoblasts in vitro decreases Foxc1 expression as does knockout of Gata4 in vivo. Functionally, GATA4 and FOXC1 interact and regulate osteoblast proteins such as RUNX2, as demonstrated by ChIP-reChIP and luciferase assays. The most enriched motif in GATA4 binding sites from ChIP-seq is for FOXC1, supporting the notion that GATA4 and FOXC1 cooperate in regulating osteoblast differentiation. Together, these data demonstrate the interactions of the transcription factors ERα, GATA4, and FOXC1 to regulate each other's expression and other osteoblast differentiation genes.

Genetics implicates overactive osteogenesis in the development of diffuse idiopathic skeletal hyperostosis

Diffuse idiopathic skeletal hyperostosis (DISH) is a condition where adjacent vertebrae become fused through formation of osteophytes. The genetic and epidemiological etiology of this condition is not well understood. Here, we implemented a machine learning algorithm to assess the prevalence and severity of the pathology in ~40,000 lateral DXA scans in the UK Biobank Imaging cohort. We find that DISH is highly prevalent, above the age of 45, ~20% of men and ~8% of women having multiple osteophytes. Surprisingly, we find strong phenotypic and genetic association of DISH with increased bone mineral density and content throughout the entire skeletal system. Genetic association analysis identified ten loci associated with DISH, including multiple genes involved in bone remodeling (RUNX2, IL11, GDF5, CCDC91, NOG, and ROR2). Overall, this study describes genetics of DISH and implicates the role of overactive osteogenesis as a key driver of the pathology.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

Prolonged FOS activity disrupts a global myogenic transcriptional program by altering 3D chromatin architecture in primary muscle progenitor cells

BACKGROUND

The AP-1 transcription factor, FBJ osteosarcoma oncogene (FOS), is induced in adult muscle satellite cells (SCs) within hours following muscle damage and is required for effective stem cell activation and muscle repair. However, why FOS is rapidly downregulated before SCs enter cell cycle as progenitor cells (i.e., transiently expressed) remains unclear. Further, whether boosting FOS levels in the proliferating progeny of SCs can enhance their myogenic properties needs further evaluation.

METHODS

We established an inducible, FOS expression system to evaluate the impact of persistent FOS activity in muscle progenitor cells ex vivo. We performed various assays to measure cellular proliferation and differentiation, as well as uncover changes in RNA levels and three-dimensional (3D) chromatin interactions.

RESULTS

Persistent FOS activity in primary muscle progenitor cells severely antagonizes their ability to differentiate and form myotubes within the first 2 weeks in culture. RNA-seq analysis revealed that ectopic FOS activity in muscle progenitor cells suppressed a global pro-myogenic transcriptional program, while activating a stress-induced, mitogen-activated protein kinase (MAPK) transcriptional signature. Additionally, we observed various FOS-dependent, chromosomal re-organization events in A/B compartments, topologically associated domains (TADs), and genomic loops near FOS-regulated genes.

CONCLUSIONS

Our results suggest that elevated FOS activity in recently activated muscle progenitor cells perturbs cellular differentiation by altering the 3D chromosome organization near critical pro-myogenic genes. This work highlights the crucial importance of tightly controlling FOS expression in the muscle lineage and suggests that in states of chronic stress or disease, persistent FOS activity in muscle precursor cells may disrupt the muscle-forming process.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Shaping modern human skull through epigenetic, transcriptional and post-transcriptional regulation of the RUNX2 master bone gene

RUNX2 encodes the master bone transcription factor driving skeletal development in vertebrates, and playing a specific role in craniofacial and skull morphogenesis. The anatomically modern human (AMH) features sequence changes in the RUNX2 locus compared with archaic hominins' species. We aimed to understand how these changes may have contributed to human skull globularization occurred in recent evolution. We compared in silico AMH and archaic hominins' genomes, and used mesenchymal stromal cells isolated from skull sutures of craniosynostosis patients for in vitro functional assays. We detected 459 and 470 nucleotide changes in noncoding regions of the AMH RUNX2 locus, compared with the Neandertal and Denisovan genomes, respectively. Three nucleotide changes in the proximal promoter were predicted to alter the binding of the zinc finger protein Znf263 and long-distance interactions with other cis-regulatory regions. By surface plasmon resonance, we selected nucleotide substitutions in the 3'UTRs able to affect miRNA binding affinity. Specifically, miR-3150a-3p and miR-6785-5p expression inversely correlated with RUNX2 expression during in vitro osteogenic differentiation. The expression of two long non-coding RNAs, AL096865.1 and RUNX2-AS1, within the same locus, was modulated during in vitro osteogenic differentiation and correlated with the expression of specific RUNX2 isoforms. Our data suggest that RUNX2 may have undergone adaptive phenotypic evolution caused by epigenetic and post-transcriptional regulatory mechanisms, which may explain the delayed suture fusion leading to the present-day globular skull shape.

Skeletal and molecular findings in 51 Cleidocranial dysplasia patients from Turkey

Loss or decrease of function in runt-related transcription factor 2 encoded by RUNX2 is known to cause a rare autosomal-dominant skeletal disorder, cleidocranial dysplasia (CCD). Clinical spectrum and genetic findings in 51 CCD patients from 30 unrelated families are herein presented. In a majority of the patients, facial abnormalities, such as delayed fontanel closure (89%), parietal and frontal bossing (80%), metopic groove (77%), midface hypoplasia (94%), and abnormal mobility of shoulders (90%), were recorded following clinical examination. In approximately one-half of the subjects, wormian bone (51%), short stature (43%), bell-shaped thorax (42%), wide pubic symphysis (50%), hypoplastic iliac wing (59%), and chef's hat sign (44%) presented in available radiological examinations. Scoliosis was identified in 28% of the patients. Investigation of RUNX2 revealed small sequence alterations in 90% and gross deletions in 10% of the patients; collectively, 23 variants including 11 novel changes (c.29_30insT, c.203delAinsCG, c.423 + 2delT, c.443_454delTACCAGATGGGAinsG, c.505C > T, c.594_595delCTinsG, c.636_637insC, c.685 + 5G > A, c.1088G > T, c.1281delC, Exon 6-9 deletion) presented high allelic heterogeneity. Novel c.29_30insT is unique in affecting the P1-driven long isoform of RUNX2, which is expected to disrupt the N-terminal region of RUNX2; this was shown in two unrelated phenotypically discordant patients. The clinical findings highlighted mild intra-familial genotype-phenotype correlation in our CCD cohort.

Transcriptional Programming in Arteriosclerotic Disease: A Multifaceted Function of the Runx2 (Runt-Related Transcription Factor 2)

Despite successful therapeutic strategies in the prevention and treatment of arteriosclerosis, the cardiovascular complications remain a major clinical and societal issue worldwide. Increased vascular calcification promotes arterial stiffness and accelerates cardiovascular morbidity and mortality. Upregulation of the Runx2 (Runt-related transcription factor 2), an essential osteogenic transcription factor for bone formation, in the cardiovascular system has emerged as an important regulator for adverse cellular events that drive cardiovascular pathology. This review discusses the regulatory mechanisms that are critical for Runx2 expression and function and highlights the dynamic and complex cross talks of a wide variety of posttranslational modifications, including phosphorylation, acetylation, ubiquitination, and O-linked β-N-acetylglucosamine modification, in regulating Runx2 stability, cellular localization, and osteogenic transcriptional activity. How the activation of an array of signaling cascades by circulating and local microenvironmental factors upregulates Runx2 in vascular cells and promotes Runx2-mediated osteogenic transdifferentiation of vascular smooth muscle cells and expression of inflammatory cytokines that accelerate macrophage infiltration and vascular osteoclast formation is summarized. Furthermore, the increasing appreciation of a new role of Runx2 upregulation in promoting vascular smooth muscle cell phenotypic switch, and Runx2 modulated by O-linked β-N-acetylglucosamine modification and Runx2-dependent repression of smooth muscle cell-specific gene expression are discussed. Further exploring the regulation of this key osteogenic transcription factor and its new perspectives in the vasculature will provide novel insights into the transcriptional regulation of vascular smooth muscle cell phenotype switch, reprograming, and vascular inflammation that promote the pathogenesis of arteriosclerosis.

miR-27a-3p Targets ATF3 to Reduce Calcium Deposition in Vascular Smooth Muscle Cells

Vascular calcification, the ectopic deposition of calcium in blood vessels, develops in association with various metabolic diseases and atherosclerosis and is an independent predictor of morbidity and mortality associated with these diseases. Herein, we report that reduction of microRNA-27a-3p (miR-27a-3p) causes an increase in activating transcription factor 3 (ATF3), a novel osteogenic transcription factor, in vascular smooth muscle cells. Both microRNA (miRNA) and mRNA microarrays were performed with rat vascular smooth muscle cells, and reciprocally regulated pairs of miRNA and mRNA were selected after bioinformatics analysis. Inorganic phosphate significantly reduced the expression of miR-27a-3p in A10 cells. The transcript level was also reduced in vitamin D₃-administered mouse aortas. miR-27a-3p mimic reduced calcium deposition, whereas miR-27a-3p inhibitor increased it. The Atf3 mRNA level was upregulated in a cellular vascular calcification model, and miR-27a-3p reduced the Atf3 mRNA and protein levels. Transfection with Atf3 could recover the miR-27a-3p-induced reduction of calcium deposition. Our results suggest that reduction of miR-27a-3p may contribute to the development of vascular calcification by de-repression of ATF3.

A review of genetics of nasal development and morphological variation

The nose is central in the determination of facial esthetics. The variations in its structural characteristics greatly influence the ultimate dentoskeletal positioning at the end of an orthodontic therapy. A careful insight into its developmental etiology will greatly aid the health care professional in identifying patient's real concern about the facial appearance. This in turn will aid in the fabrication of a better treatment plan regarding the end placement goals for the teeth and jaws in all the three dimensions of space. However, this important structure is often missed as a part of the diagnostic and treatment planning regime owing to the lack of meticulous understanding of its developmental etiology by the orthodontists. The development of the nose in the embryo occurs in pre skeletal and skeletal phases by a well-coordinated and regulated interaction of multiple signaling cascades with the crucial importance of each factor in the entire mechanism. The five key factors, which control frontonasal development are sonic hedgehog (SHH), fibroblast growth factors (FGF), transforming growth factor β (TGFβ), wingless (WNT) proteins, and bone morphogenetic protein (BMP). The recent evidence suggests the association of various nasal dimensions and their related syndromes with multiple genes. The revelation of nasal genetic makeup in totality will aid in ascertaining the direction of growth, which will govern our orthodontic treatment results and will also act as a harbinger for potential genetic editing and tissue engineering. This article describes at length the morphological and genetic aspect of nasal growth and development in light of the gender and racial variability along with the emphasis on the importance of knowing these nasal features with regard to diagnosis and treatment planning in orthodontics.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Higher order genomic organization and epigenetic control maintain cellular identity and prevent breast cancer

Cells establish and sustain structural and functional integrity of the genome to support cellular identity and prevent malignant transformation. In this review, we present a strategic overview of epigenetic regulatory mechanisms including histone modifications and higher order chromatin organization (HCO) that are perturbed in breast cancer onset and progression. Implications for dysfunctions that occur in hormone regulation, cell cycle control, and mitotic bookmarking in breast cancer are considered, with an emphasis on epithelial-to-mesenchymal transition and cancer stem cell activities. The architectural organization of regulatory machinery is addressed within the contexts of translating cancer-compromised genomic organization to advances in breast cancer risk assessment, diagnosis, prognosis, and identification of novel therapeutic targets with high specificity and minimal off target effects.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

Nuclear organization mediates cancer-compromised genetic and epigenetic control

Nuclear organization is functionally linked to genetic and epigenetic regulation of gene expression for biological control and is modified in cancer. Nuclear organization supports cell growth and phenotypic properties of normal and cancer cells by facilitating physiologically responsive interactions of chromosomes, genes and regulatory complexes at dynamic three-dimensional microenvironments. We will review nuclear structure/function relationships that include: 1. Epigenetic bookmarking of genes by phenotypic transcription factors to control fidelity and plasticity of gene expression as cells enter and exit mitosis; 2. Contributions of chromatin remodeling to breast cancer nuclear morphology, metabolism and effectiveness of chemotherapy; 3. Relationships between fidelity of nuclear organization and metastasis of breast cancer to bone; 4. Dynamic modifications of higher-order inter- and intra-chromosomal interactions in breast cancer cells; 5. Coordinate control of cell growth and phenotype by tissue-specific transcription factors; 6. Oncofetal epigenetic control by bivalent histone modifications that are functionally related to sustaining the stem cell phenotype; and 7. Noncoding RNA-mediated regulation in the onset and progression of breast cancer. The discovery of components to nuclear organization that are functionally related to cancer and compromise gene expression have the potential for translation to innovative cancer diagnosis and targeted therapy.

Expression of the ectodomain-releasing protease ADAM17 is directly regulated by the osteosarcoma and bone-related transcription factor RUNX2

Osteoblast differentiation is controlled by transcription factor RUNX2 which temporally activates or represses several bone-related genes, including those encoding extracellular matrix proteins or factors that control cell-cell, and cell-matrix interactions. Cell-cell communication in the many skeletal pericellular micro-niches is critical for bone development and involves paracrine secretion of growth factors and morphogens. This paracrine signaling is in part regulated by "A Disintegrin And Metalloproteinase" (ADAM) proteins. These cell membrane-associated metalloproteinases support proteolytic release ("shedding") of protein ectodomains residing at the cell surface. We analyzed microarray and RNA-sequencing data for Adam genes and show that Adam17, Adam10, and Adam9 are stimulated during BMP2 mediated induction of osteogenic differentiation and are robustly expressed in human osteoblastic cells. ADAM17, which was initially identified as a tumor necrosis factor alpha (TNFα) converting enzyme also called (TACE), regulates TNFα-signaling pathway, which inhibits osteoblast differentiation. We demonstrate that Adam17 expression is suppressed by RUNX2 during osteoblast differentiation through the proximal Adam17 promoter region (-0.4 kb) containing two functional RUNX2 binding motifs. Adam17 downregulation during osteoblast differentiation is paralleled by increased RUNX2 expression, cytoplasmic-nuclear translocation and enhanced binding to the Adam17 proximal promoter. Forced expression of Adam17 reduces Runx2 and Alpl expression, indicating that Adam17 may negatively modulate osteoblast differentiation. These findings suggest a novel regulatory mechanism involving a reciprocal Runx2-Adam17 negative feedback loop to regulate progression through osteoblast differentiation. Our results suggest that RUNX2 may control paracrine signaling through regulation of ectodomain shedding at the cell surface of osteoblasts by directly suppressing Adam17 expression.

RUNX2 expression in thyroid and breast cancer requires the cooperation of three non-redundant enhancers under the control of BRD4 and c-JUN

Aberrant reactivation of embryonic pathways is a common feature of cancer. RUNX2 is a transcription factor crucial during embryogenesis that is aberrantly reactivated in many tumors, including thyroid and breast cancer, where it promotes aggressiveness and metastatic spreading. Currently, the mechanisms driving RUNX2 expression in cancer are still largely unknown. Here we showed that RUNX2 transcription in thyroid and breast cancer requires the cooperation of three distantly located enhancers (ENHs) brought together by chromatin three-dimensional looping. We showed that BRD4 controls RUNX2 by binding to the newly identified ENHs and we demonstrated that the anti-proliferative effects of bromodomain inhibitors (BETi) is associated with RUNX2 transcriptional repression. We demonstrated that each RUNX2 ENH is potentially controlled by a distinct set of TFs and we identified c-JUN as the principal pivot of this regulatory platform. We also observed that accumulation of genetic mutations within these elements correlates with metastatic behavior in human thyroid tumors. Finally, we identified RAINs, a novel family of ENH-associated long non-coding RNAs, transcribed from the identified RUNX2 regulatory unit. Our data provide a new model to explain how RUNX2 expression is reactivated in thyroid and breast cancer and how cancer-driving signaling pathways converge on the regulation of this gene.

GATA4 Directly Regulates Runx2 Expression and Osteoblast Differentiation

GATA4 is a zinc‐finger transcription factor that is a pioneer factor in various tissues and regulates tissue‐specific gene regulation. In vivo deletion of Gata4 using Cre‐recombinase under the control of the Col1a1 2.3 kb promoter showed significantly reduced values for trabecular bone properties by microCT analysis of femur and tibia of 14‐week‐old male and female mice, suggesting GATA4 is necessary for maintaining normal adult bone phenotype. Quantitative PCR analysis revealed higher expression of Gata4 in trabecular bone compared with cortical bone, suggesting a role for GATA4 in maintaining normal trabecular bone mass. In vivo and in vitro, reduction of Gata4 correlates with reduced Runx2 gene expression, along with reduced osteoblast mineralization. To determine if Runx2 is a direct target of GATA4, chromatin immunoprecipitation (ChIP) was performed, and it demonstrated that GATA4 is recruited to the two Runx2 promoters and an enhancer region. Furthermore, when Gata4 is knocked down, the chromatin at the Runx2 region is not open, as detected by DNase assays and ChIP with antibodies to the open chromatin marks H3K4me2 (histone 3 lysine 4 dimethylation) and H3K27ac (histone 3 lysine 27 acetylation) and the closed chromatin mark H3K27me2 (histone 3 lysine 27 trimethylation). Together, the data suggest that GATA4 binds near the Runx2 promoter and enhancer and helps maintain open chromatin to regulate Runx2 expression leading to bone mineralization. © 2017 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Genome-wide DNase hypersensitivity, and occupancy of RUNX2 and CTCF reveal a highly dynamic gene regulome during MC3T3 pre-osteoblast differentiation

The ability to discover regulatory sequences that control bone-related genes during development has been greatly improved by massively parallel sequencing methodologies. To expand our understanding of cis-regulatory regions critical to the control of gene expression during osteoblastogenesis, we probed the presence of open chromatin states across the osteoblast genome using global DNase hypersensitivity (DHS) mapping. Our profiling of MC3T3 mouse pre-osteoblasts during differentiation has identified more than 224,000 unique DHS sites. Approximately 65% of these sites are dynamic during temporal stages of osteoblastogenesis, and a majority of them are located within non-promoter (intergenic and intronic) regions. Nearly half of all DHS sites (both constitutive and dynamic) overlap binding events of the bone-essential RUNX2 and/or the chromatin-related CTCF transcription factors. This finding reinforces the role of these regulatory proteins as essential components of the bone gene regulome. We observe a reduction in chromatin accessibility throughout the genome between pre-osteoblast and early osteoblasts. Our analysis also defined a class of differentially expressed genes that harbor DHS peaks centered within 1 kb downstream of transcriptional end sites (TES). These DHSs at the 3’-flanks of genes exhibit dynamic changes during differentiation that may impact regulation of the osteoblast genome. Taken together, the distribution of DHS regions within non-promoter locations harboring osteoblast and chromatin related transcription factor binding motifs, reflect novel cis-regulatory requirements to support temporal gene expression in differentiating osteoblasts.

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.

Epigenetic regulation of bone cells

Bone is a major organ in the skeletal system that supports and protects muscles and other organs, facilitates movement and hematopoiesis, and forms a reservoir of minerals including calcium. The cells in the bone, such as osteoblasts, osteoclasts, and osteocytes, orchestrate sequential and balanced regulatory mechanisms to maintain bone and are capable of differentiating in bones. Bone development and remodeling require a precise regulation of gene expressions in bone cells, a process governed by epigenetic mechanisms such as histone modification, DNA methylation, and chromatin structure. Importantly, lineage-specific transcription factors can determine the epigenetic regulation of bone cells. Emerging data suggest that perturbation of epigenetic programs can affect the function and activity of bone cells and contributes to pathogenesis of bone diseases, including osteoporosis. Thus, understanding epigenetic regulations in bone cells would be important for early diagnosis and future therapeutic approaches.

Statin Function as an Anti-inflammation Therapy for Depression in Patients With Coronary Artery Disease by Downregulating Interleukin-1β

It is well known that inflammation contributes to the development of coronary artery disease (CAD) and depressive symptoms. Previous studies have shown that long-term application of statin reduces the occurrence of depression in patients with CAD. However, the mechanism remains unclear. We hypothesized that inflammation contributes to depression in patients with CAD and statin function as an anti-inflammation therapy for those depressive patients. Patients with confirmed CAD hospitalized in the Department of Cardiology of Tongji Hospital in Shanghai, China, were enrolled. Depression was identified as none (ND), mild (MiD), moderate (MoD), or severe (SD) on the basis of scores of the patient health questionnaire with 9 items. Inflammatory factors in peripheral blood were measured using a chemiluminescence immunoassay and Bio-plex. Luciferase expression level was detected using the Dual-Luciferase Reporter Assay System for IL-1β or NF-κB expression by transfection in human umbilical vein endothelial cells, and patient serum was added. Data obtained from 217 patients with CAD were analyzed. The IL-1β level of CAD with SD was 14.70, which was significantly higher than that of CAD with ND 7.52, MiD 7.73, or MoD 8.63. Luciferase reporter gene analysis showed that IL-1β or NF-κB expression level was upregulated by the serum of CAD and depression patients. After the addition of atorvastatin, IL-1β or NF-κB luciferase reporter expression level decreased. It suggested that depression in patients with CAD is associated with inflammation. Statin may function as an anti-inflammation therapy for depression in patients with CAD by downregulation of IL-1β.

Anabolic and Antiresorptive Modulation of Bone Homeostasis by the Epigenetic Modulator Sulforaphane, a Naturally Occurring Isothiocyanate

Bone degenerative pathologies like osteoporosis may be initiated by age-related shifts in anabolic and catabolic responses that control bone homeostasis. Here we show that sulforaphane (SFN), a naturally occurring isothiocyanate, promotes osteoblast differentiation by epigenetic mechanisms. SFN enhances active DNA demethylation viaTet1andTet2and promotes preosteoblast differentiation by enhancing extracellular matrix mineralization and the expression of osteoblastic markers (Runx2,Col1a1,Bglap2,Sp7,Atf4, andAlpl). SFN decreases the expression of the osteoclast activator receptor activator of nuclear factor-κB ligand (RANKL) in osteocytes and mouse calvarial explants and preferentially induces apoptosis in preosteoclastic cells via up-regulation of theTet1/Fas/Caspase 8 and Caspase 3/7 pathway. These mechanistic effects correlate with higher bone volume (∼20%) in both normal and ovariectomized mice treated with SFN for 5 weeks compared with untreated mice as determined by microcomputed tomography. This effect is due to a higher trabecular number in these mice. Importantly, no shifts in mineral density distribution are observed upon SFN treatment as measured by quantitative backscattered electron imaging. Our data indicate that the food-derived compound SFN epigenetically stimulates osteoblast activity and diminishes osteoclast bone resorption, shifting the balance of bone homeostasis and favoring bone acquisition and/or mitigation of bone resorptionin vivo Thus, SFN is a member of a new class of epigenetic compounds that could be considered for novel strategies to counteract osteoporosis.

Multiple levels of epigenetic control for bone biology and pathology

Multiple dimensions of epigenetic control contribute to regulation of gene expression that governs bone biology and pathology. Once confined to DNA methylation and a limited number of post-translational modifications of histone proteins, the definition of epigenetic mechanisms is expanding to include contributions of non-coding RNAs and mitotic bookmarking, a mechanism for retaining phenotype identity during cell proliferation. Together these different levels of epigenetic control of physiological processes and their perturbations that are associated with compromised gene expression during the onset and progression of disease, have contributed to an unprecedented understanding of the activities (operation) of the genomic landscape. Here, we address general concepts that explain the contribution of epigenetic control to the dynamic regulation of gene expression during eukaryotic transcription. This article is part of a Special Issue entitled Epigenetics and Bone.

Dynamic regulation of Schwann cell enhancers after peripheral nerve injury

Myelination of the peripheral nervous system is required for axonal function and long term stability. After peripheral nerve injury, Schwann cells transition from axon myelination to a demyelinated state that supports neuronal survival and ultimately remyelination of axons. Reprogramming of gene expression patterns during development and injury responses is shaped by the actions of distal regulatory elements that integrate the actions of multiple transcription factors. We used ChIP-seq to measure changes in histone H3K27 acetylation, a mark of active enhancers, to identify enhancers in myelinating rat peripheral nerve and their dynamics after demyelinating nerve injury. Analysis of injury-induced enhancers identified enriched motifs for c-Jun, a transcription factor required for Schwann cells to support nerve regeneration. We identify a c-Jun-bound enhancer in the gene for Runx2, a transcription factor induced after nerve injury, and we show that Runx2 is required for activation of other induced genes. In contrast, enhancers that lose H3K27ac after nerve injury are enriched for binding sites of the Sox10 and early growth response 2 (Egr2/Krox20) transcription factors, which are critical determinants of Schwann cell differentiation. Egr2 expression is lost after nerve injury, and many Egr2-binding sites lose H3K27ac after nerve injury. However, the majority of Egr2-bound enhancers retain H3K27ac, indicating that other transcription factors maintain active enhancer status after nerve injury. The global epigenomic changes in H3K27ac deposition pinpoint dynamic changes in enhancers that mediate the effects of transcription factors that control Schwann cell myelination and peripheral nervous system responses to nerve injury.