JME
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Journal of Molecular Endocrinology (2004) 33, 585-608    DOI: 10.1677/jme.1.01554
© 2004 Society for Endocrinology

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary figures
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (28)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Triqueneaux, G.
Right arrow Articles by Laudet, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Triqueneaux, G.
Right arrow Articles by Laudet, V.

The orphan receptor Rev-erb{alpha} gene is a target of the circadian clock pacemaker

Gérard Triqueneaux*, Sandrine Thenot*, Tomoko Kakizawa, Marina P Antoch1, Rachid Safi, Joseph S Takahashi1, Franck Delaunay2 and Vincent Laudet

Laboratoire de Biologie Moléculaire et Cellulaire, CNRS UMR 5161, Ecole Normale Supérieur de Lyon, 46 allée d’Italie, 69364 Lyon cedex, France
1 Howard Hughes Medical Institute, Northwestern University, Department of Neurobiology and Physiology, 2205 Tech Drive, Evanston, IL 60208, USA
2 Laboratoire de Physiologie des Membranes Cellulaires, CNRS UMR 6078, Université de Nice-Sophia Antipolis, Chemin du Lazaret, 06 238 Villefranche-sur-mer, France

(Requests for offprints should be addressed to V Laudet; Email: Vincent.Laudet{at}ens-lyon.fr)

* (G Triqueneaux and S Thenot contributed equally to this work) Back


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Rev-erb{alpha} is a ubiquitously expressed orphan nuclear receptor which functions as a constitutive transcriptional repressor and is expressed in vertebrates according to a robust circadian rhythm. We report here that two Rev-erb{alpha} mRNA isoforms, namely Rev-erb{alpha}1 and Rev-erb{alpha} 2, are generated through alternative promoter usage and that both show a circadian expression pattern in an in vitro system using serum-shocked fibroblasts. Both promoter regions P1 (Rev-erb{alpha}1) and P2 (Rev-erb{alpha}2) contain several E-box DNA sequences which function as response elements for the core circadian-clock components: CLOCK and BMAL1. The CLOCK–BMAL1 heterodimer stimulates the activity of both P1 and P2 promoters in transient transfection assay by 3–6-fold. This activation was inhibited by the overexpression of CRY1, a component of the negative limb of the circadian transcriptional loop. Critical E-box elements were mapped within both promoters. This regulation is conserved in vertebrates since we found that the CLOCK–BMAL1 heterodimer also regulates the zebrafish Rev-erb{alpha} gene. In line with these data Rev-erb{alpha} circadian expression was strongly impaired in the livers of Clock mutant mice and in the pineal glands of zebrafish embryos treated with Clock and Bmal1 antisense oligonucleotides. Together these data demonstrate that CLOCK is a critical regulator of Rev-erb{alpha} circadian gene expression in evolutionarily distant vertebrates and suggest a role for Rev-erb{alpha} in the circadian clock output.


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Circadian rhythms in physiology and behavior are observed throughout the animal and plant kingdoms as well as in fungi and bacteria. They are believed to enable organisms to anticipate and adapt to rhythmic changes in their environment (Lowrey & Takahashi 2000, Young & Kay 2001, Roenneberg & Merrow 2003). These daily oscillations persist under constant conditions; that is, in the absence of external time cues. In mammals, circadian rhythms influence many if not most aspects of physiology and behavior, including sleep/wake cycles, energy metabolism, heartbeat, blood pressure and body temperature.

Recently, genetic and biochemical approaches have identified genes that contribute to the generation of circadian rhythms in mammals as well as in Drosophila, zebrafish, Arabidopsis and Neurospora (Young & Kay 2001). The current model states that, in vertebrates, a transcriptional-translational feedback loop is established by the isochronal action of transcriptional regulators such as Period 1–3 (PER 1–3) and Cryptochrome 1–2 (CRY 1–2). These proteins enter the nucleus at a specific time of the day to block the interaction between the CLOCK–BMAL1 and the MOP4–BMAL1 heterodimers and their cognate E-box element, thereby interfering with the transcription of clock-regulated genes possibly by modifying histone acetylation (Chong et al. 2000, Lowrey & Takahashi 2000, Etchegaray et al. 2003). Several recent reports revealed the importance of the regulation of the nuclear entry of PER and CRY proteins. Indeed, it has been shown that phosphorylation by casein kinase I{delta} and I{varepsilon} plays an important role in regulating the subcellular location and stability of these proteins (Lowrey & Takahashi 2000, Shearman et al. 2000, Lee et al. 2001, Vielhaber et al. 2001, Akashi et al. 2002, Yagita et al. 2002).

In mammals, the ‘master clock’ controlling circadian rhythms resides in the suprachiasmatic nuclei (SCN) of the hypothalamus (reviewed in Ripperger and Schibler 2001). To keep pace with the solar day/night cycle, the master clock can be entrained by light received through photoreceptors in the retina (Foster 1998, Abe et al. 1999). Recent evidence, however, suggests that other endogenous oscillators independent of the central pacemaker in the SCN exist and are widespread in peripheral tissues in mammals, Drosophila and zebrafish (Tosini and Menaker 1996, Plautz et al. 1997, Whitmore et al. 1998, Yamakazi et al. 2000; reviewed by Brown and Schibler 1999). The ultimate proof of the existence of peripheral oscillators has been the demonstration that circadian rhythms of clock gene expression can be induced by a serum shock administered to serum-starved immortalized rat fibroblasts in culture (Balsalobre et al. 1998, Yagita et al. 2001). The existence of these peripheral clocks can be extended to embryos since early zebrafish embryos contain a circadian clock that can be reset by light (Delaunay et al. 2000). All these data suggest that most cells, if not all, may have a circadian clock. It is believed that the SCN clock entrains the phase of peripheral clocks via chemical cues, such as rhythmically secreted hormones. Furthermore, peripheral clocks can be regulated independently from the master clock since in rodents restricted feeding during the day can completely reverse the phase of circadian oscillators in the liver but has apparently no effect on the central oscillator in the SCN (Damiola et al. 2000, Stokkan et al. 2001).

Although the core molecular pacemaker generating circadian rhythm has been defined both in the SCN and peripheral organs, the molecular outputs that ultimately regulate circadian control of cellular physiology, organ function and behavior are poorly understood (Jin et al. 1999). Specifically, the link between circadian transcriptional output and physiology under circadian control is missing. To decipher how the circadian clock is able to control output pathways it is important to identify clock-controlled genes, i.e. genes that are under the direct transcriptional control of the CLOCK–BMAL1 heterodimer.

Rev-erb{alpha} (NR1D1; Nuclear Receptors Nomenclature Committee 1999), a gene that encodes an orphan member of the nuclear receptor superfamily (Laudet & Gronnemeyer 2002) is a potential candidate as a clock-controlled gene. Indeed, we and others have shown that its expression exhibits a robust circadian rhythm in shocked fibroblasts or human hepatic cells, in rodent liver and also in zebrafish (Balsalobre et al. 1998, Delaunay et al. 2000, Torra et al. 2000, Grundschober et al. 2001). Rev-erb{alpha} appears particularly interesting as a putative clock-controlled gene because it was already known that its expression is under tight transcriptional control by a number of factors. The characterization of its promoter in humans lead to the demonstration that the expression of this gene is down-regulated by its own product, which behaves as a potent transcriptional repressor (Harding & Lazar 1995, Adelmant et al. 1996). In addition, human Rev-erb{alpha} promoter activity is inhibited by glucocorticoids (Torra et al. 2000) that are known to play an important role in the resetting of circadian time in peripheral tissues (Balsalobre et al. 2000). Finally, Rev-erb{alpha} expression is regulated by fibrates, hypolipidemic drugs that are ligands of the nuclear hormone receptor peroxisome proliferator-activated receptor {alpha} (PPAR{alpha}; Fruchart et al. 1999). Interestingly, PPAR{alpha} itself is expressed according to circadian rhythm (Lemberger et al. 1996), suggesting that the Rev-erb{alpha} gene can integrate several levels of regulation, both at the circadian and physiological levels. The link between Rev-erb{alpha} and circadian rhythm has been reinforced recently by the observation that Rev-erb{alpha}-deficient mice exhibit a circadian phenotype and that Rev-erb{alpha} controls the cyclic expression of Bmal1 (Preitner et al. 2002).

In this paper we report that the two promoters governing the expression of the Rev-erb{alpha} gene in mammals contain E-box DNA elements and generate circadian transcripts. We show that both promoters are activated by CLOCK–BMAL1 heterodimers in transient transfections and that CLOCK–BMAL1 heterodimers bind to E-box sequences and activate transcription through these E-boxes. Rev-erb{alpha} expression is strongly reduced in the livers of Clock mutant mice suggesting that the activation observed in transient transfections also occurs in vivo. In addition, we show that this regulation is evolutionarily conserved since CLOCK–BMAL1 also regulates the expression of the zebrafish Rev-erb{alpha} gene. Taken together, these data demonstrate clearly that the Rev-erb{alpha} gene is a clock-controlled gene.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Plasmid constructs and vectors

