Bioscience Reports

Original Paper

HSG cells, a model in the submandibular clock

Yoshiaki Onishi

Abstract

Circadian rhythm of vital processes is essential to health, and various tissues show unique peripheral rhythms. HSG is the human submandibular gland cell line that has been used for analysing the effects of steroids and growth factors. In the present study, we analysed the transcriptional regulation of the BMAL1 gene, a critical component of the mammalian clock system to investigate the possibility of using HSG cells as a model system of the submandibular clock. The BMAL1 gene was expressed with circadian oscillation after stimulation with dexamethasone, and its regulatory region contained two recognition motifs for ROR (retinoic acid-receptor-related orphan receptor) and ROREs [RORα (ROR α-subunit)-binding elements] in hypomethylated CpG islands with an open chromatin structure. REV-ERBα was expressed with circadian oscillation, and knockdown experiments suggested that REV-ERBα is involved in circadian transcription of the BMAL1 gene in HSG cells. These results are similar to those in NIH 3T3 cells, a standard model for the circadian system, whereas RORα required for REV-ERBα antagonism was expressed very little in HSG cells. These findings show that in the salivary gland cell line HSG there is a rhythm in the core oscillator components BMAL1 and REV-ERBα, indicating that circadian-based transcriptional regulation can be modelled in this peripheral cell type.

  • Bmal1
  • circadian rhythm
  • HSG cell
  • Rev-erbα
  • salivary gland
  • transcription

INTRODUCTION

Circadian rhythms in behaviour and physiology are regulated by the master clock located in the SCN (suprachiasmatic nuclei) of the hypothalamus [1]. In turn, peripheral clocks directly regulate many local rhythms that probably feed back to the SCN through hypothalamic integration [2]. Circadian rhythms of vital processes are essential to health, and disturbed circadian timing comprises an increasingly relevant factor in major systemic illnesses [3] and central nervous system diseases [4].

The molecular mechanism of the circadian oscillator consists of autoregulatory transcriptional and translational feedback loops that have both positive and negative elements [1]. The key transcription factors, CLOCK and BMAL1, form heterodimers that bind to E-box enhancer sequences and activate the transcription of the Per and Cry genes. The PER and CRY proteins subsequently repress transcription at their own promoters through negative feedback caused by acting on the CLOCK–BMAL1 complex [5]. This feedback loop system controls the central clock in the SCN and the peripheral clocks in most peripheral tissues.

Since the activity of Bmal1−/− mice immediately becomes completely arrhythmic under constant darkness, BMAL1 is apparently an essential and non-redundant component of the mammalian clock [6]. The level of Bmal1 transcripts robustly oscillates in the SCN and in peripheral clock cells, and the circadian regulation of Bmal1 transcription contributes to the formation of interconnected feedback loops [5]. The Bmal1 promoter contains two recognition motifs for ROR (retinoic acid-receptor-related orphan receptor) and ROREs [RORα (RORα-subunit)-binding elements] (REV-ERB orphan nuclear receptors) [7]. REV-ERBα represses Bmal1 expression [8] and RORα activates Bmal1 transcription [9]. The opposing activities of these receptors are important in the maintenance of circadian clock function [10]. In addition, we found that the unique chromatin structure around ROREs is required for intracellular circadian transcriptional regulation of the Bmal1 gene in NIH 3T3 cells [11].

Non-invasive salivary sampling facilitates dynamic tests of hormone function and assessment of biological rhythms. This is because the amount of secreted saliva and its content of numerous enzymes, hormones, growth factors, immunoglobulins and other bioactive substances have circadian rhythms. The circadian expression of clock genes in salivary glands suggests that a peripheral clock system is active in these glands [12]. Despite the likelihood of easily obtaining information about circadian timing from saliva, the circadian clock mechanism in salivary glands has not yet been elucidated in detail.

Here, we characterized BMAL1 transcriptional regulation in HSG cells [13] and showed that there is a rhythm in the core oscillator components BMAL1 and REV-ERBα in this peripheral cell type.

MATERIALS AND METHODS

Cell culture

NIH 3T3 and HSG cells were cultured in DMEM (Dulbecco's modified Eagle's medium) supplemented with 10% (v/v) FBS (fetal bovine serum) and a mixture of 100 units/ml penicillin and 100 μg/ml streptomycin, in a humidified incubator at 37°C under a 5% CO2 atmosphere.

Real-time quantitative RT–PCR (reverse transcription–PCR)

