|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Department of Veterinary Biosciences,3
Neuroscience Program,4 University of Illinois at Urbana-Champaign, Urbana, Illinois 61802
ABSTRACT
A molecular device that measures time on a daily, or circadian, scale is a nearly ubiquitous feature of eukaryotic organisms. A core group of clock genes, whose coordinated function is required for this timekeeping, is expressed both in the central clock and within numerous peripheral organs. We examined expression of clock genes in the rat ovary. Transcripts for core oscillator elements (Arntl,Clock,Per1,Per2, andCry1) were present in the ovary as indicated by quantitative real-time RT-PCR. Rhythmic expression patterns ofArntlandPer2transcripts and protein products were out of phase with respect to the central oscillator and in complete antiphase to each other. Expression ofArntlwas significantly elevated after the LH surge on the day of proestrus. Finally, hCG treatment induced cyclic expression of bothArntlandPer2gene products in hypophysectomized, immature rats primed with eCG. Collectively, these data suggest that the core underpinnings of the transcriptional/translational feedback loop that drives circadian rhythmicity is present in the rat ovary. Furthermore, the study identifies LH as a potential regulator of circadian clock gene rhythms in the ovary.
circadian rhythm, gene regulation, luteinizing hormone, mechanisms of hormone action, ovary
INTRODUCTION
Endogenous circadian rhythms measure time on a scale of
24 h and orchestrate myriad physiological processes to allow organismic synchrony with the external environment. In mammals, entrainment to alternating periods of environmental light and darkness occurs at the level of the suprachiasmatic nucleus (SCN) of the basal hypothalamus primarily through glutamatergic neurotransmission via the retinohypothalamic tract [17]. Cellular circadian rhythms are generated by a highly conserved molecular mechanism comprised of a core group of clock genes and their protein products. These molecular components form interlocking transcriptional/translational autoregulatory feedback loops that oscillate with a near 24-h periodicity [811]. The transcription factor, aryl hydrocarbon receptor nuclear translocator-like (ARNTL; also frequently referred to as BMAL1 or MOP3), provides positive drive to the system by forming a heterodimer with CLOCK through protein-protein interactions between their period-arnt-single minded (PAS) domains. This ARNTL/CLOCK heterodimer drives transcription of cryptochrome (Cry1 and Cry2), Period (Pers 1, 2, and 3), and Timeless genes through binding to E-box elements in their promoter regions [811]. PER and CRY proteins undergo posttranslational modifications, form multimeric complexes and feed back to inhibit their own transcription [12]. In addition, PER2 acts as a positive regulator for Arntl [13, 14].
Circadian oscillations are present in virtually all tissues. Microarray analysis has revealed that 2%10% of all genes exhibit circadian oscillations, but that rhythmicity of individual genes is tissue-specific, with less than 5% overlap between tissues [15]. Clock genes are rhythmically expressed in various peripheral tissues, including liver, muscle, kidney, and heart [1624]. Regulation of molecular circadian rhythmicity in peripheral tissues remains enigmatic. The SCN likely provides an important source of the timing signal. Peak expression of Per2 and Arntl in liver, muscle, and kidney occurs in late subjective day and late subjective night, respectively, contrasting significantly with their corresponding peaks in the SCN. Furthermore, SCN ablation destroys rhythms of Per1, Per2, and Arntl in the liver [22, 25]. Despite recent evidence that peripheral organs may retain rhythmic features in vitro without the SCN [26], the overwhelming majority of data suggests that the SCN is required to synchronize cellular oscillations within a tissue.
The SCN plays an important role in regulating timing of the LH surge within the female reproductive axis. SCN-lesioned female rats lack the ovulatory LH surge and fail to ovulate [27]. Although it occurs only once every 4 or 5 days in the intact rat, the LH surge is precisely timed to occur on the afternoon of proestrus [28, 29]; daily LH surges are present in ovariectomized, estrogen-treated rats [28, 29]. Furthermore, external stimuli that alter the timing of the circadian activity rhythm likewise have the same effect on the timing of the LH surge [30]. This study examines events downstream of the LH surge to begin to determine whether the circadianregulated LH surge may participate in regulation of circadian clock gene rhythms in the ovary.
Investigation of the role of clock genes in reproductive function remains in its infancy. Several studies in the mouse testis failed to detect rhythmic expression of Arntl, Per1, and Cry1 [3134]; other studies revealed clear oscillations of testicular Per1, Per3, and Arntl [18, 35, 36]. Within the female reproductive system, Arntl and Per2 transcripts are rhythmically expressed in the rat oviduct [37]. However, clock gene expression in the ovary has not been investigated. This study begins to characterize expression and regulation of core clock genes in the ovary. We hypothesized that oscillations in core elements of the molecular clock are present in the ovary and proposed that the LH surge may regulate expression of one or more clock genes. Temporal expression patterns of Arntl and Per2 transcripts over a 24-h period and on 2 disparate days of the rat estrous cycle revealed daily oscillatory patterns of expression. Clock gene rhythms were absent in hypophysectomized, immature rats treated only with eCG. However, rhythmicity was induced in these animals upon treatment with LH (hCG). Arntl and Per2 expression patterns are consistent with the presence of functional molecular clockwork machinery in the ovary, and suggest that LH may play a role in regulating ovarian clock gene rhythms through actions on Arntl and Per2.
MATERIALS AND METHODS
Intact Adult Experiments
Eight-week-old female Sprague-Dawley rats (Harlan, Indianapolis, IN) were maintained in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals, and all protocols were approved by the University of Illinois Animal Care and Use Committee. Food and water were provided ad libitum. Rats were housed under a strict 12L:12D lighting schedule. Under these lighting conditions, denoted as LD, zeitgeber time 0 (ZT 0) represents the time of lights-on in the colony. The corresponding time of lights-off is denoted as ZT 12. Animals were allowed to entrain to this lighting schedule for 2 wk before experimentation. Animals were killed by decapitation every 4 h over a 24-h period (n = 4 per group) starting at ZT 0. Ovaries were dissected from the animals and snap frozen in liquid nitrogen for immunoblot or mRNA analysis, or placed in 10% neutral buffered formalin for immunohistochemistry (IHC) and in situ hybridization. For analysis of mRNA and protein from SCN, hypothalamic brain slices (500 µM) containing the paired SCN were trimmed under a dissecting microscope to preserve only the SCN and underlying optic chiasm as described previously [38]. These reduced slices retain the electrophysiological and biochemical properties observed in full-size SCN slices [39]. Samples were frozen immediately in liquid nitrogen and stored at 80°C until processed for quantitative real-time RT-PCR (qPCR) or immunoblot analysis.
Estrous cycle stage was assessed from vaginal smears obtained every morning at ZT 4. Diestrus I was defined as vaginal smears containing abundant leukocytes mixed with cornified cells that appeared the day following a smear comprised entirely of cornified cells. Proestrus was defined as vaginal smears containing mostly small, nucleated epithelial cells, singly or in sheets, and the absence of leukocytes. Animals were monitored for at least two full estrous cycles before proceeding with the study, and those with irregular estrous cycles were omitted. After determination of cycle stage, animals were killed every 4 h over a 24-h period beginning at ZT 6 (2 h after staging) and ending at ZT 2 the following day (n = 3 per treatment group).
Hypophysectomized Immature Female Experiments
To examine the effects of gonadotropins on ovarian clock components, juvenile hypophysectomized animals were used. Female Sprague-Dawley rats were hypophysectomized at 20 days of age at Charles River Laboratories (Wilmington, MA). Animals arrived at 27 days of age, and were allowed to acclimate to a strict 12L:12D schedule for 48 h before treatment. Glucose (5%) was provided in the drinking water and along with food ad libitum. Animals were randomly assigned to hCG treatment and control groups (eCG followed by vehicle). On Day 30, both groups were administered 10 IU eCG in PBS (s.c.; Sigma, St. Louis, MO) to induce follicular development [40]. After 52 h, the treatment group was administered 10 IU hCG in PBS and controls received vehicle. Animals were killed by decapitation 0, 4, 8, 12, and 24 h after hCG treatment. Postmortem assessment was performed to confirm complete hypophysectomy (n = 35 per group).
Quantitative Real-time RT-PCR
After disruption and homogenization of the tissue, total ovarian or SCN RNA was extracted using TRIzol reagent according to the manufacturer's protocol (Invitrogen, Carlsbad, CA). Concentration, purity, and quality of the RNA were determined using the NanoDrop ND-1000 UV spectrophotometer at 260 nm and 280 nm and by gel electrophoresis. Total RNA (1 µg) was reverse transcribed using 200 ng random hexamers, 200 U SSII reverse transcriptase, 10 mM DTT, and 1x First-Strand Buffer (Invitrogen) at 43°C for 1 h in a volume of 20 µl (Mycycler; Bio-Rad, Hercules, CA). Negative controls included omission of reverse transcriptase and omission of template. Primers were chosen using Primer Express 2.0 software (Applied Biosystems, Foster City, CA) for optimum use in qPCR. A BLASTN search was performed in GenBank to ensure that all primers were unique to the gene of interest. To avoid amplification from genomic DNA contamination, all primer sets spanned a large exon-intron-exon junction (see Table 1). All qPCR reactions were performed in the ABI prism 7700 (Applied Biosystems) in triplicate with 12.5 µl SYBR Green Master Mix (Applied Biosystems), 5 µl of 5x cDNA, and primers in a total volume of 25 µl. Initial denaturation occurred at 95°C for 10 min, followed by 40 cycles of 15 s at 95°C and 15 s at 60°C. A heat dissociation was performed at the end of every run to assure specificity. Optimum primer concentration was determined for each primer pair by comparing the
Rn from a matrix of forward and reverse primer concentrations to the Ct (the fractional cycle number at which the fluorescence passes the threshold) for each pair. Optimum primer concentrations are listed in Table 2. An arbitrary threshold of fluorescence was set within the exponential phase of amplification for each gene product and kept constant for all assays. The cycle at which amplification of the product exceeded this threshold was determined and designated as the Ct value. Dissociation curves, gel analysis, and sequencing of certain PCR products confirmed gene-specific product amplification.
|
Standard Curves and Calculations
Relative standard curves were created with serially diluted total RNA samples from the ovaries and SCN used in the experiments (modified from [41]). Three sets of pooled RNA samples were serially diluted by a factor of two starting at 2 µg and ending at 125 ng. This range was determined to be well within the detection level and sensitivity of the qPCR assay. RNA samples used for standard curves are representative of 36 animals. After reverse transcription, the resulting cDNA was PCR amplified in the ABI 7700 in duplicate. Average Ct values were plotted on a logarithmic scale against the log2 of the total RNA concentration. All amplicons had a good doubling efficiency and demonstrated a linear relationship with total RNA concentration (Table 2). Resulting standard curves were used to calculate the relative amount of Arntl and Per2 in the ovary. Relative amount was converted to percent of the maximum amount for each experiment. Percentages were averaged, and SEM was calculated for each group. Data are presented as the mean ± SEM.
The standard curve for Arntl was: y = 1.2171x + 33.517 (R2 = 0.9924; E = 1.77); and for Per2, y = 1.1343x + 32.717 (R2 = 0.9924, E = 1.84); where E = efficiency of amplicon doubling with each cycle 2(1/slope). Resulting standard curves were used to calculate relative amounts of Arntl and Per2 according to the equation: relative amount = 2([Ct Intercept]/[slope]).
Histological Analysis of mRNA by Digoxigenin-Labeled Probes
Probes The Per2 and Arntl cDNA fragment-containing vectors were linearized with restriction enzymes and then used as templates for sense or antisense cRNA probes. Probes were generated using digoxigenin (DIG)-11-UTP (Roche) with standard protocols for cRNA synthesis. The Per2 cDNA fragment-containing vector was generously donated by Hitoshi Okamura (Kobe University School of Medicine, Kobe, Japan), and the Arntl cDNA fragment-containing vector was generously donated by Alena Sumova (Institute of Physiology, Academy of Sciences of the Czech Republic, Prague, Czech Republic).
In situ hybridization Tissues were fixed, paraffin embedded, and sectioned the same as for IHC. Briefly, after xylenes and hydration in serial ethanol, tissues were treated for 20 min at room temperature in 0.2 N HCl. Sections were treated in 2x saline sodium citrate (SSC; 1x SSC = 150 mM NaCl, 15 mM sodium citrate; pH 7.0) at 70°C, washed in 1x PBS with 0.3% Triton-X (1x PBS = 140 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4; pH 7.2), digested in 2 µg/ml Proteinase K (0.1 M Tris buffer, 50 mM EDTA; pH 8.0) for 15 min at 37°C, washed in 1x PBS with 100 mM Glycine, postfixed in 4% paraformaldehyde, treated in 0.25% acetic anhydride in 0.1 M triethanolamine for 10 min, and then washed in 2x SSC. After an incubation step in 4x SSC/50% deionized formamide for 1 h at 42°C, sections were incubated in hybridization solution (Sigma) containing either DIG-labeled Per2 or Arntl cRNA for 16 h at 42°C. Posthybridization rinses were performed in 2x SSC/50% formamide at 42°C, and then with decreasing concentrations of SSC. Sections were prepared for alkaline phosphatase staining by several washes in wash buffer (100 mM Tris-HCl, 150 mM NaCl; pH 7.5) and then in blocking buffer (1% heat-inactivated normal sheep serum in wash buffer). They were then incubated in the alkaline phosphatase-conjugated antibody (Roche) diluted 1:500 in blocking buffer for 2 h at room temperature. After several washes in wash buffer, sections were pretreated for 10 min in detection buffer (100 mM Tris-HCl, 100 mM NaCl, 50 mM MgCl2; pH 9.5). Finally, the sections were incubated in a solution containing 5-bromo-4-chloro-3-indolyl phosphate toluidium salt (0.18 mg/ml) and nitroblue tetrazolium salt (0.45 mg/ml) (Roche) in detection buffer for 30 min. The coloring reaction was stopped by immersing slides in 1xPBS. The use of sense Per2 and Arntl probes revealed no specific hybridization signals in the ovary sections. Additional controls included omission of the cRNA probe, which also revealed no hybridization signals.
Immunoblot
Ovarian and SCN protein was isolated from the final phenol-ethanol fraction using the standard TRIzol protocol (Invitrogen). The SCN sample, used as a positive control for this study, consisted of pooled SCN reduced slices (n = 12) obtained at different times of day. Resuspended protein was quantified using the Micro BCA protein Assay (Pierce, Rockford, IL). Samples were heated at 100°C for 5 min with an equal volume of Laemmli sample buffer/2-mercaptoethanol (Bio-Rad) and then centrifuged at 13 000 x g for 5 min before separation by sodium dodecyl sulfate-PAGE on an 8% gel (25 µg protein/lane). After transfer (30 min at 100 V), nitrocellulose membranes were stained with 1x Ponceau (Sigma) and imaged to assess uniform loading and transfer. Milk (5%; 1 h at room temperature) was used as a blocking agent. Blots were incubated with primary antibody solution overnight at 4°C (PER2 affinity purified 5 µg/ml; ADI, San Antonio, TX; ARNTL affinity purified 5 µg/ml; ABR, Golden, CO). Blots were washed and incubated in secondary antibody (donkey anti-rabbit:horseradish peroxidase [HRP], 1:5 000; 1 h at room temperature). After a final wash, membranes were developed using chemiluminescent substrate according to the manufacturer's protocol (Pierce). Densitometry using the digital imaging software Scion Image for Windows 4.0.3 (Scion Corporation) provided quantitative analysis of data. The relative density of ARNTL and PER2 immunoreactivity was normalized against alpha-tubulin. Data were plotted as fold change from the lowest point.
Immunohistochemistry
Ovaries were fixed in 10% neutral buffered formalin for 48 h at room temperature. After paraffin embedding, 4 µm sections were cut. After deparaffinization, antigen retrieval, and blocking in normal goat serum, slides were incubated overnight in primary antibody, PER21-A (ADI) or ARNTL (ABR). PER2 antibody was diluted to 20 µg/ml and ARNTL antiserum was diluted to 10 µg/ml. Sections were washed and incubated with an HRP-labeled secondary antibody for 1 h (Vector Laboratories, Burlingame, CA). The peroxidase antibody complex was visualized using aminoethyl carbazole substrate (Zymed, San Francisco, CA). Control experiments included omission of the primary antibody and preabsorption with the peptide specific to the primary antibody. Sections were counterstained with hematoxylin.
Statistical Analysis
All statistical analyses were performed in consultation with a statistician (Dr. David Schaeffer, University of Illinois). Least squares ANOVA was used to analyze all data except Arntl data from the cycle-dependence study. Tukey post hoc analysis was used when ANOVA returned a value of P
0.05. Because the Arntl qPCR data in the cycle-dependence study did not meet the criteria for homogeneity of variance in the least squares ANOVA, these data were analyzed using generalized estimated equations (GEE) model as implemented with SUDAAN version 9.0 (RTI International).
RESULTS
Clock Gene Transcript Expression in the Ovary
Clock gene products known to comprise the core oscillator (Arntl, Clock, Per1, Per2, Cry1), as well as a mammalian ortholog of the Drosophila canonical clock gene, Timeless, were expressed in the rat ovary (Fig. 1). When samples were obtained every 4 h over a 24-h period, distinct patterns of expression emerged for several of the transcripts. Estrous cycle stage was not determined in these animals. Both Arntl (7-fold; P < 0.001) and Per2 (2-fold; P < 0.05) displayed sinusoidal-like expression patterns with highly statistically significant changes over the 24-h sampling period. Ovarian expression of Arntl peaked at ZT 0/24 (Fig. 2A) and was lowest at ZT 12. Ovarian peak expression for Arntl was out of phase with the SCN (Fig. 2A, dotted line); peak expression of Arntl was observed approximately 4 h later in the ovary. Ovarian expression of Per2 peaked at ZT 16, with a nadir at ZT 4 (Fig. 2B). Similar to Arntl, Per2 expression was out of phase with the SCN (dotted line), with the ovarian peak occurring approximately 46 h later (Fig. 2B). Arntl and Per2 expression are in complete antiphase to each other in the ovary, which indicates a potential for the molecular clockworks to function within the ovary.
|
|
To examine which ovarian cell types may contribute to oscillatory expression of clock genes, we performed in situ hybridization for Arntl and Per2 in the ovary. ZT 0 and ZT 15 were selected to represent those times when the two transcripts were expected to be high and low, as indicated by qPCR data (Fig. 2). Estrous cycle stage was not assessed in these animals. Arntl transcripts were present in granulosa and theca layers of growing and antral follicles, as well as corpora lutea and stromal fibroblasts (Fig 3, AD). When comparisons were made between ZT 0 and ZT 15, it was apparent that Arntl transcripts were elevated in granulosa cells, and particularly in theca layers of growing and antral follicles at ZT 0. Levels of Arntl transcripts did not change in corpora lutea or in the ovarian stroma. Similarly, Per2 transcripts were present in granulosa and theca layers of both growing and antral follicles (Fig. 3, EH). Staining for Per2 was also observed in ovarian stroma and corpora lutea (not shown). Per2 transcripts were higher in granulosa and theca cells of growing follicles at ZT 15 (Fig. 3, E and F). Although less apparent, granulosa and theca cells of antral follicles also appeared elevated at ZT 15 (Fig. 3, G and H). Positive controls from rat SCN demonstrate the specificity of the probes. Negative controls, including no probe and sense probe, showed no staining (data not shown). These data support the qPCR results, demonstrating that Arntl and Per2 transcripts are expressed in an out-of-phase relationship in the ovary. Furthermore, these data indicate that Arntl and Per2 are expressed in the same ovarian cell types.
|
Clock Gene Protein Expression in the Ovary
Immunoblot analysis was performed to determine the protein profile for the clock genes ARNTL and PER2 in the rat ovary. Estrous cycle phase was not determined in these animals. Analysis of protein extracts isolated from the same samples as used for transcript analysis revealed significant 24-h variation in the protein products ARNTL and PER2 (Fig. 4). Peak protein levels for ARNTL (Fig. 4, A and B) were observed at ZT 24/0 (ANOVA, P < 0.01 compared to ZT 12) with minimal levels occurring at ZT 1216. PER2 levels peaked at ZT 1216 (ANOVA, P < 0.01 compared to ZT 24/0), with a nadir at ZT 24/0 (Fig. 4, A and C). As previously demonstrated for Arntl and Per2 transcripts, peak expression of ARNTL and PER2 proteins was in complete antiphase.
|
IHC was performed to confirm that cellular localization of ARNTL and PER2 proteins in the rat ovary corresponds to that previously observed for the transcripts (Fig. 3). ZT 4 was selected because both proteins should be detectable at this time at an intermediate level. Estrous cycle phase was not determined in these animals. Immunoreactivity indicated that both proteins were highly expressed in specific ovarian structures (Fig. 5, A and B). Staining was absent or severely reduced upon omission of primary antibody (data not shown) or preabsorption of primary antibody against its peptide (Fig. 5, C and D). ARNTL immunoreactivity was prominent in ovarian stromal fibroblasts (Fig. 5A). However, strong immunoreactivity for ARNTL was also observed in granulosa and theca cells of growing and antral follicles, as well as in corpora lutea. The staining pattern for PER2 was similar. PER2 immunoreactivity was also observed in granulosa and theca cells, and corpora lutea (Fig. 5B). In addition, ARNTL and, especially, PER2 immunoreactivity was strong in oocytes.
|
Estrous Cycle Changes in Clock Gene Expression
Arntl transcript expression patterns were distinctive between two estrous cycle stages. As in Figure 2, Arntl was expressed in a rhythmic pattern, with peak expression occurring near the end of night. On the day of proestrus, Arntl transcript levels were significantly elevated at ZT 18 (ANOVA with Tukey post hoc, P < 0.01) when compared to same time of day on the day of Diestrus I (Fig. 6A). Notably, ZT 18 is approximately 810 h after the expected ovulatory LH surge on the day of proestrus. In contrast, Per2 gene expression was not significantly different between the 2 days of the rat estrous cycle, although the pattern of expression on both days displayed a characteristic daily oscillation, with a peak in the first half of the night (Fig. 6B).
|
Effects of Gonadotropins on Clock Gene Expression
Increased expression of Arntl at ZT 18 on the day of proestrus suggests that LH may play a role in regulating expression of this transcript. We used the hypophysectomized, immature rat model to test this hypothesis. Expression patterns and levels of Arntl and Per2 gene products were altered by hCG treatment in juvenile, hypophysectomized rats primed with eCG (control) (Fig. 7). In control animals treated with eCG alone, Arntl transcripts did not change with time of day over the 24 h course of data collection. Arntl displayed highly significant differences between the treatment and control groups across the 24-h time course (P < 0.01, ANOVA). The pattern is very different in the hCG treatment group. As early as 4 h after treatment, Arntl levels show an increasing trend. Arntl transcript levels were significantly elevated compared to the initial time point in the hCG group at 8 and 12 h after treatment (P < 0.05). Similarly, the pattern of Per2 transcript expression was also altered by hCG treatment. Per2 expression in control ovaries (Fig. 7B) showed no evidence of cycling over the 24-h time course of data collection. In contrast, Per2 transcript levels were elevated only 4 h after hCG treatment, and returned to basal levels by 24 h later (ANOVA with Tukey post hoc, P < 0.05). In situ hybridization indicated hCG treatment specifically increases Arntl and Per2 transcripts within thecal layers of large follicles (Fig. 7, CF). These data strongly indicate that LH can induce cyclic expression of circadian clock genes in the hypophysectomized, immature rat model.
|
DISCUSSION
Production of fertilizable ova is a meticulously choreographed affair that requires precise timing of events between the ovary and the brain. Coordination of organismic timing is controlled by a circadian timing system, constructed from the dynamic molecular interactions of a core group of "clock" genes. As an initial step toward determining whether a circadian clock has a role in ovarian function, we examined clock gene expression patterns in the rat ovary. The results presented here provide the first report of rhythmic clock gene expression in the ovary, although a few studies have examined clock genes in other components of the female reproductive tract [37, 42]. Arntl and Per2 transcripts and proteins are expressed in appropriate antiphase relationships in the ovary, suggesting that the molecular underpinnings of the clockworks are present. Furthermore, Arntl transcript expression is different on 2 days of the estrous cyclenamely, diestrus I and proestrus. Finally, Arntl and Per2 are upregulated following hCG treatment in juvenile, hypophysectomized animals, suggesting that the circadian-regulated LH surge, required for ovulation, may provide a circadian signal to the ovary.
The ovary expresses each of the genes considered part of the core molecular clock machinery. Similar to other peripheral tissues [18, 4346] and the SCN, gene products for Arntl and Per2 cycle in the ovary with peaks in expression that are out of phase with each other. In the SCN, increased levels of ARNTL lead to formation of heterodimers with constitutively expressed CLOCK to drive transcription of the Per, Cry, and Tim genes [8]. Accumulating PER and CRY proteins form multimeric complexes that feed back to impede the function of the ARNTL/CLOCK heterodimer, thereby inhibiting their own transcription. In addition, PER proteins act as positive regulators of Arntl transcription [13, 14]. Thus, in the functional mammalian molecular clock first described in the SCN [8], the cycling "positive element," Arntl is expressed in complete antiphase to the "negative elements," Per and Cry; levels of Arntl transcript are elevated when Per and Cry are low, and vice versa. Our data indicate that antiphase expression patterns of positive and negative elements are preserved in the ovary. Furthermore, in accordance with patterns in other peripheral tissues [18, 44, 46], peak expression of these genes in the ovary is shifted by several hours compared to the SCN.
Precise coordination of timing within the hypothalamic-pituitary-ovarian axis is essential to female fertility. Because the circadian clock is fundamental for governance of organismic timing, it seems plausible that female reproductive function may be regulated by central timing mechanisms. Analysis of the circadian influences on female reproductive function has, to date, been restricted to SCN regulation of the LH surge. Required for ovulation, the LH surge occurs on the afternoon of proestrus, when a precisely timed circadian signal manifests concomitant with an appropriate steroidal environment [28, 29]. The circadian nature of this phenomenon is revealed by daily expression of an LH surge in ovariectomized animals with chronic estrogen replacement [28]. Thus, gonadotropin releasing hormone-directed production of LH offers an interesting candidate for circadian regulation of potential ovarian rhythms.
Examination of Arntl and Per2 expression over 2 days of the estrous cycle provides the first evidence that LH may be an important regulator of circadian clock gene rhythms in the ovary. Because our studies were focused on LH as a potential ovarian entrainment factor, we focused on the day of proestrus, and made comparison to another day (diestrus I) when no LH surge occurs. Times of peak expression of Arntl and Per2 were not different between diestrus I and proestrus. However, the Arntl gene product was significantly elevated at ZT 18 on the day of proestrus, approximately 810 h after the endogenous LH surge. This finding led to the hypothesis that LH plays a role in regulating expression of this gene product.
We examined this hypothesis directly in the juvenile, hypophysectomized rat. In contrast to the adult, cycling rat, hypophysectomized, juvenile rats, which received eCG to stimulate follicle development and LH responsiveness, showed no evidence of circadian expression patterns for the canonical clock genes, Arntl and Per2. Per2 transcripts were tonically expressed at intermediate levels, whereas, under these same conditions, Arntl transcripts were also arrhythmic, albeit at constitutively low levels. These results are consistent with a lack of circadian clock function, and are reminiscent of the Per2 and Arntl levels observed in the completely arrhythmic mice lacking both Cry1 and Cry2 [13, 47, 48]. Thus, in hypophysectomized, juvenile animals, which have never experienced an LH surge, there is no apparent circadian rhythm in the ovary, and levels of Arntl and Per2 expression mimic those observed in the SCN of known circadian clock mutants. In contrast, hCG treatment induced a circadian pattern of expression for both genes. Furthermore, peak expression of Per2 and Arntl were out of phase with each other. Similar to the estrous cycle study, Arntl was significantly elevated 812 h after hCG treatment, suggesting that LH may upregulate expression of the Arntl transcript.
The expression profile for the Per2 transcript was also altered by hCG stimulation in juvenile, hypophysectomized animals. Per2 was elevated 4 h after hCG treatment. Although these results may initially appear contradictory with the findings in the estrous cycle study, it is notable that hCG was administered when Per2 transcript levels were low. In contrast, the endogenous LH surge occurred when Per2 levels were already approaching their maximum; thus, further stimulation of this gene product beyond peak levels may not have been possible. Similarly, in the SCN, light pulses will stimulate Per2 transcripts only when Per2 is endogenously low [49].
Generally, the data from the juvenile, hypophysectomized animals suggest that LH may provide a sufficient signal to initiate clock gene rhythmicity in the ovary. In this model, Per2 was elevated first, followed by an appropriately out-of-phase increase in Arntl after several hours. In the current interlocking feedback loop model of clock gene function, Per2 can drive transcription of Arntl. Whether LH acts through Per2 or a different mechanism to drive ovarian rhythms requires further investigation. Our data highlight the usefulnees of the juvenile, hypophysectomized animal as a model for endocrine induction of circadian rhythmicity in a peripheral tissue.
Induction of Arntl and Per2 transcripts by hCG was particularly evident in the theca layers of growing and antral follicles. Thus, the theca cells may be important for receiving the circadian signal and initiating rhythmicity in developing follicles. If LH is the agent that carries the circadian signal to the ovary, the theca layer is primed to receive this message because LH receptor expression is restricted to the thecal cell compartment in developing follicles [50, 51]. Furthermore, FSH can increase LH receptor expression in granulosa cells [51]. Further experimentation is required to determine the mechanisms underlying increased expression of Per2 and Arntl in response to LH treatment.
Many questions remain regarding the functional significance of circadian clock genes within the ovary. Mice bearing deletions of the negative arm of the circadian transcriptional/translational feedback loop (the Cry and Per genes) have no apparent deficits in reproductive fitness. In contrast, disruption of positive elements (Arntl and Clock) is detrimental to reproduction. Arntl-null mice have irregular estrous cycles, decreased ovulation rates, and fail to maintain pregnancy to term [42]. Clock mutant mice also have impaired fertility [5254]. Although superficial histological examination of Clock mutant ovaries has revealed no significant differences between numbers of corpora lutea and preovulatory follicles with wild-type tissues [52], comprehensive analysis of the ovary has not been completed in any circadian clock mutant animals. IHC has revealed expression of ARNTL and PER2 in ovarian follicles throughout the latter stages of follicular development and in corpora lutea, as well as in oocytes (Fig. 4). It is clear that successful follicular development requires coordinated communication between the oocyte and the surrounding somatic cells [55]. Thus, it is attractive to hypothesize that clock gene expression may play a role in directing follicular development as the oocyte itself orchestrates the rate of ovarian follicular development by conducting the growth and maturation of the surrounding granulosa and theca layers [56].
Notably, another basic helix-loop-helix, PAS domain protein similar to circadian clock proteins, the aryl hydrocarbon receptor (AhR), is important for normal ovarian follicular development [57]. AhR protein expression patterns in the ovary are similar to ARNTL and PER2 [58]. AhR has been linked to the rate of apoptosis in the developing ovarian germ line, the rate of follicle development, and, ultimately, the rate of ovulation. Likewise, another member of the basic helix-loop-helix (bHLH) proteins, SHARP2, is expressed in granulosa and theca interna layers of the developing follicles and is also stimulated by hCG [59]. It is not unreasonable to speculate that clock genes could be involved in follicular growth and maturation through potential interactions with AhR, SHARP2, and other bHLH containing transcription factors within the ovarian follicle. Finally, interactions between clock genes and the cell cycle may be reflected in expression of these genes in granulosa and theca cells of growing follicles [6063]. The results of this study suggest that clock genes have the potential to influence ovarian follicular function, or, at the very least, may assure maintenance of the delicate balance of timing within the female reproductive axis.
Physiologic circadian rhythmicity has long been considered dependent upon the pacemaker activity of the master clock located in the SCN [3, 4, 7], although this concept has recently been the subject of debate [26, 64]. SCN lesion disrupts behavioral and endocrine circadian rhythms [65, 66], as well as oscillations in peripheral clocks [22, 25]. Transplant studies indicate that the SCN determines the circadian period in the host [67]. We hypothesize that an SCN-derived LH surge provides timing cues to the ovary. Investigation of the effects of LH on other reproductive tissues, as well as nonreproductive tissues, may provide additional insight into the effects of LH on clock gene expression. However, alternate timing cues should also be considered. Transneuronal viral-tracing techniques have revealed that the SCN has direct neuronal connections with the liver and pancreas through the autonomic nervous system [68]. Circadian timing cues could also reach the ovary through direct neuronal connections from the SCN via the autonomic nervous system [69]. Continuing research efforts are aimed at determining the role of circadian clock genes in the ovary, as well as exploring whether the SCN ultimately provides timing cues for synchronization of ovarian function within the female reproductive axis.
|
ACKNOWLEDGMENTS
The authors are grateful to M. Nakai for technical assistance with IHC, D. Schaeffer for assistance with statistical analysis, J. Hickok for development of qPCR techniques, and J. Jorgensen and K. Bottum for comments on the manuscript.
FOOTNOTES
1 Supported by National Institutes of Health grant ES012948 and the Illinois Governor's Venture Technology Fund (to S.A.T.). ![]()
2 Correspondence: Shelley A. Tischkau, Department of Veterinary Biosciences, University of Illinois at Urbana-Champaign, 3840 VMBSB, 2001 S. Lincoln Ave., Urbana, IL 61802. FAX: 217-244-1652; tischkau{at}uiuc.edu ![]()
Received: 4 January 2006.
First decision: 27 January 2006.
Accepted: 26 June 2006.
REFERENCES
This article has been cited by other articles:
![]() |
M. E. Andersen, H. J. Clewell III, E. Bermudez, G. A. Willson, and R. S. Thomas Genomic Signatures and Dose-Dependent Transitions in Nasal Epithelial Responses to Inhaled Formaldehyde in the Rat Toxicol. Sci., October 1, 2008; 105(2): 368 - 383. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mukai, T.-M. Lin, R. E. Peterson, P. S. Cooke, and S. A. Tischkau Behavioral Rhythmicity of Mice Lacking AhR and Attenuation of Light-Induced Phase Shift by 2,3,7,8-Tetrachlorodibenzo-p-Dioxin J Biol Rhythms, June 1, 2008; 23(3): 200 - 210. [Abstract] [PDF] |
||||
![]() |
N. Nakao, S. Yasuo, A. Nishimura, T. Yamamura, T. Watanabe, T. Anraku, T. Okano, Y. Fukada, P. J. Sharp, S. Ebihara, et al. Circadian Clock Gene Regulation of Steroidogenic Acute Regulatory Protein Gene Expression in Preovulatory Ovarian Follicles Endocrinology, July 1, 2007; 148(7): 3031 - 3038. [Abstract] [Full Text] [PDF] |
||||
![]() |
P.-J. He, M. Hirata, N. Yamauchi, S. Hashimoto, and M.-a. Hattori The disruption of circadian clockwork in differentiating cells from rat reproductive tissues as identified by in vitro real-time monitoring system J. Endocrinol., June 1, 2007; 193(3): 413 - 420. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mukai and S. A. Tischkau Effects of Tryptophan Photoproducts in the Circadian Timing System: Searching for a Physiological Role for Aryl Hydrocarbon Receptor Toxicol. Sci., January 1, 2007; 95(1): 172 - 181. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||