Comparative epigenomic analysis of murine and human adipogenesis

A compendium of genome-wide hematopoietic transcription factor maps supports the identification of gene regulatory control mechanisms

OBJECTIVE

Key regulators of blood stem cell differentiation into the various mature hematopoietic lineages are commonly encoded by transcription factor genes. Elucidation of transcriptional regulatory mechanisms therefore holds great promise in advancing our understanding of both normal and malignant hematopoiesis. Recent technological advances have enabled the generation of genome-wide transcription factor binding maps using chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq). However, transcription factors operate in a combinatorial fashion suggesting that integrated analysis of genome-wide maps for multiple transcription factors will be essential to fully exploit these new genome-scale data sets.

MATERIALS AND METHODS

Here we have generated a compendium that integrates 53 ChIP-Seq studies covering 30 factors across all major hematopoietic lineages with a total of 754,380 binding peaks. We also used transgenic mouse assays to validate a newly predicted transcriptional enhancer.

RESULTS

Integrated analysis of all 53 ChIP-Seq studies demonstrated that cell-type identity exerts a larger influence on global transcription factor binding patterns than the nature of the individual transcription factors. Furthermore, regions highlighted by multifactor binding within specific gene loci overlap with known regulatory elements and also provide a useful guide for identifying novel elements, as demonstrated by transgenic analysis of a previously unrecognized enhancer in the Maml3 gene locus.

CONCLUSIONS

The ChIP-Seq compendium described here provides a valuable resource for the wider research community by accelerating the discovery of transcriptional mechanisms operating in the hematopoietic system.

Increased SRF transcriptional activity in human and mouse skeletal muscle is a signature of insulin resistance

Insulin resistance in skeletal muscle is a key phenotype associated with type 2 diabetes (T2D) for which the molecular mediators remain unclear. We therefore conducted an expression analysis of human muscle biopsies from patients with T2D; normoglycemic but insulin-resistant subjects with a parental family history (FH(+)) of T2D; and family history-negative control individuals (FH(–)). Actin cytoskeleton genes regulated by serum response factor (SRF) and its coactivator megakaryoblastic leukemia 1 (MKL1) had increased expression in T2D and FH(+) groups. Furthermore, striated muscle activator of Rho signaling (STARS), an activator of SRF, was upregulated in T2D and FH(+) and was inversely correlated with insulin sensitivity. Skeletal muscle from insulin-resistant mice recapitulated this gene expression pattern and showed reduced G-actin and increased nuclear localization of MKL1, each of which regulates SRF activity. Overexpression of MKL1 or reduction in G-actin decreased insulin-stimulated Akt phosphorylation, whereas reduction of STARS expression increased insulin signaling and glucose uptake. Pharmacological SRF inhibition by CCG-1423 reduced nuclear MKL1 and improved glucose uptake and tolerance in insulin-resistant mice in vivo. Thus, SRF pathway alterations are linked to insulin resistance, may contribute to T2D pathogenesis, and could represent therapeutic targets.

High conservation of transcription factor binding and evidence for combinatorial regulation across six Drosophila species

The binding of some transcription factors has been shown to diverge substantially between closely related species. Here we show that the binding of the developmental transcription factor Twist is highly conserved across six Drosophila species, revealing strong functional constraints at its enhancers. Conserved binding correlates with sequence motifs for Twist and its partners, permitting the de novo discovery of their combinatorial binding. It also includes over 10,000 low-occupancy sites near the detection limit, which tend to mark enhancers of later developmental stages. These results suggest that developmental enhancers can be highly evolutionarily constrained, presumably because of their complex combinatorial nature.

Prediction of regulatory interactions from genome sequences using a biophysical model for the Arabidopsis LEAFY transcription factor