Real-time quantitative RT–PCR was performed using LightCycler (Roche) with a LightCycler-FastStart DNA Master SYBR Green I kit (Roche) as described in [11]. The primer sequences used are as follows: Bmal1, 5′-GGCCGAATGATTGCTGAGGAAATCATGG-3′ and 5′-TTACAGCGGCCATGGCAAGTCACTAAAG-3′; Gapdh (glyceraldehyde-3-phosphate dehydrogenase), 5′-ACCACAGTCCATGCCATCAC-3′ and 5′-TCCACCACCCTGTTGCTGTA-3′; Reverbα, 5′-TGGTAGAGTTTGCCAAACACATC-3′ and 5′-TCACTGGGCGTCCACCCGGAAGGAC-3′; Rorα, 5′-TGACTCTCAGAACAACACCGTG-3′ and 5′-CGCACAATGTCTGGGTATATTGC-3′; Rorγ, 5′-CCTCATATTCCAACAACTTGGCC-3′ and 5′- TGTCTCCCTGTAGGACTTGCAGAC-3′. We constructed an authentic template comprising PCR products cloned into the pGEM-T Easy vector (Promega). Relative expression levels were evaluated using LightCycler software, version 3.5.

Reporter gene assays

Reporter gene assays were performed as described in [11]. Reporter plasmids were constructed based on pGL3-dLuc [7]. Cells were transfected using HilyMax (Dojindo) and then stimulated with dexamethasone (100 nM) for 2 h. Transient luciferase activities were assayed using a Dual Luciferase Reporter Assay System (Promega). Transcriptional activities were normalized relative to Renilla luciferase activities. Cells for real-time assays were incubated with DMEM containing 0.1 mM luciferin (Promega), 25 mM Hepes (pH 7.2) and 10% FBS, and bioluminescence was measured and integrated for 1 min at intervals of 10 min using a Kronos AB-2500 luminometer (ATTO).

CpG methylation analysis

We identified CpG islands in the BMAL1 promoter using the algorithm at http://www.urogene.org/methprimer [14]. Methylation was analysed as a bisulfite modification using EpiTect bisulfite (Qiagen) according to the manufacturer's instructions, followed by PCR cloning and sequence analysis.

ChIP (chromatin immunoprecipitation) assay

The ChIP assay was performed as described in [11]. Briefly, HSG cells expressing Myc-tagged REV-ERBα were fixed with 1% formaldehyde, and chromatin was digested with MNase (micrococcal nuclease; 1000 units/ml) for 15 min at 37°C. The Myc-tagged REV–ERBα–DNA complex was purified using the SimpleChIP enzymatic chromatin immunoprecipitation kit (Cell Signaling Technology) and anti-Myc antibody (Roche Diagnostics). Purified DNA was analysed by PCR using the primers 5′-TGGCGGGAAGAGGCAGGTC-3′ and 5′-ACACTCACCGTGGCTCGCCGGGAGC-3′.

Indirect end-labelling

The chromatin structure in vivo was determined by indirect end-labelling as described in [11]. Nuclei were isolated from HSG cells and resuspended at a DNA concentration of 1 mg/ml in 15 mM Tris (pH 7.5) containing 15 mM NaCl, 60 mM KCl, 1 mM CaCl2, 0.34 M sucrose, 15 mM 2-mercaptoethanol and 0.5 mM spermidine. MNase (1000 units/ml) was added and the reaction allowed to proceed for 10 min at 25°C. Genomic DNA was isolated and analysed by Southern blotting.

LM-PCR (ligation-mediated PCR)

LM-PCR was performed as described in [11]. Nuclei were isolated and digested with MNase as described above for indirect end-labelling. First-strand synthesis extended up to cleaved sites of genomic DNA using the first synthesis primer specific for the BMAL1 promoter (5′-GCGCCCGCACTCGGATCCCGCGG-3′). The synthesized DNA was ligated to a double-stranded linker (sense, 5′-GAATTCAGATC-3′; antisense, 5′-GCGGTGACCCGGGAGATCTGAATTC-3′) and amplified by PCR using linker antisense oligonucleotides and forward primers (5′-GAGTCCGGCGCGGGTAAACAG-3′). The amplified fragments were linearly amplified by PCR primer extension using a 32P-end-labelled extension primer (5′-GGGTAAACAGGCACCTCCGTC-3′) and the amplified fragments were resolved by electrophoresis on 8% (w/v) denaturing polyacrylamide gels.

siRNA (small interfering RNA)

Oligonucleotides of the siRNA for REV-ERBα [5′-r(CUGAAUCCCUCUAUAGUGA)d(TT)-3′ and 5′-r(UCACUAUAGAGGGAUUCAG)d(TT)-3′] and a negative control (target sequence, 5′-ATCCGCGCGATAGTACGTA-3′) (B-Bridge) were introduced into HSG cells using HilyMax (Dojindo).

