Genome-scale techniques highlight the epigenome and redefine fundamental principles of gene regulation

Molecular endocrinology of vitamin D on the epigenome level

The molecular endocrinology of vitamin D is based on the facts that i) its metabolite 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) is the high affinity ligand of the nuclear receptor vitamin D receptor (VDR) and ii) the transcription factor VDR is the unique target of 1,25(OH)₂D₃ in the nucleus. Short-term alterations of the epigenome are primarily changes in the post-translational modification status of nucleosome-forming histone proteins, the consequences of which are i) a local increase or decrease in chromatin accessibility and ii) the activation or repression of gene transcription. Vitamin D has via VDR a direct effect on the expression of several hundred primary target genes implying numerous effects on the epigenome. Next-generation sequencing methods, such as ChIP-seq and FAIRE-seq, were applied to cellular model systems of vitamin D signaling, such as THP-1 human monocytes, and provided data for a chromatin model of vitamin D signaling. Key points of this model are that i) in the absence of ligand VDR binds to a limited number of loci within accessible chromatin, ii) a stimulation with ligand increases the number of DNA-bound VDR molecules, iii) VDR's access to genomic DNA is supported by pioneer factors, such as PU.1 in monocytes, iv) VDR binding leads to local opening of chromatin and v) the binding strength of topologically associating domain anchor forming CCCTC-binding factor sites upstream and downstream of prominent VDR binding sites is changing in response to ligand stimulation. This model provides the present basis of the molecular endocrinology of vitamin D and will be in future refined by the integration of vitamin D-sensitive chromatin markers and other genome-wide data, such as the 1,25(OH)₂D₃-sensitive binding of co-factors, chromatin modifying enzymes and chromatin remodeling proteins.

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.

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.

Osteocytes, not Osteoblasts or Lining Cells, are the Main Source of the RANKL Required for Osteoclast Formation in Remodeling Bone

The cytokine receptor activator of nuclear factor kappa B ligand (RANKL), encoded by the Tnfsf11 gene, is essential for osteoclastogenesis and previous studies have shown that deletion of the Tnfsf11 gene using a Dmp1-Cre transgene reduces osteoclast formation in cancellous bone by more than 70%. However, the Dmp1-Cre transgene used in those studies leads to recombination in osteocytes, osteoblasts, and lining cells making it unclear whether one or more of these cell types produce the RANKL required for osteoclast formation in cancellous bone. Because osteoblasts, osteocytes, and lining cells have distinct locations and functions, distinguishing which of these cell types are sources of RANKL is essential for understanding the orchestration of bone remodeling. To distinguish between these possibilities, we have now created transgenic mice expressing the Cre recombinase under the control of regulatory elements of the Sost gene, which is expressed in osteocytes but not osteoblasts or lining cells in murine bone. Activity of the Sost-Cre transgene in osteocytes, but not osteoblast or lining cells, was confirmed by crossing Sost-Cre transgenic mice with tdTomato and R26R Cre-reporter mice, which express tdTomato fluorescent protein or LacZ, respectively, only in cells expressing the Cre recombinase or their descendants. Deletion of the Tnfsf11 gene in Sost-Cre mice led to a threefold decrease in osteoclast number in cancellous bone and increased cancellous bone mass, mimicking the skeletal phenotype of mice in which the Tnfsf11 gene was deleted using the Dmp1-Cre transgene. These results demonstrate that osteocytes, not osteoblasts or lining cells, are the main source of the RANKL required for osteoclast formation in remodeling cancellous bone.

Fundamentals of vitamin D hormone-regulated gene expression

Initial research focused upon several known genetic targets provided early insight into the mechanism of action of the vitamin D hormone (1,25-dihydroxyvitamin D3 (1,25(OH)2D3)). Recently, however, a series of technical advances involving the coupling of chromatin immunoprecipitation (ChIP) to unbiased methodologies that initially involved tiled DNA microarrays (ChIP-chip analysis) and now Next Generation DNA Sequencing techniques (ChIP-seq analysis) has opened new avenues of research into the mechanisms through which 1,25(OH)2D3 regulates gene expression. In this review, we summarize briefly the results of this early work and then focus on more recent studies in which ChIP-chip and ChIP-seq analyses have been used to explore the mechanisms of 1,25(OH)2D3 action on a genome-wide scale providing specific target genes as examples. The results of this work have advanced our understanding of the mechanisms involved at both genetic and epigenetic levels and have revealed a series of new principles through which the vitamin D hormone functions to control the expression of genes. This article is part of a Special Issue entitled '16th Vitamin D Workshop'.

