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

BACKGROUND: The AP-1 transcription factor, FBJ osteosarcoma oncogene (FOS), is induced in adult muscle satellite cells (SCs) within hours following muscle damage and is required for effective stem cell activation and muscle repair. However, why FOS is rapidly downregulated before SCs enter cell cycle as progenitor cells (i.e., transiently expressed) remains unclear. Further, whether boosting FOS levels in the proliferating progeny of SCs can enhance their myogenic properties needs further evaluation. METHODS: We established an inducible, FOS expression system to evaluate the impact of persistent FOS activity in muscle progenitor cells ex vivo. We performed various assays to measure cellular proliferation and differentiation, as well as uncover changes in RNA levels and three-dimensional (3D) chromatin interactions. RESULTS: Persistent FOS activity in primary muscle progenitor cells severely antagonizes their ability to differentiate and form myotubes within the first 2 weeks in culture. RNA-seq analysis revealed that ectopic FOS activity in muscle progenitor cells suppressed a global pro-myogenic transcriptional program, while activating a stress-induced, mitogen-activated protein kinase (MAPK) transcriptional signature. Additionally, we observed various FOS-dependent, chromosomal re-organization events in A/B compartments, topologically associated domains (TADs), and genomic loops near FOS-regulated genes. CONCLUSIONS: Our results suggest that elevated FOS activity in recently activated muscle progenitor cells perturbs cellular differentiation by altering the 3D chromosome organization near critical pro-myogenic genes. This work highlights the crucial importance of tightly controlling FOS expression in the muscle lineage and suggests that in states of chronic stress or disease, persistent FOS activity in muscle precursor cells may disrupt the muscle-forming process.

Reorganization of inter-chromosomal interactions in the 2q37-deletion syndrome.

Chromosomes occupy distinct interphase territories in the three-dimensional nucleus. However, how these chromosome territories are arranged relative to one another is poorly understood. Here, we investigated the inter-chromosomal interactions between chromosomes 2q, 12, and 17 in human mesenchymal stem cells (MSCs) and MSC-derived cell types by DNA-FISH We compared our findings in normal karyotypes with a three-generation family harboring a 2q37-deletion syndrome, featuring a heterozygous partial deletion of histone deacetylase 4 (HDAC4) on chr2q37. In normal karyotypes, we detected stable, recurring arrangements and interactions between the three chromosomal territories with a tissue-specific interaction bias at certain loci. These inter-chromosomal interactions were confirmed by Hi-C. Interestingly, the disease-related HDAC4 deletion resulted in displaced inter-chromosomal arrangements and altered interactions between the deletion-affected chromosome 2 and chromosome 12 and/or 17 in 2q37-deletion syndrome patients. Our findings provide evidence for a direct link between a structural chromosomal aberration and altered interphase architecture that results in a nuclear configuration, supporting a possible molecular pathogenesis.

The connection between BRG1, CTCF and topoisomerases at TAD boundaries.

The eukaryotic genome is partitioned into topologically associating domains (TADs). Despite recent advances characterizing TADs and TAD boundaries, the organization of these structures is an important dimension of genome architecture and function that is not well understood. Recently, we demonstrated that knockdown of BRG1, an ATPase driving the chromatin remodeling activity of mammalian SWI/SNF enzymes, globally alters long-range genomic interactions and results in a reduction of TAD boundary strength. We provided evidence suggesting that this effect may be due to BRG1 affecting nucleosome occupancy around CTCF sites present at TAD boundaries. In this review, we elaborate on our findings and speculate that BRG1 may contribute to the regulation of the structural and functional properties of chromatin at TAD boundaries by affecting the function or the recruitment of CTCF and DNA topoisomerase complexes.

Identifying Nuclear Matrix-Attached DNA Across the Genome.

Experimental approaches to define the relationship between gene expression and nuclear matrix attachment regions (MARs) have given contrasting and method-specific results. We have developed a next generation sequencing strategy to identify MARs across the human genome (MAR-Seq). The method is based on crosslinking chromatin to its nuclear matrix attachment sites to minimize changes during biochemical processing. We used this method to compare nuclear matrix organization in MCF-10A mammary epithelial-like cells and MDA-MB-231 breast cancer cells and evaluated the results in the context of global gene expression (array analysis) and positional enrichment of gene-regulatory histone modifications (ChIP-Seq). In the normal-like cells, nuclear matrix-attached DNA was enriched in expressed genes, while in the breast cancer cells, it was enriched in non-expressed genes. In both cell lines, the chromatin modifications that mark transcriptional activation or repression were appropriately associated with gene expression. Using this new MAR-Seq approach, we provide the first genome-wide characterization of nuclear matrix attachment in mammalian cells and reveal that the nuclear matrix-associated genome is highly cell-context dependent. J. Cell. Physiol. 232: 1295-1305, 2017. © 2016 Wiley Periodicals, Inc.

In Vivo Characterization of Linc-p21 Reveals Functional cis-Regulatory DNA Elements.

