Abstract
The mammalian circadian clock uses interlocked negative feedback loops in which the heterodimeric basic helix-loop-helix transcription factor BMAL1/CLOCK is a master regulator. While there is prominent control of liver functions by the circadian clock, the detailed links between circadian regulators and downstream targets are poorly known. Using chromatin immunoprecipitation combined with deep sequencing we obtained a time-resolved and genome-wide map of BMAL1 binding in mouse liver, which allowed us to identify over 2,000 binding sites, with peak binding narrowly centered around Zeitgeber time 6. Annotation of BMAL1 targets confirms carbohydrate and lipid metabolism as the major output of the circadian clock in mouse liver. Moreover, transcription regulators are largely overrepresented, several of which also exhibit circadian activity. Genes of the core circadian oscillator stand out as strongly bound, often at promoter and distal sites. Genomic sequence analysis of the sites identified E-boxes and tandem E1-E2 consensus elements. Electromobility shift assays showed that E1-E2 sites are bound by a dimer of BMAL1/CLOCK heterodimers with a spacing-dependent cooperative interaction, a finding that was further validated in transactivation assays. BMAL1 target genes showed cyclic mRNA expression profiles with a phase distribution centered at Zeitgeber time 10. Importantly, sites with E1-E2 elements showed tighter phases both in binding and mRNA accumulation. Finally, analyzing the temporal profiles of BMAL1 binding, precursor mRNA and mature mRNA levels showed how transcriptional and post-transcriptional regulation contribute differentially to circadian expression phase. Together, our analysis of a dynamic protein-DNA interactome uncovered how genes of the core circadian oscillator crosstalk and drive phase-specific circadian output programs in a complex tissue.