All genomic Rev-erb{alpha} regions were subcloned in pGL2 luciferase reporter plasmids after PCR with oligonucleotides containing enzymatic restriction sites. The oligonucleotides were as follows. P1, including the E1A exon as a 1350 bp fragment: 1–24 MluI, 5'-ATTAACGCGTCTGCAGGGAG CAGACCCCCCTCTA-3', top strand; 1331–1350 XhoI, 5'-ATTACTCGAGCTGTGTTGTTGTT GGAGTCTAG-3', bottom strand. P2, including exon E1B as a 1852 bp fragment: 1353–1375 KpnI, 5'-ATTAGGTACCGTACTGAGATTCT TATCTTTGCT-3', top strand; 3920–3942 XhoI, 5'-ATTACTCGAGCTAGGAAGAGAGCACGA GGGGAG-3', bottom strand. P1+P2 as a 3142 bp fragment: 1–24 MluI and 3920–3942 XhoI.

Deletion constructs in P1 were generated by PCR using the following oligonucleotides. D400-P1-Luc: 458–478 KpnI, 5'-ATTAAGGTACCCC AGGAATTCACATGCCCTTG-3', top strand, and the 1331–1350 XhoI oligonucleotide for the bottom strand. Mutations of a1 and a2 sites in the P1 promoter were produced by PCR using the following oligonucleotides. a1*D400-P1-Luc: a1*458–478 KpnI, 5'-ATTAAGGTACCCCAGG AATTCACGGGCCCTTG-3', top strand, and the 1331–1350 XhoI oligonucleotide for the bottom strand. a2*-P1-Luc: 1–24 MluI as top strand with a2*709–732 KpnI, 5'-ATTAGGTACCATTGAA TTCCAGGGAGCG-3', bottom strand. Subsequently, PCR products were ligated into P1 at the KpnI site. Each construct was sequenced on both strands.

Three deletion constructs of the P2 promoter were generated: DistalP2-Luc, D800-P2-Luc and ProximalP2-Luc. The DistalP2-Luc region was obtained by PCR with the 1353–1375 KpnI oligonucleotide, top strand, and the 2215–2190 HindIII oligonucleotide, bottom strand, 5'-ATT AAAGCTTTGCCCCTCGCACGTGGCACC-3'. The D800-P2-Luc region was obtained by PCR amplification using a 2196–3218 KpnI, 5'-ATTA AGGTACCGATCCTCGTTGGGGTGCCACG T-3'oligonucleotide, top strand, and a 3920–3942 XhoI oligonucleotide as bottom strand. The resulting PCR fragments were cloned in pGL2 vector and sequenced. The proximal P2 region was generated by enzymatic digestion of the P2-Luc vector with SmaI and PvuII before ligation.

The positive control vector (Per1E)3-Luc was generated by ligation in triplicate of the following oligonucleotide: 5'-AGATCCAAGTCCACGTG CAGGGC-3' harboring the mPer1 E-box sequence, upstream of a synthetic TATA-box oligonucleotide (5'-AGATCTGCATCGGGTATA TAATAA-3'). The whole fragment was cloned as a XhoI–HindIII fragment in the pGL2 vector. An E-box-mutated version of this oligonucleotide (5'-AGATCCAAGTCCAATTGCAGGGC-3') was used to generate the negative control vector (Per1Em)3-Luc.

Expression vectors were designed as follows: the coding region of mouse Cry1 was obtained by reverse transcriptase (RT)-PCR using the following oligonucleotides: top strand, 5'-GATCAAGCTT ACCATGGACTACAAGGACGAC-3'; bottom strand, 5'-CCAGCCTCCTTGGCCATCTTC AT-3'. The PCR product was cloned in the pCDNA3 vector and sequenced.

Expression vectors for the mouse Clock was provided by J S T and the human BMAL1 generously provided by M Ikeda (Waseda University, Tokyo, Japan). For all vectors, PCR amplifications of cDNA were performed to replace the natural ATG by a consensus Kozak sequence to ensure translation with high yield. For this, oligonucleotides containing an efficient Kozak sequence (ACC) upstream of the ATG were used as 5' primers and a 20 bp antisense oligonucleotide overlapping the stop codon as the 3' primer. Oligonucleotides are available from G T upon request. The ß-gal expression vector was constructed from the pCMV-SPORT-ßgal vector (Life Technologies, Cergy Pontoise, France) by cloning an ApaI–KpnI fragment containing the entire open reading frame into the pCDNA3 vector.

Zebrafish Rev-erb{alpha} promoter fragments were obtained by PCR amplification on a genomic DNA clone (isolated from a Danio rerio DNA library) and subcloned into a pGL3 vector. The oligonucleotides were as follows. For the 3.2 kb promoter long fragment A: top strand, 5'-GAGAGCTCGC GGCCGCGAGCTC-3'; bottom strand, 5'-TCG AAGCTTCGCACCAAATACGTGCGC-3'. The resulting PCR product was cut by SalI and HindIII restriction enzymes and subcloned into XhoI–HindIII-digested pGL3 vector. For the 1.4 kb promoter long fragment B and the 0.4 kb promoter short fragment C: top strand (B fragment), 5'-GATGAGCTCCGAGGTAAATATCACCAC-3'; top strand (C fragment), 5'-GATGAGCTCC TCTTTGACTTCGACTAC-3'; bottom strand (B and C fragments), 5'-TCGGGATCCCGCACC AAATACGTGCGC-3'. The resulting PCR products were digested by SacI and BamHI restriction enzymes and subcloned into SacI–BglII-digested pGL3 vector.

Mutated sites of E-boxes were generated using two successive PCRs with the following oligonucleotides. Top strand E-box (ZFg1), 5'-TCAGGTG GACAATTGCGCGGGGGT-3'; bottom strand E-box(ZFg1), 5'-ACCCCCGCGCAATTGTCCA CCTGA-3'; top strand E-box(Zfa), 5'-AGTCGGG TCCATAGGACACATTT-3'; bottom strand E-box(Zfa), 5'-AAATGTGTCCTATGGACCCG ACT-3'.

Accession numbers for Rev-erb{alpha} genomic sequences are: AY336123 [GenBank] , AY336124 [GenBank] , AY336125 [GenBank] and AY336126 [GenBank] for human, mouse, rat and zebrafish species respectively.

Transient transfections and reporter assays in mammalian cells

Cos1, Ros17.2/8, 3Y1 and Rat-1 cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Invitrogen) supplemented with 5% fetal calf serum and 100 U/ml penicillin/streptomycin. Typically, 40 ng of the reporter vectors were cotransfected with 100 ng of each expression vector in 12-well plates. When necessary, the final DNA concentration was adjusted to 240 ng with the pCDNA3-ßgal vector. Transfection was achieved using 2.5 x 104 cells in 0.5% fetal calf serum with 2 µl Exgen 500 transfection reagent (Euromedex, Souffelweyersheim, France) mixed with plasmids in 100 µl DMEM, according to the manufacturer’s protocol. Then, 48 h after transfection, cells were lysed in 200 µl 1 x harvest buffer (50 mM Tris, pH 7.8, 1 mM dithiothreitol, 0.1% Triton X-100) on ice. Cell lysates were vortexed briefly, and cellular debris was pelleted by centrifugation. Then 50 µl lysate were mixed with 100 µl luciferase-assay cocktail containing 1 mM luciferin. Luciferase activity was measured in a Monolight 2010 luminometer (Analytical Luminescence Laboratory, Sparks, MD, USA).

To study CRY1 inhibition, 100 ng Clock and 100 ng Bmal1 expression vectors were cotransfected with pCDNA3-Cry1 vector and, when necessary, the pCDNA3-ßgal vector, to adjust the DNA total amount.

Northern blots

Total RNA was extracted from the livers of Clock mutant or wild-type mice housed in constant darkness and killed at various circadian times. Ten µg were loaded and migrated on a 1% agarose gel. Hybridization was performed overnight in 50% formamide, 5 x SSPE, 1x Denhardt’s solution, 0.1% SDS and 0.1 mg/ml denatured salmon-sperm DNA at 63 °C. The Rev-erb{alpha} probe used was the 230 bp fragment used in the RNase protection assay. Washing was carried out in 0.5 x SSPE/0.1% SDS at 65 °C for 30 min.

Semi-quantitative RT-PCR analysis

