CIC-DUX4 Chromatin Profiling Reveals New Epigenetic Dependencies and Actionable Therapeutic Targets in CIC-Rearranged Sarcomas.

CIC-DUX4-rearranged sarcoma (CDS) is a rare and aggressive soft tissue tumor that occurs most frequently in young adults. The key oncogenic driver of this disease is the expression of the CIC-DUX4 fusion protein as a result of chromosomal rearrangements. CIC-DUX4 displays chromatin binding properties, and is therefore believed to function as an aberrant transcription factor. However, the chromatin remodeling events induced by CIC-DUX4 are not well understood, limiting our ability to identify new mechanism-based therapeutic strategies for these patients. Here, we generated a genome-wide profile of CIC-DUX4 DNA occupancy and associated chromatin states in human CDS cell models and primary tumors. Combining chromatin profiling, proximity ligation assays, as well as genetic and pharmacological perturbations, we show that CIC-DUX4 operates as a potent transcriptional activator at its binding sites. This property is in contrast with the repressive function of the wild-type CIC protein, and is mainly mediated through the direct interaction of CIC-DUX4 with the acetyltransferase p300. In keeping with this, we show p300 to be essential for CDS tumor cell proliferation; additionally, we find its pharmacological inhibition to significantly impact tumor growth in vitro and in vivo. Taken together, our study elucidates the mechanisms underpinning CIC-DUX4-mediated transcriptional regulation.

ATAC-clock: An aging clock based on chromatin accessibility.

The establishment of aging clocks highlighted the strong link between changes in DNA methylation and aging. Yet, it is not known if other epigenetic features could be used to predict age accurately. Furthermore, previous studies have observed a lack of effect of age-related changes in DNA methylation on gene expression, putting the interpretability of DNA methylation-based aging clocks into question. In this study, we explore the use of chromatin accessibility to construct aging clocks. We collected blood from 159 human donors and generated chromatin accessibility, transcriptomic, and cell composition data. We investigated how chromatin accessibility changes during aging and constructed a novel aging clock with a median absolute error of 5.27 years. The changes in chromatin accessibility used by the clock were strongly related to transcriptomic alterations, aiding clock interpretation. We additionally show that our chromatin accessibility clock performs significantly better than a transcriptomic clock trained on matched samples. In conclusion, we demonstrate that the clock relies on cell-intrinsic chromatin accessibility alterations rather than changes in cell composition. Further, we present a new approach to construct epigenetic aging clocks based on chromatin accessibility, which bear a direct link to age-related transcriptional alterations, but which allow for more accurate age predictions than transcriptomic clocks.

Contrasting effects of whole-body and hepatocyte-specific deletion of the RNA polymerase III repressor Maf1 in the mouse.

MAF1 is a nutrient-sensitive, TORC1-regulated repressor of RNA polymerase III (Pol III). MAF1 downregulation leads to increased lipogenesis in Drosophila melanogaster, Caenorhabditis elegans, and mice. However, Maf1 -/- mice are lean as increased lipogenesis is counterbalanced by futile pre-tRNA synthesis and degradation, resulting in increased energy expenditure. We compared Chow-fed Maf1 -/- mice with Chow- or High Fat (HF)-fed Maf1 hep-/- mice that lack MAF1 specifically in hepatocytes. Unlike Maf1 -/- mice, Maf1 hep-/- mice become heavier and fattier than control mice with old age and much earlier under a HF diet. Liver ChIPseq, RNAseq and proteomics analyses indicate increased Pol III occupancy at Pol III genes, very few differences in mRNA accumulation, and protein accumulation changes consistent with increased lipogenesis. Futile pre-tRNA synthesis and degradation in the liver, as likely occurs in Maf1 hep-/- mice, thus seems insufficient to counteract increased lipogenesis. Indeed, RNAseq and metabolite profiling indicate that liver phenotypes of Maf1 -/- mice are strongly influenced by systemic inter-organ communication. Among common changes in the three phenotypically distinct cohorts, Angiogenin downregulation is likely linked to increased Pol III occupancy of tRNA genes in the Angiogenin promoter.

