Mammalian Period represses and de-represses transcription by displacing CLOCK–BMAL1 from promoters in a Cryptochrome-dependent manner

Yi-Ying Chiou, Yanyan Yang, Naim Rashid, Rui Ye, Christopher P. Selby, and Aziz Sancar
PNAS October 11, 2016 113 (41) E6072-E6079; first published September 29, 2016 https://doi.org/10.1073/pnas.1612917113

  1. Contributed by Aziz Sancar, August 17, 2016 (sent for review July 15, 2016; reviewed by Carl Hirschie Johnson and Andrew C. Liu)

See related content:

Significance

The mammalian circadian clock is controlled by a transcription-translation feedback loop consisting of transcriptional activators circadian locomotor output cycles kaput (CLOCK)–brain and muscle Arnt-like protein-1 (BMAL1), which function as a complex at E/E'-box elements, and repressors Cryptochrome 1 (CRY1)/CRY2 and PER1/PER2. CRYs repress upon binding as CRY–CLOCK–BMAL1–E-box complexes. Period proteins (PERs) repress by removing the heterotrimeric complexes from the E-box. We report here that in the Cry1 promoter, the CRY1–CLOCK–BMAL1–E-box complex represses a transcriptional activator acting in cis, and removal of the heterotrimeric complex by PER2 de-represses the transcriptional activator. ChIP-seq and RNA-seq experiments identified other genes also de-repressed by PER2. These data clarify the role of PER2 and reveal the level of complexity in regulation of Cry1 and other circadian-controlled genes.

Abstract

The mammalian circadian clock is based on a transcription-translation feedback loop (TTFL) consolidated by secondary loops. In the primary TTFL, the circadian locomotor output cycles kaput (CLOCK)–brain and muscle Arnt-like protein-1 (BMAL1) heterodimer acts as the transcriptional activator, and Cryptochrome (CRY) and Period (PER) proteins function as repressors. PER represses by displacing CLOCK–BMAL1 from promoters in a CRY-dependent manner. Interestingly, genes with complex promoters may either be repressed or de-repressed by PER, depending on the particular promoter regulatory elements. Here, using mouse cell lines with defined knockout mutations in clock genes, RNA-seq, ChIP-seq, and reporter gene assays coupled with measurements of DNA–protein interactions in nuclear extracts, we elucidate the dual functions of PER as repressor and de-repressor in a context-dependent manner.

The circadian clock is the molecular system that confers a ∼24-h periodicity to many biological processes. In mammals, CLOCK, BMAL1, Cryptochrome (CRY), and Period (PER) proteins and their paralogs are essential for generating rhythmicity (1⇓⇓–4). Rhythmicity is the product of a transcription-translation feedback loop (TTFL): The circadian locomotor output cycles kaput (CLOCK)–brain and muscle Arnt-like protein-1 (BMAL1) heterodimer binds to E-boxes in the promoters of Cry1, Cry2, Per1, and Per2 genes and activates their transcription. The CRY and PER proteins, after a time lag, enter the nucleus and inhibit their own transcription (core clock circuit) (5⇓⇓⇓–9), as well as transcription of other clock-controlled genes to maintain rhythmicity at the cellular and organismal levels. This core circuit is consolidated and stabilized by the nuclear receptors NR1D1 and NR1D2 (REV-ERBα/β), which bind to the retinoic acid response elements (RREs) in some of the clock gene promoters (10⇓⇓⇓–14). Although the canonical model has provided a useful conceptual framework, it has been insufficient in constructing a mechanistic model for the clock. This is in part due to the multiple interactions among clock proteins and to the lack of an in vitro system for analyzing the clock using purified proteins and their target promoter elements along with an appropriate readout.

To simplify the analysis of clock protein functions, previously we reported the construction of cell lines that lack CRYs, PERs, and both CRYs and PERs and that express ectopic CRY or PER that can be targeted to the nucleus in a controllable manner (15). Using this system initially, we analyzed the canonical clock model by restricting our experiments to the Nr1d1 and D-site albumin promoter binding protein (Dbp) genes, which are controlled solely by the binding of CLOCK–BMAL1 to E-box regulatory elements. These experiments confirmed the repressor functions of CRY and PER but also revealed some features of their mechanisms: CRY alone can bind to the CLOCK–BMAL1–E-box complex and inhibit the transcription of cognate genes by forming a stable CRY–CLOCK–BMAL1 ternary complex in the promoter: “blocking-type” repression. In contrast, PER represses CLOCK–BMAL1–E-box-mediated transcription by removing the CLOCK–BMAL1 complex from the E-box in a CRY-dependent manner: “displacement-type” repression (15). Although these findings clarified the roles of PERs and CRYs in regulation of simple promoters exclusively controlled by E-box binding of transcription factors, they did not address the issue of the control of key clock genes with multiple cis-elements.

Here, we report the construction of cell lines with clock gene knockout mutations and expression systems to analyze the mechanisms of control of more complex promoters by PER2 protein. We show that the removal of CLOCK–BMAL1 from promoters affects gene expression in three different ways, depending on the type of regulatory elements present in the target promoters. First, in the case of a simple promoter controlled exclusively by an E-box, such as in the Nr1d1 gene, CLOCK–BMAL1 removal represses transcription. Second, in genes such as Bmal1, which is controlled exclusively by an RRE, CLOCK–BMAL1 removal from the Nr1d1 promoter indirectly activates transcription of Bmal1 by down-regulation of NR1D1/2, which represses by binding to the RRE in the Bmal1 promoter. Third, in the case of the Cry1 gene, removal of CLOCK–BMAL1 facilitates transactivation by other transcription factors that bind to the Cry1 promoter independently of NR1D1/2.

Results

Experimental System to Analyze the Role of PER in the TTFL.

In Fig. 1 we present a highly idealized model for the mammalian circadian clock: The CLOCK–BMAL1 heterodimer activates transcription of Cry1/2, Per1/2, and Nr1d1/2 by binding to E-box (CACGTG/T) sequences in the promoters of these genes. The CRY and PER proteins inhibit transcriptional activation by CLOCK–BMAL1 to close the feedback loop, CRY by binding to the CLOCK–BMAL1–E-box and interfering with transactivation, and PER by causing the displacement of CLOCK–BMAL1 from E-boxes in a CRY-dependent manner. This primary circuit is consolidated by a stabilizing loop: The Nr1d1/2 nuclear receptor genes are also controlled by CLOCK–BMAL1, and NR1D1/2 proteins in turn repress Bmal1 and Cry1 transcription by binding to the RRE elements in the promoters of these genes (12⇓–14, 16, 17). Some observations, such as the repressor activity of CLOCK–BMAL1 (18, 19) and the apparent transcriptional up-regulation by PER (20⇓⇓–23), do not readily fit into this model and require further experimentation.

Fig. 1.
Fig. 1.

TTFL model for the mammalian circadian clock. In the TTFL model, CLOCK–BMAL1 heterodimers bind to the E-boxes of Cry1/2, Per1/2, and Nr1d1/2 promoters and activate the transcription of these genes. Increased CRY protein then inhibits transcription through binding to the CLOCK–BMAL1 complex. PER protein inhibits transcription in a CRY-dependent manner by reducing CLOCK–BMAL1 binding to the E-box. This core circuit (solid lines) is stabilized by a secondary feedback loop (dashed lines) in which CLOCK–BMAL1-controlled NR1D1/2 inhibits the transcription of Bmal1 and Cry1 through binding to the RRE of their promoters.