Identification of an enhancer that increases miR-200b~200a~429 gene expression in breast cancer cells

The miR-200b~200a~429 gene cluster is a key regulator of EMT and cancer metastasis, however the transcription-based mechanisms controlling its expression during this process are not well understood. We have analyzed the miR-200b~200a~429 locus for epigenetic modifications in breast epithelial and mesenchymal cell lines using chromatin immunoprecipitation assays and DNA methylation analysis. We discovered a novel enhancer located approximately 5.1kb upstream of the miR-200b~200a~429 transcriptional start site. This region was associated with the active enhancer chromatin signature comprising H3K4me1, H3K27ac, RNA polymerase II and CpG dinucleotide hypomethylation. Luciferase reporter assays revealed the upstream enhancer stimulated the transcription of the miR-200b~200a~429 minimal promoter region approximately 27-fold in breast epithelial cells. Furthermore, we found that a region of the enhancer was transcribed, producing a short, GC-rich, mainly nuclear, non-polyadenylated RNA transcript designated miR-200b eRNA. Over-expression of miR-200b eRNA had little effect on miR-200b~200a~429 promoter activity and its production did not correlate with miR-200b~200a~429 gene expression. While additional investigations of miR-200b eRNA function will be necessary, it is possible that miR-200b eRNA may be involved in the regulation of miR-200b~200a~429 gene expression and silencing. Taken together, these findings reveal the presence of a novel enhancer, which contributes to miR-200b~200a~429 transcriptional regulation in epithelial cells.

A straw man's neogenome

The neogenome has indeed changed how to understand the relationship between genotype and phenotype. However, this does not imply a paradigm shift, but simply a normal development of a young science. Charney creates a straw man out of the myth of an immutable genetics, and conveys the wrong idea that heritability studies and gene association studies are no longer valid.

Twist2 contributes to termination of limb bud outgrowth and patterning through direct regulation of Grem1

Twist1 has been demonstrated to play critical roles in the early development of neural crest and mesodermally derived tissues including the limb. Twist2 has been less well characterised but its relatively late onset of expression suggests specific roles in the development of a number of organs. Expression of Twist2 within the developing limbs begins after formation of the limb bud and persists within the peripheral mesenchyme until digital rays condense. We have used RCAS-mediated overexpression in chick to investigate the function of Twist2 in limb development. Viral misexpression following injection into the lateral plate mesoderm results in a spectrum of hypoplastic limb phenotypes. These include generalized shortening of the entire limb, fusion of the autopod skeletal elements, loss of individual digits or distal truncation resulting in complete loss of the autopod. These phenotypes appear to result from a premature termination of limb outgrowth and manifest as defective growth in both the proximal-distal and anterior-posterior axes. In situ hybridisation analysis demonstrates that many components of the Shh/Grem1/Fgf regulatory loop that controls early limb growth and patterning are downregulated by Twist2 overexpression. Grem1 has a complementary expression pattern to Twist2 within the limb primordia and co-expression of both Grem1 and Twist2 results in a rescue of the Twist2 overexpression phenotype. We demonstrate that Twist proteins directly repress Grem1 expression via a regulatory element downstream of the open reading frame. These data indicate that Twist2 regulates early limb morphogenesis through a role in terminating the Shh/Grem1/Fgf autoregulatory loop.

The need for mouse models in osteoporosis genetics research

Osteoporosis, the progressive loss of bone mass resulting in fragility fractures, affects ∼75 million people in the United States, Europe and Japan. Bone mineral density (BMD) correlates with fracture risk and is widely used in clinical settings to predict fracture. Numerous studies have demonstrated that peak bone mass is highly heritable and consequently a number of genome-wide association studies (GWASs) have been conducted to identify the genes that regulate BMD. Traditional intercross mapping in the mouse has met with limited successes in the field of skeletal biology. With the advent of human GWAS, questions have arisen about the continued need for mouse models in genetics research. However, significant advances have been made in the field of mouse genetics, including new genetics resource populations and loci mapping techniques, which enable gene-level mapping resolution. In this review, we discuss the need for mouse models to help understand the skeletal biology underlying novel human GWAS findings, how loci discovered in the mouse can be used to complement GWAS analysis and highlight the recent advances made in the field of skeletal biology from the use of these new and developing resources. We conclude this paper with a discussion of the need for systems-level approaches in the skeletal biology field, with an emphasis on the need for pathway and network analyses.

Regulation of target gene expression by the vitamin D receptor - an update on mechanisms

