Discrete functions of nuclear receptor Rev-erbα couple metabolism to the clock

See allHide authors and affiliations

Science  26 Jun 2015:
Vol. 348, Issue 6242, pp. 1488-1492
DOI: 10.1126/science.aab3021

Multitasking around the clock

Chronic disruption of our circadian rhythms—for example, through shift work—may increase the risk of metabolic disease. Zhang et al. found that a multitasking transcription factor called Rev-erb-α regulates expression of both clock and metabolic genes through distinct mechanisms. At clock genes, it binds directly to a specific DNA sequence, displacing a competing transcription factor. At metabolic genes, it interacts not with DNA but with other transcription factors that regulate metabolic gene expression in a tissue-specific manner.

Science, this issue p. 1488


Circadian and metabolic physiology are intricately intertwined, as illustrated by Rev-erbα, a transcription factor (TF) that functions both as a core repressive component of the cell-autonomous clock and as a regulator of metabolic genes. Here, we show that Rev-erbα modulates the clock and metabolism by different genomic mechanisms. Clock control requires Rev-erbα to bind directly to the genome at its cognate sites, where it competes with activating ROR TFs. By contrast, Rev-erbα regulates metabolic genes primarily by recruiting the HDAC3 co-repressor to sites to which it is tethered by cell type–specific transcription factors. Thus, direct competition between Rev-erbα and ROR TFs provides a universal mechanism for self-sustained control of the molecular clock across all tissues, whereas Rev-erbα uses lineage-determining factors to convey a tissue-specific epigenomic rhythm that regulates metabolism tailored to the specific need of that tissue.

Circadian rhythmicity is a common feature of nearly all physiological processes (14). Each cell of the body contains a molecular clock composed of transcription factors that act on one another in interlocking feedback loops that generate near–24-hour oscillations (3, 5). A core component of the molecular clock, the nuclear receptor Rev-erbα, is expressed with a circadian rhythm (6) and represses BMAL1, a positive regulator of clock output genes (7). Rev-erbα represses many genes, often to regulate metabolism in a circadian and tissue-dependent manner (811). Thus, Rev-erbα is central to complex interactions between the core clock and metabolism.

Because Rev-erbα is a core clock component but also has tissue-specific functions, we were interested in comparing its cistromes in different mouse tissues, including liver, brain, and epididymal adipose tissue. The majority of Rev-erbα binding sites were tissue-specific (Fig. 1A), and gene ontology analyses were consistent with specialized functions of Rev-erbα (fig. S1). However, a common Rev-erbα cistrome included binding at clock genes in all tissues, consistent with its universal function in the core clock (fig. S1) (6, 7, 12).

Fig. 1 Rev-erbα represses clock genes by competing with RORα at its cognate sites.

(A) Overlap of Rev-erbα cistromes in liver (9), brain, and epididymal adipose tissue (eWAT). Most significantly enriched known motifs (abundance >10%) in common and tissue-specific cistromes are shown. (B) Heat map showing expression fold changes of genes deactivated by RORα/γ double KO (DKO/WT < –1.3, P < 0.01) and derepressed by Rev-erbα KO (αKO/WT > 1.3, P < 0.01). (C) Mean relative global run-on sequencing (GRO-seq) transcription (left) throughout 24-hour light–dark cycle, as well as oscillation amplitudes (right) of RORs/Rev-erbα common targets (red) and Rev-erbα–specific targets (blue). Time points were duplicated for clearer visualization. (D) RORα binding at clock and control genes promoters at ZT10, in WT, Rev-erbα KO (αKO), Rev-erbβ knockdown (βKD), and αKO/βKD mice liver, analyzed by ChIP-PCR. Data are expressed as mean ± SEM (* Student’s t test, P < 0.05, n = 4 mice). (E) RORα binding at clock and control gene promoters at ZT10 or ZT22 in Rev-erbα overexpression (OE) mouse liver, analyzed by ChIP-PCR. Data are expressed as mean ± SEM (* Student’s t test, P < 0.05, n = 6 or 7 mice). (F) Circadian binding of RORα at sites overlapped or not overlapped with Rev-erbα cistrome (N.S., not significant; Student’s t test, P > 0.05). (G) Percentage of common or Rev-erbα–specific target genes containing high-confidence oscillating RORα binding sites [ZT22 > 2 reads per million (rpm), ZT22/ZT10 > 1.5] within 50 kb of transcription start sites (P value from hypergeometric test).