Despite great advances in sequencing technologies, generating functional information for nonmodel organisms remains a challenge. One solution lies in an improved ability to predict genetic circuits based on primary DNA sequence in combination with detailed knowledge of regulatory proteins that have been characterized in model species. Here, we focus on the LEAFY (LFY) transcription factor, a conserved master regulator of floral development. Starting with biochemical and structural information, we built a biophysical model describing LFY DNA binding specificity in vitro that accurately predicts in vivo LFY binding sites in the Arabidopsis thaliana genome. Applying the model to other plant species, we could follow the evolution of the regulatory relationship between LFY and the AGAMOUS (AG) subfamily of MADS box genes and show that this link predates the divergence between monocots and eudicots. Remarkably, our model succeeds in detecting the connection between LFY and AG homologs despite extensive variation in binding sites. This demonstrates that the cis-element fluidity recently observed in animals also exists in plants, but the challenges it poses can be overcome with predictions grounded in a biophysical model. Therefore, our work opens new avenues to deduce the structure of regulatory networks from mere inspection of genomic sequences.

Adipogenic histone mark regulation by matrix metalloproteinase 14 in collagen-rich microenvironments

Adipogenesis is directed by both transcriptional network and posttranslational modification of chromatin structure. Although adipogenesis in vivo proceeds in collagen-rich extracellular matrix (ECM) environments, the impact of ECM proteins and their modifying enzymes on the epigenetic regulation of adipogenesis has been largely unknown. We aimed to define the role of fibrillar type I collagen and its modifying enzymes in regulating adipogenic chromatin signatures and gene regulation in the in vivo-like settings. Adipogenic cocktail induces a robust increase in the level of protranscriptional acetylated histone H3 at lysine 9 (H3K9ac) within 24 h. When cultured atop fibrillar type I collagen gel, however, H3K9ac levels in differentiating 3T3-L1 cells are substantially reduced. The suppression of adipogenic histone mark in differentiating 3T3-L1 cells is type I collagen density dependent and released by heat denaturing of the subjacent collagen substratum, pointing to the critical role played by the triple-helical structure of type I collagen. By probing adipogenic collagenolysis with a series of proteinase inhibitors, matrix metalloproteinase (MMP) family members are found to be responsible for adipogenic collagenolysis. At the same time, MMP inhibitor specifically blocked the adipogenic induction of H3K9ac. By targeting individual MMP using small interfering RNA oligos, MMP14 was identified as the major adipogenic MMP critical for H3K9 acetylation. Consistently, MMP14-null adipose tissues display diminished protranscriptional histone mark H3K9ac while maintaining repressive histone mark tri-methylated histone H3 at lysine 9 (H3K9me3). Taken together, MMP14-dependent collagenolysis plays the major role in regulating adipogenic histone marks by releasing the epigenetic constraints imposed by fibrillar type I collagen.

Extensive chromatin remodelling and establishment of transcription factor 'hotspots' during early adipogenesis

Adipogenesis is tightly controlled by a complex network of transcription factors acting at different stages of differentiation. Peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer-binding protein (C/EBP) family members are key regulators of this process. We have employed DNase I hypersensitive site analysis to investigate the genome-wide changes in chromatin structure that accompany the binding of adipogenic transcription factors. These analyses revealed a dramatic and dynamic modulation of the chromatin landscape during the first hours of adipocyte differentiation that coincides with cooperative binding of multiple early transcription factors (including glucocorticoid receptor, retinoid X receptor, Stat5a, C/EBPβ and -δ) to transcription factor 'hotspots'. Our results demonstrate that C/EBPβ marks a large number of these transcription factor 'hotspots' before induction of differentiation and chromatin remodelling and is required for their establishment. Furthermore, a subset of early remodelled C/EBP-binding sites persists throughout differentiation and is later occupied by PPARγ, indicating that early C/EBP family members, in addition to their well-established role in activation of PPARγ transcription, may act as pioneering factors for PPARγ binding.

Cross species comparison of C/EBPα and PPARγ profiles in mouse and human adipocytes reveals interdependent retention of binding sites

Background

The transcription factors peroxisome proliferator activated receptor γ (PPARγ) and CCAAT/enhancer binding protein α (C/EBPα) are key transcriptional regulators of adipocyte differentiation and function. We and others have previously shown that binding sites of these two transcription factors show a high degree of overlap and are associated with the majority of genes upregulated during differentiation of murine 3T3-L1 adipocytes.

Results

