Circadian clocks regulate many important aspects of physiology, and their disturbance leads to various medical conditions. Circadian variations have been found in immune system variables, including daily rhythms in circulating WBC numbers and serum concentration of cytokines. However, control of immune functional responses by the circadian clock has remained relatively unexplored. In this study, we show that mouse lymph nodes exhibit rhythmic clock gene expression. T cells from lymph nodes collected over 24 h show a circadian variation in proliferation after stimulation via the TCR, which is blunted in Clock gene mutant mice. The tyrosine kinase ZAP70, which is just downstream of the TCR in the T cell activation pathway and crucial for T cell function, exhibits rhythmic protein expression. Lastly, mice immunized with OVA peptide-loaded dendritic cells in the day show a stronger specific T cell response than mice immunized at night. These data reveal circadian control of the Ag-specific immune response and a novel regulatory mode of T cell proliferation, and may provide clues for more efficient vaccination strategies.

Numerous physiological processes show cyclical variations with a period close to 24 h, known as circadian rhythms, even in the absence of environmental time cues (1). These rhythms are generated by circadian clocks located in various organs and cell types in the CNS (2) and in most peripheral tissues (3). Despite this distributed nature of the circadian system, a master clock located in the suprachiasmatic nucleus (SCN) of the hypothalamus receives input from the retina and coordinates rhythms in other tissues, thereby mediating their synchronization to the environmental cycles.

At the molecular level, circadian clocks consist of a set of clock genes (1, 4). The basic mechanism relies on the activation of the clock genes Period (Per) 1 and 2 and Cryptochrome (Cry) 1 and 2 by the CLOCK/BMAL1 dimer, and the negative feedback of PER and CRY proteins on CLOCK/BMAL1 activity, and consequently their own expression. CLOCK/BMAL1 also activate genes encoding nuclear receptors such as REV-ERB and ROR factors, which regulate Bmal1 gene transcription. These feedback loops result in rhythmic expression of RNAs encoded by clock genes as well as clock-controlled genes. Up to 10% of genes in a given tissue can be clock controlled, and the set of clock-controlled genes is tissue specific (57). Circadian clocks are present in most tissues, and, often, in organs that mediate the effects of environmental influences. Tissues that play important roles in metabolism, arousal, and stress response contain circadian oscillators.

The function of the immune system is to orchestrate a highly regulated and specific cell-mediated defense against foreign particles and pathogens. T lymphocytes are key in the adaptive response to infection. In the lymph nodes, these cells encounter and recognize, through their variable TCR, peptidic fragment of the Ag in association with self-molecules of the MHC on APCs. This recognition elicits downstream signaling that induces two crucial events, as follows: T cell proliferation, which dramatically increases the frequency of Ag-specific cells, and the differentiation into effector cells that eliminate the rapidly dividing pathogen (8, 9).

The diurnal regulation of immune processes and responses has been known for a long time (10). In rodents, susceptibility to lethal doses of the bacterial endotoxin LPS presents a diurnal rhythm (11, 12) and is exacerbated upon circadian disruption (13). Similarly, mortality following pneumococcal infection varies according to time of day (14). In humans, cases of asthma symptoms and rheumatoid arthritis occur more at certain times of the day (15, 16). However, the basis for this rhythmic regulation remains unclear. Various immune cell populations and cytokines vary over the day in the blood of humans and rodents (10, 1721). One study presented persistent rhythms in lymphocyte counts when the SCN central clock was inactivated by constant light, leading the authors to propose the existence of immune oscillators (20).

The existence of such oscillators has been supported by the expression of clock genes in blood cells of both humans (2224) and rodents (13, 25, 26) as well as in lymphoid organs such as bone marrow, thymus, spleen, and lymph nodes (LNs) (13, 25, 2731). What is the functional consequence of these immune clocks? The response of NK cells to stimulating signals has been reported to exhibit a rhythm, which depends on clock genes within these cells (18). Also, the phagocytic function of macrophages (32) and the cytokine release by these cells in response to LPS (30) present a circadian rhythm, which was proposed to depend on a clock in these cells. Many of the studies published to date have focused on the innate or non-Ag–specific immune response. Much less is known about the circadian regulation of events occurring in secondary lymphoid organs such as the LNs, and more generally of the adaptive immune response.

In this study, we show that the T cell response to a TCR trigger varies in a circadian manner. This is not based on the daily variation in the number of T cell subsets in the LN nor on the expression level of the TCR. Instead, we demonstrate a circadian variation in the protein levels of the TCR-associated kinase ZAP70. Finally, we show that the in vivo Ag-specific response of T cells is influenced by the timing of immunization.

Clock mutant mice (33) were obtained from The Jackson Laboratory and bred in house. The Clock mutant mice used in the T cell proliferation experiments had a pure C57BL/6J background, whereas those used for RNA extraction and quantitative PCR experiments had a 50:50 C57BL/6J:BALBc background. All other experiments were done with C57BL/6 mice (Charles River Laboratories). In all cases, young adult males were entrained to a 12-h light:12-h dark (LD) cycle for at least 2 wk before experimentation. For RNA, cell subpopulation, and proliferation experiments, animals were placed in constant darkness (DD) after LD entrainment and sacrificed by cervical dislocation under dim red light at the indicated times on the second day of DD. For dendritic cell (DC)–OVA immunizations, mice were kept and used under LD conditions. When animals are in LD, Zeitgeber time (ZT)0 is defined as the time of lights on and ZT12 is the time of lights off. When animals are in DD, circadian time (CT)12 is defined as the time of activity onset, equivalent to the time of lights off under the previous LD cycle. Thus, CT0-12 is subjective day and CT12-24 is subjective night. All procedures involving animals were carried out in accordance with guidelines of the Canadian Council on Animal Care and approved by the Animal Care Committee of the Douglas Mental Health University Institute.

On the second day in DD, superficial inguinal, brachial, and axillary LNs were collected every 6 h over 24 h at the time points indicated for each experiment and frozen on dry ice. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA), according to the manufacturer’s instructions. Reverse transcription was done using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA). Quantitative PCR was performed using an Applied Biosystems Real-Time Cycler ABI Prism 7500 with TaqMan gene expression assays (Per2, Mm00478113_m1; Rev-erbα, Mm00520708_m1; Zap70, Mm0049255_m1; Applied Biosystems), according to manufacturer's instructions. β-actin was used as an endogenous control (4352341E; Applied Biosystems), and relative quantification of mRNA levels was evaluated using the 2−ΔΔCT method.

For each clock gene studied by RNase protection assays (RPAs), the cDNA was generated from mouse RNAs by RT-PCR and cloned into pBluescript (Stratagene). The cDNAs cover the following sequences: Per1, nt 226–562 of GenBank entry NM_011065; Per2, nt 456–585 of AF036893; Cry1, nt 1394–1909 of AF156986; Dec1, nt 83–491 of NM_016665; Npas2, nt 1897–2344 of U77969; Bmal1, nt 360–619 of AB015203; Rev-erbα, nt 1025–1326 of NM_145434; Rev-erbβ, nt 1094–1367 of U12142; Rorγ, nt 421–850 of U43508. Antisense riboprobes were prepared by in vitro transcription of linearized templates (Riboprobe system; Promega). RPA was carried out on 6–10 μg RNA using the RPAIII kit (Ambion), according to the manufacturer’s instructions, and run on denaturing acrylamide gel. Signals were quantified by densitometry and normalized to β-actin, and the maximal value was set as 1.0.

The following Abs were used for surface stainings: TCRαβ FITC (H57-597), CD8 PerCP, or FITC (53-6.7) (BD Biosciences, Mississauga, ON, Canada); CD4 allophycocyanin RM4-5 (BioLegend, San Diego, CA); CD44 allophycocyanin (IM7) and CD62L FITC (MEL-14) (Cedarlane, Burlington, ON, Canada); and Foxp3 PE (FJK-16s) and CD25 FITC (eBiosciences, San Diego, CA). Labeled cell samples were analyzed using a BD Biosciences FACSCalibur or FACSCanto. Raw data were analyzed with FlowJo 8.7.3 (Tree Star, Ashland, OR).

On the second day in DD, inguinal, brachial, and axillary LNs were collected, mechanically dissociated, and stored on ice until manipulation. LN cells were stored on ice for a maximum of 24 h. To check that time spent on ice did not have an effect on our observations, experiments were performed with inverted experimental time points (e.g., CT6 after CT14 versus CT14 after CT6). LN cells (107/ml in PBS) were labeled with 0.5 μM CFSE (Molecular Probes, Eugene, OR) for 10 min at 37°C and washed three times in RPMI 1640 supplemented with 10% FBS. CFSE-labeled cells (2 × 106) were stimulated in 24-well plates coated with either 1 μg/ml (Figs. 2, 3, Supplemental Fig. 2) or 0.1 μg/ml (Fig. 4) anti-CD3 Ab (145-2C11; Cedarlane). In one series of experiments (Fig. 4), cells were costimulated with soluble anti-CD28 Ab (clone 37.51; Biolegend): a range of concentrations of 2–5 μg/ml anti-CD28 Ab was used in each experiment, and, in each case, we used data from the highest concentration that did not saturate the proliferative response. After 52–56 h, cells were harvested and stained with anti-CD4 and anti-CD8 Abs, followed by flow cytometry analysis. We chose to let the cells proliferate for 52–56 h because this duration consistently allows sufficient proliferation without reaching saturation. Within each experiment, cells from each time point were stimulated at the same time of day and left to proliferate for the same amount of time. Relative levels of proliferation were calculated as a ratio of the percentage of cells having undergone three divisions or more to the percentage of undivided cells. These ratios were then expressed as the percentage of mean across time points (Fig. 2) and percentage of mean across time points and conditions (Figs. 3, 4, Supplemental Fig. 2).