RevDR2 and retinoid-related orphan receptor (ROR)–response element (RORE) were the most enriched motifs at Rev-erbα binding sites shared among tissues (Fig. 1A). This is consistent with earlier reports that the function of Rev-erbα as a repressive component of the molecular clock involves binding to two RORE motifs that function in the transcriptional regulation of the Bmal1 gene (7, 13). Rev-erbα recruits the nuclear receptor co-repressor–histone deacetylase 3 (NCoR-HDAC3) complex to actively repress Bmal1 transcription (13), and liver-specific deletion of HDAC3 induced Bmal1 expression (fig. S2B) at Zeitgeber time (ZT) 10, consistent with previous reports (9, 13). However, the loss of HDAC3 did not dampen circadian rhythmicity of Bmal1 or other clock components as much as the loss of Rev-erbα itself, suggesting an additional mechanism (fig. S2, A and B).

Another, non–mutually exclusive mechanism posits competition with the activating nuclear receptor ROR for the DNA binding site, which contains RevDR2/RORE motifs bound by both receptors (1417). The α and γ isoforms of ROR are most abundant in liver (18) and are expressed in a circadian manner with a peak at ZT18, antiphase to Rev-erbα (19), although the circadian variation of RORα is modest and of unclear biological importance (fig. S3, A and B) (19). Liver-specific deletion of RORα and γ markedly dampened the circadian oscillation of core clock genes in liver (fig. S3, C to E), consistent with previous reports (19). To determine target genes common to RORs and Rev-erbα, we compared gene expression in livers depleted of RORα and γ, at their peak time of expression, with gene expression from Rev-erbα knockout (KO) livers (9) (Fig. 1B). Intriguingly, genes regulated both by Rev-erbα and the RORs included clock genes such as Bmal1, Npas2, Cry1, and E4bp4 (Fig. 1B and fig. S3F) and were expressed with large circadian amplitudes, consistent with the model that Rev-erbα and RORs are both critical regulators of the clock (Fig. 1C). By contrast, Rev-erbα–specific genes had modest circadian rhythms and were enriched for liver metabolic processes (Fig. 1C and fig. S3F).

Although RORα expression was similar at ZT10 and ZT22, there was a marked difference between RORα binding to ROREs at the clock genes Bmal1 and Npas2 at these times (fig. S4A). Deletion of Rev-erbα enhanced RORα recruitment to these sites at ZT10, and this was potentiated by loss of Rev-erbβ (Fig. 1D), consistent with lower binding of RORα at ZT10 being due to competition with Rev-erbs. Conversely, hepatic overexpression of Rev-erbα reduced RORα recruitment to Bmal1 and Npas2 sites at ZT22 (Fig. 1E). Genome-wide, ~44% of RORα binding sites overlapped with Rev-erbα, and these were more likely to be circadian than were RORα-specific sites (Fig. 1F). In addition, sites of increased RORα binding at ZT22 were enriched for the RevDR2/RORE motifs bound by both Rev-erbα and RORα (fig. S4, B and C). Moreover, oscillating RORα binding sites were enriched near common target genes of RORs and Rev-erbα (Fig. 1G), further suggesting that RORα and Rev-erbα compete for binding at highly circadian genes, including core components of the molecular clock. In contrast, consistent with its expression, RORγ had a circadian binding pattern at overlapped and nonoverlapped sites (fig. S4D).

To understand why Rev-erbα and ROR tended to compete near clock genes but not Rev-erbα–specific genes, we performed chromatin immunoprecipitation (ChIP)–exonuclease followed by high-throughput sequencing (ChIP-exo) (20) in mouse liver at ZT10 to better resolve Rev-erbα binding (fig. S5A). At clock genes regulated by Rev-erbα and RORs, exemplified by Bmal1 and Cry1, the RevDR2/RORE motif was detected at ChIP-exo peaks (Fig. 2A, left). However, Rev-erbα ChIP-exo peaks were most commonly enriched for the motif bound by liver-lineage determining TF hepatocyte nuclear factor 6 (HNF6) (Fig. 2B). As exemplified by Cyp2b13 and Slc45a3, these Rev-erbα binding sites colocalized with HNF6 in mouse liver (Fig. 2A, right). Overall, the HNF6 motif was found at 1108 Rev-erbα ChIP-exo sites (Fig. 2C), the vast majority of which were also detected by HNF6 ChIP-exo in liver (21) yet did not have an RORE motif nearby (fig. S5B). The genes located nearest to these Rev-erbα/HNF6 binding sites (“Rev-erbα/HNF6-exo sites”) were enriched for lipid metabolic processes (fig. S5C), similar to Rev-erbα–specific gene regulation. Indeed, enhancer RNAs (eRNAs) at these sites bound by Rev-erbα and HNF6 had a robust circadian expression pattern (Fig. 2D) and were markedly up-regulated in livers depleted of Rev-erbα, indicating active repression of enhancer function at these sites (Fig. 2E) (22).

