Rev-erb-alpha modulates skeletal muscle oxidative capacity by regulating mitochondrial biogenesis and autophagy

The nuclear receptor Rev-erb-α modulates hepatic lipid and glucose metabolism, adipogenesis and the inflammatory response in macrophages. We show here that Rev-erb-α is highly expressed in oxidative skeletal muscle and plays a role in mitochondrial biogenesis and oxidative function, in gain- and loss-of function studies. Rev-erb-α-deficiency in skeletal muscle leads to reduced mitochondrial content and oxidative function, resulting in compromised exercise capacity. This phenotype was recapitulated in isolated fibers and in muscle cells upon Rev-erbα knock-down, while Rev-erb-α over-expression increased the number of mitochondria with improved respiratory capacity. Rev-erb-α-deficiency resulted in deactivation of the Stk11–Ampk–Sirt1–Ppargc1-α signaling pathway, whereas autophagy was up-regulated, resulting in both impaired mitochondrial biogenesis and increased clearance. Muscle over-expression or pharmacological activation of Rev-erb-α increased respiration and exercise capacity. This study identifies Rev-erb-α as a pharmacological target which improves muscle oxidative function by modulating gene networks controlling mitochondrial number and function

FOXA1 repression is associated with loss of BRCA1 and increased promoter methylation and chromatin silencing in breast cancer

FOXA1 expression correlates with the breast cancer luminal subtype and patient survival. RNA and protein analysis of a panel of breast cancer cell lines revealed that BRCA1 deficiency is associated with the downregulation of FOXA1 expression. Knockdown of BRCA1 resulted in the downregulation of FOXA1 expression and enhancement of FOXA1 promoter methylation in MCF-7 breast cancer cells, whereas the reconstitution of BRCA1 in Brca1-deficent mouse mammary epithelial cells (MMECs) promoted Foxa1 expression and methylation. These data suggest that BRCA1 suppresses FOXA1 hypermethylation and silencing. Consistently, the treatment of MMECs with the DNA methylation inhibitor 5-aza-2'-deoxycitydine induced Foxa1 mRNA expression. Furthermore, treatment with GSK126, an inhibitor of EZH2 methyltransferase activity, induced FOXA1 expression in BRCA1-deficient but not in BRCA1-reconstituted MMECs. Likewise, the depletion of EZH2 by small interfering RNA enhanced FOXA1 mRNA expression. Chromatin immunoprecipitation (ChIP) analysis demonstrated that BRCA1, EZH2, DNA methyltransferases (DNMT)1/3a/3b and H3K27me3 are recruited to the endogenous FOXA1 promoter, further supporting the hypothesis that these proteins interact to modulate FOXA1 methylation and repression. Further co-immunoprecipitation and ChIP analysis showed that both BRCA1 and DNMT3b form complexes with EZH2 but not with each other, consistent with the notion that BRCA1 binds to EZH2 and negatively regulates its methyltransferase activity. We also found that EZH2 promotes and BRCA1 impairs the deposit of the gene silencing histone mark H3K27me3 on the FOXA1 promoter. These associations were validated in a familial breast cancer patient cohort. Integrated analysis of the global gene methylation and expression profiles of a set of 33 familial breast tumours revealed that FOXA1 promoter methylation is inversely correlated with the transcriptional expression of FOXA1 and that BRCA1 mutation breast cancer is significantly associated with FOXA1 methylation and downregulation of FOXA1 expression, providing physiological evidence to our findings that FOXA1 expression is regulated by methylation and chromatin silencing and that BRCA1 maintains FOXA1 expression through suppressing FOXA1 gene methylation in breast cancer.Oncogene advance online publication, 22 December 2014; doi:10.1038/onc.2014.421.published_or_final_versio

Expression of FOXA1 and GATA-3 in breast cancer: the prognostic significance in hormone receptor-negative tumours

The expression of additional genes, other than oestrogen receptor (ER), may be important to the hormone-responsive phenotype of breast cancer. Microarray analyses have revealed that forkhead box A1 (FOXA1) and GATA binding protein 3 (GATA-3) are expressed in close association with ERalpha, both encoding for transcription factors with a potential involvement in the ERalpha-mediated action in breast cancer. The purpose of this study was to explore if the expression of FOXA1 and GATA-3 may provide an opportunity to stratify subsets of patients that could have better outcome, among the ERalpha-negative/poor prognosis breast cancer group.The present study was supported by a research grant (SFRH/BD/15316/ 2005 to AA) financed by the Portuguese Science and Technology Foundation (FCT). The authors thank Prof. Raquel Seruca ( coordinator from the Cancer Genetics group at IPATIMUP) for scientific assistance, Dr Jose Luis Costa (postdoctorate at IPATIMUP) for critically reading the manuscript before submission, and Dr Nuno Marcos ( PhD student at IPATIMUP) for artwork assistance

Discovery of directional and nondirectional pioneer transcription factors by modeling DNase profile magnitude and shape