Because there is no in vitro system to analyze the roles of different clock proteins, individually or in various combinations, on transcription of cognate genes, we developed cell lines with various clock gene knockout mutations (Table S1) and with the controllable PER2 nuclear entry system shown in Fig. 2A for carrying out “in cellulo biochemistry” for better insight into the clock mechanism. These cell lines were isolated from knockout mouse strains or from mouse embryonic fibroblasts by using transcription activator-like effector nuclease (TALEN) and CRISPR technologies. The cell lines also ectopically express PER2-ER*, which is composed of PER2 protein fused to a mutant form of the ligand-binding domain of the estrogen receptor (ER*). PER2-ER* is trapped in the cytoplasm, and it enters the nucleus upon 4-hydroxytamoxifen (4-OHT) treatment. Thus, the effect of PER2 on gene expression can be analyzed in a controllable manner as shown in Fig. 2B. As previously reported (15), PER2 and CRY1 enter the nucleus together in this system.

Fig. 2.
Fig. 2.

Cell-based system for study of mammalian clock regulation. (A) Schematic of the experimental system. PER2 fused to a mutant form of the estrogen receptor ligand-binding domain (ER*) is constrained to the cytosol when bound to endogenous heat-shock protein (HSP). Addition of 4-OHT to the cell culture displaces the HSP, and PER2-ER* then enters the nucleus. (B) Immunoblot analysis of clock proteins upon 4-OHT addition. Total and nuclear amounts of clock and control proteins (MEK2 and LaminA, which is a nuclear protein) up to 6 h following addition of 4-OHT to Per1/2−/−;PER2-ER* cells are shown. The blot shows that the total NR1D1 protein level (and thus nuclear NR1D1) decreases following incubation with 4-OHT. An increased level of nuclear CRY1 but not total CRY1 is associated with the increased nuclear PER2. This supports the idea that CRY and PER enter the nucleus together and that CRY-dependent PER2 activity is through the formation of a CRY–PER complex. NR1D1.s and NR1D1.L are short and long exposures of the same blot.

View this table:
Table S1.

Cell lines used in this study

PER2 as a Repressor and Antirepressor.

We addressed the regulatory roles of PER using a PER2-ER* ectopic expression system in cells with various clock gene knockout mutations. We analyzed transcription of three genes selected for their different promoter complexities, and reported differences in their regulation by PER2 (9⇓⇓⇓–13). The plots in Fig. 3A show the levels of mRNA expressed from the Nr1d1, Bmal1, and Cry1 genes in Per1/2−/− cells following 0–6 h of treatment with 4-OHT, which induces translocation of PER2-ER* to the nucleus. The Nr1d1 promoter is controlled solely by an E-box, the Bmal1 promoter is controlled solely by the RRE element, and the Cry1 promoter is controlled by an E-box, an RRE element, and a D-box (14). As expected (15), the data in Fig. 3A (PER2-ER* enters the nucleus of Per1/2−/− cells) show that PER2 inhibits transcription of Nr1d1 (by removing CLOCK–BMAL1 from the E-box; Fig. 3A, Left, red); however, PER2 also induces transcription of Bmal1 and Cry1 (Fig. 3A, Middle and Right, red). Because PER2 inhibits the expression of Nr1d1 and because both Bmal1 and Cry1 have RRE elements in their promoters to which NR1D1/2 binds (12⇓–14, 16, 24), it has been suggested that PER2 de-represses genes by down-regulating the NR1D1 repressor. To test this prediction, we investigated the effect of Nr1d1/2 deletions on PER2-ER*–mediated up-regulation of Bmal1 and Cry1 transcription. In the absence of NR1D1/2, there is no up-regulation of Bmal1 transcription by PER2-ER* (Fig. 3A, Middle, blue), supporting the view that PER2 up-regulates Bmal1 by down-regulating NR1D1/2. In contrast, and surprisingly, the NR1D1/2 deletions have no effect on PER2-mediated up-regulation of Cry1 (Fig. 3A, Right, blue). This agrees with the observation that NR1D1/2 is dispensable for the rhythmic expression of Cry1 (12) and indicates that PER2 up-regulates Cry1 transcription by a mechanism independent of NR1D1/2.

Fig. 3.
Fig. 3.

Effect of PER2 nuclear entry on transcription of Nr1d1, Bmal1, and Cry1 in Per1/2−/−, Per1/2−/−;Nr1d1/2−/−, Per1/2−/−;Cry1/2−/−, and Bmal1−/− cells. (A) Nuclear entry of PER2-ER* following addition of 4-OHT (red) represses Nr1d1 transcription (Left) but induces Bmal1 (Middle) and Cry1 (Right) transcription analyzed by RT-quantitative PCR. The effect of PER2 on these three genes is CRY-dependent (green). NR1D1/2 is required for PER2-induced Bmal1 transcription but is not required for the repression of Nr1d1 or for the induction of Cry1 transcription (blue). In Bmal1−/− cells (orange), nuclear entry of PER2-ER* has no effect on Nr1d1, Bmal1, or Cry1 transcription. Endogenous PER2 was not detected in Bmal1−/− cells (4). (B) PER2-dependent, CRY-dependent, NR1D1/2-independent displacement of CLOCK–BMAL1 from Nr1d1 and Cry1 E-boxes is analyzed by BMAL1–ChIP–quantitative PCR. Bmal1, which has no E-box and is not a primary target of CLOCK–BMAL1, is a negative control for ChIP. SE values are displayed in A and B.

PER2 Repressor and Antirepressor Activities Are BMAL1- and CRY-Dependent.

Based on indirect evidence, it has been previously reported that PER2 can up-regulate gene transcription by interacting with another transcription factor (25) or by removing CRY from CLOCK–BMAL1 (20). We have shown that PER2 represses transcription by removing the CLOCK–BMAL1 complex from cognate promoters in a CRY-dependent manner (15). We wished to know if the transcriptional up-regulation effect of PER was also mediated through removal of the CLOCK–BMAL1–CRY complex. To this end, we tested the effect of PER2-ER* nuclear entry on transcription of the three sentinel genes in cells lacking BMAL1. In these Bmal1 knockout cells, endogenous PER2 is undetectable (4), and thus any effect that might be seen upon 4-OHT addition can be ascribed to PER2-ER*. The results (Fig. 3A, orange) show that in the absence of BMAL1 there is no repression or de-repression of these genes upon nuclear entry of PER2-ER*. This excluded the possibility that PER2 facilitates a transcriptional activator in a BMAL1-independent manner to activate Cry1 transcription. To examine whether PER2-induced up-regulation is dependent on CRY proteins, we tested the effect of PER2-ER* on expression of all three sentinel genes in cells lacking CRY1/2. As seen in Fig. 3A, the effect of PER2 on all these three genes (red) is abolished in the Cry1/2 knockout background (green), indicating that PER2 exerts its up-regulation activity only through a CRY-mediated mechanism as it does its repressor activity.