EWSR1-ATF1 dependent 3D connectivity regulates oncogenic and differentiation programs in Clear Cell Sarcoma.

Oncogenic fusion proteins generated by chromosomal translocations play major roles in cancer. Among them, fusions between EWSR1 and transcription factors generate oncogenes with powerful chromatin regulatory activities, capable of establishing complex gene expression programs in permissive precursor cells. Here we define the epigenetic and 3D connectivity landscape of Clear Cell Sarcoma, an aggressive cancer driven by the EWSR1-ATF1 fusion gene. We find that EWSR1-ATF1 displays a distinct DNA binding pattern that requires the EWSR1 domain and promotes ATF1 retargeting to new distal sites, leading to chromatin activation and the establishment of a 3D network that controls oncogenic and differentiation signatures observed in primary CCS tumors. Conversely, EWSR1-ATF1 depletion results in a marked reconfiguration of 3D connectivity, including the emergence of regulatory circuits that promote neural crest-related developmental programs. Taken together, our study elucidates the epigenetic mechanisms utilized by EWSR1-ATF1 to establish regulatory networks in CCS, and points to precursor cells in the neural crest lineage as candidate cells of origin for these tumors.

TWIST1 expression is associated with high-risk neuroblastoma and promotes primary and metastatic tumor growth.

The embryonic transcription factors TWIST1/2 are frequently overexpressed in cancer, acting as multifunctional oncogenes. Here we investigate their role in neuroblastoma (NB), a heterogeneous childhood malignancy ranging from spontaneous regression to dismal outcomes despite multimodal therapy. We first reveal the association of TWIST1 expression with poor survival and metastasis in primary NB, while TWIST2 correlates with good prognosis. Secondly, suppression of TWIST1 by CRISPR/Cas9 results in a reduction of tumor growth and metastasis colonization in immunocompromised mice. Moreover, TWIST1 knockout tumors display a less aggressive cellular morphology and a reduced disruption of the extracellular matrix (ECM) reticulin network. Additionally, we identify a TWIST1-mediated transcriptional program associated with dismal outcome in NB and involved in the control of pathways mainly linked to the signaling, migration, adhesion, the organization of the ECM, and the tumor cells versus tumor stroma crosstalk. Taken together, our findings confirm TWIST1 as promising therapeutic target in NB.

MAF1 is a chronic repressor of RNA polymerase III transcription in the mouse.

Maf1-/- mice are lean, obesity-resistant and metabolically inefficient. Their increased energy expenditure is thought to be driven by a futile RNA cycle that reprograms metabolism to meet an increased demand for nucleotides stemming from the deregulation of RNA polymerase (pol) III transcription. Metabolic changes consistent with this model have been reported in both fasted and refed mice, however the impact of the fasting-refeeding-cycle on pol III function has not been examined. Here we show that changes in pol III occupancy in the liver of fasted versus refed wild-type mice are largely confined to low and intermediate occupancy genes; high occupancy genes are unchanged. However, in Maf1-/- mice, pol III occupancy of the vast majority of active loci in liver and the levels of specific precursor tRNAs in this tissue and other organs are higher than wild-type in both fasted and refed conditions. Thus, MAF1 functions as a chronic repressor of active pol III loci and can modulate transcription under different conditions. Our findings support the futile RNA cycle hypothesis, elaborate the mechanism of pol III repression by MAF1 and demonstrate a modest effect of MAF1 on global translation via reduced mRNA levels and translation efficiencies for several ribosomal proteins.

THAP11F80L cobalamin disorder-associated mutation reveals normal and pathogenic THAP11 functions in gene expression and cell proliferation.

