The PPARgamma locus makes long-range chromatin interactions with selected tissue-specific gene loci during adipocyte differentiation in a protein kinase A dependent manner

Differentiation signaling results in reprogramming of cellular gene expression that leads to morphological changes and functional specialization of a precursor cell. This global change in gene expression involves temporal regulation of differentiation-specific genes that are located throughout the genome, raising the idea that genome structure may also be re-organized during cell differentiation to facilitate regulated gene expression. Using in vitro adipocyte differentiation as a model, we explored whether gene organization within the nucleus is altered upon exposure of precursor cells to signaling molecules that induce adipogenesis. The peroxisome proliferator-activated receptor gamma (PPARgamma) nuclear hormone receptor is a master determinant of adipogenesis and is required for adipose differentiation. We utilized the chromosome conformation capture (3C) assay to determine whether the position of the PPARgamma locus relative to other adipogenic genes is changed during differentiation. We report that the PPARgamma2 promoter is transiently positioned in proximity to the promoters of genes encoding adipokines and lipid droplet associated proteins at 6 hours post-differentiation, a time that precedes expression of any of these genes. In contrast, the PPARgamma2 promoter was not in proximity to the EF1alpha promoter, which drives expression of a constitutively active, housekeeping gene that encodes a translation elongation factor, nor was the PPARgamma2 promoter in proximity to the promoter driving the expression of the C/EBPalpha regulatory protein. The formation of the long-range, intergenic interactions involving the PPARgamma2 promoter required the regulatory factor C/EBPbeta, elevated cyclic AMP (cAMP) levels, and protein kinase A (PKA) signaling. We conclude that genome organization is dynamically remodeled in response to adipogenic signaling, and we speculate that these transient inter-genic interactions may be formed for the purposes of selecting some of the transcriptionally silent tissue-specific loci for subsequent transcriptional activation

A TAD boundary is preserved upon deletion of the CTCF-rich Firre locus

The binding of the transcriptional regulator CTCF to the genome has been implicated in the formation of topologically associated domains (TADs). However, the general mechanisms of folding the genome into TADs are not fully understood. Here we test the effects of deleting a CTCF-rich locus on TAD boundary formation. Using genome-wide chromosome conformation capture (Hi-C), we focus on one TAD boundary on chromosome X harboring ~ 15 CTCF binding sites and located at the long non-coding RNA (lncRNA) locus Firre. Specifically, this TAD boundary is invariant across evolution, tissues, and temporal dynamics of X-chromosome inactivation. We demonstrate that neither the deletion of this locus nor the ectopic insertion of Firre cDNA or its ectopic expression are sufficient to alter TADs in a sex-specific or allele-specific manner. In contrast, Firre’s deletion disrupts the chromatin super-loop formation of the inactive X-chromosome. Collectively, our findings suggest that apart from CTCF binding, additional mechanisms may play roles in establishing TAD boundary formation

Epigenetic control of cell cycle-dependent histone gene expression is a principal component of the abbreviated pluripotent cell cycle

Self-renewal of human pluripotent embryonic stem cells proceeds via an abbreviated cell cycle with a shortened G(1) phase. We examined which genes are modulated in this abbreviated period and the epigenetic mechanisms that control their expression. Accelerated upregulation of genes encoding histone proteins that support DNA replication is the most prominent gene regulatory program at the G(1)/S-phase transition in pluripotent cells. Expedited expression of histone genes is mediated by a unique chromatin architecture reflected by major nuclease hypersensitive sites, atypical distribution of epigenetic histone marks, and a region devoid of histone octamers. We observed remarkable differences in chromatin structure--hypersensitivity and histone protein modifications--between human embryonic stem (hES) and normal diploid cells. Cell cycle-dependent transcription factor binding permits dynamic three-dimensional interactions between transcript initiating and processing factors at 5\u27 and 3\u27 regions of the gene. Thus, progression through the abbreviated G(1) phase involves cell cycle stage-specific chromatin-remodeling events and rapid assembly of subnuclear microenvironments that activate histone gene transcription to promote nucleosomal packaging of newly replicated DNA during stem cell renewal

SMARCA4 regulates gene expression and higher-order chromatin structure in proliferating mammary epithelial cells