Fig. 2 Rev-erbα binds to the genome using both DBD-dependent and DBD-independent mechanisms.

(A) Genome browser view of Rev-erbα and HNF6 ChIP-exo signals under Rev-erbα ChIP-seq peaks near clock and metabolic genes. Blue bars indicate locations of RevDR2/RORE and HNF6 motif. (B) Highly enriched known motifs found in Rev-erbα ChIP-exo peak pairs 22 to 26 base pairs apart. (C) Heat map showing 5′-end tag densities of Rev-erbα ChIP-exo centered at HNF6 motifs within 1108 peak pairs. Red and blue indicate tag density on the plus and minus strand, respectively. (D) Mean relative eRNA transcription (22) at Rev-erbα/HNF6-exo sites throughout 24-hour light–dark cycles. Data were double plotted for clearer visualization. (E) eRNA tag density (22) centered at Rev-erbα/HNF6-exo sites near Rev-erbα target genes, in Rev-erbα KO and WT liver. (F) Heat map showing Rev-erbα ChIP-seq tag densities (at ZT10) in WT, DBD mutant (DBDm), and Rev-erbα KO (αKO) mice, at DBD-dependent and -independent sites identified among 5792 high-confidence Rev-erbα peaks (peak height > 1 rpm, WT/Rev-erbα KO > 3). (G) Sequential ChIP of Rev-erbα followed by either HNF6 or immunoglobulin G ChIP in WT and DBDm mouse liver at ZT10. Data are expressed as mean ± SEM (* Student’s t test, P < 0.05, n = 3 or 4 mice per group). (H) Binding of HDAC3 at DBD-dependent and -independent sites in WT and DBDm liver (N.S., Student’s t test, P > 0.05).

To test whether the binding of Rev-erbα to the genome can be indirect, we used a mouse model with a conditional deletion of the Rev-erbα DNA binding domain (DBD). These mice have been previously studied as a model of Rev-erbα deletion (12, 23), but the targeting strategy is predicted to lead to in-frame deletion of the DBD (fig. S6A), and Rev-erbα immunoblot of mouse liver after Cre-recombination revealed an abundant species at the approximate molecular weight of the protein lacking the DBD (fig. S6B). The identity of this protein as full-length Rev-erbα lacking its DBD was confirmed by mass spectrometric analysis of Rev-erbα immunoprecipitates from recombined liver extracts (fig. S6C). Thus, this model is actually a knock-in of a DBD mutation, rather than a complete knockout of the Rev-erbα protein. We studied the function of this Rev-erbα DBD mutant in mice whose livers were also depleted of Rev-erbβ to eliminate its compensatory effects (6, 12).

ChIP-seq analysis of Rev-erbα in livers expressing only the Rev-erbα DBD mutant (“DBDm”) revealed a comparable level of binding at a subset of wild-type (WT) sites (“DBD-independent sites”), whereas binding was markedly reduced at many other sites (“DBD-dependent sites”) (Fig. 2E and fig. S7, A and B). HNF6 ChIP-seq signals (24) were more enriched at Rev-erbα DBD-independent sites than at the DBD-dependent sites (fig. S7C), suggesting that HNF6 might tether Rev-erbα to the DNA even in the absence of Rev-erbα DBD domain.