Gene Up-Regulation by PER2 Is Associated with Removal of CLOCK–BMAL1.

We then used ChIP to directly examine the binding of BMAL1 to the sentinel gene promoters in response to PER2. In Fig. 3B, Middle, there is no BMAL1 binding to the Bmal1 promoter, as it lacks an E-box (1, 2, 4). Fig. 3B, Left shows that, as expected, nuclear entry of PER2 following 4-OHT treatment causes displacement of BMAL1 from the (PER2-repressed) Nr1d1 promoter (red), and this displacement is NR1D1/2-independent (blue) but CRY-dependent (green). Importantly, the positive regulation of the (PER2-activated) Cry1 promoter shown in Fig. 3A, Right coincides with NR1D1/2-independent but CRY-dependent displacement of BMAL1 (Fig. 3B, Right). Taken together, these data indicate that BMAL1 and CRY are required for both PER-mediated transcriptional repression and up-regulation in a context-dependent manner. Indeed, the idea of repression by CRY in a CLOCK–BMAL1-dependent manner was suggested in a previous study, although there was no direct experimental evidence (19). We suggest that in the Cry1 promoter, CLOCK–BMAL1 binding to an E-box represses another activator, which is the predominant transcription activator for this gene. Here, we show that the PER2-mediated up-regulation of Cry1 gene transcription is accompanied with the removal of the CLOCK–BMAL1–CRY complex from the Cry1 promoter and thus present direct evidence for these counterintuitive roles of CLOCK–BMAL1 in gene repression and PER2 (with CRY) in gene de-repression.

Genome-Wide Effect of PER on Regulation and Promoter Binding by the CLOCK–BMAL1 Complex.

The model that emerges from the results presented here as well as our previous studies is that PER displaces CLOCK–BMAL1 from cognate gene promoters and in doing so, depending on the complement of regulatory elements in the target promoter, it either inhibits or up-regulates (de-represses) transcription in a CRY-dependent manner (15). Because there are several thousands of CLOCK–BMAL1 binding sites in the genome (26, 27), we wished to find out whether PER-mediated CLOCK–BMAL1 removal and both the repression and de-repression functions of PER are observed on a global scale. To this end, we performed ChIP-seq experiments using an anti-BMAL1 antibody to probe CLOCK–BMAL1 binding sites and RNA-seq experiments to analyze RNA levels before and after nuclear entry of PER2 in Per1/2−/−;Nr1d1/2−/−;PER2-ER* cells so as to avoid complications arising from second-order clock-controlled genes regulated by NR1D1/2 (such as Bmal1).

ChIP-seq experiments were conducted in duplicate with cells not treated with 4-OHT. Specific BMAL1 binding sites were identified based on 10–15 million input and precipitated DNA reads (Materials and Methods). There were 4,789 and 8,615 binding sites in the two experiments and 4,740 common binding sites seen in both experiments. Induction of PER2 with 4-OHT reduced the number of common BMAL1 binding sites to 483, and PER2 induction produced no novel BMAL1 binding sites. Thus, PER2 eliminated detectable binding of BMAL1 to most of the regulatory regions to which it binds. To further assess BMAL1 binding, we analyzed reads mapped to the 4,740 common BMAL1 binding sites as relative BMAL1 binding strength, expressed as “IP/input.” We observed that relative BMAL1 binding strength was negatively associated with PER2 induction and was highly correlated between the two independent experiments (Fig. 4A, Left). BMAL1 binding was stronger without 4-OHT than with 4-OHT on almost all regions (Fig. 4A, Right). To determine whether the observed decrease in BMAL1 binding was due to experiment-specific technical factors such as sequencing depth, we compared binding differences from these BMAL1 binding sites to reads randomly selected from the genome. Relative to the BMAL1 binding sites, there is no difference of BMAL1 binding between without 4-OHT and with 4-OHT in the random areas (Fig. 4B). Results of motif analysis of the BMAL1 binding sites showed high enrichment of the E-box sequence (Fig. 4C), which further supports the finding that PER2 removes CLOCK–BMAL1 from the E-box sequence of the promoters.

Fig. 4.
Fig. 4.

Effect of PER2 nuclear entry on BMAL1 promoter binding by BMAL1–ChIP-Seq in Per1/2−/−;Nr1d1/2−/− cells. (A) Strength of BMAL1 binding (ChIP/Input) to 4,740 common binding sites is shown as a scatter plot. (Left) Results with (green) and without (red) 4-OHT show high correlation between the two experiments. The results obtained with 4-OHT are closest to the origin, indicating weak or no binding of BMAL1 following PER2 induction. (Right) The strength of BMAL1 binding in the condition with (x axis) and without (y axis) 4-OHT is plotted. The distributions largely overlap in the two experiments and are above the line representing y = x, indicating reduced BMAL1 binding following PER2 induction. (B) Distribution of relative strength of BMAL1 binding –4-OHT/+4-OHT was plotted by violin and box plot. Stronger binding of BMAL1 without 4-OHT was detected on the BMAL1 binding sites (red) but not with reads randomly selected from the genome (blue). (C) Result from motif analysis of BMAL1 binding sites shows strong enrichment of the E-box sequence.

To analyze the effect of PER2 on CLOCK–BMAL1-controlled gene expression, we first defined CLOCK–BMAL1-controlled genes as genes with transcription start sites within ±5 kb from BMAL1 binding sites. Based on our ChIP-seq data, 3,317 genes were selected for the analysis. As apparent from the volcano plot in Fig. 5, nuclear entry of PER2-ER* changes the expression pattern, with 53 transcripts going down, 28 going up, and most remaining unchanged (Table S2). Interestingly, Cry1 and Cry2 are increased and Nr1d1 and Nr1d2 are decreased when PER2 enters the nucleus. To ascertain the quality of the data from BMAL1 ChIP-seq and RNA-seq, in Fig. 6 we focus on Cry1 and Nr1d1 loci, on chromosomes 10 and 11, respectively: In agreement with the RT-quantitative PCR and ChIP-quantitative PCR data in Fig. 3, nuclear entry of PER2 up-regulates Cry1 expression and represses the Nr1d1 transcript level and both are associated with a decrease in BMAL1 binding on the promoter. Similarly, increased Cry2, Dec1 (Bhlhe40), and decreased Nr1d2 transcription (Table S2) are confirmed by RT-quantitative PCR (Fig. S1).

Fig. 5.
Fig. 5.

Effect of PER2 nuclear entry on transcriptome by mRNA-Seq in Per1/2−/−;Nr1d1/2−/− cells. Effect of PER2 nuclear entry on transcription of individual genes is shown by Volcano plot with the x axis representing level of difference (+4-OHT vs. –4-OHT) and the y axis representing the level of statistical significance. Genes with a statistically significant difference between with and without 4-OHT are shown in red.

Fig. 6.
Fig. 6.