FIGURE 2.

LN T cell proliferative response to anti-CD3 or PMA/ionomycin stimulation. LNs were taken from mice in DD at six CTs. Cells from four mice per time point were pooled, labeled with CFSE, and plated on anti-CD3–coated wells (1 μg/ml). Cells were collected after 52–56 h in culture, and stained with anti-CD4 and anti-CD8 Abs, followed by flow cytometry analysis. A, CFSE fluorescence intensity profiles gated on CD4+ or CD8+ T cells. The experiment was done three times with similar results, and a representative experiment is shown. B and C, Averaged data of three independent experiments (including the one shown in A) for CD4+ (B) and CD8+ (C) T cells. A proliferation index was calculated as the percentage of cells that have divided three times or more, divided by the percentage of undivided cells. Data are expressed as a percentage of the mean of this index over six CTs (left panels, one-way ANOVA; CD4+, p < 0.05; CD8+, p > 0.05) or after grouping in two time bins (right panels, t test; *p < 0.05). D and E, Dual harmonic regression fit of the percentage of the mean of the data represented in B and C (CD4+; cosinor 24 h, p < 0.002, cosinor 12 h, p < 0.02; CD8+; cosinor 24 h, p < 0.01; cosinor 12 h, p < 0.05). F and G, Data for the proliferation of T cells (compilation of the same three experiments as in B) in response to PMA/ionomycin; CD4+ (F), CD8+ (G). Data were analyzed by one-way ANOVA; CD4+, p = 0.03; CD8+, p > 0.05. H and I, Dual harmonic regression fit of the percentage of the mean of the data represented in F and G (CD4+; cosinor 24 h, p < 0.01; cosinor 12 h, p > 0.05; CD8+; cosinor 24 h, p > 0.05; cosinor 12 h, p > 0.05).

FIGURE 2.

LN T cell proliferative response to anti-CD3 or PMA/ionomycin stimulation. LNs were taken from mice in DD at six CTs. Cells from four mice per time point were pooled, labeled with CFSE, and plated on anti-CD3–coated wells (1 μg/ml). Cells were collected after 52–56 h in culture, and stained with anti-CD4 and anti-CD8 Abs, followed by flow cytometry analysis. A, CFSE fluorescence intensity profiles gated on CD4+ or CD8+ T cells. The experiment was done three times with similar results, and a representative experiment is shown. B and C, Averaged data of three independent experiments (including the one shown in A) for CD4+ (B) and CD8+ (C) T cells. A proliferation index was calculated as the percentage of cells that have divided three times or more, divided by the percentage of undivided cells. Data are expressed as a percentage of the mean of this index over six CTs (left panels, one-way ANOVA; CD4+, p < 0.05; CD8+, p > 0.05) or after grouping in two time bins (right panels, t test; *p < 0.05). D and E, Dual harmonic regression fit of the percentage of the mean of the data represented in B and C (CD4+; cosinor 24 h, p < 0.002, cosinor 12 h, p < 0.02; CD8+; cosinor 24 h, p < 0.01; cosinor 12 h, p < 0.05). F and G, Data for the proliferation of T cells (compilation of the same three experiments as in B) in response to PMA/ionomycin; CD4+ (F), CD8+ (G). Data were analyzed by one-way ANOVA; CD4+, p = 0.03; CD8+, p > 0.05. H and I, Dual harmonic regression fit of the percentage of the mean of the data represented in F and G (CD4+; cosinor 24 h, p < 0.01; cosinor 12 h, p > 0.05; CD8+; cosinor 24 h, p > 0.05; cosinor 12 h, p > 0.05).

Close modal
FIGURE 3.

Levels of clock gene RNAs, T cell proliferation, and T cell subsets in LNs of Clock mutant mice. A, LNs were taken from four WT (+/+, ▪) and four Clock mutant mice (c/c, □) in DD at the indicated CTs. RNA was extracted for each mouse and pooled using equal amounts from each mouse within each time point, and quantitative PCR was performed for Rev-Erbα and Per2, relative to β-actin. In WT LN RNA, rhythmic expression of Rev-Erbα and Per2 was significant in a 24-h cosinor analysis (p < 0.002 and p < 0.001, respectively). Cosinor analyses (gray lines) for Clock mutant mice LN RNA were nonsignificant for both Rev-Erbα and Per2. B, LNs were taken from WT (+/+) and Clock mutant (c/c) mice in DD at the indicated CTs. Cells from four mice per time point were labeled with CFSE and plated in anti-CD3 Ab-coated wells (1 μg/ml). Cells were collected after 54 h in culture and stained for CD4 and CD8, followed by flow cytometry analysis. The experiment was done three times, with similar results; data were expressed as a percentage of the mean, and a composite of the three experiments is shown (mean ± SEM). Data were analyzed by two-way ANOVA with Bonferroni post hoc tests: day–night difference for WT/CD4+, *p < 0.01; and WT/CD8+, *p < 0.002; Clock mutant mice, both cell types, p > 0.05; comparison between genotypes, p > 0.05. C, For each CT and genotype, CD4/CD8 profile (left), CD44/CD62L profile gated on CD4+ or CD8+ T cells (middle), and CD25 profile gated on CD4+ T cells (right) are shown. The percentage of naive T cells (CD44lowCD62Lhigh), activated/memory T cells (CD44highCD62Llow), and CD25+ regulatory T cells is shown. The experiment was done three times, with similar results, and a representative experiment is shown.

FIGURE 3.

Levels of clock gene RNAs, T cell proliferation, and T cell subsets in LNs of Clock mutant mice. A, LNs were taken from four WT (+/+, ▪) and four Clock mutant mice (c/c, □) in DD at the indicated CTs. RNA was extracted for each mouse and pooled using equal amounts from each mouse within each time point, and quantitative PCR was performed for Rev-Erbα and Per2, relative to β-actin. In WT LN RNA, rhythmic expression of Rev-Erbα and Per2 was significant in a 24-h cosinor analysis (p < 0.002 and p < 0.001, respectively). Cosinor analyses (gray lines) for Clock mutant mice LN RNA were nonsignificant for both Rev-Erbα and Per2. B, LNs were taken from WT (+/+) and Clock mutant (c/c) mice in DD at the indicated CTs. Cells from four mice per time point were labeled with CFSE and plated in anti-CD3 Ab-coated wells (1 μg/ml). Cells were collected after 54 h in culture and stained for CD4 and CD8, followed by flow cytometry analysis. The experiment was done three times, with similar results; data were expressed as a percentage of the mean, and a composite of the three experiments is shown (mean ± SEM). Data were analyzed by two-way ANOVA with Bonferroni post hoc tests: day–night difference for WT/CD4+, *p < 0.01; and WT/CD8+, *p < 0.002; Clock mutant mice, both cell types, p > 0.05; comparison between genotypes, p > 0.05. C, For each CT and genotype, CD4/CD8 profile (left), CD44/CD62L profile gated on CD4+ or CD8+ T cells (middle), and CD25 profile gated on CD4+ T cells (right) are shown. The percentage of naive T cells (CD44lowCD62Lhigh), activated/memory T cells (CD44highCD62Llow), and CD25+ regulatory T cells is shown. The experiment was done three times, with similar results, and a representative experiment is shown.

Close modal
FIGURE 4.

Enhancing the level of costimulation in vitro with anti-CD28 has no effect on the day–night difference in T cell proliferation. LNs were taken from mice in DD at the indicated CTs. Cells from four mice per time point were pooled, labeled with CFSE, and plated in anti-CD3 Ab-coated wells (0.1 μg/ml). Cells were costimulated with soluble anti-CD28 Ab. A range of conditions from 2 to 5 μg/ml was used in every experiment, and in each case we used data from the highest concentration that did not saturate the proliferative response. Cells were collected after 56 h in culture and stained for CD4 (left panel) and CD8 (right panel), followed by flow cytometry analysis. The experiment was done three times, with similar results; data were expressed as a percentage of the mean, and a composite of the three experiments is shown (mean ± SEM). Data were analyzed by two-way ANOVA: day–night difference for CD4+ CT6 versus CT14 comparison, *p < 0.02; group difference, p > 0.05; group–time interaction, p > 0.05; day–night difference for CD8+ CT6 versus CT14 comparison, *p < 0.03; group difference, **p < 0.01; group–time interaction, **p < 0.01.

FIGURE 4.