The packaging of DNA into chromatin plays an important role in transcriptional regulation and nuclear processes. Brahma-related gene-1 SMARCA4 (also known as BRG1), the essential ATPase subunit of the mammalian SWI/SNF chromatin remodeling complex, uses the energy from ATP hydrolysis to disrupt nucleosomes at target regions. Although the transcriptional role of SMARCA4 at gene promoters is well-studied, less is known about its role in higher-order genome organization. SMARCA4 knockdown in human mammary epithelial MCF-10A cells resulted in 176 up-regulated genes, including many related to lipid and calcium metabolism, and 1292 down-regulated genes, some of which encode extracellular matrix (ECM) components that can exert mechanical forces and affect nuclear structure. ChIP-seq analysis of SMARCA4 localization and SMARCA4-bound super-enhancers demonstrated extensive binding at intergenic regions. Furthermore, Hi-C analysis showed extensive SMARCA4-mediated alterations in higher-order genome organization at multiple resolutions. First, SMARCA4 knockdown resulted in clustering of intra- and inter-subtelomeric regions, demonstrating a novel role for SMARCA4 in telomere organization. SMARCA4 binding was enriched at topologically associating domain (TAD) boundaries, and SMARCA4 knockdown resulted in weakening of TAD boundary strength. Taken together, these findings provide a dynamic view of SMARCA4-dependent changes in higher-order chromatin organization and gene expression, identifying SMARCA4 as a novel component of chromatin organization

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

Three-dimensional organization of chromatin is fundamental for transcriptional regulation. Tissue-specific transcriptional programs are orchestrated by transcription factors and epigenetic regulators. The RUNX2 transcription factor is required for differentiation of precursor cells into mature osteoblasts. Although organization and control of the bone-specific Runx2-P1 promoter have been studied extensively, long-range regulation has not been explored. In this study, we investigated higher-order organization of the Runx2-P1 promoter during osteoblast differentiation. Mining the ENCODE database revealed interactions between Runx2-P1 and Supt3h promoters in several non-mesenchymal human cell lines. Supt3h is a ubiquitously expressed gene located within the first intron of Runx2. These two genes show shared synteny across species from humans to sponges. Chromosome conformation capture analysis in the murine pre-osteoblastic MC3T3-E1 cell line revealed increased contact frequency between Runx2-P1 and Supt3h promoters during differentiation. This increase was accompanied by enhanced DNaseI hypersensitivity along with RUNX2 and CTCF binding at the Supt3h promoter. Furthermore, interplasmid-3C and luciferase reporter assays showed that the Supt3h promoter can modulate Runx2-P1 activity via direct association. Taken together, our data demonstrate physical proximity between Runx2-P1 and Supt3h promoters, consistent with their syntenic nature. Importantly, we identify the Supt3h promoter as a potential regulator of the bone-specific Runx2-P1 promoter. Acids Research

Enhancers in the Peril lincRNA locus regulate distant but not local genes

Abstract Background Recently, it has become clear that some promoters function as long-range regulators of gene expression. However, direct and quantitative assessment of enhancer activity at long intergenic noncoding RNA (lincRNA) or mRNA gene bodies has not been performed. To unbiasedly assess the enhancer capacity across lincRNA and mRNA loci, we performed a massively parallel reporter assay (MPRA) on six lincRNA loci and their closest protein-coding neighbors. Results For both gene classes, we find significantly more MPRA activity in promoter regions than in gene bodies. However, three lincRNA loci, Lincp21, LincEnc1, and Peril, and one mRNA locus, Morc2a, display significant enhancer activity within their gene bodies. We hypothesize that such peaks may mark long-range enhancers, and test this in vivo using RNA sequencing from a knockout mouse model and high-throughput chromosome conformation capture (Hi-C). We find that ablation of a high-activity MPRA peak in the Peril gene body leads to consistent dysregulation of Mccc1 and Exosc9 in the neighboring topologically associated domain (TAD). This occurs irrespective of Peril lincRNA expression, demonstrating this regulation is DNA-dependent. Hi-C confirms long-range contacts with the neighboring TAD, and these interactions are altered upon Peril knockout. Surprisingly, we do not observe consistent regulation of genes within the local TAD. Together, these data suggest a long-range enhancer-like function for the Peril gene body. Conclusions A multi-faceted approach combining high-throughput enhancer discovery with genetic models can connect enhancers to their gene targets and provides evidence of inter-TAD gene regulation