Virtually all of the known biological actions of the hormonal ligand 1,25-dihydroxyvitamin D(3) (1,25(OH)(2)D(3)) are mediated by the vitamin D receptor (VDR). Following binding and activation by the ligand, the VDR localizes in the nucleus to the regulatory regions of target genes and recruits chromatin-active coregulatory complexes which, in turn, modulate transcriptional output. The failure of the VDR to function due to crippling mutations results in total hereditary resistance to 1,25(OH)(2)D(3) in both mice and humans. In this review, we summarize the structural and functional properties of the VDR and the role of 1,25(OH)(2)D(3) in receptor activation, and then describe the results of recent studies using genome-wide analyses that define the overarching principles through which the VDR modulates genes expression. We also focus on the recent analysis of a specific 1,25(OH)(2)D(3) regulated gene that provides confirmation of the principles identified through these genome-wide methodologies. Taken together, these studies suggest an unanticipated increase in the complexity of the molecular processes that govern gene regulation by hormones and other factors.

Regulation of mouse Cyp24a1 expression via promoter-proximal and downstream-distal enhancers highlights new concepts of 1,25-dihydroxyvitamin D(3) action

CYP24A1 functions in vitamin D target tissues to degrade 1,25-dihydroxyvitamin D(3) (1,25(OH)(2)D(3)). Thus, the concentration of this enzyme and the regulation of its expression is a primary determinant of the overall biological activity of 1,25(OH)(2)D(3) within cells. The principle regulator of CYP24A1 expression is 1,25(OH)(2)D(3) itself, which functions through the vitamin D receptor to upregulate the transcriptional activity of the Cyp24a1 gene. In this report, we explore the mechanism of this regulation using recently developed ChIP-chip and ChIP-seq techniques that permit an unbiased search for enhancer elements that participate in this transcriptional control. Our studies both confirm a regulatory region defined earlier and located proximal to the transcriptional start site (TSS) of mouse Cyp24a1 (-160 and -265nt) and identify a novel intergenic region located downstream of the transcription unit that contains two enhancers (+35 and +37kb) that facilitate 1,25(OH)(2)D(3)-dependent upregulation of Cyp24a1 expression. Interestingly, while C/EBPβ also binds under basal conditions to a site located immediately upstream of the Cyp24a1 promoter (-345nt), occupancy by this factor is strikingly increased following 1,25(OH)(2)D(3) treatment. The locations and activities of these regulatory regions that mediate 1,25(OH)(2)D(3) actions were confirmed in mice in vivo. We conclude that the mechanism through which 1,25(OH)(2)D(3) induces the CYP24A1 enzyme, thereby autoregulating its own destruction, involves both promoter-proximal as well as downstream-distal enhancers. These findings highlight new concepts regarding the molecular mechanism of action of 1,25(OH)(2)D(3) and other hormonal regulators.

Genome-wide principles of gene regulation by the vitamin D receptor and its activating ligand

The vitamin D receptor (VDR) mediates virtually all of the known biological actions of the hormonal ligand 1,25-dihydroxyvitamin D(3) (1,25(OH)(2)D(3)). These actions are directed toward the nucleus, where the VDR binds to the regulatory regions of target genes and modulates their transcriptional output. Recent technological advances have enabled the study of transcription factor binding on a genome-wide scale in cells and tissues that are major targets of vitamin D action. In this review, the results of several of these studies are discussed wherein overarching principles of gene regulation by the vitamin D hormone are beginning to emerge. In addition, several specific genes that are regulated by 1,25(OH)(2)D(3) and which provide new insight into the increasingly complex mechanism whereby the receptor functions to modulate gene expression are considered. These studies suggest that while many of the principles that are now accepted regarding the regulation of gene expression by hormones and other regulatory factors are well grounded, others require extensive modification.

PI3Kγ deletion reduces variability in the in vivo osteolytic response induced by orthopaedic wear particles

Orthopedic wear particles activate a number of intracellular signaling pathways associated with inflammation in macrophages and we have previously shown that the phosphoinositol-3-kinase (PI3K)/Akt pathway is one of the signal transduction pathways that mediates the in vitro activation of macrophages by orthopedic wear particles. Since PI3Kγ is primarily responsible for PI3K activity during inflammation, we hypothesized that PI3Kγ mediates particle-induced osteolysis in vivo. Our results do not strongly support the hypothesis that PI3Kγ regulates the overall amount of particle-induced osteolysis in the murine calvarial model. However, our results strongly support the conclusion that variability in the amount of particle-induced osteolysis between individual mice is reduced in the PI3Kγ(-/-) mice. These results suggest that PI3Kγ contributes to osteolysis to different degrees in individual mice and that the mice, and patients, that are most susceptible to osteolysis may be so, in part, due to an increased contribution from PI3Kγ.