Twelve human THAP proteins share the THAP domain, an evolutionary conserved zinc-finger DNA-binding domain. Studies of different THAP proteins have indicated roles in gene transcription, cell proliferation and development. We have analyzed this protein family, focusing on THAP7 and THAP11. We show that human THAP proteins possess differing homo- and heterodimer formation properties and interaction abilities with the transcriptional co-regulator HCF-1. HEK-293 cells lacking THAP7 were viable but proliferated more slowly. In contrast, HEK-293 cells were very sensitive to THAP11 alteration. Nevertheless, HEK-293 cells bearing a THAP11 mutation identified in a patient suffering from cobalamin disorder (THAP11F80L) were viable although proliferated more slowly. Cobalamin disorder is an inborn vitamin deficiency characterized by neurodevelopmental abnormalities, most often owing to biallelic mutations in the MMACHC gene, whose gene product MMACHC is a key enzyme in the cobalamin (vitamin B12) metabolic pathway. We show that THAP11F80L selectively affected promoter binding by THAP11, having more deleterious effects on a subset of THAP11 targets, and resulting in altered patterns of gene expression. In particular, THAP11F80L exhibited a strong effect on association with the MMACHC promoter and led to a decrease in MMACHC gene transcription, suggesting that the THAP11F80L mutation is directly responsible for the observed cobalamin disorder.

HCF-2 inhibits cell proliferation and activates differentiation-gene expression programs.

HCF-2 is a member of the host-cell-factor protein family, which arose in early vertebrate evolution as a result of gene duplication. Whereas its paralog, HCF-1, is known to act as a versatile chromatin-associated protein required for cell proliferation and differentiation, much less is known about HCF-2. Here, we show that HCF-2 is broadly present in human and mouse cells, and possesses activities distinct from HCF-1. Unlike HCF-1, which is excluded from nucleoli, HCF-2 is nucleolar-an activity conferred by one and a half C-terminal Fibronectin type 3 repeats and inhibited by the HCF-1 nuclear localization signal. Elevated HCF-2 synthesis in HEK-293 cells results in phenotypes reminiscent of HCF-1-depleted cells, including inhibition of cell proliferation and mitotic defects. Furthermore, increased HCF-2 levels in HEK-293 cells lead to inhibition of cell proliferation and metabolism gene-expression programs with parallel activation of differentiation and morphogenesis gene-expression programs. Thus, the HCF ancestor appears to have evolved into a small two-member protein family possessing contrasting nuclear versus nucleolar localization, and cell proliferation and differentiation functions.

Differential regulation of RNA polymerase III genes during liver regeneration.

Mouse liver regeneration after partial hepatectomy involves cells in the remaining tissue synchronously entering the cell division cycle. We have used this system and H3K4me3, Pol II and Pol III profiling to characterize adaptations in Pol III transcription. Our results broadly define a class of genes close to H3K4me3 and Pol II peaks, whose Pol III occupancy is high and stable, and another class, distant from Pol II peaks, whose Pol III occupancy strongly increases after partial hepatectomy. Pol III regulation in the liver thus entails both highly expressed housekeeping genes and genes whose expression can adapt to increased demand.

Rapid Recapitulation of Nonalcoholic Steatohepatitis upon Loss of Host Cell Factor 1 Function in Mouse Hepatocytes

Host-cell factor 1 (HCF-1), encoded by the ubiquitously expressed X-linked gene Hcfc1, is an epigenetic coregulator important for mouse development and cell proliferation, including during liver regeneration. We used a hepatocyte-specific inducible Hcfc1 knock-out allele (called Hcfc1hepKO), to induce HCF-1 loss in hepatocytes of hemizygous Hcfc1hepKO/Y males by four days. In heterozygous Hcfc1hepKO/+ females, owing to random X-chromosome inactivation, upon Hcfc1hepKO allele induction, a 50/50 mix of HCF-1 positive and negative hepatocyte clusters is engineered. The livers with Hcfc1hepKO/Y hepatocytes displayed a 21-24-day terminal non-alcoholic fatty liver (NAFL) followed by non-alcoholic steatohepatitis (NASH) disease progression typical of severe NAFL disease (NAFLD). In contrast, in livers with heterozygous Hcfc1hepKO/+ hepatocytes, HCF-1-positive hepatocytes replaced HCF-1-negative hepatocytes and revealed only mild-NAFL development. Loss of HCF-1 led to loss of PGC1α protein, probably owing to its destabilization, and deregulation of gene expression particularly of genes involved in mitochondrial structure and function, likely explaining the severe Hcfc1 hepKO/Y liver pathology. Thus, HCF-1 is essential for hepatocyte function, likely playing both transcriptional and non-transcriptional roles. These genetically-engineered loss-of-HCF-1 mice can be used to study NASH as well as NAFLD resolution.