DBD-dependent sites were enriched in RORE as well as dimeric RevDR2 motifs, in agreement with direct DNA binding (fig. S7D). These motifs are also recognized by ROR, and indeed the binding of RORα at these sites decreased markedly at ZT10, when Rev-erbα competition is strongest (fig. S7E). RevDR2 motifs were depleted in DBD-independent sites, whereas ROREs still exist in a minority of sites (fig. S7D), suggesting that the RORE may facilitate, but is not required for, DBD-independent binding. The HNF6 motif was markedly enriched at DBD-independent sites (fig. S7D), and Rev-erbα binding at Rev-erbα/HNF6-exo sites was comparable between WT and DBDm mice, following the same pattern as that seen at DBD-independent sites (fig. S7F). The simultaneous binding of Rev-erbα and HNF6 at these sites was confirmed by ChIP-re-ChIP experiments in WT liver, whereas HNF6 and Rev-erbα were not colocalized at DBD-dependent sites (Fig. 2G and fig. S7G). Enhancer RNA transcription showed circadian oscillation in phase ZT22, at DBD-dependent and -independent sites (fig. S7H), suggesting active repression of Rev-erbα in both cases. In agreement with Rev-erbα functioning by recruiting the co-repressor complex, HDAC3 binding at Rev-erbα sites was reduced in Rev-erb–depleted livers (fig. S7I). In addition, the HDAC3 ChIP-seq signal in the DBDm was reduced at DBD-dependent sites but not at sites that are DBD-independent (Fig. 2H), suggesting active repression by Rev-erbα via recruitment of HDAC3.

To determine whether HNF6 is required for Rev-erbα DBD-independent binding, we performed ChIP-seq for Rev-erbα in the liver of 129S1/SvlmJ mice and compared this result with that obtained in the C57BL/6J mice. The two strains differ by ~5.4 million single-nucleotide polymorphisms (SNPs), and SNPs were predicted to cripple the HNF6 motif at 107 Rev-erb binding sites in C57BL/6J mice and 71 Rev-erb binding sites in 129S1/SvlmJ mice. Notably, Rev-erbα binding was markedly diminished at the sites in which the SNPs disrupt either the RevDR2/RORE or the HNF6 motif in C57BL/6J or 129S1/SvlmJ mice, whereas Rev-erbα binding at random SNPs tended to be unaffected (Fig. 3, A and B, and fig. S8A). Specific examples of strain-dependent binding of Rev-erbα at HNF6 sites are shown in Fig. 3, A and B, and differential binding in the two mouse strains was confirmed by ChIP–polymerase chain reaction (PCR) for HNF6 and Rev-erbα (fig. S8, B and C). Interestingly, although the HNF6 motif was most significantly associated with strain specific Rev-erbα binding by the same analysis, we also found significant association with motifs for several other TFs that play important roles in liver function, suggesting involvement of other partners in Rev-erbα binding in the absence of RORE and RevDR2 (Fig. 3C).

Fig. 3 SNP-associated strain-specific occupancy suggests HNF6-mediated binding of Rev-erbα.

(A) Heat map showing log2 fold changes of Rev-erbα binding in 129 and B6 mice. The left column in the heat map contains 141 Rev-erbα peaks where RevDR2/RORE motif scores are higher in 129 mice and lower in B6 mice, owing to the SNPs (illustrated in the left panel). Similarly, the right column contains 101 SNP-bearing Rev-erbα peaks with better RevDR2/RORE in the B6 genome. P value was calculated with Student’s t test. (B) Same analysis as in (A), focusing on SNPs disrupting HNF6 motif under Rev-erbα peaks. (C) Heat map showing –log10 P values for other motifs that are enriched in DBD-independent Rev-erbα peaks.

The preserved binding of HDAC3, mediated by Rev-erbα DBD mutant, at metabolic genes suggested that hepatic expression of these genes might also be intact relative to livers of mice in which the Rev-erbα protein is deleted. To test this hypothesis, we compared the gene expression changes in mice lacking Rev-erbα in liver with published results using the DBDm mouse model used here, in which Rev-erbα is converted to the DBD mutant and Rev-erbβ is also deleted (12). Circadian clock genes were derepressed in both situations, demonstrating that the regulation of these genes required direct binding at RevDR2/RORE sites by Rev-erbα. Overall, only ~25% of Rev-erbα target genes that were derepressed in Rev-erbα KO mice were also derepressed in the DBDm mice (“DBD-dependent Rev-erbα targets”) (Fig. 4A). Genes derepressed specifically in Rev-erbα KO mice (“DBD-independent Rev-erbα targets”) showed circadian expression peaking at ZT22 (fig. S9A) and were enriched for lipid metabolic functions (fig. S9B), suggesting that Rev-erbα regulates circadian lipid metabolic genes independent of its DBD. In support of this, DBD-independent Rev-erbα/HDAC3 sites were enriched near DBD-independent Rev-erbα targets, whereas DBD-dependent Rev-erbα/HDAC3 sites, where Rev-erbα and HDAC3 binding was markedly reduced in DBDm mice, were more enriched near DBD-dependent Rev-erbα targets (Fig. 4A). Examples of the deletion-specific regulation of metabolic genes are shown in Fig. 4, B and C. Consistent with preserved metabolic gene expression in the DBDm mice livers, these mice did not display hepatosteatosis as is characteristic of the mice with complete deletion of Rev-erbα (Fig. 4, D and E, and fig. S9C) (6, 9).