Visualization of BMAL1–ChIP-seq and RNA-seq on Nr1d1 (Top) and Cry1 (Bottom) gene loci. (Top) An 11 kbp region of chromosome 11 containing the Nr1d1 gene. (Bottom) A 57 kbp region of chromosome 10 containing the Cry1 gene. Gene structures are illustrated along the x axis beneath each panel, with E/E’-box regulatory elements indicated with green asterisks. Within each panel, the magnitude of the ChIP-seq signal (y axis) is shown in the top half, and the magnitude of the RNA-seq signal (y axis) is shown in the bottom half. Signals observed with 4-OHT are shown in red, and signals without 4-OHT are shown in blue. With both genes, lower ChIP-seq signals were observed with 4-OHT (red) than without 4-OHT (blue). However, following induction of PER2 with 4-OHT, RNA-seq data show a weaker signal in the Nr1d1 gene and a stronger signal in the Cry1 gene. Data were visualized using IGV software.

View this table:
Table S2.

Genes repressed (left) and enhanced (right) following PER2 nuclear entry after 4-OHT treatment were identified by change levels (log2FoldChange)

Fig. S1.
Fig. S1.

Nuclear entry of PER2-ER* following addition of 4-OHT (4 h) induces Cry2 (Left) and Dec1 (Middle) transcription but represses Nr1d2 (Right) transcription in Per1/2−/−;Nr1d1/2−/−;PER2-ER* cells analyzed by RT-quantitative PCR. SEs are plotted.

Decreased CLOCK–BMAL1 Binding on the Cry1 Promoter Increases Cry1 Transcription.

Finally, the experiments described so far were carried out with 4-OHT–induced nuclear entry of PER2-ER*, and it could be argued that in this model system, the observed removal of CLOCK–BMAL1 from promoters and the accompanying up-regulation of Cry1 might be due to some unique feature of our system and that the CLOCK–BMAL1 removal from the E-box upon PER2-ER* nuclear entry and Cry1 up-regulation might not be physiologically relevant. To address this issue, we conducted reporter gene assays using a 1.5 kbp promoter fragment of Cry1, which can sustain rhythmic transcription (14). This fragment possesses two D-box and two E-box sequences that partially overlap, as shown in Fig. 7A. To characterize the effect of CLOCK–BMAL1 binding on Cry1 transcription, we first compared the reporter activity of Cry1 promoters with either wild-type (WT) or mutant E-box (mutE) sequences (Fig. 7A). In a preliminary control experiment to ascertain the specificity of CLOCK–BMAL1 binding to the E-box in the Cry1 promoter, the two 1.5 kbp promoter fragments, end-labeled with biotin, were incubated with NIH 3T3 cell-free extract. The fragments were pulled down and the bound BMAL1 was visualized in the Western blot shown in Fig. 7B. BMAL1 bound to the WT promoter but only weakly to the mutant promoter.

Fig. 7.
Fig. 7.

Decreased CLOCK–BMAL1 binding increases reporter gene expression driven from the Cry1 promoter. (A) Sequences of the Cry1 promoter region including WT, mutant E-boxes, and D-boxes are shown. Nucleotides conserved in the E-box and D-box sequences are colored in green and blue, respectively. Nucleotides mutated in mutE (E-box mutation) and mutDE (D-box and E-box mutation) are colored in red. (B) CLOCK–BMAL1 binding affinity to Cry1 promoters with WT or mutE sequences. BMAL1 was pulled down from nuclear extract from NIH 3T3 cells by immobilized biotin–DNA and analyzed by Western blot. (C) Reporter gene assay of the Cry1 promoters with the WT or mutant promoter elements is shown. Reporter signals were first normalized to the signals from cotransfected control plasmid and then were normalized to the condition without Cry1 and Per2 plasmids. Note that reporter gene expression driven by the mutant Cry1 promoter is higher than expression driven by the WT Cry1 promoter. SEs are plotted in C.

Reporter assays were then performed following cotransfection of the Cry1 promoter/reporter constructs with and without PER2 and/or CRY1-expressing plasmids, which is a commonly used circadian experimental system. Fig. 7C, lane 1 shows that in the absence of exogenous PER2 or CRY1, there was much more expression from the mutE (green) than the WT (red) promoter. This increase is consistent with strong transcription from another element on the 1.5 kbp fragment that is inhibited by CLOCK–BMAL1–CRY1 binding to the WT E-box. Expression of PER2 led to substantially increased transcription from the WT promoter fragment (lanes 2 and 3). Apparently the basal level of CRY1/2 in these cells is sufficient to enable the removal of CLOCK–BMAL1 by exogenous PER2 so as to reveal strong transcription. A modest increase in transcription of the mutE promoter by exogenous PER2 is consistent with the removal of residual CLOCK–BMAL1 bound weakly to the mutant promoter (green, lane 1 vs. lanes 2 and 3). Comparing lane 4 versus lane 1 (red) shows that exogenous CRY1 inhibited transcription mediated by CLOCK–BMAL1. Increasing the amount of exogenous PER2, in the presence of exogenous CRY1 (lanes 5–8), resulted in increased transcription, as PER2 removed CLOCK–BMAL1 from the E-boxes in a CRY1-dependent manner. These data from cellular oscillators support our conclusion that up-regulation of Cry1 by PER2 is not simply due to removal of CRY1 to facilitate CLOCK–BMAL1-mediated transcription but in fact is due to the removal of the entire CRY1–CLOCK–BMAL1 “repressor complex” to facilitate the activation of transcription by another sequence element in the 1.5 kbp Cry1 promoter region.

It has been shown that Cry1 transcription could be regulated through the D-box in the promoter (14). It is possible that removal of CLOCK–BMAL1 from the E-boxes by PER facilitates the transcription of Cry1 through the nearby and overlapping D-boxes (Fig. 7A). We considered that D-box binding proteins, such as DBP, TEF, and HLF, would activate the transcription of Cry1 when not blocked by the binding of CLOCK–BMAL1–CRY to the E-boxes. To test this, we measured the reporter activity of a Cry1 promoter with mutated D-box and E-box (mutDE, Fig. 7A). As seen with the mutE promoter, mutDE promoter activity in the absence of exogenous PER2 or CRY1 was much higher than WT (lane 1, blue) and furthermore was very similar to mutE in all conditions tested. Thus, the increase in Cry1 transcription mediated by PER as a result of CLOCK–BMAL1 dissociation is not by removal of a repressor complex interfering with D-box–mediated activation. It is most likely that another transcriptional activator with a binding site in the 1.5 kbp Cry1 promoter fragment is responsible for increased Cry1 expression following removal of CLOCK–BMAL1. For purposes of discussion, this possible transcription factor will be referred to as factor T, which binds to a postulated T element in the Cry1 promoter (Fig. 8).

Fig. 8.
Fig. 8.

Model of PER-mediated regulation of Nr1d1, Bmal1, and Cry1 genes. This model represents transcription and regulation patterns of the Nr1d1, Bmal1, and Cry1 genes individually as a function of the circadian time. Gradient-colored bars represent transcription of Nr1d1 (blue), Bmal1 (orange), and Cry1 (green) genes based on nascent RNA-seq data from mouse liver (27). Above the bars are the proposed status of promoter occupation and transcription. Spacing of the T element, E-box, and RRE in this representation of the Cry1 promoter is for illustrative purposes only. The RRE element is in intron 1 of Cry1 ∼23 kbp downstream from the E-box. At the bottom of the model are the levels of NR1D1, CRY1, and PER2 in mouse liver nuclei. The levels of these and additional proteins from different times were cropped from a single image (Fig. S2). For Nr1d1, which is controlled solely by an E-box, transcription is inhibited by either binding of CRY on CLOCK–BMAL1 or removal of CRY1–CLOCK–BMAL1 from the E-box by PER. For Bmal1, which is controlled solely by an RRE, transcription is inversely correlated with the NR1D1 level. For Cry1, which is controlled by an E-box, RRE, and a putative factor T, binding of either the CRY–CLOCK–BMAL1 complex on the E-box or NR1D1 on the RRE represses the transcription of Cry1 from another element (T). PER eliminates both pathways of inhibition by removing CLOCK–BMAL1 from the E-box.