RESULTS

Oscillatory BMAL1 transcription in HSG cells

We examined expression profiles of the BMAL1 gene, which is essential for the mammalian clock, in HSG cells. The BMAL1 gene was rhythmically expressed with peaks at 20 and 42 h, and a trough at 30 h in HSG cells after stimulation with dexamethasone, and the BMAL1 expression profile in HSG cells was similar to that in NIH 3T3 cells, one of the standard models for the circadian system (Figure 1).

Figure 1 Circadian expression of the BMAL1 gene in HSG cells

Transcripts of HSG and NIH 3T3 cells stimulated with 100 nM dexamethasone were analysed by real-time quantitative RT–PCR. Levels of RNA were normalized to GAPDH or Gapdh expression, and the peak value was set at 1. Values are means±S.E.M. for triplicate assays. ■, Rev-erbα; ●, Bmal1; ▲, Rorα; △, Rorγ.

Transcriptional regulation of the BMAL1 gene was analysed using transient reporter gene assays. The shortest construct (nt –201 to +24) showed promoter activity that was suppressed and enhanced at 24 and 36 h after stimulation respectively (Figure 2). This region (nt –201 to +24) corresponds to the minimal promoter region required for rhythmic expression of the Bmal1 gene in mouse NIH 3T3 cells [11], and nucleotide homology between them was 96.9%. The construct (nt –501 to –201) did not have distinct promoter activity. These results indicate that the BMAL1 gene is transcriptionally regulated in HSG cells.

Figure 2 Functional promoter region for oscillatory transcription