Cycles of gene expression and genome response during mammalian tissue regeneration.

BACKGROUND: Compensatory liver hyperplasia-or regeneration-induced by two-thirds partial hepatectomy (PH) permits the study of synchronized activation of mammalian gene expression, particularly in relation to cell proliferation. Here, we measured genomic transcriptional responses and mRNA accumulation changes after PH and sham surgeries. RESULTS: During the first 10-20 h, the PH- and sham-surgery responses were very similar, including parallel early activation of cell-division-cycle genes. After 20 h, however, whereas post-PH livers continued with a robust and coordinate cell-division-cycle gene-expression response before returning to the resting state by 1 week, sham-surgery livers returned directly to a resting gene-expression state. Localization of RNA polymerase II (Pol II), and trimethylated histone H3 lysine 4 (H3K4me3) and 36 (H3K36me3) on genes dormant in the resting liver and activated during the PH response revealed a general de novo promoter Pol II recruitment and H3K4me3 increase during the early 10-20 h phase followed by Pol II elongation and H3K36me3 accumulation in gene bodies during the later proliferation phase. H3K36me3, generally appearing at the first internal exon, was preceded 5' by H3K36me2; 3' of the first internal exon, in about half of genes H3K36me3 predominated and in the other half H3K36me2 and H3K36me3 co-existed. Further, we observed some unusual gene profiles with abundant Pol II but little evident H3K4me3 or H3K36me3 modification, indicating that these modifications are neither universal nor essential partners to Pol II transcription. CONCLUSIONS: PH and sham surgical procedures on mice reveal striking early post-operatory gene expression similarities followed by synchronized mRNA accumulation and epigenetic histone mark changes specific to PH.

Mechanism of selective recruitment of RNA polymerases II and III to snRNA gene promoters.

RNA polymerase II (Pol II) small nuclear RNA (snRNA) promoters and type 3 Pol III promoters have highly similar structures; both contain an interchangeable enhancer and "proximal sequence element" (PSE), which recruits the SNAP complex (SNAPc). The main distinguishing feature is the presence, in the type 3 promoters only, of a TATA box, which determines Pol III specificity. To understand the mechanism by which the absence or presence of a TATA box results in specific Pol recruitment, we examined how SNAPc and general transcription factors required for Pol II or Pol III transcription of SNAPc-dependent genes (i.e., TATA-box-binding protein [TBP], TFIIB, and TFIIA for Pol II transcription and TBP and BRF2 for Pol III transcription) assemble to ensure specific Pol recruitment. TFIIB and BRF2 could each, in a mutually exclusive fashion, be recruited to SNAPc. In contrast, TBP-TFIIB and TBP-BRF2 complexes were not recruited unless a TATA box was present, which allowed selective and efficient recruitment of the TBP-BRF2 complex. Thus, TBP both prevented BRF2 recruitment to Pol II promoters and enhanced BRF2 recruitment to Pol III promoters. On Pol II promoters, TBP recruitment was separate from TFIIB recruitment and enhanced by TFIIA. Our results provide a model for specific Pol recruitment at SNAPc-dependent promoters.

Transcriptional interference by RNA polymerase III affects expression of the Polr3e gene.