Semi-quantitative PCR analyses were performed on RNA extracted from serum-shocked fibroblasts. RNAs were extracted using the Sigma GenElute Mammalian Total RNA Kit with 250 µl lysis buffer per well for cells cultivated in 12-well plates according to the manufacturer’s protocol. One-tenth of the extracted RNA was reverse-transcribed with avian myeloblastosis virus RT (Promega) at 42 °C for 1 h. 1–2 µl of the reverse-transcribed RNA was used for the PCR reaction with 200 ng primers, 2.5 mM Mg2+ and 1.5 U Taq Gold polymerase with appropriate buffer (Perkin-Elmer) in a final volume of 30 µl. For each gene analyzed, three different numbers of cycles were tested to reach linear phase in PCR reactions (see number of cycles on Fig. 2Go). PCR cycles were as follows: 94 °C for 5 min, 94 °C for 30 s, 57 °C for 40 s, 72 °C for 30 s and a final extension of 72 °C for 5 min. The oligonucleotides used for the PCR were as follows. Rev-erb{alpha}: E1A 5'-GGCTTCACTCGTCTCTCT CAGCC-3', top strand; E1B 5'-TGAGTCTTAT CTCCATATCACA-3', top strand; E2 5'-GCAC AGTGCCAAATGAGCGGGC-3', bottom strand. Per1: 5'-ATGACTGGGGCAGAGGTTGAGCC TG-3', top strand; 5'-TCATGCTTAGATCGTG AAATAGGG-3', bottom strand. Per3: 5'-ATGA CATACCAGGTGCCGGAGAGG-3', top strand; 5'-CTTCTGCAGTGTCACCAACTGAAC-3', bottom strand. Cry1: 5'-CAGCAAAATGGAGC CCCTGG-3', top strand; 5'-CACACCGCAGAG GACAAGCC-3', bottom strand. Clock: 5'-GCG AGAACTTGGCATTGAAGAG-3', top strand; 5'-TTTGCAGCTTGAGACATCGCTGGC-3', bottom strand. 28S: 5'-GTGAAAGCGGGGCCT CACGATCC-3', top strand; 5'-GTACTGAGCA GGATTACCATGGC-3', bottom strand.



View larger version (42K):
[in this window]
[in a new window]
 
Figure 2 Circadian expression of Rev-erb{alpha} transcripts in fibroblasts after 2 h horse serum-shock treatment. T0 corresponds to the serum-shock start. Expression was observed on RNA samples collected every 4 h with an additional point at 2 h during 48 h (A) or 24 h (B and C). (A) Expressions of Rev-erb{alpha}, Per1, Per3, Cry1 and Clock as a non-cyclic gene control in Rat-1 serum-shocked fibroblasts studied by semi-quantitative PCR. Each histogram is the average result of two experiments normalized on 28S analysis (not shown). For Cry1 two independent Q-PCRs were done on triplicates to refine the quantitation. (B, C) RT-PCR analysis of Rev-erb{alpha}1 (E1A) and Rev-erb{alpha}2 (E1B) transcripts in serum-shocked 3Y1 (B) and Rat-1 (C) fibroblasts. In each case the upper panel shows the result of a representative semi-quantitative RT-PCR assay whereas the lower curves are the averages from three Q-PCR experiments done on a 5700 GeneAmp Applied Biosystems apparatus. In Q-PCR experiments, a sharp analysis was done over the first 3 h. Note the difference of scales in the histograms (lower panels) and the number of PCR cycles (upper panels). All PCR products were obtained using an equal amount of reverse-transcribed mRNA amplified with primers located on two different exons hybridized with an internal specific probe. The signal obtained with 28S is a non-cycling control (upper B, C) and Q-PCR experiments were normalized on 28S values (lower B, C).

 
For the semi-quantification, half of the PCR products were loaded on 1% agarose gel. The gel was then incubated for 20 min in denaturation buffer (0.5 M NaOH/1.5 M NaCl) and transferred on to a nylon membrane (Hybond N; Amersham) in the same buffer for 3 h. Denatured DNA was fixed by ultraviolet treatment (254 nm; 0.36 J/cm2) with fluolink (Vilbert-Lourmat, Marne-La-Vallée, France). Hybridizations with labeled oligonucleotides specific to each PCR fragment were performed overnight at 37 °C in 4 x SSPE, 1 x Denhardt’s solution, 0.1% SDS and 0.1 mg/ml denatured salmon-sperm DNA buffer. Each hybridization was performed with 500 ng of specific oligonucleotide labeled at the 5' end using poly-nucleotide kinase in the presence of [{gamma}-32P]dATP. Washing was carried out in 0.5x SSPE/0.1% SDS at 30 or 37 °C for 20 min depending on the length of the oligonucleotide used. The oligonucleotides used as probes were all designed within the PCR fragment and were as follows. Rev-erb{alpha} in E2, 5'-CACCTACATTGGCTCCAGCGGATCCT CCC-3'; Per1, 5'-TGCTGAAGTAGAGCCTGAA GTTC-3'; Per3, 5'-ATGACATACCAGGTGCC GGAGAGG-3'; Cry1, 5'-GACGTGATAGGGAA GTGCAC-3'; Clock, 5'-CAGTTTTCAGCTCA GTTAGGAGCC-3'; 28S, 5'-GGGATAACTGG CTTGTGGCGGCCAAGCG-3'.

Serum-shock treatment

Circadian induction of 3Y1 or Rat-1 fibroblasts was performed on 7-day-confluent cells, maintained in starvation conditions in 12-well plates before being shocked. At T0, 50% horse serum was added to the cells for 2 h. The cells were then grown in DMEM without serum for 2 days and harvested at various circadian times.

Electrophoretic mobility shift assays (EMSAs)

CLOCK and BMAL1 proteins were synthesized in vitro using the TNT T7 Coupled Reticulocyte Lysate System (Promega) according to the manufacturer’s instructions. The two proteins were synthesized separately and mixed during the EMSA.

In other cases, either crude nuclear extracts of mouse fibroblast STO cells or nuclear extracts of STO cells transfected with Clock and Bmal1 expression vectors were used. All nuclear extracts were prepared using kit from Active Motif Europe (Rixensart, Belgium) according to the supplier’s protocol.

EMSAs were performed according to Vanacker et al.(1999) with the following modifications: 2 µl aliquots of each specific TNT reaction mixture were mixed with 20 fmol of acrylamide-purified, double-stranded DNA labeled probe in the presence or absence of the competitor and incubated for 30 min at 30 °C in incubation buffer (20 mM Tris-HCl, 50 mM NaCl, 5% glycerol, 5 mM magnesium sulfate, 1 µg poly(dIdC).poly-(dIdC) and 1 mM dithiothreitol). The mixture was then loaded on to a 5% polyacrylamide gel.

The oligonucleotides used for EMSA were as follows. mPer1 E-box: 5'-GATCCAGCACCCAA GTCCACGTGCAGGGATGTGTGA-3',top strand; 5'-GATCTCACACATCCCTGCACGTG GACTTGGGTGCTG-3', bottom strand. g1: 5'-GATCCCAGTTCTGCAATCACGTGAAGCTC TCACGTA-3', top strand; 5'-GATCTACGTG AGAGCTTCACGTGATTGCAGAACTGG-3', bottom strand. g2: 5'-GATCCCAGAGCCGGGC CCACGTGCTGCATTTGTTTA-3', top strand; 5'-GATCTAAACAAATGCAGCACGTGGGCC CGGCTCTGG-3', bottom strand. g3: 5'-GATC CTCGTTGGGGTGCCACGTGCGAGGGGCA CACA-3', top strand; 5'-GATCTGTGTGCCC CTCGCACGTGGCACCCCAACGAG-3', bottom strand. g4: 5'-GATCCGTGCGAGGGGCAC ACGTGGAGCGGGGACGTA-3', top strand; 5'-GATCTACGTCCCCGCTCCACGTGTGCCC CTCGCACG-3', bottom strand. g5: 5'-GATC CTGTCAGCTCCCACACGTGTCTGGGGATC CTA-3', top strand; 5'-GATCTAGGATCCCC AGACACGTGTGGGAGCTGACAG-3', bottom strand. g3+4: 5'-GATCCTCGTTGGGGTGCCA CGTGCGAGGGGCACACGTGGAGCGGGGA CGTGA-3', top strand; 5'-GTACTACGTCCC CGCTCCACGTGTGCCCCTCGCACGTGGC ACCCCAACGAG-3', bottom strand. a1: 5'-GA TCCTCCCCAGGAATTCACATGCCCTTGCC ATACA-3', top strand; 5'-GATCTGTATGGCAA GGGCATGTGAATTCCTGGGGAG-3', bottom strand. a2: 5'-GATCCGCTCCCTGGAATCAC ATGGTACCTGCTCCAA-3', top strand; 5'-GA TCTTGGAGCAGGTACCATGTGATTCCAGG GAGCG-3', bottom strand. a3: 5'-GATCCCCG GGAAGGGCTCACATGGCTGCAGAGCCGA -3', top strand; 5'-GATCTCGGCTCTGCAGCC ATGTGAGCCCTTCCCGGG-3', bottom strand.