Fig. 4 DBD-independent Rev-erbα sites regulate metabolic genes in liver.

(A) Top panel shows the number of DBD-dependent and -independent Rev-erbα target genes identified with microarrays in Rev-erbα KO (9) and DBDm mice (12). Bar graph shows ratios of DBD-independent and -dependent Rev-erbα/HDAC3 binding sites (29) located near two groups of Rev-erbα target genes (P value from hypergeometric test). (B) mRNA expression of lipid metabolic genes normalized to Arpp, measured by reverse transcription–quantitative polymerase chain reaction (RT-QPCR), in livers of Rev-erbα KO mice and WT mice at ZT10. (C) mRNA expression of lipid metabolic genes normalized to Arpp, measured by RT-QPCR, in livers of Rev-erbα DBDm (Rev-erbα/β double floxed mice injected with adeno-associated virus expressing cre-recombinase under a hepatocyte-specific thyroxine-binding globulin promoter (AAV-Tbg-Cre)) or control mice (floxed mice injected with AAV-Tbg–green fluorescent protein (GFP)) at ZT10. Data are expressed as mean ± SEM (* Student’s t test, P < 0.05, n = 4 mice per group). (D) Hepatic triglyceride (TG) concentrations in the same mice as in (B). (E) Hepatic TG concentrations in mice as in (C) (* Student’s t test, P < 0.05, n = 4 mice per group).

These findings demonstrate that Rev-erbα has a DBD-independent function that contributes to its regulation of liver metabolism. Other nuclear receptors, including estrogen receptor and glucocorticoid receptor, have DBD-independent activities through protein-protein interactions with other TFs, either directly or indirectly (25, 26). In liver, Rev-erbα is tethered to chromatin by hepatic lineage–determining TFs (fig. S10), and this mechanism of binding explains much of the nonoverlapping cistromes of Rev-erbα in different tissues, as well as the large proportion of binding sites without the RevDR2/RORE motif. In liver, the tethered cistrome is more enriched for genes with specialized function in hepatic metabolism, whereas the DBD-dependent cistrome is enriched for circadian clock genes and common to multiple tissues.

Circadian clocks and metabolism are tightly connected (1, 3, 4), and Rev-erbα has emerged as a transcriptional link from circadian rhythms to metabolism in multiple tissues (8). Our findings delineate a molecular hierarchy that governs how the clock is wired with metabolism. Direct competition between Rev-erb and ROR provides a universal mechanism for self-sustained control of the molecular clock across all tissues. On top of this basic landscape, circadianly expressed Rev-erb uses lineage-determination factors to convey a tissue-specific epigenomic rhythm that, through co-repressor and HDAC3, regulates metabolism tailored to the specific need of that tissue. These two modes of action may bestow on Rev-erbα the ability to stabilize the circadian oscillations of clock gene, while coupling liver metabolism to environmental and metabolic changes, perhaps through its endogenous ligand heme (27, 28). This raises the possibility that synthetic ligands that specifically affect Rev-erbα interaction with NCoR/HDAC3 without disrupting DNA binding could modulate liver metabolism with lesser effects on the integrity of the circadian clock.

Supplementary Materials

Materials and Methods

Figs. S1 to S10

Table S1

References (3039)

References and Notes

  1. Materials and methods are available as supplementary materials on Science Online.
  2. Acknowledgments: ChIP-seq and microarray data have been deposited in the Gene Expression Omnibus (GSE67973). We thank S. Sidoli and B. A. Garcia for assistance with mass spectrometry. We acknowledge the Functional Genomics Core and the Viral Vector Core of the Penn Diabetes Research Center (P30 DK19525) for next-generation sequencing and virus preparation, respectively. We thank the Penn Digestives Disease Center Morphology Core (P30 DK050306) for histology studies and the Molecular Profiling Core for microarray analysis. This work was supported by NIH R01 DK45586 (M.A.L.), K08 DK094968 (R.E.S.), R00 DK099443 (Z.S.), R01 DK098542 (D.J.S.), F32 DK102284 (S.M.A.), F30 DK104513 (M.J.E.), T32 GM0008275 (J.R.R.), and the Cox Medical Research Institute.
View Abstract

Navigate This Article