Fig. S2.
Fig. S2.

Original Western blot image used in Fig. 8 shows the changes of different circadian proteins in the nucleus of mouse liver at different Zeitgeber Times (ZT).

Discussion

The mammalian molecular circadian clock consists of a TTFL generated by four genes (core clock genes) and their paralogs. Research over the past two decades has supported the TTFL model and identified a secondary loop involving nuclear receptors Nr1d1/2. Furthermore, identification of kinases (28, 29) and ubiquitin ligases (30, 31) that modulate the activity and stability of the core clock proteins has provided insights into the mechanisms that are necessary to establish a daily rhythm of high amplitude, precise period, and flexibility for phase resetting in response to stimuli.

Despite these important developments, the negative arm of the TTFL model has remained poorly defined. In particular, the respective roles of CRY and PER in repression have been unclear. There are several reasons for this uncertainty. First, the two proteins interact strongly and aid in one another’s stability and nuclear entry (32⇓–34). Second, removal of one affects the posttranslational modification and stability of the other. Also, although CRY has been shown to be capable of repressing target genes on its own, there is no convincing evidence that PER affects the TTFL in the absence of CRY (15). Moreover, depending on the clock gene tested, PER has been reported to function as either a repressor (for Nr1d1) or an up-regulator (for Cry1 and Bmal1).

To examine the clock model, we have developed an “in cellulo biochemical” system consisting of mouse cell lines with knockout mutations in one or more of the core clock genes and that have a controllable delivery system for CRY or PER proteins (15). In this study, using the PER2 delivery system, we have obtained RT-quantitative PCR, RNA-seq, ChIP, and ChIP-seq data that have enabled us to revise the consensus mammalian clock model so as to reconcile various views regarding the roles of CRY and PER in the negative arm of the TTFL and eliminate the “internal inconsistencies” of the conventional model (3).

We previously reported that PER causes the removal of CLOCK–BMAL1 from the cognate promoters in a CRY-dependent manner and that it inhibits the transcription of certain genes with simple promoters (promoters with only an E-box control element) by this mechanism. In this study, we investigated the effect of PER2 on transcription of genes with more complex promoters. In aggregate, our work shows that PER2 exerts three types of effect on clock gene transcription, depending upon the type of promoter. Furthermore, we show that all of these effects are CRY-dependent, and moreover our data reconcile some previous seemingly contradictory reports on PER function in the circadian clock. We explain the three different effects of PER2 on circadian gene expression as follows:

  • i) Nr1d1 (repressed by PER2): This gene is predominantly regulated by an E-box. Binding of CRY1 (or CRY2) to the CLOCK–BMAL1–E-box complex represses transcription. When PER2 is abundant enough, it removes the CRY1–CLOCK–BMAL1 complex from the promoter as part of the dual repression mechanism of CRY and PER (15).

  • ii) Bmal1 (de-repressed by PER2): This gene is predominantly regulated by NR1D1/NR1D2 nuclear receptors; CRY1/CRY2 de-repress its transcription. We show that de-repression is through the effects of CRY1/2 and PER1/2 on NR1D1/NR1D2 expression. When CRY and PER are at high enough concentrations, they inhibit NR1D1/NR1D2 expression through the E-box and thus down-regulate the repressors of Bmal1 causing up-regulation of Bmal1 transcription. Knockout of Nr1d1/Nr1d2 eliminates the effect of CRYs and PERs on Bmal1 transcription, which is mediated by NR1D1/NR1D2.

  • iii) Cry1 (de-repressed by PER2): This gene has several transcriptional elements, three of which have been studied previously—D-box (day element), E-box (morning element), and RRE (evening element)—which are the targets of transcription factors DBP, CLOCK–BMAL1, and NR1D1/NR1D2, respectively (14). Up-regulation (de-repression) of the Cry1 promoter by PER2 seen in our experiments is not consistent with the accepted role of PER2 as a repressor. Kondratov et al. (19) and Liu et al. (12) observed related, unexpected results, namely the up-regulation of Cry1 in Bmal1 knockout cells. Kondratov et al. (19) suggested that with some genes, such as Cry1, there is a high “basal” level of transcription that is not substantially increased by CLOCK–BMAL1 binding to the promoter, and binding of CRY1 to the CLOCK–BMAL1 complex inhibits both CLOCK–BMAL1 activated and basal transcription. In this vein, our findings are consistent with the existence of an additional transcriptional activator (factor T) that binds to a regulatory element (T element) in the Cry1 promoter (12) and is inhibited by the CRY1–CLOCK–BMAL1 complex. In this scenario, removal of CRY1-CLOCK-BMAL1 from the Cry1 promoter by PER2 restores transcription via factor T.

Regarding the elevated Cry1 level in the Bmal1 knockout cells, Liu et al. (12) hypothesized that the absence of BMAL1 led to less NR1D1/NR1D2, which normally binds to the RRE element and represses transcription of Cry1. In fact, in their experiments, overexpression of NR1D1 in Bmal1 knockout cells did inhibit Cry1 expression. Thus, their interpretation of BMAL1 indirectly inhibiting Cry1 expression through NR1D1/NR1D2 appears appropriate to their defined experimental system, and it is reasonable to expect that NR1D1/NR1D2 influence Cry1 expression in vivo. However, our cell-based (Fig. 3) and reporter assay (Fig. 7) data show that PER2-mediated expression of Cry1 is independent of NR1D1 (our reporter lacked the NR1D1 binding site). These seemingly conflicting results are reconciled and incorporated in the revised clock model shown in Fig. 8, which also shows levels of NR1D1, CRY1, and PER2 proteins in mouse liver at various circadian times. At ZT10, NR1D1 is high and it likely represses Cry1 by binding to the RRE element, whereas at ZT18 NR1D1 is low, and this likely contributes to initiate Cry1 expression. Also, at ZT18, PER2 is high and PER2 at this time likely helps initiate Cry1 transcription by removing the inhibitory CRY1–CLOCK–BMAL1 complex. Thus, in this model, PER2 has dual roles in regulation of Cry1. First, through displacement of CRY–CLOCK–BMAL1 from the Nr1d1/Nr1d2 promoters, it represses NR1D1/NR1D2 and thus indirectly activates Cry1 expression. Second, through displacement of CRY–CLOCK–BMAL1 from the Cry1 promoter, it directly de-represses Cry1 expression.