Mutated versions of the E-box elements were obtained by substitution of two bases inside the core sequence (CACGTG to CAATTG, or CAC ATG to CCATAG) in all the relevant oligonucleotides. a2 g and Per1a substitutions were obtained with oligonucleotides bearing a CACGTG and a CACATG respectively. ZFg1, 5'-TCAGGTGGA CACGTGCGCGGGGGT-3' E-box (G), top strand; ZFg1, 5'-ACCCCCGCGCACGTGTCCA CCTGA-3' E-box (G), bottom strand. ZFg1*, 5'-TCAGGTGGACAATTGCGCGGGGGT-3' E-box (mutant), top strand; ZFg1*, 5'-ACCCCC GCGCAATTGTCCACCTGA-3' E-box (mutant), bottom strand. ZFa, 5'-AGTCGGGTCACATG GACACATTT-3' E-box (A), top strand; ZFa, 5'-AAATGTGTCCATGTGACCCGACT-3' E-box (A), bottom strand. ZFa*, 5'-AGTCGGG TCCATAGGACACATTT-3' E-box (A), top strand; ZFa*, 5'-AAATGTGTCCTATGGACCC GACT-3' E-box (A), bottom strand.

Quantitative (Q)-PCR experiments

To verify the semi-quantitative PCR and to complete the analysis of Rev-erb{alpha} isoforms, Q-PCR was performed with the ABI Prism SYBR Green Reagents (Applied Biosystems, Courtaboeuf, France). cDNAs were synthesized by reverse transcription as below. Samples contained 1x SYBR Green Master Mix, 0.5 µM primers and 1/40 synthesized cDNA in a 25 µl final volume. PCR conditions were as follows: 10 min at 95 °C, then 50 cycles of 15 s at 95 °C, 60 s at 60 °C and a final elongation cycle of 1 min at 60 °C. Absolute abundance of cDNA was calculated using the standard curve obtained with the two ranges done by serial dilutions from pure to 1/64. For the 28S RT products were diluted 1/300 before use.

Oligonucleotides used for rat samples were as follows. E1A top strand, 5'-CACGGGGCGAGA GAGGGCACC-3'; E1B top strand, 5'-TGAG TCTTATCTCCATATCACA-3'; E2 bottom strand, 5'-GAAGGGGAGCTATCATCACTG-3'; 28S top strand, 5'-GTGAAAGCGGGGCCTCA CGATCC-3'; 28S bottom strand, 5'-GTACT GAGCAGGATTACCATGGC-3'; Per1 top strand, 5'-ATGACTGGGGCAGAGGTTGAGC CTG-3'; Per1 bottom strand, 5'-TCATGCTT AGATCGTGAAATAGGG-3'. Cry1, 5'-CAGCA AAATGGAGCCCCTGG-3', top strand; 5'-CAC ACCGCAGAGGACAAGCC-3', bottom strand. Oligonucleotides used for mouse samples were as follows. E1A top strand, 5'-CAGGGGGCGAGA GAGGCCATCAC-3'; E1B top strand, 5'-AGG AAGGGGAATTCCTAAATCCC-3'; E2 bottom strand, 5'-GGCGCACTGCCAAAAGAGCGG GC-3'; 28S oligonucleotides were the same as for rat experiments.

RNA in situ hybridization

Zebrafish (Danio rerio) were kept at 28.5 °C in a 14-h light/10-h dark cycle. Embryos were collected after spawning to perform the morpholino injection. To prevent pigmentation, 0.2 mM 1-phenyl-2-thiourea (Sigma) was added to the water at 12 h post-fertilization. At 48 h post-fertilization, embryos were fixed in 4% paraformaldehyde in PBS overnight at 4 °C. Whole-mount in situ hybridization was perfomed using antisense digoxigenin-labeled Rev-erb{alpha} probe and embryos were incubated at 70 °C in a 50% formamide hybridization solution. Probes were detected using alkaline phosphatase-conjugated antibodies and visualized by 4-Nitro Blue Tetrazolium and 5-bromo-4-chloro-3-indolyl-phosphate staining (for a protocol reference, see Thisse et al. 1993).

Morpholinos

We obtained morpholinos from GeneTools (Philomath, OR, USA). The morpholinos sequences were as follows: clock1, 5'-CATCCCGG TCTATGCTGGAGGTCAT-3'; clock2, 5'-GAT AACTCGGTCTCATGGATCAGTC-3'; bmal1, 5'-TATCCATTCTTTGGTCTGCCATTAG-3'; bmal2, 5'-CAGATTTCATTTCCAGGTTGTC CAT-3'. For the control, we used standard control morpholino provided by GeneTools: 5'-CCUCUU ACCUCAGUUACAAUUUAUA-3'. We injected wild-type embryos at the one–two cell stage with 1–2 µl morpholino in 1 x Danio buffer at 0.25 or 0.5 mg/ml.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Circadian expression of the two Rev-erb{alpha} isoforms

The human and rat Rev-erb{alpha} genes generate two mRNA isoforms, Rev-erb{alpha}1 and Rev-erb{alpha}2, with different 5' regions (Lazar et al. 1989, Miyajima et al. 1989, Laudet et al. 1991, Dumas et al. 1994, G T Triqueneaux, B Staels & V L Laudet, unpublished observations; see Fig. 1Go). The resulting proteins differ only in their N-terminal A/B domain. The long isoform, Rev-erb{alpha}1, is generated from the P1 promoter located upstream of exon E1A which contains a start codon (Adelmant et al. 1996). The short isoform Rev-erb{alpha}2 is expressed at lower levels and can be generated from a newly identified promoter called P2 located upstream of the non-coding exon E1B (G T Triqueneaux, B Staels & V L Laudet, unpublished observations; see supplementary Fig. 1Go at http://jme.endocrinology-journals.org/content/vol33/issue3/). This exon is spliced to exon 2, which contains an in-frame start codon. Of note, the transcripts generated at P1 can encode the two protein isoforms by alternative start-codon usage. The relative importance of the two mechanisms (alternative promoter usage and splicing or alternative start-codon selection) used to generate Rev-erb{alpha}1 and Rev-erb{alpha}2 is, at present, unknown.



View larger version (14K):
[in this window]
[in a new window]
 