HSG cells were transfected with indicated BMAL1 promoter constructs and internal control plasmid and stimulated with dexamethasone; then bioluminescence was measured (open and closed bars indicate at 24 and 36 h after stimulation respectively). The arrow and open boxes in the map indicate the transcription start site (+1) and two recognition motifs for ROR and ROREs [−77 to −67 (5′-AAAGTAGGTTA-3′) and −41 to −31 (5′-AAAGTAGGTCA-3′)]. Analysed regions were −401 to +24, −313 to +24, −201 to +24 and −501 to −201. The value of the construct (−401 to +24) at 24 h after stimulation was set as 1. Values are means±S.E.M. for triplicate assays. *P<0.05 (Student's t test).

Stimulation with dexamethasone resulted in circadian oscillation of the reporter gene in HSG cells according to real-time reporter gene assays using the construct (nt −201 to +24) containing two ROREs (Figure 3). The initial response to dexamethasone on Bmal1 transcription was different between HSG and NIH 3T3 cells and the profiles of the transcription at 24 h after stimulation with dexamethasone showed a trough and a peak in HSG and NIH 3T3 cells respectively. However, the subsequent profiles after 36 h were similar between these cell lines, e.g. circadian phase and period. These lines of evidence suggested that a circadian system of BMAL1 transcriptional regulation exists in HSG cells, which is slightly different from the typical cell-autonomous clock system in NIH 3T3 cells.

Figure 3 Circadian oscillation of the BMAL1 promoter in HSG cells

HSG and NIH 3T3 cells were transfected with the BMAL1 promoter construct (−201 to +24), stimulated with dexamethasone and analysed using Kronos. Detrended results of independent experiments [dots, n=4 (HSG) or n=3 (NIH 3T3)] and the representative data (line) are shown.

Chromatin structure of the BMAL1 promoter

The methylation status of the BMAL1 promoter is involved in BMAL1 transcription [15], and the results of bisulfite genomic sequencing from four individual clones indicated that the BMAL1 promoter in HSG cells exists in the extreme hypomethylated state in CpG islands (Figure 4). The BMAL1 promoter region around the ROREs was hypersensitive to MNase, suggesting that the region comprises an open chromatin structure (bracket in Figure 5A). ChIP assays of HSG cells expressing Myc-tagged REV-REBα showed that the transcription factor interacted with the ROREs (Figure 5B). Taken together, these results suggested that the structure around the ROREs in the BMAL1 promoter of HSG cells is open, which allows transcription factor access. Our analysis of local chromatin structure around the ROREs using LM-PCR showed a footprint profile at 36 h that was similar to that at 0 h, but not that at 24 h (the region between −80 and −50 in Figure 6). These data imply that the footprint profile around the ROREs is related to the BMAL1 expression level (Figure 1). The bases marked with circles are protected more at 24 h than at 0 and 36 h. On the other hand, the base marked by a triangle (at −72) is protected at 0 and 36 h than at 24 h. These results suggest that transcription factors such as REV-ERB or ROR bind to the RORE at certain times of the circadian day. In addition, there are digested bands indicated by a triangle around the −88 region at 24 h, which was similar to those observed previously as the result of the circadian changes of the chromatin structure around the ROREs in NIH 3T3 cells [11]. These analyses of the chromatin structure indicated that the mechanism of BMAL1 transcription in HSG cells is similar to that in NIH 3T3 cells. The correlation between BMAL1 expression and the altered chromatin structure around the ROREs, particularly the region adjacent to the upstream RORE, suggests that the ROREs participate in the regulation of circadian BMAL1 expression in HSG cells through chromatin alterations.

Figure 4 Tandem ROREs in the hypomethylated CpG island in HSG cells

Examination of BMAL1 promoter sequence using the algorithm at http://www.urogene.org/methprimer revealed a prominent CpG island at −505 to +421 (grey line). The genomic sequence of HSG cells was analysed after bisulfite modification. Vertical lines indicate CpG sites in the BMAL1 promoter region. Filled and open circles indicate methylated and non-methylated CpG sites respectively. Results are shown from four individual clones.

Figure 5 Chromatin structure of the BMAL1 promoter in HSG cells

(A) Nuclease-hypersensitive region in the BMAL1 promoter. HSG cells were cultured for the indicated periods (0, 0 h; 24, 24 h; 36, 36 h) after stimulation with dexamethasone and, subsequently, isolated nuclei and naked genome (N) were incubated with MNase. Purified DNA was digested with XhoI, blotted and probed with DNA fragment (−851 to −599). The bracket indicates the region hypersensitive to MNase. Right panel: agarose gel used for Southern blotting. (B) Binding of the transcriptional regulator REV-ERBα to tandem ROREs in HSG cells. Myc–REV-ERBα was expressed in HSG cells and analysed by ChIP assays 24 h later using anti-Myc antibody. Cont, with control IgG; α-Myc, with anti-Myc antibody.

Figure 6 Footprinting in vivo around the ROREs in HSG cells

The numbers and ‘N’ (naked DNA sample) above the Figure indicate time (h) after stimulation and naked DNA respectively. The triangle and circles indicate the nucleotides showing higher or lower intensity at 24 h after stimulation with dexamethasone respectively. Numbers beside markers indicate the position from the transcription start site.

REV-ERBα is an important factor for circadian BMAL1 gene transcription in HSG cells

REV-ERBα mRNA expression in HSG cells showed a circadian rhythm, an antiphase to BMAL1 oscillation, and a similar relationship between REV-ERBα and BMAL1 expression in HSG cells to that in NIH 3T3 cells, suggesting that at least the circadian system functions in HSG cells (Figure 1). RT–PCR analysis revealed that the PCR product expression of RORα was very low in HSG cells, but obvious in NIH 3T3 cells. On the other hand, for another ROR member, RORγ, mRNA expression in HSG cells showed a circadian rhythm (Figure 1), whereas its expression level was very low, about a hundredth of that of REV-ERBα.

We performed knockdown experiments using siRNA for REV-ERBα to determine whether REV-ERBα is required for circadian BMAL1 gene regulation in HSG cells. After introducing REV-ERBα siRNA into HSG cells, the mRNA level of endogenous REV-ERBα decreased to approx. 40% of that in the cells harbouring control siRNA (Figure 7A). Real-time quantitative RT–PCR with siRNA showed that the amplitude of the transcriptional rhythm in HSG cells harbouring REV-ERBα siRNA was smaller than that of the cells harbouring control siRNA, whereas the oscillatory profile of cells harbouring control or REV-ERBα siRNA was observed (Figure 7). These findings suggested that REV-ERBα affects circadian expression of the BMAL1 gene.

Figure 7 REV-ERBα is important for BMAL1 gene oscillation

(A) RNA interference with REV-ERBα. HSG cells were transfected with siRNA (‘Cont’, 200 pmol of control siRNA; ‘100 and 200’, 100 and 200 pmol of REV-ERBα siRNA) and then transcripts were analysed by real-time quantitative RT–PCR. Levels of RNA were normalized to GAPDH expression, and the value of the control siRNA sample was set as 1. Values are means±S.E.M. for triplicate assays. (B) Function of REV-ERBα in circadian BMAL1 gene transcription. HSG cells were transfected with siRNA (closed and open circles, 100 pmol of control and REV-ERBα siRNA respectively) and stimulated with 100 nM dexamethasone; then transcripts were analysed by real-time quantitative RT–PCR. The levels of RNA were normalized to GAPDH expression, and the peak value was set as 1. Values are means±S.E.M. for triplicate assays.

DISCUSSION

The present study showed that HSG cells have an endogenous clock system that can be reset by stimulation with dexamethasone and that BMAL1, a critical component of the clock system, is transcriptionally regulated with circadian oscillation (Figures 1–3). These findings are consistent with those of previous reports showing that the structure and functions of HSG cells are similar to those of the salivary glands [13] and that a Clock- and Cry-dependent molecular clock system is active in salivary glands [12]. The regulatory region of the BMAL1 promoter for oscillatory transcription in HSG cells contains two ROREs in an open chromatin structure (Figure 5), which might be related to the intrinsic characteristics of the DNA sequence of the BMAL1 promoter such as CpG islands [16]. The unique chromatin structure around the ROREs that is hypomethylated to allow REV-ERBα access (Figures 4 and 5B) was altered in parallel with BMAL1 transcription (Figure 6). These common features of the chromatin structure are also found in NIH 3T3 cells [11], implying that the clock systems of HSG and NIH 3T3 cells are similar and cell-autonomous.

The ROREs in the Bmal1 promoter of NIH 3T3 cells are essential for the circadian oscillation of Bmal1 transcription [7], and the opposing activities of RORα and REV-ERBα are important in the maintenance of circadian clock function [10]. The circadian expression of Rev-erbα in NIH 3T3 cells peaked at 24 h and that of Rorα reached a nadir at 24 h, which is essentially antiphase to Rev-erbα expression (Figure 1). Surprisingly, very little RORα was detected in HSG cells, although the phase of REV-ERBα expression was similar to that in NIH 3T3 cells (Figure 1). The relationships between the circadian profiles of BMAL1 and REV-ERBα expression in HSG cells are consistent with those in mouse salivary glands [12], indicating functional similarity between HSG cells and salivary glands.

Knockdown of REV-ERBα damped the circadian oscillation of BMAL1 transcription (Figure 7) in accordance with the report that circadian Bmal1 expression in peripheral tissues is severely blunted in mice deficient in Rev-erbα [8]. These data suggest that REV-ERBα is involved in circadian transcription of the BMAL1 gene in HSG cells.

All REV-ERB (α and β) and ROR (α, β and γ) family members are crucial components of the molecular circadian clock, with functional differences among various peripheral tissues [17]. Both REV-ERBα and β are functionally redundant and are required for rhythmic Bmal1 transcription, whereas RORs contribute to Bma11 amplitude, but are dispensable for Bmal1 rhythms [18]. We found that RORγ (Figure 1) but not RORβ (results not shown) is expressed in HSG cells. However, the phase of the RORγ expression profile was similar to that of REV-ERBα oscillation, although the function of RORγ is thought to oppose that of REV-ERBα. It is possible that the co-ordinated expression of REV-ERBα and RORγ can generate BMAL1 oscillatory transcription in HSG cells by the variety of ligands to REV-ERBα and RORγ or by modulating the time delay in the function of REV-ERBα and RORγ. There is still the gap of the expression levels between REV-ERBα and RORγ, and the exact mechanism of the co-ordination for the circadian regulation of BMAL1 transcription remains obscure. Then, the involvement of other REV-ERB and ROR molecules in the salivary regulation of BMAL1 circadian transcription should be elucidated.

Here, we showed that circadian oscillation of the BMAL1 gene is induced by dexamethasone stimulation in HSG cells and is regulated through the ROREs in the promoter when the chromatin environment is altered, a process in which REV-ERBα plays a key role. In the salivary gland cell line HSG, there is a rhythm in the core oscillator components, indicating that circadian-based transcriptional regulation can be modelled in this peripheral cell type. Saliva might become a feasible clinical tool to determine physical time after the salivary clock mechanism is elucidated.

FUNDING

This work was supported by Kakenhi [grant number 18592055].

Acknowledgments

I am grateful to Professor H. Kizaki (Tokyo Dental College, Chiba, Japan) for providing HSG cells and to Dr K. Oishi (AIST) for helpful discussions.

Abbreviations: ChIP, chromatin immunoprecipitation; DMEM, Dulbecco's modified Eagle's medium; FBS, fetal bovine serum; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; LM-PCR, ligation-mediated PCR; MNase, micrococcal nuclease; ROR, retinoic acid-receptor-related orphan receptor; RORE, RORα (ROR α-subunit)-binding element; RT–PCR, reverse transcription–PCR; SCN, suprachiasmatic nuclei; siRNA, small interfering RNA

References

View Abstract