We note that our findings contradict a recent report (20) that claimed that PER2 up-regulates Cry1 transcription by removing CRY1 from the CRY1–CLOCK–BMAL1 complex at the Cry1 E-box and at the E-boxes of other genes regulated by CLOCK–BMAL1. Here, we confirm our earlier finding (15) that in fact PER2 nuclear entry causes the removal of CLOCK–BMAL1 from the Cry1 promoter in a CRY-dependent manner while at the same time up-regulating Cry1 transcription.

Because the D-box in the Cry1 promoter has been shown to play a role in Cry1 regulation, we considered transcription factors (DBP, TEF, and HLF) that bind to the D-box as potential Cry1 transcriptional activators when PER2 removes the CRY–CLOCK–BMAL1 complex from the Cry1 promoter. We tested this possibility by using reporter gene assays with a 1.5 kbp fragment carrying the Cry1 promoter including the D-box and E-box elements. As expected, in this system, CRY1 repressed and PER2 together with CRY1 de-repressed transcription of Cry1. MutE and double mutant E-box/D-box constructs, with greatly reduced CLOCK–BMAL binding, demonstrated reduced transcription stimulation following PER2 induction. Thus, DBP cannot be factor T whose activity is revealed by dissociation of the CRY1–CLOCK–BMAL1 complex. Additional studies are needed to identify the postulated T element and T factor that acts in a tonic manner (12) to regulate Cry1 transcription and whose effect is rendered circadian in the context of the entire circadian ensemble.

Materials and Methods

Cells and Antibodies.

Cell lines (Table S1) were maintained in DMEM supplemented with 10% (vol/vol) FBS in 5% (vol/vol) CO2 at 37 °C. Blasticidin (5 μg/mL) was present during cell culture and absent during experiments. To induce PER2-ER* nuclear entry, 4-hydoxytamoxifen was added to the medium to 1 μM when cells were ∼90% confluent.

The Bmal1−/−, Per1/2−/−, and Cry1/2−/−;Per1/2−/− mouse embryonic fibroblasts have been described (15, 32). The Per1/2−/−;Nr1d1/2−/- mouse embryonic fibroblast cell line was made by CRISPR technology using LentiCRISPRv2 (35) obtained from Addgene to mutate the Nr1d1/2 alleles in Per1/2−/− fibroblasts. Single colonies lacking NR1D1/2 proteins were isolated and screened by Western blot, and mutational inactivation was confirmed by genomic DNA sequencing. Mutated sequences in the isolated clone (PN-P2ER) are shown in Fig. S3. Lack of NR1D1/2 protein and movement of PER2-ER* to the nucleus are shown in Fig. S4.

Fig. S3.
Fig. S3.

Target sequences (highlighted in red) for Nr1d1 and Nr1d2 knockout are both located in Exon2 of the genes. Sequences of Nr1d1 and Nr1d2 near the target site from Per1/2−/−;Nr1d1/2−/−;PER2-ER* cells indicate frame-shifted mutation of both genes on both alleles (Nr1d1: –8 nt, –8 nt; Nr1d2: –2 nt, –8 nt). Locations of primers used to amplify genomic DNA for sequencing and location of gRNA target sequence are shown. Red arrows indicate the expected cutting sites on DNA.

Fig. S4.
Fig. S4.

Western blot data show that there is no detectable NR1D1 or NR1D2 protein in Per1/2−/−;Nr1d1/2−/−;PER2-ER* cells (Left). Nuclear PER2-ER* level is induced in these cells after 4-OHT treatment (Right).

DNA containing mPER2-ER* was subcloned from pBABE-puro-mPER2-ER* into plasmid pWZL-blast to create pWZL-blast-mPER2-ER*. Cells stably expressing P-P2ER, PN-P2ER, and B-P2ER were made by retrovirus infection using pWZL-blast-mPER2-ER* (15).

Anti-mCRY1 (IgM-type monoclonal) antibodies were described previously (18). Anti-CLOCK (Bethyl Laboratories), anti-BMAL1 (Bethyl Laboratories), anti-CRY2 (Bethyl Laboratories), anti-PER2 (Alpha Diagnostic International), anti–Rev-Erbαα (NR1D1) (Cell Signaling Technology), anti–Rev-Erbβ (NR1D2) (Santa Cruz Biotechnology), anti-MEK2 (BD Biosciences), and anti-Lamin A/C (EMD Millipore) antibodies were obtained from commercial sources.

ChIP and mRNA Real-Time PCR.

ChIP and mRNA real-time PCR were performed as previous described with minimal modifications (15). Protein G Dynabeads (Thermo Fisher Scientific) were used for purification of immune complexes, and eluted DNA was purified using phenol-chloroform extraction and ethanol precipitation. Real-time PCR assays were performed using an ABI 7500 system (Applied Biosystems) using primers shown in Tables S3 and S4.

View this table:
Table S3.

Primer sets used for RT-quantitative PCR

View this table:
Table S4.

Primer sets used for ChIP-quantitative PCR

Real-time PCR assays measured expression of RNA following addition of 4-OHT. For each gene, in each cell line, the values for RNA expressed after 4-OHT addition were normalized to give the percent expression relative to the level expressed at zero time. It should be noted that at zero time, the RNA levels of some clock RNAs varied between cell lines, compared with expression of control RNA (GAPDH). Thus, in Fig. 3A, Right, although all cell lines showed a similar level of Cry1 RNA expression at zero time, in Bmal1−/− cells, Cry1 RNA was expressed constitutively at a high level.

Library Generation and Next Generation Sequencing.

ChIP-seq.

DNA libraries of input or BMAL1-ChIP from Per1/2−/−;Nr1d1/2−/−;PER2-ER* cells with or without 4-OHT treatment (4 h) were made using ThruPLEX DNA-seq Kit (RUBICON GENOMICS). Libraries were sequenced using Illumina HiSEq. 2000 (1 × 50). Two independent experiments were performed. The number of reads obtained from each sample was around 15–20 million.

RNA-seq.

Total RNA was extracted from PN-P2ER cells with or without 4-OHT treatment (4 h) using TRIzol RNA extraction (Thermo Fisher Scientific). After phase separation, RNA was purified using the PureLink RNA Mini Kit (Thermo Fisher Scientific). Libraries were made using TruSeq stranded mRNA preparation kit (Illumina) and sequenced using Illumina HiSEq. 2000 (1 × 50 × 2 lanes). The number of reads obtained from each sample was around 80∼90 million.

Next Generation Sequencing Data Analysis.

The sequences were trimmed to remove adaptor sequences using BBDuk tool (bbmap/35.82) (36) using parameters “ktrim = r k = 23 hdist = 1 minlen = 50.”

ChIP-Seq analysis.

Sequencing reads from each experiment were mapped to the mouse genome (mm10) using Bowtie2/2.2.8 (37) using default parameters. Low-quality reads were filtered using Samtools/1.3 (38) with option “-q 10.” Enrichment regions (peaks) relative to the matching input control for each experiment were determined using Callpeaks function in MACS/2015–04-20 (39) with options “-f BAM -g mm -B -q 0.01.” Given the peak regions called in each experiment, we next determined the common peak regions using Intersect function in Bedtools/2.25.0 (40). From two experiments, 4,740 common peak regions were determined in the condition without 4-OHT treatment, and 483 regions were determined in the condition with 4-OHT treatment. Because the regions determined in the condition with 4-OHT treatment were included in the regions determined in the condition without 4-OHT, we used these 4,740 regions for further quantitative analysis. To verify the specificity of our data, we performed motif analysis of these regions by MEME-ChIP (41) and the result showed high enrichment of the E-box sequence (CACGTG). To compare with random regions from the genome, we generated random regions with the same length distribution using the ShuffleBed function in Bedtools/2.25.0.