Enhancing the level of costimulation in vitro with anti-CD28 has no effect on the day–night difference in T cell proliferation. LNs were taken from mice in DD at the indicated CTs. Cells from four mice per time point were pooled, labeled with CFSE, and plated in anti-CD3 Ab-coated wells (0.1 μg/ml). Cells were costimulated with soluble anti-CD28 Ab. A range of conditions from 2 to 5 μg/ml was used in every experiment, and in each case we used data from the highest concentration that did not saturate the proliferative response. Cells were collected after 56 h in culture and stained for CD4 (left panel) and CD8 (right panel), followed by flow cytometry analysis. The experiment was done three times, with similar results; data were expressed as a percentage of the mean, and a composite of the three experiments is shown (mean ± SEM). Data were analyzed by two-way ANOVA: day–night difference for CD4+ CT6 versus CT14 comparison, *p < 0.02; group difference, p > 0.05; group–time interaction, p > 0.05; day–night difference for CD8+ CT6 versus CT14 comparison, *p < 0.03; group difference, **p < 0.01; group–time interaction, **p < 0.01.

Close modal

On the second day in DD, superficial inguinal, brachial, and axillary LNs were collected every 6 h over 24 h, and 107 cells were washed in PBS, pelleted, frozen on dry ice, and stored at −80°C. Total proteins were extracted in lysis buffer (10 mM HEPES [pH 8], 0.1 mM EDTA [pH 8], 50 mM NaCl, 50 mM KCl, 5 mM MgCl2, 4 mM spermidine, 0.7 mM DTT, 100 μg/ml BSA, 17% v/v glycerol, 0.5 mM PMSF, 0.1% Nonidet P-40, and Complete Mini protease inhibitor [Roche]). Extracts were run on SDS-PAGE and immunoblotted with anti–β-actin (Sigma-Aldrich), anti-protein kinase C (PKC) A-9 (Santa Cruz Biotechnology, Santa Cruz, CA), anti–lymphocyte-specific protein tyrosine kinase (LCK), clone 3A5 (Santa Cruz Biotechnology), or anti-ZAP70 clone 1E7.2 (Upstate, Lake Placid, NY), followed by anti-mouse HRP-conjugated Ab (Sigma-Aldrich), and finally ECL detection reagent (Perkin Elmer, Waltham, MA). ZAP70 levels were quantified by densitometry, normalized to actin, and calculated as percentage of mean across 24 h, and the data from four independent experiments were then averaged.

Bone marrow cells (2 × 105 cells/ml) were plated on 6-well plates in complete RPMI 1640, and GM-CSF (500 U/ml; Invitrogen) and IL-4 (250 U/ml; made in house) were added on days 0, 2, 3, and 6 of culture. To induce maturation of DCs, 1 μg/ml LPS was added on day 6. Also on day 6, cells were pulsed overnight with 2 μg/ml OVA257–264 (SIINFEKL) peptide (Midwest Biotech, Fishers, IN). Cultures were underlaid on 14.7% Histodenz (Sigma-Aldrich) and centrifuged at 1200 × g for 20 min. DC-OVAs were collected from the interface and washed three times in PBS before injection. Three mice per time point were injected i.v. with either 1.25 × 106 DC-OVAs or sterile PBS. Injections at ZT18 were performed in dim red light. Cytotoxic T cell response to i.v. immunization is initiated in the spleen. Accordingly, following this type of immunization, we cannot detect any T cell expansion in the LNs (data not shown). Thus, after 7 d, the mice were sacrificed at the times indicated, and spleens were harvested and stored on ice until time of staining. Total spleen cell number and percentage of CD8+ T cells in spleen did not vary between the day and night time points (data not shown). Splenocytes were stained with Kb-OVA tetramers (made in-house) and anti-CD3, anti-CD8, and anti-CD44 Abs, followed by flow cytometry analysis.

For restimulation and assessment of IFN-γ production, splenocytes from each mouse were cultured with 1 μg/ml OVA peptide in complete RPMI 1640 for 6 h at 37°C. For the last 2 h, 10 μg brefeldin A (Sigma-Aldrich) per ml of cells was added. Cells were fixed and then stained with anti–IFN-γ Ab (clone XMG1.2; Invitrogen) or isotype control diluted in 0.5% saponin (Sigma-Aldrich). Cells were washed, cell surface stained with anti-CD8 and anti-CD44 Abs, and analyzed by flow cytometry.

All 24 h time courses (Figs. 1A, 2, 5A, 5C, 5D, Supplemental Fig. 1C) were analyzed by one-way ANOVA (effect of time) with Bonferroni post hoc test where applicable. Proliferation assays at two time points with two conditions (Figs. 3B, 4, Supplemental Fig. 2) and clock gene expression in Clock mutant and wild-type (WT) mice were analyzed by two-way ANOVA (for effect of time and group) with Bonferroni post hoc tests. Data of DC–OVA immunizations (Fig. 6) were analyzed by two-tailed paired t test. Differences were considered to be significant if p < 0.05.

FIGURE 1.

Clock genes are rhythmically expressed in LNs and show similar expression in inguinal and brachial LNs. A, LNs were taken from mice in DD at the indicated CT. Rev-Erbα and Per2 RNA levels were analyzed in samples from individual mice by quantitative PCR, relative to β-actin. n = 4 mice per CT. Analysis by one-way ANOVA revealed a significant effect of time for Rev-Erbα (p < 0.01; Bonferroni post hoc tests for CT2 versus CT8 and CT8 versus CT14, p < 0.05; CT8 versus CT20, p < 0.001) and for Per2 (p < 0.01; Bonferroni post hoc tests for CT14 versus all other CTs, p < 0.05). Cosinor analysis (gray lines) revealed a significant 24-h rhythm for both Rev-Erbα (p < 0.0001) and Per2 (p < 0.0001). The experiment was done twice, with similar results. B, Inguinal (▪) and brachial (△) LNs were collected, as in A, and used to prepare RNA. Equal amounts of RNA from six mice were pooled for each time point and analyzed by RPA for the indicated clock genes. RNA curves were quantified by densitometry and normalized to β-actin.

FIGURE 1.

Clock genes are rhythmically expressed in LNs and show similar expression in inguinal and brachial LNs. A, LNs were taken from mice in DD at the indicated CT. Rev-Erbα and Per2 RNA levels were analyzed in samples from individual mice by quantitative PCR, relative to β-actin. n = 4 mice per CT. Analysis by one-way ANOVA revealed a significant effect of time for Rev-Erbα (p < 0.01; Bonferroni post hoc tests for CT2 versus CT8 and CT8 versus CT14, p < 0.05; CT8 versus CT20, p < 0.001) and for Per2 (p < 0.01; Bonferroni post hoc tests for CT14 versus all other CTs, p < 0.05). Cosinor analysis (gray lines) revealed a significant 24-h rhythm for both Rev-Erbα (p < 0.0001) and Per2 (p < 0.0001). The experiment was done twice, with similar results. B, Inguinal (▪) and brachial (△) LNs were collected, as in A, and used to prepare RNA. Equal amounts of RNA from six mice were pooled for each time point and analyzed by RPA for the indicated clock genes. RNA curves were quantified by densitometry and normalized to β-actin.

Close modal
FIGURE 5.

Expression of TCR and downstream signaling molecules over the day in LNs. LNs were taken from WT mice in DD at the indicated CTs, and then processed for surface staining, Western blotting, or quantitative PCR. A, Cells from four individual mice per time point were stained using anti-TCRαβ Ab, followed by flow cytometry. The mean fluorescence intensity for lymphocytes is plotted with each point representing a single mouse. Analysis by one-way ANOVA showed no significant daily variation of the expression level of TCRαβ (p > 0.05). B, Protein extracts pooled from four mice were resolved by SDS-PAGE and immunoblotted with Abs for LCK, PKC, ZAP70, and ACTIN. The experiment was done four times, with similar results, and a representative blot is shown. C, Quantification of ZAP70 protein expression from four experiments (the one shown in B and three other experiments), with similar results; data were expressed as a percentage of the mean, and a composite of the four experiments is shown. Analysis by one-way ANOVA revealed a significant effect of time (p < 0.05; Bonferroni post hoc test for CT8 versus CT20, p < 0.05). Cosinor analysis (gray line) revealed a significant 24-h rhythm (p < 0.0001). D, Quantitative PCR was done for Zap70, relative to β-actin, on RNA extracted from LNs of four individual mice per CT. One-way ANOVA showed no significant effect of time (p = 0.17). Cosinor analysis (gray line) revealed a significant 24-h rhythm (p < 0.0001). The experiment was done twice, with similar results (but with significant effect of time in ANOVA for the other experiment).

FIGURE 5.