Here we have mapped all binding sites of C/EBPα and PPARγ in human SGBS adipocytes and compared these with the genome-wide profiles from mouse adipocytes to systematically investigate what biological features correlate with retention of sites in orthologous regions between mouse and human. Despite a limited interspecies retention of binding sites, several biological features make sites more likely to be retained. First, co-binding of PPARγ and C/EBPα in mouse is the most powerful predictor of retention of the corresponding binding sites in human. Second, vicinity to genes highly upregulated during adipogenesis significantly increases retention. Third, the presence of C/EBPα consensus sites correlate with retention of both factors, indicating that C/EBPα facilitates recruitment of PPARγ. Fourth, retention correlates with overall sequence conservation within the binding regions independent of C/EBPα and PPARγ sequence patterns, indicating that other transcription factors work cooperatively with these two key transcription factors.

Conclusions

This study provides a comprehensive and systematic analysis of what biological features impact on retention of binding sites between human and mouse. Specifically, we show that the binding of C/EBPα and PPARγ in adipocytes have evolved in a highly interdependent manner, indicating a significant cooperativity between these two transcription factors.

Species-specific strategies underlying conserved functions of metabolic transcription factors

The winged helix protein FOXA2 and the nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ) are highly conserved, regionally expressed transcription factors (TFs) that regulate networks of genes controlling complex metabolic functions. Cistrome analysis for Foxa2 in mouse liver and PPARγ in mouse adipocytes has previously produced consensus-binding sites that are nearly identical to those used by the corresponding TFs in human cells. We report here that, despite the conservation of the canonical binding motif, the great majority of binding regions for FOXA2 in human liver and for PPARγ in human adipocytes are not in the orthologous locations corresponding to the mouse genome, and vice versa. Of note, TF binding can be absent in one species despite sequence conservation, including motifs that do support binding in the other species, demonstrating a major limitation of in silico binding site prediction. Whereas only approximately 10% of binding sites are conserved, gene-centric analysis reveals that about 50% of genes with nearby TF occupancy are shared across species for both hepatic FOXA2 and adipocyte PPARγ. Remarkably, for both TFs, many of the shared genes function in tissue-specific metabolic pathways, whereas species-unique genes fail to show enrichment for these pathways. Nonetheless, the species-unique genes, like the shared genes, showed the expected transcriptional regulation by the TFs in loss-of-function experiments. Thus, species-specific strategies underlie the biological functions of metabolic TFs that are highly conserved across mammalian species. Analysis of factor binding in multiple species may be necessary to distinguish apparent species-unique noise and reveal functionally relevant information.

miR-130 suppresses adipogenesis by inhibiting peroxisome proliferator-activated receptor gamma expression

Adipose tissue development is tightly regulated by altering gene expression. MicroRNAs are strong posttranscriptional regulators of mammalian differentiation. We hypothesized that microRNAs might influence human adipogenesis by targeting specific adipogenic factors. We identified microRNAs that showed varying abundance during the differentiation of human preadipocytes into adipocytes. Among them, miR-130 strongly affected adipocyte differentiation, as overexpressing miR-130 impaired adipogenesis and reducing miR-130 enhanced adipogenesis. A key effector of miR-130 actions was the protein peroxisome proliferator-activated receptor γ (PPARγ), a major regulator of adipogenesis. Interestingly, miR-130 potently repressed PPARγ expression by targeting both the PPARγ mRNA coding and 3' untranslated regions. Adipose tissue from obese women contained significantly lower miR-130 and higher PPARγ mRNA levels than that from nonobese women. Our findings reveal that miR-130 reduces adipogenesis by repressing PPARγ biosynthesis and suggest that perturbations in this regulation is linked to human obesity.

Molecular basis for gene-specific transactivation by nuclear receptors

Nuclear receptors (NRs) are key transcriptional regulators of metazoan physiology and metabolism. Different NRs bind to similar or even identical core response elements; however, they regulate transcription in a highly receptor- and gene-specific manner. These differences in gene activation can most likely be accounted for by mechanisms involving receptor-specific interactions with DNA as well as receptor-specific interactions with protein complexes binding to adjacent and distant DNA sequences. Here, we review key molecular aspects of transactivation by NRs with special emphasis on the recent advances in the molecular mechanisms responsible for receptor- and gene-specific transcriptional activation. This article is part of a Special Issue entitled: Translating nuclear receptors from health to disease.