To quantify the relative binding strength of BMAL1, we first tabulated the number of unique mapping reads per region using BamtoBed and Intersect functions (with option “-c”) in Bedtools/2.25.0. Then, the numbers were normalized to number of reads per 10 million reads. We then determined the relative peak binding strength, defined as (number of reads in ChIP + c)/(number of reads in Input + c) in these regions. The constant “c” was used to smooth fluctuations in fold change due to very small denominator values, shrinking the fold change toward 1 for small counts. We chose c = 5 in this case as a conservative measure.

RNA-Seq.

Sequencing reads from each experiment were mapped to the mouse genome (mm10) using bbmap/35.82 (36) with the option “maxindel = 100000 intronlen = 10 ambig = random qin = 33.” Mapped reads were quantified using featureCounts (42) with option “-T 4 -t exon.” Gene expression analysis was performed using Deseq2 (43). To analyze CLOCK–BMAL1 controlled genes, the transcription start sites of annotated genes were compared with the BMAL1 binding sites (extended to 10 kb) identified in our ChIP-seq using Intersect functions in Bedtools/2.25.0. Data of annotated genes with a transcription start site within 5 kb from a BMAL1 binding site were selected from the Deseq2 analysis and were shown by Volcano blot.

Reporter Gene Assay and Biotin-DNA Pull Down.

The mouse Cry1 promoter (1.5 kb from –1208 to +328) was amplified by PCR from genomic DNA of PN-P2ER* cells using the primer sequences provided in a previous publication (14). The amplified DNA fragment was cloned into pGL4.16 plasmid (Promega) using NheI and XhoI. Two E-boxes identified in this region (AACGTG, CACGTG) were mutated to AAGCTG and CAGCTG (mutE). D-boxes were mutated using the published sequence (14) and plasmid mutE as template to generate the mCry1 promoter with mutDE. Mutagenesis was performed by Q5 mutagenesis kit (NEB). Reporter gene assay was performed by transfection of pGL4.16 (250 ng), pBind (50 ng), pcDNA3-mCry1 (50 ng), pcDNA3-mPER2 (100, 200, 400, and 800 ng), and pcDNA4-myc-his (to final 1,150 ng) into NIH 3T3 cells in 24-well plates using Lipofetamine 3000 (Thermo Fisher Scientific). After 24 h, reporter gene expression was analyzed using the Dual-Luciferase Reporter Assay System (Promega).

For biotin–DNA pull-down analysis, the 1.5 kbp NheI–XhoI WT and mutant Cry1 promoter fragments were first purified after NheI and EcoRV (next to XhoI) digestion of the reporter plasmids. Biotin-11-dUTP was incorporated into the DNA fragments by fill-in of the 5′-overhangs of the NheI-digested sites using exo− Klenow (NEB). DNA was purified by microSpin G50 column (GE Healthcare Life Sciences) and phenol-chloroform/ethanol precipitation to remove residual Biotin-dUTP. Binding of CLOCK–BMAL1 was performed in 30 μL reactions with 600 ng of poly-dI:C, 30 μg of NIH 3T3 nuclear extract, 0.33 M urea, 0.33% Nonidet P-40, 100 mM NaCl, and 300 fmol biotin–DNA. After 30 min incubation on ice, 20 μL of each reaction was mixed with 5 μL Dynabeads M280 streptavidin (Thermo Fisher Scientific) and incubated for another 30 min with rotation. After washing with TE buffer, BMAL1 was eluted with SDS sample buffer and detected by Western blot.

Footnotes

  • ↵1Y.-Y.C., Y.Y., and N.R. contributed equally to this work.

  • ↵2To whom correspondence should be addressed. Email: aziz_sancar{at}med.unc.edu.

  • Author contributions: Y.-Y.C., R.Y., and A.S. designed research; Y.-Y.C., Y.Y., and R.Y. performed research; Y.-Y.C., Y.Y., and N.R. analyzed data; and C.P.S. and A.S. wrote the paper.

  • Reviewers: C.H.J., Vanderbilt University; and A.C.L., University of Memphis.

  • The authors declare no conflict of interest.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1612917113/-/DCSupplemental.