Overlapping gene arrangements can potentially contribute to gene expression regulation. A mammalian interspersed repeat (MIR) nested in antisense orientation within the first intron of the Polr3e gene, encoding an RNA polymerase III (Pol III) subunit, is conserved in mammals and highly occupied by Pol III. Using a fluorescence assay, CRISPR/Cas9-mediated deletion of the MIR in mouse embryonic stem cells, and chromatin immunoprecipitation assays, we show that the MIR affects Polr3e expression through transcriptional interference. Our study reveals a mechanism by which a Pol II gene can be regulated at the transcription elongation level by transcription of an embedded antisense Pol III gene.

Diurnal regulation of RNA polymerase III transcription is under the control of both the feeding-fasting response and the circadian clock.

RNA polymerase III (Pol III) synthesizes short noncoding RNAs, many of which are essential for translation. Accordingly, Pol III activity is tightly regulated with cell growth and proliferation by factors such as MYC, RB1, TRP53, and MAF1. MAF1 is a repressor of Pol III transcription whose activity is controlled by phosphorylation; in particular, it is inactivated through phosphorylation by the TORC1 kinase complex, a sensor of nutrient availability. Pol III regulation is thus sensitive to environmental cues, yet a diurnal profile of Pol III transcription activity is so far lacking. Here, we first use gene expression arrays to measure mRNA accumulation during the diurnal cycle in the livers of (1) wild-type mice, (2) arrhythmic Arntl knockout mice, (3) mice fed at regular intervals during both night and day, and (4) mice lacking the Maf1 gene, and so provide a comprehensive view of the changes in cyclic mRNA accumulation occurring in these different systems. We then show that Pol III occupancy of its target genes rises before the onset of the night, stays high during the night, when mice normally ingest food and when translation is known to be increased, and decreases in daytime. Whereas higher Pol III occupancy during the night reflects a MAF1-dependent response to feeding, the rise of Pol III occupancy before the onset of the night reflects a circadian clock-dependent response. Thus, Pol III transcription during the diurnal cycle is regulated both in response to nutrients and by the circadian clock, which allows anticipatory Pol III transcription.

Human MAF1 targets and represses active RNA polymerase III genes by preventing recruitment rather than inducing long-term transcriptional arrest.

RNA polymerase III (Pol III) is tightly controlled in response to environmental cues, yet a genomic-scale picture of Pol III regulation and the role played by its repressor MAF1 is lacking. Here, we describe genome-wide studies in human fibroblasts that reveal a dynamic and gene-specific adaptation of Pol III recruitment to extracellular signals in an mTORC1-dependent manner. Repression of Pol III recruitment and transcription are tightly linked to MAF1, which selectively localizes at Pol III loci, even under serum-replete conditions, and increasingly targets transcribing Pol III in response to serum starvation. Combining Pol III binding profiles with EU-labeling and high-throughput sequencing of newly synthesized small RNAs, we show that Pol III occupancy closely reflects ongoing transcription. Our results exclude the long-term, unproductive arrest of Pol III on the DNA as a major regulatory mechanism and identify previously uncharacterized, differential coordination in Pol III binding and transcription under different growth conditions.

Loss of the RNA polymerase III repressor MAF1 confers obesity resistance.

MAF1 is a global repressor of RNA polymerase III transcription that regulates the expression of highly abundant noncoding RNAs in response to nutrient availability and cellular stress. Thus, MAF1 function is thought to be important for metabolic economy. Here we show that a whole-body knockout of Maf1 in mice confers resistance to diet-induced obesity and nonalcoholic fatty liver disease by reducing food intake and increasing metabolic inefficiency. Energy expenditure in Maf1(-/-) mice is increased by several mechanisms. Precursor tRNA synthesis was increased in multiple tissues without significant effects on mature tRNA levels, implying increased turnover in a futile tRNA cycle. Elevated futile cycling of hepatic lipids was also observed. Metabolite profiling of the liver and skeletal muscle revealed elevated levels of many amino acids and spermidine, which links the induction of autophagy in Maf1(-/-) mice with their extended life span. The increase in spermidine was accompanied by reduced levels of nicotinamide N-methyltransferase, which promotes polyamine synthesis, enables nicotinamide salvage to regenerate NAD(+), and is associated with obesity resistance. Consistent with this, NAD(+) levels were increased in muscle. The importance of MAF1 for metabolic economy reveals the potential for MAF1 modulators to protect against obesity and its harmful consequences.