Intranuclear and higher-order chromatin organization of the major histone gene cluster in breast cancer

Alterations in nuclear morphology are common in cancer progression. However, the degree to which gross morphological abnormalities translate into compromised higher-order chromatin organization is poorly understood. To explore the functional links between gene expression and chromatin structure in breast cancer, we performed RNA-seq gene expression analysis on the basal breast cancer progression model based on human MCF10A cells. Positional gene enrichment identified the major histone gene cluster at chromosome 6p22 as one of the most significantly upregulated (and not amplified) clusters of genes from the normal-like MCF10A to premalignant MCF10AT1 and metastatic MCF10CA1a cells. This cluster is subdivided into three sub-clusters of histone genes that are organized into hierarchical topologically associating domains (TADs). Interestingly, the sub-clusters of histone genes are located at TAD boundaries and interact more frequently with each other than the regions in-between them, suggesting that the histone sub-clusters form an active chromatin hub. The anchor sites of loops within this hub are occupied by CTCF, a known chromatin organizer. These histone genes are transcribed and processed at a specific sub-nuclear microenvironment termed the major histone locus body (HLB). While the overall chromatin structure of the major HLB is maintained across breast cancer progression, we detected alterations in its structure that may relate to gene expression. Importantly, breast tumor specimens also exhibit a coordinate pattern of upregulation across the major histone gene cluster. Our results provide a novel insight into the connection between the higher-order chromatin organization of the major HLB and its regulation during breast cancer progression

The PPARγ2 promoter maximally interacts with adipokine and perilipin gene promoters in 3T3-L1 cells 6 hours post-differentiation.

<p>(A) Schematic map of chromosomal locations of loci of interest. Lep: leptin; Plin1: perilipin 1; Plin2: perilipin 2; Eef1α1: EF1α; AdipoQ; adiponectin. 3C assays were performed at the indicated times to measure long-range, intergenic interactions between the PPARγ2 promoter and (B) the adiponectin promoter, (C) the leptin promoter, (D) the perilipin 1 promoter, (E) the perilipin 2/ADRP promoter, (F) the C/EBPα promoter, or (G) the EF1α promoter. Interaction frequencies are presented relative to the intragenic interaction frequency observed between different regions of the TFIIH locus, which served as a normalization control. The interaction frequency at 0 h was set to 1; the interaction frequencies for the remaining time points are presented relative to that value. Results indicate the average of three independent experiments +/− standard deviation. Values for the timepoint of maximum interaction frequency (6 h) were compared to values at time 0 by one-tailed t-test. ***<i>p</i><0.001.</p

Elevation of cAMP levels is required for the formation of long-range, inter-genic interactions between the PPARγ2 promoter and the indicated promoters.

<p>The EF1α promoter is shown as a control. 3C assays were performed on 3T3-L1 cells at 0 or at 6 h post-differentiation in the presence of complete differentiation cocktail or in cocktail lacking the indicated component. Interaction frequencies were quantified as described in the legends for <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0086140#pone-0086140-g001" target="_blank">Figures 1</a>–<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0086140#pone-0086140-g002" target="_blank">2</a>. Results indicate the average of three independent experiments +/− standard deviation. Values for the 6 h samples differentiated in complete cocktail were compared to values for the samples differentiated in cocktail lacking IBMX by one-tailed t-test. <i>p</i> values are indicated on the graph. ins; insulin. dex; dexamethasone.</p

PKA inhibitors prevent the formation of long-range, inter-genic interactions between the PPARγ2 promoter and (A) the adiponectin promoter, (B) the leptin promoter, (C) the perilipin 1 promoter, (D) the perilipin 2/ADRP promoter.

<p>The EF1α promoter (E) is shown as a control. Interaction frequencies were quantified as described in the legends for <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0086140#pone-0086140-g001" target="_blank">Figures 1</a>–<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0086140#pone-0086140-g002" target="_blank">2</a>. Results indicate the average of three independent experiments +/− standard deviation. Values for the 6 h samples differentiated in complete cocktail were compared to values for the samples differentiated for 6 h in complete cocktail containing the PKI or H89 inhibitors by one-tailed t-test. ***<i>p</i><0.001.</p