Comparative epigenomic analysis of murine and human adipogenesis

Genome-wide map of nuclear protein degradation shows NCoR1 turnover as a key to mitochondrial gene regulation

Transcription factor activity and turnover are functionally linked, but the global patterns by which DNA-bound regulators are eliminated remain poorly understood. We established an assay to define the chromosomal location of DNA-associated proteins that are slated for degradation by the ubiquitin-proteasome system. The genome-wide map described here ties proteolysis in mammalian cells to active enhancers and to promoters of specific gene families. Nuclear-encoded mitochondrial genes in particular correlate with protein elimination, which positively affects their transcription. We show that the nuclear receptor corepressor NCoR1 is a key target of proteolysis and physically interacts with the transcription factor CREB. Proteasome inhibition stabilizes NCoR1 in a site-specific manner and restrains mitochondrial activity by repressing CREB-sensitive genes. In conclusion, this functional map of nuclear proteolysis links chromatin architecture with local protein stability and identifies proteolytic derepression as highly dynamic in regulating the transcription of genes involved in energy metabolism.

In Vivo and in vitro dynamics of undifferentiated embryonic cell transcription factor 1

Pluripotent stem cells retain the ability to differentiate into the three germ layers and germline. As a result, there is a major interest in characterizing regulators that establish and maintain pluripotency. The network of transcription factors continues to expand in complexity, and one factor, undifferentiated embryonic cell transcription factor 1 (UTF1), has recently moved more into the limelight. To facilitate the study of UTF1, we report the generation and characterization of two reporter lines that enable efficient tracking, mapping, and purification of endogenous UTF1. In particular, we include a built-in biotinylation system in our targeted locus that allows efficient and reliable pulldown. We also use this reporter to show the dynamic regulation of Utf1 in distinct stem cell conditions and demonstrate its utility for reprogramming studies. The multipurpose design of the reporter lines enables many directions of future study and should lead to a better understanding of UTF1's diverse roles.

Integrated analysis of transcript-level regulation of metabolism reveals disease-relevant nodes of the human metabolic network

Metabolic diseases and comorbidities represent an ever-growing epidemic where multiple cell types impact tissue homeostasis. Here, the link between the metabolic and gene regulatory networks was studied through experimental and computational analysis. Integrating gene regulation data with a human metabolic network prompted the establishment of an open-sourced web portal, IDARE (Integrated Data Nodes of Regulation), for visualizing various gene-related data in context of metabolic pathways. Motivated by increasing availability of deep sequencing studies, we obtained ChIP-seq data from widely studied human umbilical vein endothelial cells. Interestingly, we found that association of metabolic genes with multiple transcription factors (TFs) enriched disease-associated genes. To demonstrate further extensions enabled by examining these networks together, constraint-based modeling was applied to data from human preadipocyte differentiation. In parallel, data on gene expression, genome-wide ChIP-seq profiles for peroxisome proliferator-activated receptor (PPAR) γ, CCAAT/enhancer binding protein (CEBP) α, liver X receptor (LXR) and H3K4me3 and microRNA target identification for miR-27a, miR-29a and miR-222 were collected. Disease-relevant key nodes, including mitochondrial glycerol-3-phosphate acyltransferase (GPAM), were exposed from metabolic pathways predicted to change activity by focusing on association with multiple regulators. In both cell types, our analysis reveals the convergence of microRNAs and TFs within the branched chain amino acid (BCAA) metabolic pathway, possibly providing an explanation for its downregulation in obese and diabetic conditions.

Combinatorial regulation of lipoprotein lipase by microRNAs during mouse adipogenesis

MicroRNAs (miRNAs) regulate gene expression directly through base pairing to their targets or indirectly through participating in multi-scale regulatory networks. Often miRNAs take part in feed-forward motifs where a miRNA and a transcription factor act on shared targets to achieve accurate regulation of processes such as cell differentiation. Here we show that the expression levels of miR-27a and miR-29a inversely correlate with the mRNA levels of lipoprotein lipase (Lpl), their predicted combinatorial target, and its key transcriptional regulator peroxisome proliferator-activated receptor gamma (Pparg) during 3T3-L1 adipocyte differentiation. More importantly, we show that Lpl, a key lipogenic enzyme, can be negatively regulated by the two miRNA families in a combinatorial fashion on the mRNA and functional level in maturing adipocytes. This regulation is mediated through the Lpl 3'UTR as confirmed by reporter gene assays. In addition, a small mathematical model captures the dynamics of this feed-forward motif and predicts the changes in Lpl mRNA levels upon network perturbations. The obtained results might offer an explanation to the dysregulation of LPL in diabetic conditions and could be extended to quantitative modeling of regulation of other metabolic genes under similar regulatory network motifs.

Epigenetics in adipose tissue, obesity, weight loss, and diabetes

Given the role that diet and other environmental factors play in the development of obesity and type 2 diabetes, the implication of different epigenetic processes is being investigated. Although it is well known that external factors can cause cell type-dependent epigenetic changes, including DNA methylation, histone tail modifications, and chromatin remodeling, the regulation of these processes, the magnitude of the changes and the cell types in which they occur, the individuals more predisposed, and the more crucial stages of life remain to be elucidated. There is evidence that obese and diabetic people have a pattern of epigenetic marks different from nonobese and nondiabetic individuals. The main long-term goals in this field are the identification and understanding of the role of epigenetic marks that could be used as early predictors of metabolic risk and the development of drugs or diet-related treatments able to delay these epigenetic changes and even reverse them. But weight gain and insulin resistance/diabetes are influenced not only by epigenetic factors; different epigenetic biomarkers have also been identified as early predictors of weight loss and the maintenance of body weight after weight loss. The characterization of all the factors that are able to modify the epigenetic signatures and the determination of their real importance are hindered by the following factors: the magnitude of change produced by dietary and environmental factors is small and cumulative; there are great differences among cell types; and there are many factors involved, including age, with multiple interactions between them.