Quantifying ChIP-seq data: a spiking method providing an internal reference for sample-to-sample normalization.

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments are widely used to determine, within entire genomes, the occupancy sites of any protein of interest, including, for example, transcription factors, RNA polymerases, or histones with or without various modifications. In addition to allowing the determination of occupancy sites within one cell type and under one condition, this method allows, in principle, the establishment and comparison of occupancy maps in various cell types, tissues, and conditions. Such comparisons require, however, that samples be normalized. Widely used normalization methods that include a quantile normalization step perform well when factor occupancy varies at a subset of sites, but may miss uniform genome-wide increases or decreases in site occupancy. We describe a spike adjustment procedure (SAP) that, unlike commonly used normalization methods intervening at the analysis stage, entails an experimental step prior to immunoprecipitation. A constant, low amount from a single batch of chromatin of a foreign genome is added to the experimental chromatin. This "spike" chromatin then serves as an internal control to which the experimental signals can be adjusted. We show that the method improves similarity between replicates and reveals biological differences including global and largely uniform changes.

Gene duplication and neofunctionalization: POLR3G and POLR3GL.

RNA polymerase III (Pol III) occurs in two versions, one containing the POLR3G subunit and the other the closely related POLR3GL subunit. It is not clear whether these two Pol III forms have the same function, in particular whether they recognize the same target genes. We show that the POLR3G and POLR3GL genes arose from a DNA-based gene duplication, probably in a common ancestor of vertebrates. POLR3G- as well as POLR3GL-containing Pol III are present in cultured cell lines and in normal mouse liver, although the relative amounts of the two forms vary, with the POLR3G-containing Pol III relatively more abundant in dividing cells. Genome-wide chromatin immunoprecipitations followed by high-throughput sequencing (ChIP-seq) reveal that both forms of Pol III occupy the same target genes, in very constant proportions within one cell line, suggesting that the two forms of Pol III have a similar function with regard to specificity for target genes. In contrast, the POLR3G promoter--not the POLR3GL promoter--binds the transcription factor MYC, as do all other promoters of genes encoding Pol III subunits. Thus, the POLR3G/POLR3GL duplication did not lead to neo-functionalization of the gene product (at least with regard to target gene specificity) but rather to neo-functionalization of the transcription units, which acquired different mechanisms of regulation, thus likely affording greater regulation potential to the cell.

HCFC1 is a common component of active human CpG-island promoters and coincides with ZNF143, THAP11, YY1, and GABP transcription factor occupancy.

In human transcriptional regulation, DNA-sequence-specific factors can associate with intermediaries that orchestrate interactions with a diverse set of chromatin-modifying enzymes. One such intermediary is HCFC1 (also known as HCF-1). HCFC1, first identified in herpes simplex virus transcription, has a poorly defined role in cellular transcriptional regulation. We show here that, in HeLa cells, HCFC1 is observed bound to 5400 generally active CpG-island promoters. Examination of the DNA sequences underlying the HCFC1-binding sites revealed three sequence motifs associated with the binding of (1) ZNF143 and THAP11 (also known as Ronin), (2) GABP, and (3) YY1 sequence-specific transcription factors. Subsequent analysis revealed colocalization of HCFC1 with these four transcription factors at ∼90% of the 5400 HCFC1-bound promoters. These studies suggest that a relatively small number of transcription factors play a major role in HeLa-cell transcriptional regulation in association with HCFC1.