Figure 1 The Rev-erb{alpha} gene generates two N-terminal distinct isoforms, Rev-erb{alpha}1 and Rev-erb{alpha}2, from two different promoters, P1 and P2. Transcription start sites are indicated by arrows. Untranslated sequences are shown as dotted lines whereas translated ones are represented as thick lines. The start codons are indicated by circles showing that the Rev-erb{alpha}1 mRNA can generate both the Rev-erb{alpha}1 and Rev-erb{alpha}2 protein isoforms by alternative intitiator codon usage. The A/B, C and E domains are indicated. The 5' and 3' primers used to detect specifically Rev-erb{alpha}1 and Rev-erb{alpha}2 transcripts are indicated by arrowheads on E1A, E1B and E2 exons respectively.

 
Rev-erb{alpha} gene expression has been shown to be under circadian regulation in vivo and in serum-shocked tissue-culture cells (Balsalobre et al. 1998, Delaunay et al. 2000, Torra et al. 2000, Grundschober et al. 2001). We observed this rhythmic expression in serum-shocked Rat-1 or 3Y1 rat fibroblasts using a strategy that did not discriminate between the two Rev-erb{alpha} isoforms (Fig. 2AGo and results not shown). To determine whether expression of both Rev-erb{alpha}1 and Rev-erb{alpha}2 was circadian, we designed a Q-PCR-based assay using isoform-specific primers for exon 1 (Figs 1Go, 2B and 2CGo). Analysis of serum-shocked fibroblastic 3Y1 and Rat-1 cells showed that both isoforms exhibited a circadian expression. Interestingly, this experiment showed that whereas circadian expression of both transcripts was induced by serum, the initial serum response was different. Transcripts initiated at P1 were repressed whereas those transcribed from P2 were induced. This suggests that transcription of these two isoforms is differentially regulated by serum. It is important to note that, in both cell lines, there are P1-generated transcripts much more than those formed at P2. Indeed we detected in those cells only the Rev-erb{alpha}1 isoform using an antibody directed against the ligand binding domain (LBD) (from N Preitner, Department of Molecular Biology, NCCR Frontiers of Genetics, Sciences II, University of Geneva, Switzerland; see supplementary Fig. 2AGo at http://jme.endocrinology-journals.org/content/vol33/issue3/). In every Q-PCR experiment, five additional cycles were needed for P2 transcript detection compared with transcripts initiated at P1. Taken together, these data suggest that (1) the activity of both promoters is circadian, (2) the two promoters display different responses to serum and (3) the respective phase of Per transcripts and Rev-erb{alpha} transcripts is compatible with a down-regulation of Rev-erb{alpha} by the negative components of the circadian clock, PER and CRY.

Rev-erb{alpha} P1 and P2 promoters contain circadian clock-response elements

The circadian regulation of P1 and P2 prompted us to search within these promoter regions for E-box sequences, which are known response elements for the circadian transcriptional activating heterodimer CLOCK–BMAL1 (Gekakis et al. 1998). Comparative analysis of the human, rat (Rattus norvegicus) and mouse (Mus musculus) P1 promoters revealed one canonical E-box element (CACGTG; g1 in Fig. 3AGo) and two divergent E-boxes (CACATG; a1 and a2). We also found three conserved canonical E-boxes (g2, g3 and g4 in Fig. 3AGo) and one divergent E-box (a3) in the human and rat P2 region. These E-boxes were also conserved in the mouse P2 promoter (results not shown). An additional canonical E-box (g5) was found only in the rodent promoters instead of a divergent one in human (a4).



View larger version (34K):
[in this window]
[in a new window]
 
Figure 3 (A) Genomic structure of Rev-erb{alpha} promoters in rat (Rattus norvegicus) and human. P1 and P2 are the two promoters. E1A and E1B are the first alternative exons spliced to exon 2. The numbered E-boxes, g1–g5 (black boxes), correspond to classical E-boxes (CACGTG) and a1–a4 (grey boxes) to the divergent ones (CACATG). Transcription start sites are indicated by arrows. (B) Alignment of the E-boxes found in P1 and P2 rat and human promoters with E-boxes found in various circadian promoters. Numbers on the right indicates the E-boxes position in each promoter. The alignment was obtained by Seqvu software analysis (The Garvan Institute of Medical Research, Sydney, Australia). Boxed bases show conserved positions and the arrow indicates the base similarity at position -1. At -3 and +10 positions were found with a higher frequency of G and C respectively. A consensus sequence derived from this comparison is boxed on the bottom of the Figure.

 
Alignment of these E-box sequences with those shown to bind CLOCK–BMAL1 heterodimers and to confer circadian regulation of other promoters (see Fig. 3BGo) showed that Rev-erb{alpha} E-boxes contain flanking sequences that are reminiscent of those ‘circadian’ E-boxes with a pyrimidine base at position –1 and, a G at –3 and often a C at position +10 (Fig. 3BGo). Thus, Rev-erb{alpha} promoters contain all the necessary DNA elements required for a regulation by CLOCK–BMAL1 and these elements are conserved in mammals.

The CLOCK–BMAL1 heterodimer regulates Rev-erb{alpha} promoters

We next explored whether the CLOCK–BMAL1 heterodimer was able to activate P1, P2 or a construct containing P1+P2 in transient co-transfection assays. Since CLOCK–BMAL1 activity was shown to be tissue-specific in certain cases (Chen & Baler 2000), we performed these experiments in three different cell lines: Cos-1 cells, human osteosarcoma Ros17.2/8 cells and rat fibroblastic 3Y1 cells. As a positive control we used a construct in which we cloned three canonical E-boxes derived from the mouse Per1 promoter upstream of a TATA box and the luciferase gene to generate the (Per1E)3-Luc construct. The mutated version of these E-boxes from this construct give rise to (Per1Em)3-Luc, which was used as a negative control. In all three cell lines the CLOCK–BMAL1 heterodimer activated the (Per1E)3-Luc construct between 4- and 6-fold and was inactive on the mutated version (Fig. 4AGo). CLOCK alone or BMAL1 alone didn’t affect the promoter activity (supplementary Fig. 3Go). Interestingly, CLOCK–BMAL1 was able to activate 3–5-fold the activity of P1, P2 or P1+P2 in the three different cell lines (Fig. 4AGo). The association of P1 and P2 in the same vector did not confer an increased sensitivity to CLOCK–BMAL1.



View larger version (33K):
[in this window]
[in a new window]
 
Figure 4 CLOCK–BMAL1 heterodimers activate both the P1 and P2 promoters. (A) Transactivation assays were performed on Rev-erb{alpha} promoters in Cos-1, Ros17.2/8 and 3Y1 growing cells. Results are shown in fold activation and are the means±SEM from at least three independent experiments. 40 ng pGL2-luciferase reporter plasmid were used and mixed with 100 ng pCDNA3-Clock and 100 ng pCDNA3-Bmal1 or 200 ng pCDNA3-ßgal as a neutral vector. 48 h later cells were harvested and luciferase activity was analyzed. (Per1E)3-Luc and (Per1Em)3-Luc are positive and negative controls, containing wild-type or mutated E-boxes respectively in front of a minimal promoter. (B) Cry1 inhibits CLOCK–BMAL1-induced activation of the Rev-erb{alpha} promoters. The left panel shows the effect of Cry1 on the positive (Per1E)3-Luc and the negative (Per1Em)3-Luc reporter vectors. The right-hand panel shows the effect of Cry1 on P1+P2, P1 and P2 promoters. This experiment was performed on Ros17.2/8 cells but similar results were obtained on Cos-1 cells. 20 ng of pCDNA3-Cry1 were used and the total amount of DNA was adjusted to 300 ng; CB, CLOCK–BMAL1.

 
We then tested whether CRY1, a repressor of the master molecular oscillator (Kume et al. 1999), was able to inhibit this activation. The addition of low amounts of Cry1 expression vector in co-transfection experiments virtually abolished CLOCK–BMAL1-induced activation of P1+P2-, P1- or P2-driven reporter gene (Fig. 4BGo). This effect was specific as the control expression vector (ß-gal) induced only a marginal effect on CLOCK–BMAL1-mediated activation of both promoters. This specific inhibitory effect was comparable to that observed using the positive-control reporter vector (Per1E)3-Luc. These data suggest that CLOCK and BMAL1 are directly involved in the circadian regulation of both Rev-erb{alpha} isoforms.

The CLOCK–BMAL1 heterodimer binds to and transactivates Rev-erb{alpha} E-box elements

To map the regions involved in CLOCK–BMAL1 binding and trans-activation, we designed a series of P1 and P2 deletion constructs that were all tested for their ability to be activated by CLOCK–BMAL1 in transient assay in Cos-1 cells.

Deletion of the canonical E-box (g1) from P1 sequence (D400-P1-Luc construct) decreased but did not abolish its ability to be activated by CLOCK–BMAL1 (Fig. 5AGo), most likely due to the presence of two divergent CACATG E-boxes (a1 and a2). Trans-activation was not affected by introducing the mutation into a1 E-boxes (a1*D400-P1-Luc construct), strongly suggesting a critical role for the a2 site in P1 circadian regulation. This site also corresponds to the major transcriptional start site of the human P1 promoter (Adelmant et al. 1996) but to a minor one in rat (results not shown). Mutation of this site (a2*-P1-Luc construct) in the context of the complete P1 promoter totally abolished CLOCK–BMAL1 activation but not the basal activity of the promoter (Fig. 5BGo). These data indicate that CLOCK–BMAL1 activates the P1 promoter through the a2 site.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 5 (A) Mapping of the E-box mediating CLOCK–BMAL1 activation on the Rev-erb{alpha} promoters. Trans-activation assays were performed on mutated Rev-erb{alpha} promoter fragments in Cos-1 cells. Results are shown in fold activation and are the means±SEM from at least three independent experiments. 40 ng pGL2-luciferase reporter plasmid were used and mixed with 100 ng pCDNA3-Clock and 100 ng pCDNA3-Bmal1 or 200 ng pCDNA3-ß gal as a neutral vector. (Per1E)3-Luc and (Per1Em)3-Luc are positive and negative controls containing wild-type or mutated E-boxes respectively in front of a minimal promoter. Cells were harvested 48 h after transfection and luciferase activity was analyzed with a luminometer. The E-boxes discussed in the text are shown as black boxes (classical E-box CACGTG) or green boxes (divergent E-box CACATG). An asterisk indicates a mutated E-box. (B) Absolute activities of P1-Luc and a2*P1-Luc (mutated a2) promoters (grey bars) and their CLOCK–BMAL1 activation responses (blue bars) indicated the specificity of the a2 E-box element.

 
For P2, we identified a cluster of closely spaced E-boxes (three canonical, one divergent) in the distal part of the promoter and one downstream E-box in the proximal region. Deletion of the distal promoter region led to a construct (ProximalP2-Luc) containing only the latter isolated site (g5 in the rat promoter). This construct still exhibited a strong constitutive activity of the promoter (results not shown) but was not activated by CLOCK–BMAL1 (Fig. 5AGo). This suggests that the distal region of the P2 promoter is the target of CLOCK–BMAL1. This was confirmed using an internal deletion mutant lacking the P2 proximal region (DistalP2-Luc) which was fully activated by CLOCK–BMAL1. Using a new construct (D800-P2-Luc) adding g3 and g4 to the P2 proximal region we restored a partial activation (about 2-fold), indicating that they were active.

To determine whether the CLOCK and BMAL1 proteins were able to bind the identified Rev-erb{alpha} E-boxes, we set up a gel-shift assay using in vitro-translated proteins and the mouse Per1 proximal promoter E-box sequence as a probe together with competitor oligonucleotides corresponding to each of the Rev-erb{alpha} E-boxes. As expected, the CLOCK–BMAL1 heterodimer bound strongly to the mPer1 E-box probe (Fig. 6AGo). This binding is reduced by adding a molar excess of unlabeled mPer1 E-box but not by the mutated version of the Per1 E-box (Per1*). Surprisingly, g1 site, which was not sufficient for CLOCK–BMAL1 activation of the P1 promoter in our transient transfection experiments, competed strongly with CLOCK–BMAL1 binding to the mPer1 E-box. When using the a1 site-specific oligonucleotide as a probe, we observed very weak, if any, binding of CLOCK–BMAL1 consistent with the fact that mutation of this site does not modify CLOCK–BMAL1 activation of the Rev-erb{alpha} gene (Fig. 6BGo, left-hand panel). In contrast, we detected a specific, albeit weak, binding of CLOCK–BMAL1 using the a2 probe (Fig. 6BGo, middle panel). We then performed similar gel-shift assay (Fig. 6CGo) using nuclear extracts of mouse STO fibroblasts in which the CLOCK and BMAL1 proteins were transiently expressed. If BMAL1 alone only marginally affected the amount of the shifted complex, the overexpression of CLOCK increased the binding strongly, suggesting that endogenous CLOCK is a limiting factor in the crude nuclear extract. As shown in Fig. 6CGo, CLOCK–BMAL1 binding to an a2-specific probe was strongly reduced by addition of an a2 unlabeled competitor but not by an unrelated oligonucleotide or by the mutated a2 site (a2*; results not shown).



View larger version (37K):
[in this window]
[in a new window]
 
Figure 6 CLOCK–BMAL1 heterodimers bind to the E-boxes of the Rev-erb{alpha} promoters. EMSAs were performed on E-box from the mouse Per1 promoter (A) or E-boxes from the Rev-erb{alpha} promoter (B and C) as probes. For A and B, the CLOCK and BMAL1 proteins were produced independently in reticulocytes lysates whereas for C the EMSA was performed using nuclear extracts from mouse STO fibroblasts in which Clock and Bmal1 expression vectors were transiently transfected. In each case, specific complexes are indicated by black arrows whereas non-specific ones are shown by empty circles. (A) The classical E-boxes of the Rev-erb{alpha} promoters compete with the binding of CLOCK–BMAL1 to the mPer1 sequence. Lane 1, C, CLOCK protein alone; lane 2, B, BMAL1 protein alone; lanes 3–19, CB, CLOCK–BMAL1 heterodimer. Competition experiments were done with a molar excess of 10- and 100-fold of the relevant E-box as indicated by the triangles. g3+4 indicates the competition by an oligonucleotide containing both the g3 and g4 E-boxes that are adjacent in the P2 promoter. (B) Direct binding of CLOCK–BMAL1 to the divergent E-boxes. The a1, a2 and a3 E-boxes were used as probes in these EMSA. Lane 1, U, untranslated reticulocyte lysates; lanes 3 and 4, CB, CLOCK–BMAL1 heterodimer; the triangle indicates increasing amount of competitor (10- and 100-fold) which in each case is the mPer1 E-box. (C) EMSA performed with nuclear-extracts of STO fibroblasts producing BMAL1 (B, lane 2), CLOCK (CLK, lane 3) or CLOCK–BMAL1 (lanes 4–8). Nuclear extracts of cells producing the ß-galactosidase protein were used as controls (crude, lane 1). The probe used was the a2 element; the competitors added with a molar excess (10- or 100-fold, indicated by triangles) are indicated above each lane. Bottom panel, quantitation of the specific binding of the a2 probe.

 
Using a similar competition gel-shift assay we observed that all E-boxes of the P2 promoter, except the g5 site located in its distal region, were able to bind the CLOCK–BMAL1 heterodimer (Fig. 6AGo). This is consistent with our transient transfection data, which showed the importance of the distal region for CLOCK–BMAL1 activation. Since the g3 and g4 sites were adjacent, we tested these two sites together with a single oligonucleotide in a competition experiment and detected a strong competition ability of this DNA fragment in good agreement with their activity in transient transfection. As for the P1 promoter, the divergent a3 site competed with the binding of CLOCK–BMAL1 to the mPer1 E-box probe and displayed direct specific binding activity when labeled as a probe (Fig. 6BGo, right-hand panel). All these data suggest that the four E-boxes of the P2 promoter distal region are able to bind CLOCK–BMAL1 heterodimer to mediate its transcriptional activation.

Since in P1 the main E-box responsible for CLOCK–BMAL1 activation is the divergent CACATG E-box (a2), we decided to study in more detail the importance of the core sequence of the E-box compared to the adjacent sequences. As shown in Fig. 7Go, using crude nuclear extracts of STO fibroblasts we obtained a specific retarded complex at the same level as the bona fide CLOCK–BMAL1 complex with the divergent a2 E-box probe (Fig. 6CGo), as well as the canonical CACGTG E-box from Per1 promoter. The complex was competed out by an excess of specific competitors but not by unrelated sequences (see Fig. 7Go, lanes 14 and 15) or by the mutated version of the probes (a2* and Per1*; see Fig. 7Go, lanes 6, 7, 12 and 13). In all cases, it is clear that the a2 unlabeled competitor was more efficient at decreasing the binding than the Per1 unlabeled competitor (compare lanes 2 and 3 with lanes 8 and 9 on Fig. 7A and BGo). Interestingly, when we replaced the core sequence CACATG of the a2 site with a canonical CACGTG (a2 g element) we observed that the resulting sequence competed equally well for the binding to either the a2 or the Per1 probe (see Fig. 7AGo, lanes 4 and 5 and Fig. 7BGo, lanes 10 and 11). This was not the case when the Per1 canonical sequence was mutated to a CACATG-type element (Per1a) since this sequence competed poorly for the binding to the a2 probe (Fig. 7AGo, lanes 10 and 11) and only moderately for the binding to the Per1 probe (Fig. 7BGo, lanes 4 and 5). All these data clearly suggest that the context of the E-box is important for CLOCK–BMAL1 binding: the adjacent positions allow the divergent a2 element to bind CLOCK–BMAL1 whereas the context is less important for canonical E-boxes.



View larger version (39K):
[in this window]
[in a new window]
 
Figure 7 Influence of the E-box core sequence for CLOCK–BMAL1 binding. EMSAs were performed on the divergent a2 E-box from the rat P1 promoter (A) or the canonical E-box from the Per1 promoter (B) as probes with crude nuclear extracts from STO fibroblasts. In each case, competition experiments were performed with two different molar excess (10- or 100-fold, indicated by triangles) of the Per1 E-box sequence (Per1), the Per1 E-box in which the core sequence has been changed to CACATG (Per1a), the mutated Per1 E-box (CAATTG; Per1*), the rat a2 E-box (CACATG; a2), the a2 element that has been changed toward CACGTG (a2 g), the mutated a2 E-box (CCATAG; a2*) or an unrelated oligonucleotide (unr).

 
Rev-erb{alpha} is a target gene of Clock in vivo

If Rev-erb{alpha} is a target of the circadian pacemaker as suggested by our transfection and gel-shift data, then an alteration of the circadian clock function should impair Rev-erb{alpha} gene expression. To address this question, we first compared Rev-erb{alpha} circadian gene expression profiles in the livers of wild-type and Clock mutant mice. Clock mutant mice express a dominant-negative version of CLOCK, which is defective in transactivation (King et al. 1997, Gekakis et al. 1998). At the behavioral level, they exhibit lengthening of an endogenous period followed by complete arhythmicity after long exposure to constant darkness conditions (Vitaterna et al. 1994). Northern blot analysis of Rev-erb{alpha} expression levels in wild-type and Clock mutant mice kept under conditions of constant darkness was performed using total RNA samples extracted from livers collected at 4-h intervals during the third cycle of constant darkness. In the livers of wild-type mice Rev-erb{alpha} transcript demonstrates robust oscillation with a peak of expression occuring at 66 h of constant darkness, which corresponds to circadian time (CT) 6 (Fig. 8AGo). In contrast, mutant mice maintained under the same conditions showed a marked decrease of Rev-erb{alpha} expression level with significantly reduced amplitude and the phase of expression delayed by approximately 8 h. Therefore, this experiment indicates that Rev-erb{alpha} is a target gene of CLOCK in vivo.



View larger version (39K):
[in this window]
[in a new window]
 
Figure 8 (A) Northern blot analysis of Rev-erb{alpha} expression in Clock mutant mice. (A, B) The scale above the Figure indicates subjective day (hatched) or night (black), respectively. Time is indicated either as the number of hours after constant darkness (above the scale), or circadian time (CT below the scale): CT0 is the time when light would normally be switched on. Each lane was loaded with 10 µg total liver RNA. Animals were maintained in dark/dark conditions and three animals were killed at each indicated time. (A) Membranes were first probed with a Rev-erb{alpha} probe (detecting the two isoforms), exposed, washed and rehybridized with 28S probe to normalize the signal. Normalized Rev-erb{alpha} expression levels are shown as histograms at the bottom of the Figure. These histograms are the results of two independent Northern blot experiments. (B) Differential analysis of transcripts emanating from transcripts starting from P1 (E1A, upper panel) or P2 (E1B, lower panel) in wild-type or Clock mutant mice by Q-PCR experiments. Note the different scales of the histograms.

 
Since our transient transfection and gel-shift assays suggest that the two Rev-erb{alpha} promoters are direct targets of CLOCK–BMAL1, we checked whether the circadian expression of the transcripts emanating from each promoter was altered in Clock mice. Using Q-PCR analysis we observed that the transcripts starting at P1 as well as those from P2 are altered in the Clock mice suggesting that the two promoters are regulated independently by CLOCK–BMAL1 (Fig. 8BGo). These data clearly indicate that, in vivo, Rev-erb{alpha} expression is under the control of the Clock gene and that the Rev-erb{alpha} gene is a Clock-controlled gene.

The zebrafish Rev-erb{alpha} gene is also a Clock-controlled gene

We observed previously that Rev-erb{alpha} expression is circadian in zebrafish, as it is in mammals (Delaunay et al. 2000). To determine if Rev-erb{alpha} is regulated by the CLOCK–BMAL1 heterodimer in zebrafish, we isolated 3.2 kb of the 5' flanking region of the zebrafish Rev-erb{alpha} gene. Interestingly, just upstream of exon 1, we noticed a region harboring 61% sequence identity with rat P1 promoter over 165 bp (Fig. 9AGo). Of note, this region contains two E-boxes, one of which appears to be homologous to the a2 E-box. The other one is a canonical E-box that is not present in rat or human. In addition, we observed that three canonical E-boxes are present in the 5' region of the 3.2 kb P1 region (Fig. 9AGo). Between exon 1 and 2 we did not find any exon reminiscent of exon 1B that may be under the control of a specific promoter.



View larger version (45K):
[in this window]
[in a new window]
 
Figure 9 (A) Genomic structure of Rev-erb{alpha} promoter in zebrafish. The two first coding exons E1 and E2 are indicated. The numbered E-boxes ZFg1–ZFg4 (black boxes) correspond to classical E-boxes (CACGTG) and ZFa (grey boxes) to the unique divergent one (CACATG) which is homologous to the rat and human a2 element. The Rev-DR2 element is indicated by an oval. The sequence of the zebrafish proximal promoter region compared with the rat proximal P1 sequence is shown under the genomic structure. E-boxes are boxed whereas the Rev-DR2 element is underlined. The start of the rat exon 1A is emboldened. The various constructs A, B and C used in transient transfections (see panel B) are indicated. (B) CLOCK–BMAL1 heterodimers activates the zebrafish Rev-erb{alpha} promoter in COS-1 cells. The construct used, A, B, C or C mutated in the two E-boxes ZFg1 and ZFa (C2 m) are indicated below the histogram. Results are the means±SEM from at least three independent experiments. 40 ng pGL2-luciferase reporter plasmid were used and mixed with 100 ng pCDNA3-Clock and 100 ng pCDNA3-Bmal1 or 200 ng pCDNA3 as a neutral vector. Cells were harvested after 48 h and luciferase activity was analyzed. (C) CLOCK–BMAL1 binds to the zebrafish E-boxes. EMSAs were performed on the a2 E-box from the rat P1 promoter as probe with crude nuclear extract from STO fibroblasts. In each case, competition experiments were performed with two different molar excess (10- or 100-fold, indicated by triangles) of either the a2 E-box (lanes 3 and 4), the ZFg1 E-box (lanes 5 and 6), the mutated version of the ZFg1 E-box (CAATTG; ZFg1*; lanes 7 and 8), the divergent ZFa E-box (lanes 9 and 10), the mutated version of the ZFa E-box Per1 E-box sequence (CCATAG; ZFa*; lanes 11 and 12) or an unrelated oligonucleotide (unr; lanes 13 and 14). Specific complexes are indicated by an arrowhead whereas non-specific ones are shown by an empty circle. (D) Whole-mount in situ hybridization analysis of Rev-erb{alpha} expression in the zebrafish epiphysis at 48 h post-fertilization (hpf). All analyses were performed on at least 25 embryos. Wild-type and injected embryos were shown as follows: random morpholino (Mo; 0.5 mg/ml), CLOCK1+CLOCK2 morpholinos (0.25 mg/ml+0.25 mg/ml) and BMAL1+BMAL2 morpholinos (0.25 mg/ml+ 0.25 mg/ml). Rev-erb{alpha} probe was labeled with UTP-digoxigenin for overnight at 70 °C. Higher magnification at the pineal gland level is depicted below each stage in the left-hand inset. In each case the right-hand inset shows the expression of Otx5, a non-circadian gene expressed in the pineal gland, the parapineal and the retina (results not shown). For both probes, all the staining was performed under precisely identical conditions.

 
When tested in COS-1 cells, the zebrafish Rev-erb{alpha} upstream sequence (construct A) clearly exhibited a promoter activity and was activated 3–4-fold by co-transfection of the Clock and Bmal1 expression vectors (Fig. 9BGo). When we tested deletion constructs (B and C) we observed that if the three canonical E-boxes clearly play a role in CLOCK–BMAL1 activation the two downstream E-boxes were also important. Indeed, mutation of these two E-boxes (C2 m construct) abolished promoter activation (Fig. 9BGo) whereas mutation of each of them is not sufficient (results not shown). In accordance with these findings, these two sequences compete efficiently for the binding of nuclear extracts containing CLOCK–BMAL1 to the rat a2 probe (Fig. 9CGo). The three upstream E-boxes are also able to bind the CLOCK–BMAL1 heterodimer (results not shown).

To provide an independent confirmation of the relevance of these data in vivo, we used morpholinos (Nasevicius & Ekker 2000) to knock-down the expression of either CLOCK or BMAL1 proteins during zebrafish embryogenesis. Due to extensive gene duplication, the zebrafish genome contains at least two Clock genes and two Bmal1 genes (F Delaunay & V Laudet, unpublished observations), we decided to co-inject two morpholinos able to block the protein synthesis of both genes (results not shown). We checked that these morpholinos effectively block the translation of the relevant protein in an in vitro expression system. In zebrafish, Rev-erb{alpha} is expressed specifically in the pineal gland at 48 h of development (Fig. 9DGo; Delaunay et al. 2000). Interestingly, the injection of a control morpholino at the one-cell stage did not modify this expression (Fig. 9DGo). In contrast, the injection of morpholinos that blocked either Clock or Bmal1 mRNA translation strongly decreased the expression of the endogenous Rev-erb{alpha} in the pineal gland at 48 h post-fertilization. The effect of these morpholinos was specific since they did not affect the expression of Otx5, a gene that is not circadian in the pineal gland (Fig. 9DGo; Gamse et al. 2002). Taken together these data are in accordance with those obtained using the Clock mutant mice and suggest that in vivo the Clock and Bmal1 genes regulate Rev-erb{alpha} expression effectively.


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Two Rev-erb{alpha} isoforms are under circadian regulation

Rev-erb{alpha} is an orphan nuclear hormone receptor that was cloned more than 10 years ago in several mammalian species and for which no well-defined physiological role has been found (reviewed in Laudet & Gronemeyer 2000). Importantly, this receptor lacks the C-terminal AF2 domain that is required for nuclear receptor ligand-dependent transcriptional activation and consequently behaves as a constitutive repressor. This suggests that Rev-erb{alpha} may be a true orphan receptor, the activity of which might be regulated through mechanisms other than ligand (Renaud et al. 2000). In this study, we show that the Rev-erb{alpha} gene encodes two transcripts, namely Rev-erb{alpha}1 and Rev-erb{alpha}2, that differ only in their 5' region and which are both regulated in a circadian manner. Our functional data suggest that the short Rev-erb{alpha}2 isoform is generated from an internal promoter located in the first intron. This finding suggests that these two Rev-erb{alpha} isoforms could repress the transcription of target genes in a manner dependent on the time of day and independent of ligand. N-terminus isoforms are commonly found within the nuclear receptor superfamilly and they have been shown, in some instances, to result in proteins with different target gene specificity (for examples see Mulac-Jericevic et al. 2000 and Laudet & Gronemeyer 2002).

The circadian regulation of Rev-erb{alpha} could be controlled at either the transcriptional or post-transcriptional levels. However, since this rhythmic expression was observed in serum-shocked fibro-blasts, an in vitro system, which mimicks in vivo free-running conditions (Balsalobre et al. 1998, Grundschober et al. 2001), we hypothesize that it is under the transcriptional control of the circadian clock pacemaker. In addition, analysis of the rat and human Rev-erb{alpha} 5' flanking regions revealed several E-box DNA motifs, which are circadian system-response elements. Our paper presents two main lines of evidence strongly suggesting that the Rev-erb{alpha} gene is effectively under direct control of CLOCK and BMAL1: (1) functional data showing that the Rev-erb{alpha} promoters are activated by the CLOCK–BMAL1 heterodimer through binding to specific E-boxes in mammals as well as in zebrafish and (2) genetic evidence showing that Rev-erb{alpha} expression is decreased and the amplitude of its rhythmic expression severely affected in Clock mutant mice, as well as zebrafish embryos injected with CLOCK or BMAL1 morpholinos.

Rev-erb{alpha} is a target of the circadian pacemaker

The phase of Rev-erb{alpha} circadian expression both in vivo and in serum-shocked fibroblasts is consistent with control of this gene by the circadian oscillator. CLOCK and BMAL1 protein expression cycle with peak values at CT21–CT3 and a minimal value at CT6–CT12 in mouse liver (Lee et al. 2001). This expression pattern is exactly antiphasic to PER1, PER2, CRY1 and CRY2 expression as well as that of Rev-erb{alpha}, which peaks in the liver at CT8 (Balsalobre et al. 1998, Oishi et al. 1998, Torra et al. 2000). Interestingly, the pattern of Per2, Per3, Rev-erb{alpha} and Clock rhythmic expression is inverted in zebrafish when compared with nocturnal rodents but the relative phases are identical (Whitmore et al. 1998, Delaunay et al. 2000, 2003). Together, these expression data are strongly suggestive of an activation of Rev-erb{alpha} expression by the positive limb of the clock (CLOCK–BMAL1 heterodimer) and a repression by the negative limb including the PER and CRY proteins. Indeed, it has recently been shown that the phase of Rev-erb{alpha} mRNA accumulation is considerably advanced in Per2Brdm1 mutant mice (Preitner et al. 2002). This observation is further supported by the near-complete loss of Rev-erb{alpha} circadian expression in Clock mice livers and in ‘anti-CLOCK’ morpholino-treated and serum-shocked fibroblasts. These experiments provide the genetic evidence that Rev-erb{alpha} is an output of the circadian clock system. Because Rev-erb{alpha} is a transcriptional regulator, it is likely to act as a molecular link between the circadian oscillator and downstream target genes. For instance, Rev-erb{alpha} has been shown to bind to the promoter of the cellular retinol-binding protein I (CRBP1) for which mRNA circadian oscillation in the mouse liver has been described (Harding et al. 1995, Zheng et al. 2001). In addition, a direct role of Rev-erb{alpha} in the circadian pacemaker, at least in part through repression of Bmal1 gene expression, has been proposed recently (Preitner et al. 2002).

Our functional data strongly suggest that the CLOCK–BMAL1 heterodimer controls Rev-erb{alpha} expression directly. Both promoters were activated by CLOCK–BMAL1 and this activation was abolished by CRY protein. Interestingly, deletion and mutation analysis pointed to a divergent CACATG E-box in the P1 promoter, which is located in one of the transcription initiation sites of this promoter (Adelmant et al. 1996). In the P2 promoter, a complex of four E-boxes have been found (three canonical and one divergent) to respond to CLOCK–BMAL1. Gel-shift assays clearly showed that all these elements, including the divergent E-boxes, were bound by CLOCK–BMAL1. Interestingly, such divergent E-boxes have already been identified in the regulatory region of the liver transcription factor albumin-D site-binding protein (DBP) and they were shown to be recognized by CLOCK–BMAL1 (Ripperger et al. 2000). Our gel-shift experiments led us to propose that the context of these divergent E-boxes, and not only the core sequence itself, is important in determining their efficacy for CLOCK–BMAL1 binding. Indeed we show that E-boxes are less tolerant to variations in the core sequence than divergent E-boxes. These data suggest that CLOCK–BMAL1 can regulates genes that are devoid of canonical E-boxes and suggest that the promoter sequences of putative target genes of CLOCK–BMAL1 should be studied for divergent CACATG E-boxes with adjacent bases favorable for CLOCK–BMAL1 binding. In that respect, we note that all the divergent E-box of the Rev-erb{alpha} contains a T at position –1. This TCACATG motif would be interesting to use in a bioinformatic search for putative CLOCK–BMAL1 targets.

All the identified elements in the rat Rev-erb{alpha} promoter are conserved in the human and mouse sequences, stressing their functional relevance. In addition, we provide evidence suggesting that the structure of the P1 promoter and its regulation by CLOCK–BMAL1 is also conserved in the zebrafish Rev-erb{alpha} gene, suggesting that Rev-erb{alpha} may be under the control of the master oscillator in vertebrate species that have been diverged for more than 300 million years. This conservation suggests that the role of Rev-erb{alpha} in circadian regulation is under strong positive selection and thus really important for the organism’s fitness. This is in accordance with the circadian phenotype of the Rev-erb{alpha}-deficient mice. Experiments are under way in our laboratory to better delineate the role played by Rev-erb{alpha} and its paralogue Rev-erbß in circadian clocks in both mammals and zebrafish.

Although a similar response to CLOCK–BMAL1 was observed for both Rev-erb{alpha} promoters, in agreement with the rhythmicity of both Rev-erb{alpha} transcripts, the serum-shock experiments suggest that these two promoters are differentially serum-regulated. It is likely that the transient and rapid differential regulation of P1 and P2 activities found after a serum shock is not the result of a circadian regulation but rather the action of serum-induced factors such as AP1. The basis for these different early effects of serum is unknown and may result from a promoter-specific response to growth factors. Of note, both P1 and P2 contain AP1 sites conserved between human and rodent promoters and it is possible that these sites do not respond identically to the serum treatment.

Multiple interlocked loops in circadian rhythm

Interestingly, a recent report shows that Rev-erb{alpha} mRNA accumulation is considerably advanced in Per2Brdm1 mutant mice, an observation that also suggests that Rev-erb{alpha} expression is under the control of the circadian pacemaker (Preitner et al. 2002). This report also demonstrates that Rev-erb{alpha} controls the cyclic expression of BMAL1, since in Rev-erb{alpha}-knockout mice BMAL1 expression remains constant. These data, together with our demonstration that Rev-erb{alpha} is a target of the molecular oscillator, suggest that Rev-erb{alpha} is part of a regulatory loop which plays an important role in generating circadian rhythm. This situation is reminiscent of the case of DBP, which is also controlled by the clock and able to regulate Clock gene expression (Ripperger et al. 2000, Yamaguchi et al. 2000; reviewed in Roenneberg & Merrow 2003). Interestingly, in the case of Rev-erb{alpha}, as for DBP, other closely related genes also appear to play a role in circadian rhythm (Fig. 10Go). The three paralogues DBP, hepatic leukemia factor (HLF) and thyrotroph embryonic factor (TEF), which are all transcriptional activators, appear to cycle with identical phase, whereas the related gene, adenovirus E4 promoter-binding protein (E4 BP4), encoding a transcriptional repressor which does not contain the proline and acidic amino acid rich (PAR) activation domain, cycles in an opposite phase (Mitsui et al. 2001). E4 BP4 has been shown to compete with DBP, HLF and TEF for binding on the same DNA target sequences. In the case of Rev-erb{alpha}, it is interesting to observe that its paralog, Rev-erb{alpha}, also displays circadian expression, whereas two of the three closely related Retinoid-related Orphan Receptor (ROR) genes, RORß and ROR{gamma}, have been shown to cycle in brain and/or liver with a phase opposed to that of Rev-erb{alpha} (Andre et al. 1998, Panda et al. 2002, Preitner et al. 2002, Storch et al. 2002). RORs are transcriptional activators that bind to the same target sequence as Rev-erbs, namely the Rev