Genome-wide profiling of transcription factor binding and epigenetic marks in adipocytes by ChIP-seq

The recent advances in high-throughput sequencing combined with various other technologies have allowed detailed and genome-wide insight into the transcriptional networks that control adipogenesis. Chromatin immunoprecipitation (ChIP) combined with high-throughput sequencing (ChIP-seq) is one of the most widely used of these technologies. Using these methods, association of transcription factors, cofactors, and epigenetic marks can be mapped to DNA in a genome-wide manner. Here, we provide a detailed protocol for performing ChIP-seq analyses in preadipocytes and adipocytes. We have focused mainly on critical points, limitations of the assay, and quality controls required in order to obtain reproducible ChIP-seq data.

Impaired leptin gene expression and release in cultured preadipocytes isolated from individuals born with low birth weight

Low birth weight (LBW) is associated with increased risk of the development of type 2 diabetes (T2D). The appetite-regulating hormone leptin is released from mature adipocytes, and its production may be decreased in immature preadipocytes from LBW individuals. We recruited 14 men born with LBW and 13 controls born with normal birth weight (NBW). Biopsy samples were obtained from subcutaneous abdominal fat depots, and preadipocytes were isolated and cultured. Gene expression of leptin and selected differentiation markers were analyzed during preadipocyte differentiation, and cell culture media were collected to analyze leptin secretion. DNA methylation of CpG sites in the leptin promoter was measured using pyrosequencing. We found that differentiating preadipocytes from LBW individuals showed reduced leptin gene expression and a corresponding reduced leptin release compared with NBW individuals. Mean DNA methylation of the proximal LEP promoter was increased in LBW compared with NBW individuals. The notion of impaired adipocyte maturation in LBW individuals was supported by a lower mRNA expression of the differentiation markers; fatty acid binding protein 4, peroxisome proliferator-activated receptor γ, and GLUT4. Our findings are consistent with impaired preadipocyte maturation, contributing to an increased risk of the development of T2D in LBW subjects.

The obesogen tributyltin

The obesogen hypothesis postulates the role of environmental chemical pollutants that disrupt homeostatic controls and adaptive mechanisms to promote adipose-dependent weight gain leading to obesity and metabolic syndrome complications. One of the most direct molecular mechanisms for coupling environmental chemical exposures to perturbed physiology invokes pollutants mimicking endogenous endocrine hormones or bioactive dietary signaling metabolites that serve as nuclear receptor ligands. The organotin pollutant tributyltin can exert toxicity through multiple mechanisms but most recently has been shown to bind, activate, and mediate RXR-PPARγ transcriptional regulation central to lipid metabolism and adipocyte biology. Data in support of long-term obesogenic effects on whole body adipose tissue are also reported. Organotins represent an important model test system for evaluating the impact and epidemiological significance of chemical insults as contributing factors for obesity and human metabolic health.

Peroxisome proliferator-activated receptor γ and C/EBPα synergistically activate key metabolic adipocyte genes by assisted loading

Peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer binding protein α (C/EBPα) are key activators of adipogenesis. They mutually induce the expression of each other and have been reported to cooperate in activation of a few adipocyte genes. Recently, genome-wide profiling revealed a high degree of overlap between PPARγ and C/EBPα binding in adipocytes, suggesting that cooperativeness could be mediated through common binding sites. To directly investigate the interplay between PPARγ and C/EBPα at shared binding sites, we established a fibroblastic model system in which PPARγ and C/EBPα can be independently expressed. Using RNA sequencing, we demonstrate that coexpression of PPARγ and C/EBPα leads to synergistic activation of many key metabolic adipocyte genes. This is associated with extensive C/EBPα-mediated reprogramming of PPARγ binding and vice versa in the vicinity of these genes, as determined by chromatin immunoprecipitation combined with deep sequencing. Our results indicate that this is at least partly mediated by assisted loading involving chromatin remodeling directed by the leading factor. In conclusion, we report a novel mechanism by which the key adipogenic transcription factors, PPARγ and C/EBPα, cooperate in activation of the adipocyte gene program.

Chromatin modifiers and remodellers: regulators of cellular differentiation

Cellular differentiation is, by definition, epigenetic. Genome-wide profiling of pluripotent cells and differentiated cells suggests global chromatin remodelling during differentiation, which results in a progressive transition from a fairly open chromatin configuration to a more compact state. Genetic studies in mouse models show major roles for a variety of histone modifiers and chromatin remodellers in key developmental transitions, such as the segregation of embryonic and extra-embryonic lineages in blastocyst stage embryos, the formation of the three germ layers during gastrulation and the differentiation of adult stem cells. Furthermore, rather than merely stabilizing the gene expression changes that are driven by developmental transcription factors, there is emerging evidence that chromatin regulators have multifaceted roles in cell fate decisions.

Understanding variation in transcription factor binding by modeling transcription factor genome-epigenome interactions