Expression of TCR and downstream signaling molecules over the day in LNs. LNs were taken from WT mice in DD at the indicated CTs, and then processed for surface staining, Western blotting, or quantitative PCR. A, Cells from four individual mice per time point were stained using anti-TCRαβ Ab, followed by flow cytometry. The mean fluorescence intensity for lymphocytes is plotted with each point representing a single mouse. Analysis by one-way ANOVA showed no significant daily variation of the expression level of TCRαβ (p > 0.05). B, Protein extracts pooled from four mice were resolved by SDS-PAGE and immunoblotted with Abs for LCK, PKC, ZAP70, and ACTIN. The experiment was done four times, with similar results, and a representative blot is shown. C, Quantification of ZAP70 protein expression from four experiments (the one shown in B and three other experiments), with similar results; data were expressed as a percentage of the mean, and a composite of the four experiments is shown. Analysis by one-way ANOVA revealed a significant effect of time (p < 0.05; Bonferroni post hoc test for CT8 versus CT20, p < 0.05). Cosinor analysis (gray line) revealed a significant 24-h rhythm (p < 0.0001). D, Quantitative PCR was done for Zap70, relative to β-actin, on RNA extracted from LNs of four individual mice per CT. One-way ANOVA showed no significant effect of time (p = 0.17). Cosinor analysis (gray line) revealed a significant 24-h rhythm (p < 0.0001). The experiment was done twice, with similar results (but with significant effect of time in ANOVA for the other experiment).

Close modal
FIGURE 6.

Day–night difference in T cell activation in response to Ag presentation in vivo. Bone marrow DCs loaded with OVA peptide were injected i.v. into three mice at either ZT6 or ZT18 (mid day and mid night, respectively). Seven days after injection, spleens were harvested and processed for surface labeling (A–C) or in vitro restimulation (D, E). A and B, Splenocytes of mice injected at ZT6 or 18, and dissected at the same respective times 1 wk later, were stained with anti-CD3, anti-CD8, anti-CD44 Abs, and Kb-OVA tetramer, and analyzed by flow cytometry. A, Shows representative FACS plots gated on CD8+CD3+, whereas B shows the data of one experiment (left; n = 3 per group and time point; paired t test, *p < 0.01), where data are represented relative to control mice injected with vehicle (PBS) and the ZT6/ZT18 ratio for each of three independent experiments (right; mean ± SEM). C, Experiment as in B, except that mice injected at ZT6 were sacrificed at ZT18, and vice versa, 7 d plus 12 h in both cases. The experiment was done twice, with similar results, and a representative experiment is shown (n = 3 per group and time point; paired t test, *p < 0.01). D and E, Splenocytes of mice injected at ZT6 or 18, and dissected at the same respective times 1 wk later, were put in culture and restimulated with OVA peptide for 6 h. Cells were then stained for CD8, CD44, and IFN-γ and analyzed by flow cytometry. D, Shows representative FACS plots, whereas E shows a representative experiment (n = 3 per group and time point; paired t test, *p < 0.01; done twice with similar results).

FIGURE 6.

Day–night difference in T cell activation in response to Ag presentation in vivo. Bone marrow DCs loaded with OVA peptide were injected i.v. into three mice at either ZT6 or ZT18 (mid day and mid night, respectively). Seven days after injection, spleens were harvested and processed for surface labeling (A–C) or in vitro restimulation (D, E). A and B, Splenocytes of mice injected at ZT6 or 18, and dissected at the same respective times 1 wk later, were stained with anti-CD3, anti-CD8, anti-CD44 Abs, and Kb-OVA tetramer, and analyzed by flow cytometry. A, Shows representative FACS plots gated on CD8+CD3+, whereas B shows the data of one experiment (left; n = 3 per group and time point; paired t test, *p < 0.01), where data are represented relative to control mice injected with vehicle (PBS) and the ZT6/ZT18 ratio for each of three independent experiments (right; mean ± SEM). C, Experiment as in B, except that mice injected at ZT6 were sacrificed at ZT18, and vice versa, 7 d plus 12 h in both cases. The experiment was done twice, with similar results, and a representative experiment is shown (n = 3 per group and time point; paired t test, *p < 0.01). D and E, Splenocytes of mice injected at ZT6 or 18, and dissected at the same respective times 1 wk later, were put in culture and restimulated with OVA peptide for 6 h. Cells were then stained for CD8, CD44, and IFN-γ and analyzed by flow cytometry. D, Shows representative FACS plots, whereas E shows a representative experiment (n = 3 per group and time point; paired t test, *p < 0.01; done twice with similar results).

Close modal

All data of 24 h time courses were also analyzed to evaluate the circadian rhythmicity using a cosinor fit. In all cases, a harmonic regression model was applied to the data using nlmixed SAS procedure (SAS Institute, Cary, NC). In the case of the proliferation data (Fig. 2), the fitting model used both a 24-h period and a 12-h period, as follows:

where Yi is the ith smoothed value of proliferation for a given T cell subset at time ti, and μ is the mean proliferation per subset. Circadian amplitude (Ai) is defined as the difference between the 24-h mean and the maximum along the cosinor curve, and ϕi is the acrophase in clock hours.

For all other time courses, the fitting model used a 24-h period only, as follows:

where Yi is the ith smoothed value of RNA or protein levels at time ti, and μ is the mean levels across the time course. Circadian amplitude (Ai) is defined as the difference between the 24-h mean and the maximum along the cosinor curve, and ϕi is the acrophase in clock hours. Residual error is assumed to be independent and to have a normal distribution εi ∼ N(0,σ2). Each parameter was tested for significance with p < 0.05. The T statistic was used to evaluate whether the Y intercept, the amplitude, and the phase of the 24-h rhythm differed from zero.

As a first step in studying circadian processes in the LNs, we evaluated clock gene expression. Mice were entrained on a LD cycle and then put in DD. Superficial LNs were collected every 6 h over 24 h on the second day in DD, and RNA was extracted. Quantitative PCR revealed that the clock genes Rev-Erbα and Per2 are rhythmically expressed in LNs (Fig. 1A), and that the phase of the rhythms is similar to that in other peripheral tissues (3, 26, 34). In addition to Rev-Erbα and Per2, RPAs also showed expression of the clock genes Per1, Cry1, Rev-Erbβ, Bmal1, Rorγ, Npas2, and Dec1, in a circadian pattern very similar for both inguinal and brachial LNs (Fig. 1B). This reveals that the molecular clockwork is present in the LNs, suggesting circadian regulation of events in these organs.

Because previous research has shown rhythms of leukocyte populations in the blood of rodents and humans, we surveyed the relative amounts of different T cell subsets and B cells in LNs sampled from WT mice under DD. There were no variations across the day in the abundance of B cells, CD4+ or CD8+, naive, effector, activated, or regulatory T cells (Table I, Supplemental Fig. 1).

Table I.
Lymphocyte subsets do not vary over the day in mouse LNs
CT2 (×105 Cells)CT8 (×105 Cells)CT14 (×105 Cells)CT20 (×105 Cells)p
CD4+ T cells Total 10.5 ± 0.3 11.7 ± 0.3 11.4 ± 0.3 10.8 ± 0.2 >0.05 
 Effector 1.2 ± 0.1 0.9 ± 0.1 1 ± 0.1 0.8 ± 0.1 >0.05 
 Naive 7.4 ± 0.1 8.7 ± 0.1 8.5 ± 0.1 7.9 ± 0.1 >0.05 
 Regulatory 0.2 ± 0.02 0.1 ± 0.01 0.2 ± 0.01 0.1 ± 0.01 >0.05 
CD8+ T cells Total 8.3 ± 0.1 9.7 ± 0.1 9.9 ± 0.1 8.7 ± 0.1 >0.05 
 Effector 0.5 ± 0.02 0.3 ± 0.01 0.4 ± 0.02 0.4 ± 0.01 >0.05 
 Naive 5.9 ± 0.1 7.0 ± 0.1 7.1 ± 0.1 6.1 ± 0.2 >0.05 
B cells  11.9 ± 0.1 11.8 ± 0.2 11.9 ± 0.2 12.1 ± 0.1 >0.05 
CT2 (×105 Cells)CT8 (×105 Cells)CT14 (×105 Cells)CT20 (×105 Cells)p
CD4+ T cells Total 10.5 ± 0.3 11.7 ± 0.3 11.4 ± 0.3 10.8 ± 0.2 >0.05 
 Effector 1.2 ± 0.1 0.9 ± 0.1 1 ± 0.1 0.8 ± 0.1 >0.05 
 Naive 7.4 ± 0.1 8.7 ± 0.1 8.5 ± 0.1 7.9 ± 0.1 >0.05 
 Regulatory 0.2 ± 0.02 0.1 ± 0.01 0.2 ± 0.01 0.1 ± 0.01 >0.05 
CD8+ T cells Total 8.3 ± 0.1 9.7 ± 0.1 9.9 ± 0.1 8.7 ± 0.1 >0.05 
 Effector 0.5 ± 0.02 0.3 ± 0.01 0.4 ± 0.02 0.4 ± 0.01 >0.05 
 Naive 5.9 ± 0.1 7.0 ± 0.1 7.1 ± 0.1 6.1 ± 0.2 >0.05 
B cells  11.9 ± 0.1 11.8 ± 0.2 11.9 ± 0.2 12.1 ± 0.1 >0.05 