References

  1. ↵
    1. Partch CL,
    2. Green CB,
    3. Takahashi JS
    (2014) Molecular architecture of the mammalian circadian clock. Trends Cell Biol 24(2):90–99.
    .
    OpenUrlCrossRefPubMed
  2. ↵
    1. Hardin PE,
    2. Panda S
    (2013) Circadian timekeeping and output mechanisms in animals. Curr Opin Neurobiol 23(5):724–731.
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Levi F,
    2. Schibler U
    (2007) Circadian rhythms: Mechanisms and therapeutic implications. Annu Rev Pharmacol Toxicol 47:593–628.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Reppert SM,
    2. Weaver DR
    (2002) Coordination of circadian timing in mammals. Nature 418(6901):935–941.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Vitaterna MH, et al.
    (1999) Differential regulation of mammalian period genes and circadian rhythmicity by cryptochromes 1 and 2. Proc Natl Acad Sci USA 96(21):12114–12119.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Shearman LP, et al.
    (2000) Interacting molecular loops in the mammalian circadian clock. Science 288(5468):1013–1019.
    .
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Kume K, et al.
    (1999) mCRY1 and mCRY2 are essential components of the negative limb of the circadian clock feedback loop. Cell 98(2):193–205.
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Zheng B, et al.
    (2001) Nonredundant roles of the mPer1 and mPer2 genes in the mammalian circadian clock. Cell 105(5):683–694.
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Zheng B, et al.
    (1999) The mPer2 gene encodes a functional component of the mammalian circadian clock. Nature 400(6740):169–173.
    .
    OpenUrlCrossRefPubMed
  10. ↵
    1. Bugge A, et al.
    (2012) Rev-erbα and Rev-erbβ coordinately protect the circadian clock and normal metabolic function. Genes Dev 26(7):657–667.
    .
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Cho H, et al.
    (2012) Regulation of circadian behaviour and metabolism by REV-ERB-α and REV-ERB-β. Nature 485(7396):123–127.
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Liu AC, et al.
    (2008) Redundant function of REV-ERBalpha and beta and non-essential role for Bmal1 cycling in transcriptional regulation of intracellular circadian rhythms. PLoS Genet 4(2):e1000023.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Preitner N, et al.
    (2002) The orphan nuclear receptor REV-ERBalpha controls circadian transcription within the positive limb of the mammalian circadian oscillator. Cell 110(2):251–260.
    .
    OpenUrlCrossRefPubMed
  14. ↵
    1. Ukai-Tadenuma M, et al.
    (2011) Delay in feedback repression by cryptochrome 1 is required for circadian clock function. Cell 144(2):268–281.
    .
    OpenUrlCrossRefPubMed
  15. ↵
    1. Ye R, et al.
    (2014) Dual modes of CLOCK:BMAL1 inhibition mediated by Cryptochrome and Period proteins in the mammalian circadian clock. Genes Dev 28(18):1989–1998.
    .
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Etchegaray JP,
    2. Lee C,
    3. Wade PA,
    4. Reppert SM
    (2003) Rhythmic histone acetylation underlies transcription in the mammalian circadian clock. Nature 421(6919):177–182.
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Xu H, et al.
    (2015) Cryptochrome 1 regulates the circadian clock through dynamic interactions with the BMAL1 C terminus. Nat Struct Mol Biol 22(6):476–484.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Nguyen KD, et al.
    (2013) Circadian gene Bmal1 regulates diurnal oscillations of Ly6C(hi) inflammatory monocytes. Science 341(6153):1483–1488.
    .
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Kondratov RV,
    2. Shamanna RK,
    3. Kondratova AA,
    4. Gorbacheva VY,
    5. Antoch MP
    (2006) Dual role of the CLOCK/BMAL1 circadian complex in transcriptional regulation. FASEB J 20(3):530–532.
    .
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Akashi M, et al.
    (2014) A positive role for PERIOD in mammalian circadian gene expression. Cell Reports 7(4):1056–1064.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Chappuis S, et al.
    (2013) Role of the circadian clock gene Per2 in adaptation to cold temperature. Mol Metab 2(3):184–193.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Hampp G, et al.
    (2008) Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood. Curr Biol 18(9):678–683.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Kaasik K,
    2. Lee CC
    (2004) Reciprocal regulation of haem biosynthesis and the circadian clock in mammals. Nature 430(6998):467–471.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Ueda HR, et al.
    (2005) System-level identification of transcriptional circuits underlying mammalian circadian clocks. Nat Genet 37(2):187–192.
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. Schmutz I,
    2. Ripperger JA,
    3. Baeriswyl-Aebischer S,
    4. Albrecht U
    (2010) The mammalian clock component PERIOD2 coordinates circadian output by interaction with nuclear receptors. Genes Dev 24(4):345–357.
    .
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Koike N, et al.
    (2012) Transcriptional architecture and chromatin landscape of the core circadian clock in mammals. Science 338(6105):349–354.
    .
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Menet JS,
    2. Rodriguez J,
    3. Abruzzi KC,
    4. Rosbash M
    (2012) Nascent-Seq reveals novel features of mouse circadian transcriptional regulation. eLife 1:e00011.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Vielhaber E,
    2. Eide E,
    3. Rivers A,
    4. Gao ZH,
    5. Virshup DM
    (2000) Nuclear entry of the circadian regulator mPER1 is controlled by mammalian casein kinase I epsilon. Mol Cell Biol 20(13):4888–4899.
    .
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Eide EJ,
    2. Vielhaber EL,
    3. Hinz WA,
    4. Virshup DM
    (2002) The circadian regulatory proteins BMAL1 and cryptochromes are substrates of casein kinase Iepsilon. J Biol Chem 277(19):17248–17254.
    .
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Siepka SM, et al.
    (2007) Circadian mutant Overtime reveals F-box protein FBXL3 regulation of cryptochrome and period gene expression. Cell 129(5):1011–1023.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Yoo SH, et al.
    (2013) Competing E3 ubiquitin ligases govern circadian periodicity by degradation of CRY in nucleus and cytoplasm. Cell 152(5):1091–1105.
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Ye R,
    2. Selby CP,
    3. Ozturk N,
    4. Annayev Y,
    5. Sancar A
    (2011) Biochemical analysis of the canonical model for the mammalian circadian clock. J Biol Chem 286(29):25891–25902.
    .
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Nangle SN, et al.
    (2014) Molecular assembly of the period-cryptochrome circadian transcriptional repressor complex. eLife 3:e03674.
    .
    OpenUrlPubMed
  34. ↵
    1. Lee C,
    2. Etchegaray JP,
    3. Cagampang FR,
    4. Loudon AS,
    5. Reppert SM
    (2001) Posttranslational mechanisms regulate the mammalian circadian clock. Cell 107(7):855–867.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Sanjana NE,
    2. Shalem O,
    3. Zhang F
    (2014) Improved vectors and genome-wide libraries for CRISPR screening. Nat Methods 11(8):783–784.
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Bushnell B
    (2014) BBMap. Available at https://sourceforge.net/projects/bbmap/.
    .
  37. ↵
    1. Langmead B,
    2. Salzberg SL
    (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9(4):357–359.
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Li H, et al., 1000 Genome Project Data Processing Subgroup
    (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25(16):2078–2079.
    .
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Zhang Y, et al.
    (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9(9):R137.
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Quinlan AR,
    2. Hall IM
    (2010) BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 26(6):841–842.
    .
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Ma W,
    2. Noble WS,
    3. Bailey TL
    (2014) Motif-based analysis of large nucleotide data sets using MEME-ChIP. Nat Protoc 9(6):1428–1450.
    .
    OpenUrlCrossRefPubMed
  42. ↵
    1. Liao Y,
    2. Smyth GK,
    3. Shi W
    (2014) featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30(7):923–930.
    .
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15(12):550.
    .
    OpenUrlCrossRefPubMed
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Molecular mechanism of mammalian circadian clock
Yi-Ying Chiou, Yanyan Yang, Naim Rashid, Rui Ye, Christopher P. Selby, Aziz Sancar
Proceedings of the National Academy of Sciences Oct 2016, 113 (41) E6072-E6079; DOI: 10.1073/pnas.1612917113
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Proceedings of the National Academy of Sciences: 116 (16)
Current Issue

Submit

Jump to section

You May Also be Interested in

PNAS Profile of NAS member Angelika Amon, winner of the 2019 Vilcek Prize in Biomedical Science
Featured Profile
PNAS Profile of NAS member Angelika Amon, winner of the 2019 Vilcek Prize in Biomedical Science
Image courtesy of Samara Vise (Koch Institute for Integrative Cancer Research, Cambridge, MA.)
Controlled burn in mixed conifer forest near Glenbrook, Lake Tahoe. Image courtesy of Alan H. Taylor.
California fires, precipitation, and jet stream
Simulations of future California climate suggest that under high CO2 conditions, by the late 21st century, increasing temperatures could lead to reduced snowpack, drier summers, and increased fire risk, independent of changes in winter precipitation.
Image courtesy of Alan H. Taylor.
Natural cast of a Eubrontes giganteus footprint from very shallow water lake sediments of the Hartford Basin at Dinosaur State Park, Rocky Hill, Connecticut. Image courtesy of Paul E. Olsen.
Mapping past Solar System dynamics
Researchers report a system for obtaining information about planetary motions approximately 200 million years ago, and the system could help constrain models of Solar System evolution and provide tests of gravitational models.
Image courtesy of Paul E. Olsen.
Prenatal testosterone exposure is tied to long-term socioeconomic traits. Image courtesy of iStock/Avril Morgan.
Twins and prenatal testosterone exposure
Compared with female co-twins not exposed to a male co-twin inside the uterus, female co-twins exposed to a male co-twin had 15.2% higher probability of dropping out of high school, 8.6% lower earnings at age 32, and 5.8% fewer children, according to a study.
Image courtesy of iStock/Avril Morgan.