Despite explosive growth in genomic datasets, the methods for studying epigenomic mechanisms of gene regulation remain primitive. Here we present a model-based approach to systematically analyze the epigenomic functions in modulating transcription factor-DNA binding. Based on the first principles of statistical mechanics, this model considers the interactions between epigenomic modifications and a cis-regulatory module, which contains multiple binding sites arranged in any configurations. We compiled a comprehensive epigenomic dataset in mouse embryonic stem (mES) cells, including DNA methylation (MeDIP-seq and MRE-seq), DNA hydroxymethylation (5-hmC-seq), and histone modifications (ChIP-seq). We discovered correlations of transcription factors (TFs) for specific combinations of epigenomic modifications, which we term epigenomic motifs. Epigenomic motifs explained why some TFs appeared to have different DNA binding motifs derived from in vivo (ChIP-seq) and in vitro experiments. Theoretical analyses suggested that the epigenome can modulate transcriptional noise and boost the cooperativity of weak TF binding sites. ChIP-seq data suggested that epigenomic boost of binding affinities in weak TF binding sites can function in mES cells. We showed in theory that the epigenome should suppress the TF binding differences on SNP-containing binding sites in two people. Using personal data, we identified strong associations between H3K4me2/H3K9ac and the degree of personal differences in NFκB binding in SNP-containing binding sites, which may explain why some SNPs introduce much smaller personal variations on TF binding than other SNPs. In summary, this model presents a powerful approach to analyze the functions of epigenomic modifications. This model was implemented into an open source program APEG (Affinity Prediction by Epigenome and Genome, http://systemsbio.ucsd.edu/apeg).

Emerging roles of zinc finger proteins in regulating adipogenesis

Proteins containing the zinc finger domain(s) are named zinc finger proteins (ZFPs), one of the largest classes of transcription factors in eukaryotic genomes. A large number of ZFPs have been studied and many of them were found to be involved in regulating normal growth and development of cells and tissues through diverse signal transduction pathways. Recent studies revealed that a small but increasing number of ZFPs could function as key transcriptional regulators involved in adipogenesis. Due to the prevalence of obesity and metabolic disorders, the investigation of molecular regulatory mechanisms of adipocyte development must be more completely understood in order to develop novel and long-term impact strategies for ameliorating obesity. In this review, we discuss recent work that has documented that ZFPs are important functional contributors to the regulation of adipogenesis. Taken together, these data lead to the conclusion that ZFPs may become promising targets to combat human obesity.

Genomic and epigenomic regulation of adipose tissue inflammation in obesity

Chronic inflammation of adipose tissue is viewed as a hallmark of obesity and contributes to the development of type 2 diabetes and cardiovascular disease. According to current models, nutrient excess causes metabolic and structural changes in adipocytes, which initiate transcriptional programs leading to the expression of inflammatory molecules and the subsequent recruitment of immune cells. Recent advances in deciphering the underlying mechanisms revealed that key regulatory events occur at the genomic and epigenomic levels. Here we review these advances because they offer a better understanding of the mechanisms behind the complex obesogenic program in adipose tissue, and because they may help in defining new therapeutic strategies that prevent, restrict, and resolve inflammation in the context of obesity.

The orphan nuclear receptors at their 25-year reunion

The nuclear receptor superfamily includes many receptors, identified based on their similarity to steroid hormone receptors but without a known ligand. The study of how these receptors are diversely regulated to interact with genomic regions to control a plethora of biological processes has provided critical insight into development, physiology, and the molecular pathology of disease. Here we provide a compendium of these so-called orphan receptors and focus on what has been learned about their modes of action, physiological functions, and therapeutic promise.

Peroxisome proliferator-activated receptor γ regulates genes involved in insulin/insulin-like growth factor signaling and lipid metabolism during adipogenesis through functionally distinct enhancer classes

The nuclear receptor peroxisome proliferator-activated receptor (PPAR) is a transcription factor whose expression is induced during adipogenesis and that is required for the acquisition and control of mature adipocyte functions. Indeed, PPAR induces the expression of genes involved in lipid synthesis and storage through enhancers activated during adipocyte differentiation. Here, we show that PPAR also binds to enhancers already active in preadipocytes as evidenced by an active chromatin state including lower DNA methylation levels despite higher CpG content. These constitutive enhancers are linked to genes involved in the insulin/insulin-like growth factor signaling pathway that are transcriptionally induced during adipogenesis but to a lower extent than lipid metabolism genes, because of stronger basal expression levels in preadipocytes. This is consistent with the sequential involvement of hormonal sensitivity and lipid handling during adipocyte maturation and correlates with the chromatin structure dynamics at constitutive and activated enhancers. Interestingly, constitutive enhancers are evolutionary conserved and can be activated in other tissues, in contrast to enhancers controlling lipid handling genes whose activation is more restricted to adipocytes. Thus, PPAR utilizes both broadly active and cell type-specific enhancers to modulate the dynamic range of activation of genes involved in the adipogenic process.

Closing the (nuclear) envelope on the genome: how nuclear lamins interact with promoters and modulate gene expression

The nuclear envelope shapes the functional organization of the nucleus. Increasing evidence indicates that one of its main components, the nuclear lamina, dynamically interacts with the genome, including the promoter region of specific genes. This seems to occur in a manner that accords developmental significance to these interactions. This essay addresses key issues raised by recent data on the association of nuclear lamins with the genome. We discuss how lamins interact with large chromatin domains and with spatially restricted regions on gene promoters. We address the relationship between these interactions, chromatin modifications and gene expression outcomes. Lamin-genome contacts are redistributed after cell division and during stem cell differentiation, with evidence of lineage specificity. Thus, we also speculate on a developmental role of lamin interactions with specific genes. Finally, we highlight how concepts arising from this recent work lay the foundations of future challenges and investigations.

PLZF expression during colorectal cancer development and in normal colorectal mucosa according to body size, as marker of colorectal cancer risk

Promyelocytic leukemia zinc finger protein (PLZF) is a protein involved in various signaling, growth regulatory, and differentiation pathways, including development/function of some T cells. Here, we aimed at the detection of PLZF during colorectal carcinogenesis, using immunofluorescence, and at the evaluation of the colocalization of PLZF with CD2 and CD56 positive cells (T, γδ, NK, and NKT cells), using confocal-microscopy, along colorectal carcinogenesis, since its earliest stages, that is, dysplastic aberrant crypt foci (ACF). Furthermore, we analyzed PLZF in the normal colonic mucosa (NM) according to anthropometric parameters of the subject. NM exhibited strong CD56 fluorescent staining. This infiltration was lost in both ACF and colorectal carcinoma (CRC), while PLZF presence increased from NM to ACF and CRC. Strong association was found between CD56+ colonic mucosa cell infiltration and body mass index. Interestingly, an increased stromal PLZF-reactivity was present in NM of obese subjects. This study shows that overexpression of PLZF and exclusion of NK cells in dysplastic microenvironment are very early events in the stepwise sequence leading to CRC and that lower levels of CD56+ cells in NM, together with increased levels of PLZF+ cells, can be a reflection of colon cancer risk due to obesity.

Beyond the ENCODE project: using genomics and epigenomics strategies to study enhancer evolution

The complex expression patterns observed for many genes are often regulated by distal transcription enhancers. Changes in the nucleotide sequences of enhancers may therefore lead to changes in gene expression, representing a central mechanism by which organisms evolve. With the development of the experimental technique of chromatin immunoprecipitation (ChIP), in which discrete regions of the genome bound by specific proteins can be identified, it is now possible to identify transcription factor binding events (putative cis-regulatory elements) in entire genomes. Comparing protein-DNA binding maps allows us, for the first time, to attempt to identify regulatory differences and infer global patterns of change in gene expression across species. Here, we review studies that used genome-wide ChIP to study the evolution of enhancers. The trend is one of high divergence of cis-regulatory elements between species, possibly compensated by extensive creation and loss of regulatory elements and rewiring of their target genes. We speculate on the meaning of the differences observed and discuss that although ChIP experiments identify the biochemical event of protein-DNA interaction, it cannot determine whether the event results in a biological function, and therefore more studies are required to establish the effect of divergence of binding events on species-specific gene expression.

Nuclear lamins: making contacts with promoters

The nuclear lamina guards the genome and in many ways contributes to regulating nuclear function. Increasing evidence indicates that the lamina dynamically interacts with chromatin mainly through large repressive domains, and recent data suggest that at least some of the lamin-genome contacts may be developmentally significant. In an attempt to provide an additional meaning to lamin-genome contacts, a recent study characterized the association of gene promoters with A-type lamins in progenitor and differentiated cells. Here, we discuss how A-type lamins interact with spatially defined promoter regions, and the relationship between these interactions, associated chromatin marks and gene expression outputs. We discuss the impact of A-type lamins on nucleus-wide and local chromatin organization. We also address how lamin-promoter interactions are redistributed during differentiation of adipocyte progenitors into adipocytes. Finally, we propose a model of lineage-specific "unlocking" of developmentally regulated loci and its significance in cellular differentiation.

Adipogenic transcriptome profiling using high throughput technologies

The recent technological innovations in the area of functional genomics, gene expression/transcriptomic profiling can provide new insights to understand the molecular basis of adipogenesis. The focus of this review is to highlight the recent advances in our understanding of the complex interplay of gene expression events and the regulatory mechanisms of adipogenesis in mouse cell lines, humans and farm animals. All these studies have employed the availability of constantly evolving high throughput 'omics' technologies including microarrays, RNA-Seq, chromatin immunoprecipitation, next generation sequencing, RNAi, miRNA profiling and quantitative PCR arrays.

Early B-cell factor-1 (EBF1) is a key regulator of metabolic and inflammatory signaling pathways in mature adipocytes

EBF1 plays a crucial role in early adipogenesis; however, despite high expression in mature adipocytes, its function in these cells is currently unknown. To identify direct and indirect EBF1 targets in fat, we undertook a combination of transcriptional profiling of EBF1-deficient adipocytes and genome-wide EBF1 location analysis. Our results indicate that many components of metabolic and inflammatory pathways are positively and directly regulated by EBF1, including PI3K/AKT, MAPK, and STAT1 signaling. Accordingly, we observed significant reduction of multiple signaling events in EBF1 knockdown cells as well as a reduction in insulin-stimulated glucose uptake and lipogenesis. Inflammatory signaling, gene expression, and secretion of inflammatory cytokines were also significantly affected by loss of EBF1 in adipocytes, although ChIP-sequencing results suggest that these actions are indirect. We also found that EBF1 occupies some 35,000 sites in adipocytes, most of which occur in enhancers. Significantly, comparison with three other published EBF1 ChIP-sequencing data sets in B-cells reveals both gene- and cell type-specific patterns of EBF1 binding. These results advance our understanding of the transcriptional mechanisms regulating signaling pathways in mature fat cells and indicate that EBF1 functions as a key integrator of signal transduction, inflammation, and metabolism.

A wavelet-based method to exploit epigenomic language in the regulatory region

MOTIVATION

Epigenetic landscapes in the regulatory regions reflect binding condition of transcription factors and their co-factors. Identifying epigenetic condition and its variation is important in understanding condition-specific gene regulation. Computational approaches to explore complex multi-dimensional landscapes are needed.

RESULTS

To study epigenomic condition for gene regulation, we developed a method, AWNFR, to classify epigenomic landscapes based on the detected epigenomic landscapes. Assuming mixture of Gaussians for a nucleosome, the proposed method captures the shape of histone modification and identifies potential regulatory regions in the wavelet domain. For accuracy estimation as well as enhanced computational speed, we developed a novel algorithm based on down-sampling operation and footprint in wavelet. We showed the algorithmic advantages of AWNFR using the simulated data. AWNFR identified regulatory regions more effectively and accurately than the previous approaches with the epigenome data in mouse embryonic stem cells and human lung fibroblast cells (IMR90). Based on the detected epigenomic landscapes, AWNFR classified epigenomic status and studied epigenomic codes. We studied co-occurring histone marks and showed that AWNFR captures the epigenomic variation across time.

AVAILABILITY AND IMPLEMENTATION

The source code and supplemental document of AWNFR are available at http://wonk.med.upenn.edu/AWNFR.

Analysis of in vitro insulin-resistance models and their physiological relevance to in vivo diet-induced adipose insulin resistance

Diet-induced obesity (DIO) predisposes individuals to insulin resistance, and adipose tissue has a major role in the disease. Insulin resistance can be induced in cultured adipocytes by a variety of treatments, but what aspects of the in vivo responses are captured by these models remains unknown. We use global RNA sequencing to investigate changes induced by TNF-α, hypoxia, dexamethasone, high insulin, and a combination of TNF-α and hypoxia, comparing the results to the changes in white adipose tissue from DIO mice. We found that different in vitro models capture distinct features of DIO adipose insulin resistance, and a combined treatment of TNF-α and hypoxia is most able to mimic the in vivo changes. Using genome-wide DNase I hypersensitivity followed by sequencing, we further examined the transcriptional regulation of TNF-α-induced insulin resistance, and we found that C/EPBβ is a potential key regulator of adipose insulin resistance.

Mapping human epigenomes

As the second dimension to the genome, the epigenome contains key information specific to every type of cells. Thousands of human epigenome maps have been produced in recent years thanks to rapid development of high throughput epigenome mapping technologies. In this review, we discuss the current epigenome mapping toolkit and utilities of epigenome maps. We focus particularly on mapping of DNA methylation, chromatin modification state, and chromatin structures, and emphasize the use of epigenome maps to delineate human gene regulatory sequences and developmental programs. We also provide a perspective on the progress of the epigenomics field and challenges ahead.

Global differences in specific histone H3 methylation are associated with overweight and type 2 diabetes

Background

Epidemiological evidence indicates yet unknown epigenetic mechanisms underlying a propensity for overweight and type 2 diabetes. We analyzed the extent of methylation at lysine 4 and lysine 9 of histone H3 in primary human adipocytes from 43 subjects using modification-specific antibodies.

Results

The level of lysine 9 dimethylation was stable, while adipocytes from type 2 diabetic and non-diabetic overweight subjects exhibited about 40% lower levels of lysine 4 dimethylation compared with cells from normal-weight subjects. In contrast, trimethylation at lysine 4 was 40% higher in adipocytes from overweight diabetic subjects compared with normal-weight and overweight non-diabetic subjects. There was no association between level of modification and age of subjects.

Conclusions

The findings define genome-wide molecular modifications of histones in adipocytes that are directly associated with overweight and diabetes, and thus suggest a molecular basis for existing epidemiological evidence of epigenetic inheritance.

DNMT1 is regulated by ATP-citrate lyase and maintains methylation patterns during adipocyte differentiation

During adipocyte differentiation, significant epigenomic changes occur in association with the implementation of the adipogenic program. We have previously shown that histone acetylation increases during differentiation in a manner dependent on acetyl coenzyme A (acetyl-CoA) production by the enzyme ATP-citrate lyase (ACL). Whether ACL regulates nuclear targets in addition to histones during differentiation is not clear. In this study, we report that DNA methyltransferase 1 (DNMT1) levels in adipocytes are controlled in part by ACL and that silencing of DNMT1 can accelerate adipocyte differentiation. DNMT1 gene expression is induced early in 3T3-L1 adipocyte differentiation during mitotic clonal expansion and is critical for maintenance of DNA and histone H3K9 methylation patterns during this period. In the absence of DNMT1, adipocyte-specific gene expression and lipid accumulation occur precociously. Later in differentiation, DNMT1 levels decline in an ACL-dependent manner. ACL-mediated suppression of DNMT1 occurs at least in part by promoting expression of microRNA 148a (miR-148a), which represses DNMT1. Ectopic expression of miR-148a accelerates differentiation under standard conditions and can partially rescue a hypermethylation-mediated differentiation block. The data suggest a role for DNMT1 in modulating the timing of differentiation and describe a novel ACL-miR-148a-dependent mechanism for regulating DNMT1 during adipogenesis.

Acute genome-wide effects of rosiglitazone on PPARγ transcriptional networks in adipocytes

Peroxisome proliferator-activated receptor γ (PPARγ) is a master regulator of adipocyte differentiation, and genome-wide studies indicate that it is involved in the induction of most adipocyte genes. Here we report, for the first time, the acute effects of the synthetic PPARγ agonist rosiglitazone on the transcriptional network of PPARγ in adipocytes. Treatment with rosiglitazone for 1 hour leads to acute transcriptional activation as well as repression of a number of genes as determined by genome-wide RNA polymerase II occupancy. Unlike what has been shown for many other nuclear receptors, agonist treatment does not lead to major changes in the occurrence of PPARγ binding sites. However, rosiglitazone promotes PPARγ occupancy at many preexisting sites, and this is paralleled by increased occupancy of the mediator subunit MED1. The increase in PPARγ and MED1 binding is correlated with an increase in transcription of nearby genes, indicating that rosiglitazone, in addition to activating the receptor, also promotes its association with DNA, and that this is causally linked to recruitment of mediator and activation of genes. Notably, both rosiglitazone-activated and -repressed genes are induced during adipogenesis. However, rosiglitazone-activated genes are markedly more associated with PPARγ than repressed genes and are highly dependent on PPARγ for expression in adipocytes. By contrast, repressed genes are associated with the other key adipocyte transcription factor CCAAT-enhancer binding proteinα (C/EBPα), and their expression is more dependent on C/EBPα. This suggests that the relative occupancies of PPARγ and C/EBPα are critical for whether genes will be induced or repressed by PPARγ agonist.

Lamin A/C-promoter interactions specify chromatin state-dependent transcription outcomes

The nuclear lamina is implicated in the organization of the eukaryotic nucleus. Association of nuclear lamins with the genome occurs through large chromatin domains including mostly, but not exclusively, repressed genes. How lamin interactions with regulatory elements modulate gene expression in different cellular contexts is unknown. We show here that in human adipose tissue stem cells, lamin A/C interacts with distinct spatially restricted subpromoter regions, both within and outside peripheral and intra-nuclear lamin-rich domains. These localized interactions are associated with distinct transcriptional outcomes in a manner dependent on local chromatin modifications. Down-regulation of lamin A/C leads to dissociation of lamin A/C from promoters and remodels repressive and permissive histone modifications by enhancing transcriptional permissiveness, but is not sufficient to elicit gene activation. Adipogenic differentiation resets a large number of lamin-genome associations globally and at subpromoter levels and redefines associated transcription outputs. We propose that lamin A/C acts as a modulator of local gene expression outcome through interaction with adjustable sites on promoters, and that these position-dependent transcriptional readouts may be reset upon differentiation.

The evolution of lineage-specific regulatory activities in the human embryonic limb

The evolution of human anatomical features likely involved changes in gene regulation during development. However, the nature and extent of human-specific developmental regulatory functions remain unknown. We obtained a genome-wide view of cis-regulatory evolution in human embryonic tissues by comparing the histone modification H3K27ac, which provides a quantitative readout of promoter and enhancer activity, during human, rhesus, and mouse limb development. Based on increased H3K27ac, we find that 13% of promoters and 11% of enhancers have gained activity on the human lineage since the human-rhesus divergence. These gains largely arose by modification of ancestral regulatory activities in the limb or potential co-option from other tissues and are likely to have heterogeneous genetic causes. Most enhancers that exhibit gain of activity in humans originated in mammals. Gains at promoters and enhancers in the human limb are associated with increased gene expression, suggesting they include molecular drivers of human morphological evolution.

The role of MyoD1 and histone modifications in the activation of muscle enhancers

MyoD1 is a key regulator that orchestrates skeletal muscle differentiation through the regulation of gene expression. Although many studies have focused on its role in transcriptional control at gene promoters, less is known regarding the role of MyoD1 in the assembly of active enhancers. Here, we discuss novel data that point to the ability of MyoD1 to mediate the assembly of active enhancers that augment the transcription of genes essential for muscle development and lineage specification. Based on genome-wide studies of epigenetic marks that typify active enhancers, we recently identified the compendium of distal regulatory elements that dictate transcriptional programs during myogenesis. Superimposition of MyoD1 binding sites upon the locations of muscle enhancers revealed its unequivocal binding to a core region of nearly a third of condition-specific muscle enhancers. Further studies exploring deposition of enhancer-related epigenetic marks in myoblasts lacking MyoD1 demonstrate the dependence of muscle enhancer assembly on the presence of MyoD1. We propose a model wherein MyoD1 mediates recruitment of Set7, H3K4me1, H3K27ac, p300, and RNAP II to MyoD1-bound enhancers to establish condition-specific activation of muscle genes. Moreover, muscle enhancers are modulated through coordinated binding of transcription factors, including c-Jun, Jdp2, Meis, and Runx1, which are recruited to muscle enhancers in a MyoD1-dependent manner. Thus, MyoD1 and enhancer-associated transcription factors function coordinately to assemble and regulate enhancers, thereby augmenting expression of muscle-related genes.

Chromatin and epigenetic features of long-range gene regulation

The precise regulation of gene transcription during metazoan development is controlled by a complex system of interactions between transcription factors, histone modifications and modifying enzymes and chromatin conformation. Developments in chromosome conformation capture technologies have revealed that interactions between regions of chromatin are pervasive and highly cell-type specific. The movement of enhancers and promoters in and out of higher-order chromatin structures within the nucleus are associated with changes in expression and histone modifications. However, the factors responsible for mediating these changes and determining enhancer:promoter specificity are still not completely known. In this review, we summarize what is known about the patterns of epigenetic and chromatin features characteristic of elements involved in long-range interactions. In addition, we review the insights into both local and global patterns of chromatin interactions that have been revealed by the latest experimental and computational methods.

Transcriptional regulation of endothelial dysfunction in atherosclerosis: an epigenetic perspective

Atherosclerosis is a progressive human pathology that encompasses several stages of development. Endothelial dysfunction represents an early sign of lesion within the vasculature. A number of risk factors for atherosclerosis, including hyperlipidemia, diabetes, and hypertension, target the vascular endothelium by re-programming its transcriptome. These profound alterations taking place on the chromatin rely on the interplay between sequence specific transcription factors and the epigenetic machinery. The epigenetic machinery, in turn, tailor individual transcription events key to atherogenesis to intrinsic and extrinsic insults dictating the development of atherosclerotic lesions. This review summarizes our current understanding of the involvement of the epigenetic machinery in endothelial injury during atherogenesis.

An aberrant transcription factor network essential for Wnt signaling and stem cell maintenance in glioblastoma

Glioblastoma (GBM) is thought to be driven by a subpopulation of cancer stem cells (CSCs) that self-renew and recapitulate tumor heterogeneity yet remain poorly understood. Here, we present a comparative analysis of chromatin state in GBM CSCs that reveals widespread activation of genes normally held in check by Polycomb repressors. These activated targets include a large set of developmental transcription factors (TFs) whose coordinated activation is unique to the CSCs. We demonstrate that a critical factor in the set, ASCL1, activates Wnt signaling by repressing the negative regulator DKK1. We show that ASCL1 is essential for the maintenance and in vivo tumorigenicity of GBM CSCs. Genome-wide binding profiles for ASCL1 and the Wnt effector LEF-1 provide mechanistic insight and suggest widespread interactions between the TF module and the signaling pathway. Our findings demonstrate regulatory connections among ASCL1, Wnt signaling, and collaborating TFs that are essential for the maintenance and tumorigenicity of GBM CSCs.

Transcriptional and epigenetic dynamics during specification of human embryonic stem cells

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

Bridging epigenomics and complex disease: the basics

The DNA sequence largely defines gene expression and phenotype. However, it is becoming increasingly clear that an additional chromatin-based regulatory network imparts both stability and plasticity to genome output, modifying phenotype independently of the genetic blueprint. Indeed, alterations in this "epigenetic" control layer underlie, at least in part, the reason for monozygotic twins being discordant for disease. Functionally, this regulatory layer comprises post-translational modifications of DNA and histones, as well as small and large noncoding RNAs. Together these regulate gene expression by changing chromatin organization and DNA accessibility. Successive technological advances over the past decade have enabled researchers to map the chromatin state with increasing accuracy and comprehensiveness, catapulting genetic research into a genome-wide era. Here, aiming particularly at the genomics/epigenomics newcomer, we review the epigenetic basis that has helped drive the technological shift and how this progress is shaping our understanding of complex disease.

Absolute quantification of transcription factors during cellular differentiation using multiplexed targeted proteomics

The cellular abundance of transcription factors (TFs) is an important determinant of their regulatory activities. Deriving TF copy numbers is therefore crucial to understanding how these proteins control gene expression. We describe a sensitive selected reaction monitoring-based mass spectrometry assay that allowed us to determine the copy numbers of up to ten proteins simultaneously. We applied this approach to profile the absolute levels of key TFs, including PPARγ and RXRα, during terminal differentiation of mouse 3T3-L1 pre-adipocytes. Our analyses revealed that individual TF abundance differs dramatically (from ∼250 to >300,000 copies per nucleus) and that their dynamic range during differentiation can vary up to fivefold. We also formulated a DNA binding model for PPARγ based on TF copy number, binding energetics and local chromatin state. This model explains the increase in PPARγ binding sites during the final differentiation stage that occurs despite a concurrent saturation in PPARγ copy number.

RNA-seq identified a super-long intergenic transcript functioning in adipogenesis

RNA transcripts are generally classified into polyA-plus and polyA-minus subgroups due to the presence or absence of a polyA tail at the 3′ end. Even though a number of physiologically and pathologically important polyA-minus RNAs have been recently identified, a systematic analysis of the expression and function of these transcripts in adipogenesis is still elusive. To study the potential function of the polyA-minus RNAs in adipogenesis, a dynamic expressional profiling was performed in the induced differentiation of 3T3-L1 cells. In addition to identifying thousands of novel intergenic transcripts, differentiation-synchronized expression was characterized for many of them. Among these, several large intergenic transcripts were found to be upregulated by more than 19-fold during differentiation. Further study demonstrated a fat tissue-specific expression pattern for these regions and identified an adipogenesis-associated long non-coding RNA. Collectively, these lines of evidence contribute to the characterization of a super-long intergenic transcript functioning in adipogenesis.

EBF2 determines and maintains brown adipocyte identity

The master transcription factor Pparγ regulates the general differentiation program of both brown and white adipocytes. However, it has been unclear whether Pparγ also controls fat lineage-specific characteristics. Here, we show that early B cell factor-2 (Ebf2) regulates Pparγ binding activity to determine brown versus white adipocyte identity. The Ebf DNA-binding motif was highly enriched within brown adipose-specific Pparγ binding sites that we identified by genome-wide ChIP-Seq. Of the Ebf isoforms, Ebf2 was selectively expressed in brown relative to white adipocytes and was bound at brown adipose-specific Pparγ target genes. When expressed in myoblasts or white preadipose cells, Ebf2 recruited Pparγ to its brown-selective binding sites and reprogrammed cells to a brown fat fate. Brown adipose cells and tissue from Ebf2-deficient mice displayed a loss of brown-specific characteristics and thermogenic capacity. Together, these results identify Ebf2 as a key transcriptional regulator of brown fat cell fate and function.

Enabling interspecies epigenomic comparison with CEpBrowser

SUMMARY

We developed the Comparative Epigenome Browser (CEpBrowser) to allow the public to perform multi-species epigenomic analysis. The web-based CEpBrowser integrates, manages and visualizes sequencing-based epigenomic datasets. Five key features were developed to maximize the efficiency of interspecies epigenomic comparisons.

AVAILABILITY

CEpBrowser is a web application implemented with PHP, MySQL, C and Apache. URL: http://www.cepbrowser.org/.

Identification and analysis of murine pancreatic islet enhancers

AIMS/HYPOTHESIS

The paucity of information on the epigenetic barriers that are blocking reprogramming protocols, and on what makes a beta cell unique, has hampered efforts to develop novel beta cell sources. Here, we aimed to identify enhancers in pancreatic islets, to understand their developmental ontologies, and to identify enhancers unique to islets to increase our understanding of islet-specific gene expression.

METHODS

We combined H3K4me1-based nucleosome predictions with pancreatic and duodenal homeobox 1 (PDX1), neurogenic differentiation 1 (NEUROD1), v-Maf musculoaponeurotic fibrosarcoma oncogene family, protein A (MAFA) and forkhead box A2 (FOXA2) occupancy data to identify enhancers in mouse islets.

RESULTS

We identified 22,223 putative enhancer loci in in vivo mouse islets. Our validation experiments suggest that nearly half of these loci are active in regulating islet gene expression, with the remaining regions probably poised for activity. We showed that these loci have at least nine developmental ontologies, and that islet enhancers predominately acquire H3K4me1 during differentiation. We next discriminated 1,799 enhancers unique to islets and showed that these islet-specific enhancers have reduced association with annotated genes, and identified a subset that are instead associated with novel islet-specific long non-coding RNAs (lncRNAs).

CONCLUSIONS/INTERPRETATIONS

Our results indicate that genes with islet-specific expression and function tend to have enhancers devoid of histone methylation marks or, less often, that are bivalent or repressed, in embryonic stem cells and liver. Further, we identify a subset of enhancers unique to islets that are associated with novel islet-specific genes and lncRNAs. We anticipate that these data will facilitate the development of novel sources of functional beta cell mass.

Oslo Epigenetics Symposium 2012. Oslo, Norway, 8-9 November 2012

The Oslo Epigenetics Symposium 2012 held in Oslo, Norway, brought together ten speakers from several European countries and the USA for an evening public lecture and a full day of presentations on emerging topics in the field of epigenetics, gene regulation and organization of the cell nucleus.

Genome-wide map of quantified epigenetic changes during in vitro chondrogenic differentiation of primary human mesenchymal stem cells

Background

For safe clinical application of engineered cartilage made from mesenchymal stem cells (MSCs), molecular mechanisms for chondrogenic differentiation must be known in detail. Changes in gene expression and extracellular matrix synthesis have been extensively studied, but the epigenomic modifications underlying these changes have not been described. To this end we performed whole-genome chromatin immunoprecipitation and deep sequencing to quantify six histone modifications, reduced representation bisulphite sequencing to quantify DNA methylation and mRNA microarrays to quantify gene expression before and after 7 days of chondrogenic differentiation of MSCs in an alginate scaffold. To add to the clinical relevance of our observations, the study is based on primary bone marrow-derived MSCs from four donors, allowing us to investigate inter-individual variations.

Results

We see two levels of relationship between epigenetic marking and gene expression. First, a large number of genes ontogenetically linked to MSC properties and the musculoskeletal system are epigenetically prepatterned by moderate changes in H3K4me3 and H3K9ac near transcription start sites. Most of these genes remain transcriptionally unaltered. Second, transcriptionally upregulated genes, more closely associated with chondrogenesis, are marked by H3K36me3 in gene bodies, highly increased H3K4me3 and H3K9ac on promoters and 5' end of genes, and increased H3K27ac and H3K4me1 marking in at least one enhancer region per upregulated gene. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression. Inter-donor variability analysis shows high level of similarity between the donors for this data set.

Conclusions

Histone modifications, rather than DNA methylation, provide the primary epigenetic control of early differentiation of MSCs towards the chondrogenic lineage.

Long noncoding RNAs regulate adipogenesis

The prevalence of obesity has led to a surge of interest in understanding the detailed mechanisms underlying adipocyte development. Many protein-coding genes, mRNAs, and microRNAs have been implicated in adipocyte development, but the global expression patterns and functional contributions of long noncoding RNA (lncRNA) during adipogenesis have not been explored. Here we profiled the transcriptome of primary brown and white adipocytes, preadipocytes, and cultured adipocytes and identified 175 lncRNAs that are specifically regulated during adipogenesis. Many lncRNAs are adipose-enriched, strongly induced during adipogenesis, and bound at their promoters by key transcription factors such as peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer-binding protein α (CEBPα). RNAi-mediated loss of function screens identified functional lncRNAs with varying impact on adipogenesis. Collectively, we have identified numerous lncRNAs that are functionally required for proper adipogenesis.

Genome-wide identification of enhancers in skeletal muscle: the role of MyoD1

To identify the compendium of distal regulatory elements that govern myogenic differentiation, we generated chromatin state maps based on histone modifications and recruitment of factors that typify enhancers in myoblasts and myotubes. We found a striking concordance between the locations of these newly defined enhancers, MyoD1-binding events, and noncoding RNA transcripts. These enhancers recruit several sequence-specific transcription factors in a spatially constrained manner around MyoD1-binding sites. Remarkably, MyoD1-null myoblasts show a wholesale loss of recruitment of these factors as well as diminished monomethylation of H3K4 (H3K4me1) and acetylation of H3K27 (H3K27ac) and reduced recruitment of Set7, an H3K4 monomethylase. Surprisingly, we found that H3K4me1, but not H3K27ac, could be restored by re-expression of MyoD1 in MyoD1(-/-) myoblasts, although re-expression of this factor in MyoD1-null myotubes restored both histone modifications. Our studies identified a role for MyoD1 in condition-specific enhancer assembly through recruitment of transcription factors and histone-modifying enzymes that shape muscle differentiation.

Dual functions of TAF7L in adipocyte differentiation

The diverse transcriptional mechanisms governing cellular differentiation and development of mammalian tissue remains poorly understood. Here we report that TAF7L, a paralogue of TFIID subunit TAF7, is enriched in adipocytes and white fat tissue (WAT) in mouse. Depletion of TAF7L reduced adipocyte-specific gene expression, compromised adipocyte differentiation, and WAT development as well. Ectopic expression of TAF7L in myoblasts reprograms these muscle precursors into adipocytes upon induction. Genome-wide mRNA-seq expression profiling and ChIP-seq binding studies confirmed that TAF7L is required for activating adipocyte-specific genes via a dual mechanism wherein it interacts with PPARγ at enhancers and TBP/Pol II at core promoters. In vitro binding studies confirmed that TAF7L forms complexes with both TBP and PPARγ. These findings suggest that TAF7L plays an integral role in adipocyte gene expression by targeting enhancers as a cofactor for PPARγ and promoters as a component of the core transcriptional machinery.DOI:http://dx.doi.org/10.7554/eLife.00170.001.

What obesity research tells us about epigenetic mechanisms

The pathophysiology of obesity is extremely complex and is associated with extensive gene expression changes in tissues throughout the body. This situation, combined with the fact that all gene expression changes are thought to have associated epigenetic changes, means that the links between obesity and epigenetics will undoubtedly be vast. Much progress in identifying epigenetic changes induced by (or inducing) obesity has already been made, with candidate and genome-wide approaches. These discoveries will aid the clinician through increasing our understanding of the inheritance, development and treatment of obesity. However, they are also of great value for epigenetic researchers, as they have revealed mechanisms of environmental interactions with epigenetics that can produce or perpetuate a disease state. Here, we will review the evidence for four mechanisms through which epigenetics contributes to obesity: as downstream effectors of environmental signals; through abnormal global epigenetic state driving obesogenic expression patterns; through facilitating developmental programming and through transgenerational epigenetic inheritance.

The epigenome and its role in diabetes

Both genetic and environmental factors play critical roles in the development of diabetes. Epidemiological evidence and data from clinical studies suggest the persistence of a "metabolic memory" of past exposures to environmental factors or glycemic control. Epigenetic mechanisms are regarded as one of the likeliest candidates underlying these phenomena. On the other hand, owing to the recent elucidation of mechanisms that erase epigenetic marks, it has gradually become recognized that epigenetic regulation is a more dynamic process than previously thought. A technological breakthrough in epigenome research in the past decade was the development of high-throughput sequencing. This new technology lets us investigate the epigenome in a global and comprehensive manner, and provides previously unrecognized findings and insights. This review presents an overview of the recent progress in our understanding of epigenetic regulation in type 1 and type 2 diabetes research.