Superficial LNs were taken from mice in DD at CTs. Cells of individual mice were stained for CD4, CD8, CD62L, and CD44; or CD4 and CD25; or IgM and B220. Naive T cells are those that stain CD4+ or CD8+ and CD44lowCD62Lhigh; effector T cells are CD4+ or CD8+ and CD44highCD62Llow; regulatory T cells are CD4+CD25+; and mature B cells are IgM+B220high. The experiment was done three times, except for staining for CD62L and CD44 (twice), with three to four mice each time. All numbers shown are mean ± SEM for the two to three experiments grouped, i.e., for 7–11 mice per time point (individual data points are shown in Supplemental Fig. 1C). The values in the table are the absolute number of cells in three LNs. All variations across CTs are nonsignificant by one-way ANOVA (rightmost column).

Given that the LN is the main site of T cell stimulation, we next investigated T cell proliferation potential over the 24-h cycle. LNs were collected from mice every 4 h over 24 h, under DD conditions. Cells were stained with the fluorescent cytosol marker CFSE and stimulated in vitro using an Ab directed against a TCR-associated polypeptide, CD3ε. This gives direct assessment of the ability of T cells to proliferate in response to Ag because anti-CD3 stimulation mimics Ag recognition by T cells. CD4+ and CD8+ T cells from LNs sampled from animals in the late day or in the night, and stimulated with anti-CD3, divided much more quickly than when LNs were taken in the early or mid day (Fig. 2A–E). As expected, CD8+ T cells proliferated more than CD4+ T cells. This rhythm is not due to the order of sampling, as it was also observed when day sampling was done after night sampling (Supplemental Fig. 2). This rhythm cannot be explained by a variation of T cells in the samples, because, as described above, the abundance of each T cell subset does not vary over the day in LNs.

In parallel to anti-CD3 stimulation, we also stimulated cells with PMA and ionomycin. There was no rhythm in the proliferation of CD8+ T cells stimulated by PMA/ionomycin (Fig. 2G, 2I), but interestingly, there was a low-amplitude rhythm in CD4+ T cells after PMA/ionomycin stimulation with a different profile than the one observed with anti-CD3 stimulation (Fig. 2F, 2H). An absence of rhythms in CD8+ T cells in response to PMA/ionomycin and a weak rhythm with a profile different from anti-CD3 stimulation in CD4+ T cells together suggest that the proliferation rhythm after triggering of the TCR is primarily due to mechanisms upstream of cell cycle regulation. This rhythmic control could be at the level of signaling events proximal to the intracellular tail of the TCR complex.

The experiments described above were done with mice housed in DD, showing the endogenous nature of the rhythm. To confirm that the rhythmic T cell response is dependent on the circadian clock, we triggered proliferation of cells from mice mutant for the gene Clock (33). We first confirmed that rhythms of clock gene expression in the LNs of Clock mutant mice are blunted (Fig. 3A). Whereas LN cells sampled from WT mice at CT14 (early subjective night) are more responsive to mitogenic stimulation than cells sampled at CT6 (mid subjective day), cells from Clock homozygous mutant littermates did not show such a day–night difference (Fig. 3B). This is not due to a difference in the relative amounts of cell subsets between WT and Clock mutant mice at CT6 or CT14 (Fig. 3C).

Both T cells and APCs are involved in T cell stimulation. Although we use a TCR-specific Ab in the above in vitro experiments, T cells receive the required costimulatory signals from APCs present in the culture. To shed light on the contribution of T cells in rhythmic proliferation, we enhanced the costimulatory signal using soluble anti-CD28. LN cells sampled from mice at CT14 are more responsive to mitogenic stimulation than cells sampled at CT6 regardless of whether the costimulatory signal comes only from the APCs included in the culture or with the enhancement of costimulation using soluble anti-CD28 Ab (Fig. 4). These results indicate that the effect of costimulation has less of an influence on rhythmic proliferation than do factors inherent to the T cells themselves.

We then examined T cell proliferation-relevant molecules that could be under circadian control. Our in vitro T cell proliferation experiments indicated that the rhythm of T cell proliferative response may be based on changes in the TCR-proximal signaling (Fig. 2). Thus, we first investigated cell surface expression of the TCR itself, and showed that it does not present a circadian rhythm in LNs (Fig. 5A). Thus, the rhythm in proliferation is not due to a variation in the number of TCRs over the day.

We next examined key molecules downstream of the TCR. After the engagement of the TCR, LCK is activated and phosphorylates the ζ-chains of the TCR-associated CD3 molecule. The cytosolic tyrosine kinase ZAP70 is then recruited and phosphorylated, and plays a pivotal role in initiating downstream signaling, including PKCθ activation (9). LNs were collected every 6 h over 24 h under DD conditions, and protein extracts were immunoblotted for LCK, PKCθ, and ZAP70. Although LCK and PKCθ did not present a circadian rhythm (Fig. 5B), ZAP70 protein expression was found to be rhythmic with peak levels at CT8 (Fig. 5B, 5C). Notably, Zap70 mRNA is also rhythmically expressed (Fig. 5D), with a profile similar to that of the protein, which indicates that this molecule might be controlled at the transcriptional level by a T cell circadian clock. A ZAP70 rhythm might underlie rhythmic T cell function. However, the high levels of proliferation potential of LN T cells persist throughout the late day and night, and present a bimodal profile (Fig. 2D, 2E), implying that additional molecular mechanisms are probably involved.

We next investigated T cell activation in response to Ag presentation in vivo. Bone marrow-derived DCs loaded with the class I-restricted OVA peptide (DC-OVA) were injected in mice either in the middle of the day (ZT6) or in the middle of the night (ZT18). Seven days later, mice were sacrificed at the same time point as injection. OVA-specific CD8+ T cells were quantified in the spleen using the Kb-OVA tetramer. Mice injected at ZT6 showed a higher proportion of OVA-specific CD8+ T cells than mice injected at ZT18 (Fig. 6A, 6B). To distinguish whether group differences were due to the injection or dissection time, we also injected mice at ZT6 and sacrificed them 7 d later, but at ZT18, and vice versa. Regardless of dissection time, mice injected at ZT6 showed a markedly higher response to immunization than mice injected at ZT18 (Fig. 6C). To see whether the T cells activated by DC-OVA acquired effector functions, splenocytes were restimulated with the OVA peptide for 6 h, followed by intracellular staining to detect IFN-γ. There were more IFN-γ–producing cells when animals were initially injected in the day than in the night (Fig. 6D, 6E), paralleling the results obtained with tetramer staining. The fold difference is proportional for these two markers of T cell expansion and activation. Thus, we can also conclude that the T cell activation day–night difference is not due to a more efficient staining of splenic T cells with tetramer at ZT6. These data indicate that T cells respond differentially to Ag presentation depending on time of day in vivo.

The results presented in this study reveal clear variations over the day of T cell function both in vitro and in vivo. When stimulated with a mitogenic signal in culture, T cells exhibit a circadian rhythm of proliferation levels. Notably, this is not due to variation of lymphocyte subsets in the LN or of CD28 costimulation from APCs, but might depend on the circadian control of TCR-dependent signaling proteins such as ZAP70, whose levels present a circadian rhythm. As shown with DC–OVA immunizations, T cell response to Ag presentation also presents a time dependency in vivo.

Previous work has uncovered a rhythm of proliferation of blood T cells in humans, with maximal levels in the night (35). Proliferation of rat LN cells in response to Con A, a lectin that activates T cells by bridging their TCRs, was also shown to present a daily variation, but with variable phases and amplitudes across conditions and studies (3638). However, in the latter case, rats were first treated with Freund's adjuvant, which activates the innate response, and the amplitude of the proliferation rhythm was greatly reduced without this prior stimulation. The same authors found rhythms of lymphocyte subsets in rat submaxillary LNs (36, 37) (but not in other studies; see Ref. 39 for example), which contrasts our data showing constant numbers of T cell subsets and B cells in mouse superficial LNs. Melatonin can alter the balance of T cell subsets (40). Because rats produce a high-amplitude melatonin rhythm (41), whereas C57BL/6 mice produce very low levels of this hormone (42), an effect of melatonin may explain the difference observed between the studies.

Notably, these previous studies were all conducted under alternation of light and dark, that is, in conditions that make it impossible to distinguish between endogenous control (e.g., by the clock) and environment-driven daily variations. Our in vitro T cell proliferation experiments were conducted with mice housed in DD, thus eliminating the possibility that the phenomenon results from light-driven processes. Consistent with a control by the circadian system, LN T cells from Clock mutant mice did not display the day–night difference in T cell proliferation. These mutant mice have a mutation in the gene encoding the protein CLOCK, leading to the skipping of one exon and to a dominant-negative form of this transcription factor (43, 44). These mice display a long free-running period in DD, and eventually arrhythmicity after 2–3 wk (33). Moreover, in these mutants, the expression of the genes controlled by the CLOCK/BMAL1 dimer is blunted (4547), and the amplitude of the corticosterone rhythm is reduced (48) or abolished (46) depending on the genetic background. These mutant mice also display shorter total sleep duration and reduced recovery sleep after sleep deprivation (49). Both corticosterone (21, 5052) and sleep (10, 53) are known to have an impact on the immune response. Yet, a detailed analysis of immune responses of Clock mutant mice has not been reported. The amplitude of leukocyte count rhythms is reduced in the peripheral blood of Clock mutant mice (48, 54). Notably, changes in cell numbers are unlikely to explain the proliferation described in this work, as we found no difference in levels of lymphocyte subsets in LNs of Clock mutant versus WT mice.