Tiling histone H3 lysine 4 and 27 methylation in zebrafish using high-density microarrays

Background

Uncovering epigenetic states by chromatin immunoprecipitation and microarray hybridization (ChIP-chip) has significantly contributed to the understanding of gene regulation at the genome-scale level. Many studies have been carried out in mice and humans; however limited high-resolution information exists to date for non-mammalian vertebrate species.

Principal Findings

We report a 2.1-million feature high-resolution Nimblegen tiling microarray for ChIP-chip interrogations of epigenetic states in zebrafish (Danio rerio). The array covers 251 megabases of the genome at 92 base-pair resolution. It includes ∼15 kb of upstream regulatory sequences encompassing all RefSeq promoters, and over 5 kb in the 5′ end of coding regions. We identify with high reproducibility, in a fibroblast cell line, promoters enriched in H3K4me3, H3K27me3 or co-enriched in both modifications. ChIP-qPCR and sequential ChIP experiments validate the ChIP-chip data and support the co-enrichment of trimethylated H3K4 and H3K27 on a subset of genes. H3K4me3- and/or H3K27me3-enriched genes are associated with distinct transcriptional status and are linked to distinct functional categories.

Conclusions

We have designed and validated for the scientific community a comprehensive high-resolution tiling microarray for investigations of epigenetic states in zebrafish, a widely used developmental and disease model organism.

Genome-wide analysis of glucocorticoid receptor binding regions in adipocytes reveal gene network involved in triglyceride homeostasis

Glucocorticoids play important roles in the regulation of distinct aspects of adipocyte biology. Excess glucocorticoids in adipocytes are associated with metabolic disorders, including central obesity, insulin resistance and dyslipidemia. To understand the mechanisms underlying the glucocorticoid action in adipocytes, we used chromatin immunoprecipitation sequencing to isolate genome-wide glucocorticoid receptor (GR) binding regions (GBRs) in 3T3-L1 adipocytes. Furthermore, gene expression analyses were used to identify genes that were regulated by glucocorticoids. Overall, 274 glucocorticoid-regulated genes contain or locate nearby GBR. We found that many GBRs were located in or nearby genes involved in triglyceride (TG) synthesis (Scd-1, 2, 3, GPAT3, GPAT4, Agpat2, Lpin1), lipolysis (Lipe, Mgll), lipid transport (Cd36, Lrp-1, Vldlr, Slc27a2) and storage (S3-12). Gene expression analysis showed that except for Scd-3, the other 13 genes were induced in mouse inguinal fat upon 4-day glucocorticoid treatment. Reporter gene assays showed that except Agpat2, the other 12 glucocorticoid-regulated genes contain at least one GBR that can mediate hormone response. In agreement with the fact that glucocorticoids activated genes in both TG biosynthetic and lipolytic pathways, we confirmed that 4-day glucocorticoid treatment increased TG synthesis and lipolysis concomitantly in inguinal fat. Notably, we found that 9 of these 12 genes were induced in transgenic mice that have constant elevated plasma glucocorticoid levels. These results suggested that a similar mechanism was used to regulate TG homeostasis during chronic glucocorticoid treatment. In summary, our studies have identified molecular components in a glucocorticoid-controlled gene network involved in the regulation of TG homeostasis in adipocytes. Understanding the regulation of this gene network should provide important insight for future therapeutic developments for metabolic diseases.

Genome-wide interplay of nuclear receptors with the epigenome

The nuclear receptor superfamily consists of DNA binding transcription factors that are involved in regulating a wide variety of processes such as metabolism, development, reproduction, and immune responses. Upon binding, nuclear receptors modulate transcription through affecting the local chromatin environment via recruitment of various coregulatory proteins. The recent development of new high-throughput sequencing methods allowed for the first time the comprehensive examination of nuclear receptor action in the context of the epigenome. Here, we discuss how recent genome-wide analyses have provided important new insights on the interplay of nuclear receptors and the epigenome in health and disease. This article is part of a Special Issue entitled: Translating nuclear receptors from health to disease.