We describe protein interaction quantitation (PIQ), a computational method for modeling the magnitude and shape of genome-wide DNase I hypersensitivity profiles to identify transcription factor (TF) binding sites. Through the use of machine-learning techniques, PIQ identified binding sites for >700 TFs from one DNase I hypersensitivity analysis followed by sequencing (DNase-seq) experiment with accuracy comparable to that of chromatin immunoprecipitation followed by sequencing (ChIP-seq). We applied PIQ to analyze DNase-seq data from mouse embryonic stem cells differentiating into prepancreatic and intestinal endoderm. We identified 120 and experimentally validated eight 'pioneer' TF families that dynamically open chromatin. Four pioneer TF families only opened chromatin in one direction from their motifs. Furthermore, we identified 'settler' TFs whose genomic binding is principally governed by proximity to open chromatin. Our results support a model of hierarchical TF binding in which directional and nondirectional pioneer activity shapes the chromatin landscape for population by settler TFs.National Institutes of Health (U.S.) (Common Fund 5UL1DE019581)National Institutes of Health (U.S.) (Common Fund RL1DE019021)National Institutes of Health (U.S.) (Common Fund 5TL1EB008540)National Institutes of Health (U.S.) (Grant 1U01HG007037)National Institutes of Health (U.S.) (Grant 5P01NS055923

Unraveling the Regulatory Mechanisms Underlying Tissue-Dependent Genetic Variation of Gene Expression

It is known that genetic variants can affect gene expression, but it is not yet completely clear through what mechanisms genetic variation mediate this expression. We therefore compared the cis-effect of single nucleotide polymorphisms (SNPs) on gene expression between blood samples from 1,240 human subjects and four primary non-blood tissues (liver, subcutaneous, and visceral adipose tissue and skeletal muscle) from 85 subjects. We characterized four different mechanisms for 2,072 probes that show tissue-dependent genetic regulation between blood and non-blood tissues: on average 33.2% only showed cis-regulation in non-blood tissues; 14.5% of the eQTL probes were regulated by different, independent SNPs depending on the tissue of investigation. 47.9% showed a different effect size although they were regulated by the same SNPs. Surprisingly, we observed that 4.4% were regulated by the same SNP but with opposite allelic direction. We show here that SNPs that are located in transcriptional regulatory elements are enriched for tissue-dependent regulation, including SNPs at 3′ and 5′ untranslated regions (P = 1.84×10−5 and 4.7×10−4, respectively) and SNPs that are synonymous-coding (P = 9.9×10−4). SNPs that are associated with complex traits more often exert a tissue-dependent effect on gene expression (P = 2.6×10−10). Our study yields new insights into the genetic basis of tissue-dependent expression and suggests that complex trait associated genetic variants have even more complex regulatory effects than previously anticipated

Epistasis of Transcriptomes Reveals Synergism between Transcriptional Activators Hnf1α and Hnf4α

The transcription of individual genes is determined by combinatorial interactions between DNA–binding transcription factors. The current challenge is to understand how such combinatorial interactions regulate broad genetic programs that underlie cellular functions and disease. The transcription factors Hnf1α and Hnf4α control pancreatic islet β-cell function and growth, and mutations in their genes cause closely related forms of diabetes. We have now exploited genetic epistasis to examine how Hnf1α and Hnf4α functionally interact in pancreatic islets. Expression profiling in islets from either Hnf1a+/− or pancreas-specific Hnf4a mutant mice showed that the two transcription factors regulate a strikingly similar set of genes. We integrated expression and genomic binding studies and show that the shared transcriptional phenotype of these two mutant models is linked to common direct targets, rather than to known effects of Hnf1α on Hnf4a gene transcription. Epistasis analysis with transcriptomes of single- and double-mutant islets revealed that Hnf1α and Hnf4α regulate common targets synergistically. Hnf1α binding in Hnf4a-deficient islets was decreased in selected targets, but remained unaltered in others, thus suggesting that the mechanisms for synergistic regulation are gene-specific. These findings provide an in vivo strategy to study combinatorial gene regulation and reveal how Hnf1α and Hnf4α control a common islet-cell regulatory program that is defective in human monogenic diabetes

An Integrated Pipeline for the Genome-Wide Analysis of Transcription Factor Binding Sites from ChIP-Seq

ChIP-Seq has become the standard method for genome-wide profiling DNA association of transcription factors. To simplify analyzing and interpreting ChIP-Seq data, which typically involves using multiple applications, we describe an integrated, open source, R-based analysis pipeline. The pipeline addresses data input, peak detection, sequence and motif analysis, visualization, and data export, and can readily be extended via other R and Bioconductor packages. Using a standard multicore computer, it can be used with datasets consisting of tens of thousands of enriched regions. We demonstrate its effectiveness on published human ChIP-Seq datasets for FOXA1, ER, CTCF and STAT1, where it detected co-occurring motifs that were consistent with the literature but not detected by other methods. Our pipeline provides the first complete set of Bioconductor tools for sequence and motif analysis of ChIP-Seq and ChIP-chip data

Study of FoxA Pioneer Factor at Silent Genes Reveals Rfx-Repressed Enhancer at Cdx2 and a Potential Indicator of Esophageal Adenocarcinoma Development

Understanding how silent genes can be competent for activation provides insight into development as well as cellular reprogramming and pathogenesis. We performed genomic location analysis of the pioneer transcription factor FoxA in the adult mouse liver and found that about one-third of the FoxA bound sites are near silent genes, including genes without detectable RNA polymerase II. Virtually all of the FoxA-bound silent sites are within conserved sequences, suggesting possible function. Such sites are enriched in motifs for transcriptional repressors, including for Rfx1 and type II nuclear hormone receptors. We found one such target site at a cryptic “shadow” enhancer 7 kilobases (kb) downstream of the Cdx2 gene, where Rfx1 restricts transcriptional activation by FoxA. The Cdx2 shadow enhancer exhibits a subset of regulatory properties of the upstream Cdx2 promoter region. While Cdx2 is ectopically induced in the early metaplastic condition of Barrett's esophagus, its expression is not necessarily present in progressive Barrett's with dysplasia or adenocarcinoma. By contrast, we find that Rfx1 expression in the esophageal epithelium becomes gradually extinguished during progression to cancer, i.e, expression of Rfx1 decreased markedly in dysplasia and adenocarcinoma. We propose that this decreased expression of Rfx1 could be an indicator of progression from Barrett's esophagus to adenocarcinoma and that similar analyses of other transcription factors bound to silent genes can reveal unanticipated regulatory insights into oncogenic progression and cellular reprogramming

The transcriptional activity of hepatocyte nuclear factor 4 alpha is inhibited via phosphorylation by ERK1/2

Hepatocyte nuclear factor 4 alpha (HNF4alpha) nuclear receptor is a master regulator of hepatocyte development, nutrient transport and metabolism. HNF4alpha is regulated both at the transcriptional and post-transcriptional levels by different mechanisms. Several kinases (PKA, PKC, AMPK) were shown to phosphorylate and decrease the activity of HNF4alpha. Activation of the ERK1/2 signalling pathway, inducing proliferation and survival, inhibits the expression of HNF4alpha. However, based on our previous results we hypothesized that HNF4alpha is also regulated at the post-transcriptional level by ERK1/2. Here we show that ERK1/2 is capable of directly phosphorylating HNF4alpha in vitro at several phosphorylation sites including residues previously shown to be targeted by other kinases, as well. Furthermore, we also demonstrate that phosphorylation of HNF4alpha leads to a reduced trans-activational capacity of the nuclear receptor in luciferase reporter gene assay. We confirm the functional relevance of these findings by demonstrating with ChIP-qPCR experiments that 30-minute activation of ERK1/2 leads to reduced chromatin binding of HNF4alpha. Accordingly, we have observed decreasing but not disappearing binding of HNF4alpha to the target genes. In addition, 24-hour activation of the pathway further decreased HNF4alpha chromatin binding to specific loci in ChIP-qPCR experiments, which confirms the previous reports on the decreased expression of the HNF4a gene due to ERK1/2 activation. Our data suggest that the ERK1/2 pathway plays an important role in the regulation of HNF4alpha-dependent hepatic gene expression

Genome-Wide Analysis Reveals a Major Role in Cell Fate Maintenance and an Unexpected Role in Endoreduplication for the Drosophila FoxA Gene Fork Head

Transcription factors drive organogenesis, from the initiation of cell fate decisions to the maintenance and implementation of these decisions. The Drosophila embryonic salivary gland provides an excellent platform for unraveling the underlying transcriptional networks of organ development because Drosophila is relatively unencumbered by significant genetic redundancy. The highly conserved FoxA family transcription factors are essential for various aspects of organogenesis in all animals that have been studied. Here, we explore the role of the single Drosophila FoxA protein Fork head (Fkh) in salivary gland organogenesis using two genome-wide strategies. A large-scale in situ hybridization analysis reveals a major role for Fkh in maintaining the salivary gland fate decision and controlling salivary gland physiological activity, in addition to its previously known roles in morphogenesis and survival. The majority of salivary gland genes (59%) are affected by fkh loss, mainly at later stages of salivary gland development. We show that global expression of Fkh cannot drive ectopic salivary gland formation. Thus, unlike the worm FoxA protein PHA-4, Fkh does not function to specify cell fate. In addition, Fkh only indirectly regulates many salivary gland genes, which is also distinct from the role of PHA-4 in organogenesis. Our microarray analyses reveal unexpected roles for Fkh in blocking terminal differentiation and in endoreduplication in the salivary gland and in other Fkh-expressing embryonic tissues. Overall, this study demonstrates an important role for Fkh in determining how an organ preserves its identity throughout development and provides an alternative paradigm for how FoxA proteins function in organogenesis