Because Clock mutant mice lack a robust functional clock in all cell types, our experiment does not indicate which cell types or which, if any, humoral or neuronal cues contribute to circadian control of T cell response. Our results show a clear circadian variation in T cell proliferation in whole LN extracts, containing both T cells and APCs. Hayashi et al. (32) and Keller et al. (30) showed rhythmic clock gene expression in macrophages, which are one type of APCs, and rhythms in phagocytic activity and cytokine release function, respectively. Thus, the contribution of either a T cell clock, an APC clock, or other clocks in the organism (e.g., the SCN) to T cell proliferation remains to be determined. To address the potential contribution of APCs to the in vitro T cell proliferation rhythm, we performed an experiment in which the CD28 costimulatory signal was strongly enhanced. CD28 engagement is a key signal for the priming of naive T cells. In this experiment, high CD28 engagement mimicked a strong and constant CD28 costimulation signal from APCs. Even under these conditions, T cells still exhibited a marked day–night difference in their capacity to proliferate. These data reduce the likelihood of a rhythmic contribution from APCs, and indicate that the T cells themselves have the intrinsic ability to respond differentially to TCR trigger according to the time of day.

The results we show illustrating a strong day–night difference in response to immunization in vivo provide additional insights. The DCs injected both at ZT6 and ZT18 were from the same culture of matured bone marrow cells collected at a single time point from donor mice. Given that the DC–OVA immunization approach bypasses events prior to Ag presentation to T cells, the effect we see in vivo is most likely not due to an APC-intrinsic clock. For the same reasons, this day–night difference of response to DC–OVA immunization is probably not explained by the stronger innate responses observed in the rest period (11, 12, 30, 53). Rather, the day–night difference probably relies either on a T cell-intrinsic factor, on the effect of a systemic cue acting on T cells (immunomodulatory hormones, autonomic signals, sleep/wake cycles, spleen environment), or on a combination of rhythmic T cell and rhythmic systemic cue processes.

Such a multilevel additive circadian regulation has been described for other physiological processes and tissues. For example, the high-amplitude rhythm of corticosterone synthesis and release depends on both SCN central clock function and the local adrenal gland clock (55, 56). Also, circadian regulation in the liver relies on both the hepatocyte clock and systemic cues (57, 58). A complex, organism-wide phenomenon such as the immune response is well placed for a multigated circadian control. Within this complex response, a central event like the T cell response to Ag and proliferation is a prime position for such a point of control.

As an initial step in identifying points of circadian control of T cell function, we assessed the abundance over the day of several proteins downstream of the TCR. We found that ZAP70 protein levels follow a circadian rhythm. This is paralleled by a rhythm in the mRNA, suggesting regulation at the transcriptional level by a T cell clock machinery. Interestingly, the search for putative transcription factor binding sites in the Zap70 gene revealed two regions of high similarity between mouse and human genomes that contained E-boxes, DNA elements that can be bound by the clock transcription factors CLOCK and BMAL1 (1), just upstream of the first exon and within the first intron (data not shown). Could a rhythm of ZAP70 protein levels lead to a daily variation of the capacity of T cells to proliferate? Variations in abundance of proteins modulating TCR responses were shown to impact the amount of engagement of the TCR in complexes with MHC-Ag or the strength of signaling downstream of the TCR (59). Changes in expression levels of ZAP70 itself were shown to be functionally relevant in vivo: in the thymus, ZAP70 protein levels vary throughout thymocyte development, and these changes in abundance underlie the precise timing of development of single-positive CD4+ and CD8+ T cells (60). Thus, the circadian rhythm of ZAP70 levels in LN T cells could contribute to the rhythm of proliferation potential. Nevertheless, the high levels of proliferation potential persist throughout the late day and night, implying that other mechanisms are likely to be involved. Additional levels of regulation could include the following: rhythms of other proteins downstream of the TCR, rhythms of T cell responsiveness to cytokines or to APC costimulatory signals, rhythms in the levels to which signaling pathways can be activated, or a daily variation in the spatiotemporal regulation of the TCR or TCR-associated factors.

A time-of-day preference for immune potentiation could be advantageous to optimize the response at phases of the day during which a pathogen is most likely to be encountered. Alternatively, a period of low responsiveness could prevent a hypersensitivity that could lead to deleterious complications, including allergy and autoimmunity. Beyond this, circadian control of T cell proliferation and the adaptive immune response provides a new paradigm for exploration of disease susceptibility. In particular, a higher incidence of cancer was shown in conditions in which circadian rhythms were disrupted (6163), and in particular in shift workers (6466).

Furthermore, our striking observation that the time of immunization has a strong impact on the magnitude of the adaptive immune response as measured by looking at T cell activation, proliferation, and acquisition of effector functions many days later, lends serious support to the novel idea that the efficiency of vaccination might benefit from careful selection of the time of administration. A few pilot studies with subjects receiving vaccine for hepatitis A or B or influenza have indicated that, indeed, the outcome of vaccination could depend on the time of administration (6769). Our work gives a clear and controlled demonstration of this concept and provides an easily tractable system to dissect the mechanisms at play.

We thank Aude Villemain, Nathalie Henley, Martine Dupuis, Jean-Francois Daudelin, Rachel Morin, and Shalon Liu for technical assistance; Philippe Boudreau for help with cosinor analyses; Susan Westfall for critical reading of the manuscript; and all members of the laboratories of N.C. and N.L. for continuous support and discussions.

This work was supported by Canadian Institutes for Health Research Operating Grant MOP-84441 (to N.C. and N.L.), Fonds de la Recherche en Santé du Québec salary awards (to N.C. and N.L.), a McGill University Faculty of Medicine award (to N.C.), and a Fonds de la Recherche en Santé du Québec graduate scholarship (to E.E.F.).

The online version of this article contains supplemental material.

Abbreviations used in this article:

Cry

Cryptochrome

CT

circadian time

DC

dendritic cell

DD

constant darkness

LCK

lymphocyte-specific protein tyrosine kinase

LD

light:dark

LN

lymph node

Per

Period

PKC

protein kinase C

RPA

RNase protection assay

SCN

suprachiasmatic nucleus

WT

wild-type

ZT

Zeitgeber time.