The Linc-p21 locus, encoding a long non-coding RNA, plays an important role in p53 signaling, cell-cycle regulation, and tumor suppression. However, despite extensive study, confusion exists regarding its mechanism of action: is activity driven by the transcript acting in trans, in cis, or by an underlying functional enhancer? Here, using a knockout mouse model and a massively parallel enhancer assay, we delineate the functional elements at this locus. We observe that, even in tissues with no detectable Linc-p21 transcript, deletion of the locus significantly affects local gene expression, including of the cell-cycle regulator Cdkn1a. To characterize this RNA-independent regulatory effect, we systematically interrogated the underlying DNA sequence for enhancer activity at nucleotide resolution and confirmed the existence of multiple enhancer elements. Together, these data suggest that, in vivo, the cis-regulatory effects mediated by Linc-p21, in the presence or absence of transcription, are due to DNA enhancer elements.

RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells.

RUNX1 is a transcription factor functioning both as an oncogene and a tumor suppressor in breast cancer. RUNX1 alters chromatin structure in cooperation with chromatin modifier and remodeling enzymes. In this study, we examined the relationship between RUNX1-mediated transcription and genome organization. We characterized genome-wide RUNX1 localization and performed RNA-seq and Hi-C in RUNX1-depleted and control MCF-7 breast cancer cells. RNA-seq analysis showed that RUNX1 depletion led to up-regulation of genes associated with chromatin structure and down-regulation of genes related to extracellular matrix biology, as well as NEAT1 and MALAT1 lncRNAs. Our ChIP-Seq analysis supports a prominent role for RUNX1 in transcriptional activation. About 30% of all RUNX1 binding sites were intergenic, indicating diverse roles in promoter and enhancer regulation and suggesting additional functions for RUNX1. Hi-C analysis of RUNX1-depleted cells demonstrated that overall three-dimensional genome organization is largely intact, but indicated enhanced association of RUNX1 near Topologically Associating Domain (TAD) boundaries and alterations in long-range interactions. These results suggest an architectural role for RUNX1 in fine-tuning local interactions rather than in global organization. Our results provide novel insight into RUNX1-mediated perturbations of higher-order genome organization that are functionally linked with RUNX1-dependent compromised gene expression in breast cancer cells.

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.

Chromatin interaction analysis reveals changes in small chromosome and telomere clustering between epithelial and breast cancer cells.

BACKGROUND: Higher-order chromatin structure is often perturbed in cancer and other pathological states. Although several genetic and epigenetic differences have been charted between normal and breast cancer tissues, changes in higher-order chromatin organization during tumorigenesis have not been fully explored. To probe the differences in higher-order chromatin structure between mammary epithelial and breast cancer cells, we performed Hi-C analysis on MCF-10A mammary epithelial and MCF-7 breast cancer cell lines. RESULTS: Our studies reveal that the small, gene-rich chromosomes chr16 through chr22 in the MCF-7 breast cancer genome display decreased interaction frequency with each other compared to the inter-chromosomal interaction frequency in the MCF-10A epithelial cells. Interestingly, this finding is associated with a higher occurrence of open compartments on chr16-22 in MCF-7 cells. Pathway analysis of the MCF-7 up-regulated genes located in altered compartment regions on chr16-22 reveals pathways related to repression of WNT signaling. There are also differences in intra-chromosomal interactions between the cell lines; telomeric and sub-telomeric regions in the MCF-10A cells display more frequent interactions than are observed in the MCF-7 cells. CONCLUSIONS: We show evidence of an intricate relationship between chromosomal organization and gene expression between epithelial and breast cancer cells. Importantly, this work provides a genome-wide view of higher-order chromatin dynamics and a resource for studying higher-order chromatin interactions in two cell lines commonly used to study the progression of breast cancer.

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.

Epigenetic landscape during osteoblastogenesis defines a differentiation-dependent Runx2 promoter region.

Runx2 is a developmentally regulated gene in vertebrates and is essential for bone formation and skeletal homeostasis. The induction of runx2-P1 isoform transcripts is a hallmark of early osteoblastogenesis. Although previous in vitro studies have defined a minimal Runx2-P1 promoter sequence with well-characterized functional elements, several lines of evidence suggest that transcription of the Runx2-P1 isoform relies on elements that extend beyond the previously defined P1 promoter boundaries. In this study, we examined Runx2-P1 transcriptional regulation in a cellular in vivo context during early osteoblastogenesis of MC3T3-E1 cultures and BMSCs induced towards the bone lineage by multi-layered analysis of the Runx2-P1 gene promoter using the following methodologies: 1) sequence homology among several mammalian species, 2) DNaseI hypersensitivity coupled with massively parallel sequencing (DNase-seq), and 3) chromatin immunoprecipitation of activating histone modifications coupled with massively parallel sequencing (ChIP-seq). These epigenetic features have allowed the demarcation of boundaries that redefine the minimal Runx2-P1 promoter to include a 336-bp sequence that mediates responsiveness to osteoblast differentiation. We also find that an additional level of control is contributed by a regulatory region in the 5'-UTR of Runx2-P1.

The PPARγ 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 (PPARγ) 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 PPARγ locus relative to other adipogenic genes is changed during differentiation. We report that the PPARγ2 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 PPARγ2 promoter was not in proximity to the EF1α promoter, which drives expression of a constitutively active, housekeeping gene that encodes a translation elongation factor, nor was the PPARγ2 promoter in proximity to the promoter driving the expression of the C/EBPα regulatory protein. The formation of the long-range, intergenic interactions involving the PPARγ2 promoter required the regulatory factor C/EBPß, 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.