1
Duguay
D.
,
Cermakian
N.
.
2009
.
The crosstalk between physiology and circadian clock proteins.
Chronobiol. Int.
26
:
1479
1513
.
2
Guilding
C.
,
Piggins
H. D.
.
2007
.
Challenging the omnipotence of the suprachiasmatic timekeeper: are circadian oscillators present throughout the mammalian brain?
Eur. J. Neurosci.
25
:
3195
3216
.
3
Balsalobre
A.
2002
.
Clock genes in mammalian peripheral tissues.
Cell Tissue Res.
309
:
193
199
.
4
Ko
C. H.
,
Takahashi
J. S.
.
2006
.
Molecular components of the mammalian circadian clock.
Hum. Mol. Genet.
15
(
Spec. No. 2
):
R271
R277
.
5
Duffield
G. E.
2003
.
DNA microarray analyses of circadian timing: the genomic basis of biological time.
J. Neuroendocrinol.
15
:
991
1002
.
6
Panda
S.
,
Antoch
M. P.
,
Miller
B. H.
,
Su
A. I.
,
Schook
A. B.
,
Straume
M.
,
Schultz
P. G.
,
Kay
S. A.
,
Takahashi
J. S.
,
Hogenesch
J. B.
.
2002
.
Coordinated transcription of key pathways in the mouse by the circadian clock.
Cell
109
:
307
320
.
7
Storch
K. F.
,
Lipan
O.
,
Leykin
I.
,
Viswanathan
N.
,
Davis
F. C.
,
Wong
W. H.
,
Weitz
C. J.
.
2002
.
Extensive and divergent circadian gene expression in liver and heart.
Nature
417
:
78
83
.
8
Goldrath
A. W.
,
Bevan
M. J.
.
1999
.
Selecting and maintaining a diverse T-cell repertoire.
Nature
402
:
255
262
.
9
Smith-Garvin
J. E.
,
Koretzky
G. A.
,
Jordan
M. S.
.
2009
.
T cell activation.
Annu. Rev. Immunol.
27
:
591
619
.
10
Haus
E.
,
Smolensky
M. H.
.
1999
.
Biologic rhythms in the immune system.
Chronobiol. Int.
16
:
581
622
.
11
Halberg
F.
,
Johnson
E. A.
,
Brown
B. W.
,
Bittner
J. J.
.
1960
.
Susceptibility rhythm to E. coli endotoxin and bioassay.
Proc. Soc. Exp. Biol. Med.
103
:
142
144
.
12
Marpegan
L.
,
Leone
M. J.
,
Katz
M. E.
,
Sobrero
P. M.
,
Bekinstein
T. A.
,
Golombek
D. A.
.
2009
.
Diurnal variation in endotoxin-induced mortality in mice: correlation with proinflammatory factors.
Chronobiol. Int.
26
:
1430
1442
.
13
Castanon-Cervantes
O.
,
Wu
M.
,
Ehlen
J. C.
,
Paul
K.
,
Gamble
K. L.
,
Johnson
R. L.
,
Besing
R. C.
,
Menaker
M.
,
Gewirtz
A. T.
,
Davidson
A. J.
.
2010
.
Dysregulation of inflammatory responses by chronic circadian disruption.
J. Immunol.
185
:
5796
5805
.
14
Shackelford
P. G.
,
Feigin
R. D.
.
1973
.
Periodicity of susceptibility to pneumococcal infection: influence of light and adrenocortical secretions.
Science
182
:
285
287
.
15
Cutolo
M.
,
Seriolo
B.
,
Craviotto
C.
,
Pizzorni
C.
,
Sulli
A.
.
2003
.
Circadian rhythms in RA.
Ann. Rheum. Dis.
62
:
593
596
.
16
Sutherland
E. R.
2005
.
Nocturnal asthma.
J. Allergy Clin. Immunol.
116
:
1179
1186, quiz 1187
.
17
Abo
T.
,
Kawate
T.
,
Itoh
K.
,
Kumagai
K.
.
1981
.
Studies on the bioperiodicity of the immune response. I. Circadian rhythms of human T, B, and K cell traffic in the peripheral blood.
J. Immunol.
126
:
1360
1363
.
18
Arjona
A.
,
Sarkar
D. K.
.
2006
.
The circadian gene mPer2 regulates the daily rhythm of IFN-gamma.
J. Interferon Cytokine Res.
26
:
645
649
.
19
Born
J.
,
Lange
T.
,
Hansen
K.
,
Mölle
M.
,
Fehm
H. L.
.
1997
.
Effects of sleep and circadian rhythm on human circulating immune cells.
J. Immunol.
158
:
4454
4464
.
20
Deprés-Brummer
P.
,
Bourin
P.
,
Pages
N.
,
Metzger
G.
,
Lévi
F.
.
1997
.
Persistent T lymphocyte rhythms despite suppressed circadian clock outputs in rats.
Am. J. Physiol.
273
:
R1891
R1899
.
21
Dimitrov
S.
,
Benedict
C.
,
Heutling
D.
,
Westermann
J.
,
Born
J.
,
Lange
T.
.
2009
.
Cortisol and epinephrine control opposing circadian rhythms in T cell subsets.
Blood
113
:
5134
5143
.
22
Archer
S. N.
,
Viola
A. U.
,
Kyriakopoulou
V.
,
von Schantz
M.
,
Dijk
D. J.
.
2008
.
Inter-individual differences in habitual sleep timing and entrained phase of endogenous circadian rhythms of BMAL1, PER2 and PER3 mRNA in human leukocytes.
Sleep
31
:
608
617
.
23
James
F. O.
,
Boivin
D. B.
,
Charbonneau
S.
,
Bélanger
V.
,
Cermakian
N.
.
2007
.
Expression of clock genes in human peripheral blood mononuclear cells throughout the sleep/wake and circadian cycles.
Chronobiol. Int.
24
:
1009
1034
.
24
Kusanagi
H.
,
Hida
A.
,
Satoh
K.
,
Echizenya
M.
,
Shimizu
T.
,
Pendergast
J. S.
,
Yamazaki
S.
,
Mishima
K.
.
2008
.
Expression profiles of 10 circadian clock genes in human peripheral blood mononuclear cells.
Neurosci. Res.
61
:
136
142
.
25
Arjona
A.
,
Sarkar
D. K.
.
2005
.
Circadian oscillations of clock genes, cytolytic factors, and cytokines in rat NK cells.
J. Immunol.
174
:
7618
7624
.
26
Oishi
K.
,
Sakamoto
K.
,
Okada
T.
,
Nagase
T.
,
Ishida
N.
.
1998
.
Antiphase circadian expression between BMAL1 and period homologue mRNA in the suprachiasmatic nucleus and peripheral tissues of rats.
Biochem. Biophys. Res. Commun.
253
:
199
203
.
27
Alvarez
J. D.
,
Sehgal
A.
.
2005
.
The thymus is similar to the testis in its pattern of circadian clock gene expression.
J. Biol. Rhythms
20
:
111
121
.
28
Granda
T. G.
,
Liu
X. H.
,
Smaaland
R.
,
Cermakian
N.
,
Filipski
E.
,
Sassone-Corsi
P.
,
Lévi
F.
.
2005
.
Circadian regulation of cell cycle and apoptosis proteins in mouse bone marrow and tumor.
FASEB J.
19
:
304
306
.
29
Hashiramoto
A.
,
Yamane
T.
,
Tsumiyama
K.
,
Yoshida
K.
,
Komai
K.
,
Yamada
H.
,
Yamazaki
F.
,
Doi
M.
,
Okamura
H.
,
Shiozawa
S.
.
2010
.
Mammalian clock gene Cryptochrome regulates arthritis via proinflammatory cytokine TNF-alpha.
J. Immunol.
184
:
1560
1565
.
30
Keller
M.
,
Mazuch
J.
,
Abraham
U.
,
Eom
G. D.
,
Herzog
E. D.
,
Volk
H. D.
,
Kramer
A.
,
Maier
B.
.
2009
.
A circadian clock in macrophages controls inflammatory immune responses.
Proc. Natl. Acad. Sci. USA
106
:
21407
21412
.
31
Tsinkalovsky
O.
,
Smaaland
R.
,
Rosenlund
B.
,
Sothern
R. B.
,
Hirt
A.
,
Steine
S.
,
Badiee
A.
,
Abrahamsen
J. F.
,
Eiken
H. G.
,
Laerum
O. D.
.
2007
.
Circadian variations in clock gene expression of human bone marrow CD34+ cells.
J. Biol. Rhythms
22
:
140
150
.
32
Hayashi
M.
,
Shimba
S.
,
Tezuka
M.
.
2007
.
Characterization of the molecular clock in mouse peritoneal macrophages.
Biol. Pharm. Bull.
30
:
621
626
.
33
Vitaterna
M. H.
,
King
D. P.
,
Chang
A. M.
,
Kornhauser
J. M.
,
Lowrey
P. L.
,
McDonald
J. D.
,
Dove
W. F.
,
Pinto
L. H.
,
Turek
F. W.
,
Takahashi
J. S.
.
1994
.
Mutagenesis and mapping of a mouse gene, Clock, essential for circadian behavior.
Science
264
:
719
725
.
34
Preitner
N.
,
Damiola
F.
,
Lopez-Molina
L.
,
Zakany
J.
,
Duboule
D.
,
Albrecht
U.
,
Schibler
U.
.
2002
.
The orphan nuclear receptor REV-ERBalpha controls circadian transcription within the positive limb of the mammalian circadian oscillator.
Cell
110
:
251
260
.
35
Bollinger
T.
,
Bollinger
A.
,
Skrum
L.
,
Dimitrov
S.
,
Lange
T.
,
Solbach
W.
.
2009
.
Sleep-dependent activity of T cells and regulatory T cells.
Clin. Exp. Immunol.
155
:
231
238
.
36
Bonacho
M. G.
,
Cardinali
D. P.
,
Castrillón
P.
,
Cutrera
R. A.
,
Esquifino
A. I.
.
2001
.
Aging-induced changes in 24-h rhythms of mitogenic responses, lymphocyte subset populations and neurotransmitter and amino acid content in rat submaxillary lymph nodes during Freund’s adjuvant arthritis.
Exp. Gerontol.
36
:
267
282
.
37
Esquifino
A. I.
,
Chacon
F.
,
Cano
P.
,
Marcos
A.
,
Cutrera
R. A.
,
Cardinali
D. P.
.
2004
.
Twenty-four-hour rhythms of mitogenic responses, lymphocyte subset populations and amino acid content in submaxillary lymph nodes of growing male rats subjected to calorie restriction.
J. Neuroimmunol.
156
:
66
73
.
38
Esquifino
A. I.
,
Gonzalez
M. E.
,
Lafuente
A.
.
1996
.
Possible interactions of cyclosporine and hyperprolactinemia modulating the episodic secretion of prolactin.
Proc. Soc. Exp. Biol. Med.
213
:
206
211
.
39
Jiménez-Ortega
V.
,
Fernández-Mateos
M. P.
,
Barquilla
P. C.
,
Cardinali
D. P.
,
Esquifino
A. I.
.
2011
.
Continuous versus discontinuous drinking of an ethanol liquid diet in peripubertal rats: effect on 24-h variation of lymph node and splenic mitogenic responses and lymphocyte subset populations.
Alcohol
45
:
183
192
.
40
Castrillón
P. O.
,
Esquifino
A. I.
,
Varas
A.
,
Zapata
A.
,
Cutrera
R. A.
,
Cardinali
D. P.
.
2000
.
Effect of melatonin treatment on 24-h variations in responses to mitogens and lymphocyte subset populations in rat submaxillary lymph nodes.
J. Neuroendocrinol.
12
:
758
765
.
41
Klein
D. C.
,
Weller
J. L.
,
Moore
R. Y.
.
1971
.
Melatonin metabolism: neural regulation of pineal serotonin: acetyl coenzyme A N-acetyltransferase activity.
Proc. Natl. Acad. Sci. USA
68
:
3107
3110
.
42
Vivien-Roels
B.
,
Malan
A.
,
Rettori
M. C.
,
Delagrange
P.
,
Jeanniot
J. P.
,
Pévet
P.
.
1998
.
Daily variations in pineal melatonin concentrations in inbred and outbred mice.
J. Biol. Rhythms
13
:
403
409
.
43
Gekakis
N.
,
Staknis
D.
,
Nguyen
H. B.
,
Davis
F. C.
,
Wilsbacher
L. D.
,
King
D. P.
,
Takahashi
J. S.
,
Weitz
C. J.
.
1998
.
Role of the CLOCK protein in the mammalian circadian mechanism.
Science
280
:
1564
1569
.
44
King
D. P.
,
Zhao
Y.
,
Sangoram
A. M.
,
Wilsbacher
L. D.
,
Tanaka
M.
,
Antoch
M. P.
,
Steeves
T. D.
,
Vitaterna
M. H.
,
Kornhauser
J. M.
,
Lowrey
P. L.
, et al
.
1997
.
Positional cloning of the mouse circadian clock gene.
Cell
89
:
641
653
.
45
Jin
X.
,
Shearman
L. P.
,
Weaver
D. R.
,
Zylka
M. J.
,
de Vries
G. J.
,
Reppert
S. M.
.
1999
.
A molecular mechanism regulating rhythmic output from the suprachiasmatic circadian clock.
Cell
96
:
57
68
.
46
Kennaway
D. J.
,
Owens
J. A.
,
Voultsios
A.
,
Varcoe
T. J.
.
2006
.
Functional central rhythmicity and light entrainment, but not liver and muscle rhythmicity, are Clock independent.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
291
:
R1172
R1180
.
47
Mongrain
V.
,
Ruan
X.
,
Dardente
H.
,
Fortier
E. E.
,
Cermakian
N.
.
2008
.
Clock-dependent and independent transcriptional control of the two isoforms from the mouse Rorgamma gene.
Genes Cells
13
:
1197
1210
.
48
Oishi
K.
,
Ohkura
N.
,
Kadota
K.
,
Kasamatsu
M.
,
Shibusawa
K.
,
Matsuda
J.
,
Machida
K.
,
Horie
S.
,
Ishida
N.
.
2006
.
Clock mutation affects circadian regulation of circulating blood cells.
J. Circadian Rhythms
4
:
13
.
49
Naylor
E.
,
Bergmann
B. M.
,
Krauski
K.
,
Zee
P. C.
,
Takahashi
J. S.
,
Vitaterna
M. H.
,
Turek
F. W.
.
2000
.
The circadian clock mutation alters sleep homeostasis in the mouse.
J. Neurosci.
20
:
8138
8143
.
50
Hiemke
C.
,
Brunner
R.
,
Hammes
E.
,
Müller
H.
,
Meyer zum Büschenfelde
K. H.
,
Lohse
A. W.
.
1995
.
Circadian variations in antigen-specific proliferation of human T lymphocytes and correlation to cortisol production.
Psychoneuroendocrinology
20
:
335
342
.
51
Sterzer
P.
,
Wiegers
G. J.
,
Reul
J. M.
.
2004
.
Long-term in vivo administration of glucocorticoid hormones attenuates their capacity to accelerate in vitro proliferation of rat splenic T cells.
Endocrinology
145
:
3630
3638
.
52
Wiegers
G. J.
,
Labeur
M. S.
,
Stec
I. E.
,
Klinkert
W. E.
,
Holsboer
F.
,
Reul
J. M.
.
1995
.
Glucocorticoids accelerate anti-T cell receptor-induced T cell growth.
J. Immunol.
155
:
1893
1902
.
53
Lange
T.
,
Dimitrov
S.
,
Born
J.
.
2010
.
Effects of sleep and circadian rhythm on the human immune system.
Ann. N. Y. Acad. Sci.
1193
:
48
59
.
54
Gorbacheva
V. Y.
,
Kondratov
R. V.
,
Zhang
R.
,
Cherukuri
S.
,
Gudkov
A. V.
,
Takahashi
J. S.
,
Antoch
M. P.
.
2005
.
Circadian sensitivity to the chemotherapeutic agent cyclophosphamide depends on the functional status of the CLOCK/BMAL1 transactivation complex.
Proc. Natl. Acad. Sci. USA
102
:
3407
3412
.
55
Oster
H.
,
Damerow
S.
,
Kiessling
S.
,
Jakubcakova
V.
,
Abraham
D.
,
Tian
J.
,
Hoffmann
M. W.
,
Eichele
G.
.
2006
.
The circadian rhythm of glucocorticoids is regulated by a gating mechanism residing in the adrenal cortical clock.
Cell Metab.
4
:
163
173
.
56
Son
G. H.
,
Chung
S.
,
Choe
H. K.
,
Kim
H. D.
,
Baik
S. M.
,
Lee
H.
,
Lee
H. W.
,
Choi
S.
,
Sun
W.
,
Kim
H.
, et al
.
2008
.
Adrenal peripheral clock controls the autonomous circadian rhythm of glucocorticoid by causing rhythmic steroid production.
Proc. Natl. Acad. Sci. USA
105
:
20970
20975
.
57
Kornmann
B.
,
Schaad
O.
,
Bujard
H.
,
Takahashi
J. S.
,
Schibler
U.
.
2007
.
System-driven and oscillator-dependent circadian transcription in mice with a conditionally active liver clock.
PLoS Biol.
5
:
e34
.
58
Lamia
K. A.
,
Storch
K. F.
,
Weitz
C. J.
.
2008
.
Physiological significance of a peripheral tissue circadian clock.
Proc. Natl. Acad. Sci. USA
105
:
15172
15177
.
59
Feinerman
O.
,
Veiga
J.
,
Dorfman
J. R.
,
Germain
R. N.
,
Altan-Bonnet
G.
.
2008
.
Variability and robustness in T cell activation from regulated heterogeneity in protein levels.
Science
321
:
1081
1084
.
60
Saini
M.
,
Sinclair
C.
,
Marshall
D.
,
Tolaini
M.
,
Sakaguchi
S.
,
Seddon
B.
.
2010
.
Regulation of Zap70 expression during thymocyte development enables temporal separation of CD4 and CD8 repertoire selection at different signaling thresholds.
Sci. Signal.
3
:
ra23
.
61
Filipski
E.
,
Delaunay
F.
,
King
V. M.
,
Wu
M. W.
,
Claustrat
B.
,
Gréchez-Cassiau
A.
,
Guettier
C.
,
Hastings
M. H.
,
Francis
L.
.
2004
.
Effects of chronic jet lag on tumor progression in mice.
Cancer Res.
64
:
7879
7885
.
62
Filipski
E.
,
King
V. M.
,
Li
X.
,
Granda
T. G.
,
Mormont
M. C.
,
Liu
X.
,
Claustrat
B.
,
Hastings
M. H.
,
Lévi
F.
.
2002
.
Host circadian clock as a control point in tumor progression.
J. Natl. Cancer Inst.
94
:
690
697
.
63
Fu
L.
,
Pelicano
H.
,
Liu
J.
,
Huang
P.
,
Lee
C.
.
2002
.
The circadian gene Period2 plays an important role in tumor suppression and DNA damage response in vivo.
Cell
111
:
41
50
.
64
Hansen
J.
2001
.
Light at night, shiftwork, and breast cancer risk.
J. Natl. Cancer Inst.
93
:
1513
1515
.
65
Schernhammer
E. S.
,
Laden
F.
,
Speizer
F. E.
,
Willett
W. C.
,
Hunter
D. J.
,
Kawachi
I.
,
Fuchs
C. S.
,
Colditz
G. A.
.
2003
.
Night-shift work and risk of colorectal cancer in the nurses’ health study.
J. Natl. Cancer Inst.
95
:
825
828
.
66
Straif
K.
,
Baan
R.
,
Grosse
Y.
,
Secretan
B.
,
El Ghissassi
F.
,
Bouvard
V.
,
Altieri
A.
,
Benbrahim-Tallaa
L.
,
Cogliano
V.
.
2007
.
Carcinogenicity of shift-work, painting, and fire-fighting.
Lancet Oncol.
8
:
1065
1066
.
67
Feigin
R. D.
,
Jaeger
R. F.
,
McKinney
R. W.
,
Alevizatos
A. C.
.
1967
.
Live, attenuated Venezuelan equine encephalomyelitis virus vaccine. II. Whole-blood amino-acid and fluorescent-antibody studies following immunization.
Am. J. Trop. Med. Hyg.
16
:
769
777
.
68
Langlois
P. H.
,
Smolensky
M. H.
,
Glezen
W. P.
,
Keitel
W. A.
.
1995
.
Diurnal variation in responses to influenza vaccine.
Chronobiol. Int.
12
:
28
36
.
69
Phillips
A. C.
,
Gallagher
S.
,
Carroll
D.
,
Drayson
M.
.
2008
.
Preliminary evidence that morning vaccination is associated with an enhanced antibody response in men.
Psychophysiology
45
:
663
666
.

The authors have no financial conflicts of interest.