Testing the circadian gene hypothesis in prostate cancer: a population-based case-control study

doi:10.1158/0008-5472.CAN-09-0648

Cancer Res. Author manuscript; available in PMC 2010 December 15.

Published in final edited form as:

Cancer Res. 2009 December 15; 69(24): 9315–9322.

doi: 10.1158/0008-5472.CAN-09-0648

PMCID: PMC2955869

NIHMSID: NIHMS152966

Testing the circadian gene hypothesis in prostate cancer: a population-based case-control study

Yong Zhu,¹ Richard G. Stevens,² Aaron E. Hoffman,¹ Liesel M. FitzGerald,³ Erika M. Kwon,⁴ Elaine A. Ostrander,⁴ Scott Davis,³ Tongzhang Zheng,¹ and Janet L. Stanford³

¹Department of Epidemiology and Public Health, Yale University School of Medicine, New Haven, CT

²Department of Community Medicine and Health Care, University of Connecticut Health Center, Farmington, CT

³Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA

⁴National Human Genome Research Institute, Cancer Genetics Branch, National Institutes of Health, Bethesda, MD

Correspondence to: Yong Zhu, PhD, Department of Epidemiology and Public Health, Yale University School of Medicine, New Haven, CT 06520, Phone: (203) 785-4844; Fax: (203) 737-6023; Email: yong.zhu/at/yale.edu; and Janet L. Stanford, PhD, Epidemiology Program, Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA. Phone: (206) 667-2716 Fax: (206) 667-4787 ; Email: jstanfor/at/fhcrc.org

The publisher's final edited version of this article is available free at Cancer Res

See other articles in PMC that cite the published article.

Abstract

Circadian genes are responsible for maintaining the ancient adaptation of a 24-hour circadian rhythm and influence a variety of cancer-related biological pathways, including the regulation of sex hormone levels. However, few studies have been undertaken to investigate the role of circadian genes in the development of prostate cancer, the most common cancer type among men (excluding non-melanoma skin cancer). The current genetic association study tested the circadian gene hypothesis in relation to prostate cancer by genotyping a total of 41 tagging and amino acid altering SNPs in ten circadian-related genes in a population-based case-control study of Caucasian men (N=1,308 cases and 1,266 controls). Our results showed that at least one SNP in nine core circadian genes (rs885747 and rs2289591 in PER1, rs7602358 in PER2, rs1012477 in PER3, rs1534891 in CSNK1E, rs12315175 in CRY1, rs2292912 in CRY2, rs7950226 in ARNTL, rs11133373 in CLOCK, and rs1369481, rs895521, and rs17024926 in NPAS2) was significantly associated with susceptibility to prostate cancer (either overall risk or risk of aggressive disease), and the risk estimate for four SNPs in three genes (rs885747 and rs2289591 in PER1, rs1012477 in PER3 and rs11133373 in CLOCK) varied by disease aggressiveness. Further analyses of haplotypes were consistent with these genotyping results. Findings from this candidate gene association study support the hypothesis of a link between genetic variants in circadian genes and prostate cancer risk, warranting further confirmation and mechanistic investigation of circadian biomarkers in prostate tumorigenesis.

Keywords: Circadian Genes, Prostate Cancer, Genetic Polymorphisms

Introduction

The universal 24-hour oscillation, which arises from an ancient adaptation to the rotation of the earth, is one of the fundamental components of our biological world and has a profound impact on biochemical, physiological, and behavioral processes in almost all living organisms. Disrupted circadian rhythm caused by 'light-at-night' was hypothesized over a decade ago to explain, in part, the high risk of human breast and prostate cancers in industrialized societies (1, 2). An expert panel assembled by the International Agency for Research on Cancer (IARC) recently concluded that “shift-work that involves circadian disruption is probably carcinogenic to humans” based on “limited evidence in humans” and “sufficient evidence in experimental animals for the carcinogenicity of light during the daily dark period” (3).

The role of circadian disruption in prostate tumorigenesis has been supported by observations of pilots, flight attendants, and shift workers, who are more likely to have disrupted circadian cycles due to abnormal work hours. Studies from pilot cohorts have shown an excess risk of prostate cancer (4–6) that increases with increasing number of flight hours (6). Similarly, male rotating-shift workers have a higher risk of prostate cancer as compared with day- or night-only workers (7, 8). Furthermore, studies of the blind who have a reduced or abolished sensitivity to environmental light exposures have observed a reduced risk of prostate cancer (9, 10), which provides additional evidence suggesting a role for circadian disruption in the etiology of prostate cancer.

Circadian rhythms are controlled and maintained by several core circadian genes via transcription-translation feedback loops that include positive activators, e.g., CLOCK, neuronal PAS domain protein 2 (NPAS2), and aryl hydrocarbon receptor nuclear translocator-like (ARNTL), negative effectors, e.g., two cryptochrome (CRY1 and CRY2) and three period (PER1, PER2, and PER3) genes, and modulators (casein kinase 1, ε (CSNK1E); see Kondratov et al. (11) for an extensive discussion on the molecular organization of the circadian system). Emerging data have demonstrated that circadian genes are also involved in regulating cancer-related biological pathways such as cell proliferation and apoptosis by controlling the expression of tumor suppressor genes, cell cycle genes, and genes encoding caspases and transcription factors (reviewed in (12)).

Given the potential role of circadian genes in tumorigenesis, it has been hypothesized that genetic variations in these genes could be associated with an individual’s susceptibility to cancer, especially hormone-related cancers such as prostate cancer (13). In the current study, we tested the circadian gene hypothesis by systematically genotyping genetic variants in 10 circadian genes in subjects from two population-based case-control studies of prostate cancer.

Subjects, Materials and Methods

Study Population

The study population consists of participants from two population-based case-control studies in residents of King County, Washington (Study I and Study II), which have been described previously (14, 15). Briefly, subjects diagnosed with histologically confirmed prostate cancer were ascertained from the Seattle-Puget Sound SEER cancer registry. In Study I, cases were diagnosed between January 1, 1993, and December 31, 1996 and were 40–64 years of age at diagnosis. In Study II, cases were diagnosed between January 1, 2002, and December 31, 2005 and were 35–74 years of age at diagnosis. Overall, 2,244 eligible prostate cancer patients were identified, 1,754 (78.2%) were interviewed, and blood samples yielding sufficient DNA for genotyping were drawn from 1,457 (83.1%) cases who completed a study interview. Demographic and clinical characteristics did not differ between cases that did and did not provide a DNA sample. A comparison group of controls without a self-reported physician’s diagnosis of prostate cancer was identified using one-step random digit telephone dialing. Controls were frequency matched to cases by five-year age groups and recruited evenly throughout each ascertainment period for cases. A total of 2,448 men were identified who met the eligibility criteria, 1,645 (67.2%) were interviewed and blood samples were drawn and DNA prepared from 1,352 men (82.2%) using standard protocols. Demographic characteristics did not differ between controls that did and did not provide a DNA sample. For the current analyses, only Caucasian participants with DNA were included (1,308 cases and 1,266 controls).

Data Collection

Subjects in both studies completed in-person interviews conducted by trained male interviewers using similar standardized questionnaires. Questions pertained to the time prior to reference date, which was the date of prostate cancer diagnosis for cases and for controls, was a pre-assigned random date that approximated the distribution of cases’ diagnosis dates. The questionnaire collected information on social and demographic factors, family structure and cancer history, and medical history including prostate cancer screening history. Clinical information on cases, including Gleason score, stage of disease and serum PSA level at diagnosis was obtained from the cancer registry. All study procedures were reviewed and approved by the Fred Hutchinson Cancer Research Center institutional review board and written informed consent was obtained from all study subjects before participation.

SNP Selection and Genotyping

SNPs were identified using the Tagger algorithm (16), which is implemented in the Haploview interface (17) of HapMap’s genome browser, Release 22. SNP sets were identified to capture all of the variation found within each circadian gene with r² ≥ 0.8. In addition to the tagging SNPs, seven SNPs that cause amino acid changes were included in the genotyping panel. We also examined data from the National Cancer Institute (NCI) Cancer Genetic Markers of Susceptibility (CGEMS) project (http://cgems.cancer.gov/) (18) and of the eight SNPs found to be significantly associated with prostate cancer risk in CGEMS (global P<0.05), six were included in the current study. Due to SNP compatibility limitations associated with the genotyping assay, two were unable to be included in the current study. Genotyping for all SNPs was performed at the National Human Genome Research Institute (Dr. Ostrander’s laboratory) using the SNPlex Genotyping System (Applied Biosystems, Inc., Foster City, CA) according to the manufacturer’s protocol. The details of this assay have been described previously (19, 20). Briefly, the genotyping system employs multiplex PCR technology and capillary electrophoresis in order to determine the genotypes of up to 48 different SNPs per SNPlex array. The GeneMapper software package (Applied Biosystems) was used to assign genotypes for each SNP. Replicate samples (n=141) were interspersed throughout all genotyping batches, and the concordance levels for blind duplicate samples were > 99% for all SNPs assayed. All genotyping scores, including quality control data, were re-checked by different laboratory personnel and the accuracy of each assay was confirmed.

Statistical Analysis

All statistical analyses were performed using SAS Version 9.1 (SAS Institute, Cary, NC), unless otherwise noted. A chi-square test was used to test for departures from Hardy-Weinberg equilibrium (HWE) for each SNP in the control population. SNPs not in HWE (p<0.01) were excluded from subsequent analyses. A goodness-of-fit chi-square test was used to determine whether the data were best described using a recessive, dominant or co-dominant model. Since no SNPs were significant under a recessive model only, and using the homozygous common plus heterozygous genotypes as the referent category did not dramatically alter the risk estimates for any SNP, all genotyping results are presented under both a co-dominant and dominant model. Adjusted odds ratios (ORs) with 95% confidence intervals (CI) and two-sided P-values were calculated by unconditional logistic regression to estimate the relative risk associated with each genotype. Adjustments were made to control for age and family history of prostate cancer in first-degree relatives. In addition, an indicator variable for study (I or II) was tested for confounding effects in the overall model, and no evidence of confounding was detected. Further, an interaction term for study and the combined variant genotypes for all significant SNPs in the overall population was tested, and no strong effect modification was identified (P>0.01). As such, the study covariate was not included in the final models. In addition to the overall associations, separate risk estimates were obtained for less aggressive and more aggressive tumors using multinomial logistic regression. Tests for interaction were also conducted using a Wald chi-square comparing the parameter estimates obtained for each genotype in less aggressive versus more aggressive cases. A binomial test was conducted to determine the probability of obtaining as many or more significant SNPs by chance alone, as were actually observed in our analysis.

Haplotype estimates were calculated by the PHASE program (21) using the genotyping information from all markers in each gene. Odds ratios and 95% confidence intervals for each haplotype were determined by unconditional multivariate logistic regression, using all other haplotypes as the referent category. The same covariates that were included in the main effects model were also used in the haplotype model. Separate risk estimates were not obtained for haplotypes with a frequency less than 1%, but these individuals remained in the analysis as members of the reference group.

Results

Compared to controls, prostate cancer cases reported a higher proportion of first-degree relatives with a history of prostate cancer, a greater frequency of prostate-specific antigen (PSA) screening and a higher serum PSA level (p<0.001; Table 1). Genotyping was completed for 44 SNPs in nine core circadian genes and one circadian related gene, TIMELESS. A list of these SNPs, including the minor allele frequency (MAF) detected in this study population and p-values for HWE, is presented in Table 2. The MAF of each SNP was very similar to that found in the CEPH HapMap population (www.hapmap.org/), which like this study, is population-based. Three SNPs (rs11133376 in CLOCK, rs1534891 in CSNK1E and rs2291739 in TIMELESS) did not fit HWE (P<0.01) and were thus excluded from further analysis.

Table 1

Distribution of selected characteristics of prostate cancer cases and population-based controls

Table 2

List of genetic variations in circadian genes included in the analysis

Genotypes of circadian genes and prostate cancer risk

Of the remaining 41 SNPs, 17 were significantly associated with risk of prostate cancer (either overall risk or risk of more aggressive disease), whereas ~ 6 (41 × 0.05 × 3) would be expected by chance alone (P=0.0001; Table 3). Specifically, assuming a dominant model three SNPs located in NPAS2 (rs895521: OR=0.83; 95% CI, 0.70–0.97; rs1369481: OR=0.81; 95% CI, 0.69–0.95; and rs17024926: OR=1.25; 95% CI, 1.07–1.47) and rs7950226 in ARNTL (OR=1.22: 95% CI, 1.02–1.46) were significantly associated with prostate cancer risk, as was the homozygous rare genotype of rs1534891 in CSNK1E (OR=2.65: 95% CI, 1.16–5.95). We also analyzed these SNPs for interaction with family history of prostate cancer in a first-degree relative. Only one SNP, rs17024926 in NPAS2, showed significant interaction (P for interaction=0.013 for the group with any variant allele), with a significant association observed only among men with no family history (OR=1.36: 95% CI, 1.15–1.63). In addition, since heterodimerization between the protein products of ARNTL and NPAS2 is required for proper functioning, and each of these genes harbored at least one significant SNP, we tested for interactions among each of the significant variants in these genes (rs7950226 in ARNTL combined with rs1369481, rs895521, or rs17024926 in NPAS2). In each case, elevated relative risks were observed for individuals carrying the risk allele at both loci, but the interactions were not statistically significant.

Table 3

Significant associations between genetic variants in circadian genes and prostate cancer risk

In addition to the main effects model, a multinomial regression analysis was conducted to determine whether any of the circadian gene variants were associated with a composite variable indicating aggressive disease. For these analyses, controls were compared simultaneously to prostate cancer cases stratified into two groups: 1) less aggressive (localized stage, Gleason score ≤7 (3+4), and a diagnostic PSA level <20 ng/mL); and 2) more aggressive (regional or distant stage, or Gleason score ≥7 (4+3) or a diagnostic PSA level ≥20 ng/mL). These results showed that eight SNPs were significantly associated with risk of less aggressive prostate cancer, including five assuming a dominant model: rs1012477 in PER3 (OR=1.25; 95% CI, 1.03–1.52), rs2292912 in CRY2 (OR=0.82; 95% CI, 0.69–0.99), rs7950226 in ARNTL (OR=1.22; 95% CI, 1.00–1.49), and SNPs rs17024926 and rs1369481 in NPAS2 (OR=1.28; 95% CI, 1.07–1.53; OR=0.80; 95% CI, 0.67–0.96, respectively); two assuming a recessive model: rs1534891 in CSNK1E (OR=3.09; 95% CI, 1.32–7.21) and rs12315175 in CRY1 (OR=1.55; 95% CI, 1.01–2.39); and, finally, the heterozygous genotype of rs7602358 in PER2 (OR=1.24; 95% CI, 1.03–1.50). Four SNPs were significantly associated with risk of more aggressive prostate cancer, including two assuming a dominant model: rs11133373 in CLOCK (OR=0.79; 95% CI, 0.63–0.99), and rs895521 in NPAS2 (OR=0.79; 95% CI, 0.62–0.99); and the homozygous rare genotypes of two SNPs in PER1 [rs885747 (OR=0.71; 95% CI, 0.51–0.99), and rs2289591 (OR=1.70; 95% CI, 1.08–2.66)]. In total, the effect of one or more genotypes in four SNPs varied significantly according to disease aggressiveness (P for interaction<0.05), including two SNPs in PER1 (rs2289591 and rs885747), one SNP in PER3 (rs1012477), and one SNP in CLOCK (rs11133373). Of note, the effect of all five SNPs that were significantly associated with prostate cancer risk in the full population did not vary significantly by aggressiveness. Genotyping results, including prostate cancer risk estimates for each of the non-significant SNPs, are provided in Supplementary Table 1.

Haplotypes of circadian genes and prostate cancer risk

In addition to the single-SNP analyses, haplotypes were estimated using the genotyping data from all SNPs in each gene. The haplotype analysis identified seven haplotypes in two different genes that were significantly associated with prostate cancer risk in the full population (Table 4). These included four different haplotypes in ARNTL, with a combined frequency of 13.94% in the population, as well as three haplotypes in NPAS2 with a combined frequency of 16.11%. The effect of one common haplotype in PER1 (24.82%), varied significantly according to aggressiveness, as it was significantly associated with prostate cancer risk in the more aggressive disease group (OR=1.23; 95% CI, 1.03–1.48, P for interaction<0.05) but not in the less aggressive disease group (OR=0.96; 95% CI, 0.83–1.11).

Table 4

Significant associations between haplotypes of circadian genes and prostate cancer risk and aggressiveness status

Discussion

In this hypothesis driven association study of ten candidate circadian genes, three, ARNTL, CSNK1E, and NPAS2, were significantly associated with susceptibility to prostate cancer in the full population, and SNPs in an additional six genes, CLOCK, CRY1, CRY2, PER1, PER2, and PER3, were significantly associated with a subgroup comprised of less or more aggressive cancers. These results support the hypothesis of a link between genetic variations in circadian genes and prostate cancer, and are consistent with previous genetic association findings from breast cancer (22, 23), prostate cancer (24), and non-Hodgkin’s lymphoma (NHL) (21).

It is not unexpected that alterations in these genes might have a broad impact on cancer susceptibility since as many as 10% of all genes in the mammalian genome are under some form of circadian control, including a number of transcripts with relevance to tumorigenesis (25). In addition, the molecular clockwork of the circadian system operates via interacting feedback loops among each of the core genes. As such, disturbances in one gene could result in phenotypic consequences across the entire system. As a result of this interconnectedness, if alterations to one or a few of the circadian genes were in fact causally associated with cancer risk, this phenotype could potentially be elicited by perturbing any one of the core circadian genes, with the consequence that a number of circadian gene variants could be significantly associated with disease risk without being directly causal.

Our results are also congruent with findings from a recently released NCI genome-wide association study on prostate cancer that showed some significant associations between variants in circadian genes and prostate cancer risk. The CGEMS project genotyped 550,000 SNPs in 1,182 prostate cancer cases and 1,174 controls from the Prostate, Lung, Cervical, and Ovarian Cancer Screening Trial (18). A total of 104 SNPs in circadian genes were included in the CGEMS project. Eight of these, located in four genes, NPAS2 (2 SNPs), CSNK1E (3 SNPs), CRY1 (2 SNPs), and CRY2 (1 SNP), were significantly associated with prostate cancer risk (p < 0.05 for all SNPs), including two (NPAS2 rs895521 and CRY2 rs12281674) with p ≤ 0.01. Six of these eight variants were genotyped in the current study and of these, two, CSNK1E rs1534891 and NPAS2 rs895521, were also significantly associated with prostate cancer risk. There are several possibilities as to why the results of the two studies differed. First, our study is population-based and geographically homogenous, while CGEMS is neither. Since prostate cancer rates vary widely by geography and prostate cancer is known to be a genetically heterogeneous disease, it is not surprising that different genetic variants are highlighted in these two studies. Nevertheless, both studies support a role of circadian genes in prostate cancer, which warrants further confirmation and validation in future studies.

The clock-cancer connection in prostate cancer is further supported by discoveries demonstrating that the circadian clock may function as a tumor suppressor at the systemic, cellular, and molecular levels through its involvement in cell proliferation, apoptosis, cell cycle control, and DNA damage response (26–28). For example, mice with mutant mPer2 had impaired DNA damage responses to gamma irradiation, which made them more cancer-prone (26). CSNK1E also functions in promoting cell proliferation by stabilizing β-Catenin (29). β-Catenin can interact with transcription factors of the T-cell-specific transcription factor/lymphoid enhancer factor-1 (TCF/LEF) family to regulate transcription (30) and promote tumorigenesis (31). PER1 also plays an important role in regulating growth and DNA damage control, and it interacts with proteins in the cell-cycle pathway (32). A recent study also demonstrated that the circadian gene NPAS2 affects several gene pathways associated with the DNA damage response, including cell cycle, apoptosis, and DNA repair (33). As such, polymorphisms in circadian genes may change their protein expression and activity, thereby affecting downstream cancer-related biological pathways.

Apart from their tumor suppressor role, circadian genes might affect prostate cancer risk by mediating levels of serum androgen (13). This in turn may act as a carcinogen by increasing cellular proliferation, thereby increasing the chance of random DNA copy errors. The prostate is an androgen-dependent organ, since androgens are essential for its normal growth and maintenance. Animal studies have shown that large amounts of androgens given to rodents can induce prostate cancer (34). A human study has also observed that long-term exposure to high levels of androgens can promote the development of prostate cancer (35), although a large pooled analysis of 18 prospective studies did not support this finding (36). Further, men taking the drug finasteride, which reduces androgen exposure within the prostate, compared to men on placebo were shown to have about a 25% lower prevalence of prostate cancer over a seven year period (37). Levels of testosterone, the principle circulating androgen, have also been shown to play a pivotal role in the differentiation and maintenance of prostate cancer cells (38). The impact of circadian rhythms on the expression of androgens has been previously documented in animals (39). In men, a circadian pattern in the levels of serum total testosterone has also been observed, and this circadian rhythm becomes blunted with normal aging (40). Circadian rhythms in plasma levels of cortisol, dehydroepiandrosterone, delta 4-androstenedione, testosterone and dihydrotestosterone have been detected in healthy young men (41). In addition, circadian gene variants have been correlated to alterations in serum sex hormone levels, including testosterone and sex hormone binding globulin (42). These observations demonstrate the potential impact of the circadian rhythm on androgen expression, and suggest that circadian gene variants may alter the risk of hormonally-mediated cancers by influencing the regulation of hormone pathways.

Among all ten circadian genes examined, NPAS2 gives the most robust association with prostate cancer risk. Importantly, NPAS2 has been linked to the risk of two other cancer types: hormone-related breast cancer (23) and immune-related NHL (43). The possible involvement of NPAS2 in tumorigenesis has been supported by an animal study demonstrating that ARNTL/NPAS2 heterodimers directly suppress transcription of the oncogene c-myc (26). A recent in vitro study further showed that the circadian gene NPAS2 affects several pathways associated with the DNA damage response, including the cell cycle, apoptosis, and DNA repair (33). Although the explicit mechanisms underlying these associations remain unclear, there is strong evidence indicating that NPAS2 may serve as a biomarker for an individual’s risk of cancers. As such, further study is needed in order to elucidate the functional significance of NPAS2 variants, and to further understand their role in carcinogenesis.

In conclusion, the genotype and haplotype approaches, using both tagging and amino acid altering SNPs in this population-based case-control study, provide comprehensive evidence supporting a role for the molecular clockwork pathway in prostate tumorigenesis. The strengths of this study include the relatively large sample size available for genotyping and the novel findings associated with the emerging area of circadian genetics in cancer. In addition, while previous epidemiologic studies have noted a link between circadian disruption and prostate cancer risk, no previous study has systematically examined the effects of variants in circadian genes on prostate cancer susceptibility. A major limitation is one that is shared with most genetic epidemiology studies, i.e., the lack of functional data relating the variants under study to a well-defined phenotypic effect. Nevertheless, the genetic polymorphisms identified in this study will build a molecular basis for future mechanistic studies and phenotypic assessments of circadian genes in relation to prostate cancer. Given the newly identified role of the circadian clock in tumorigenesis, unraveling the association between circadian genes and prostate cancer will add to our understanding of a fundamental aspect of cellular processes in tumorigenesis, and may facilitate the development of novel risk and prognostic biomarkers for prostate cancer.

Supplementary Material

Click here to view.^{(323K, doc)}

Acknowledgments

This work was supported by the National Institutes of Health (grants R01 CA56678, R01 CA082664, R01 CA092579, R01 CA122676, and R03 CA110937), with additional support from the Fred Hutchinson Cancer Research Center, Yale University and the National Human Genome Research Institute.

References

1. Stevens RG. Electric power use and breast cancer: a hypothesis. Am J Epidemiol. 1987;125:556–561. [PubMed]

2. Stevens RG, Davis S, Thomas DB, Anderson LE, Wilson BW. Electric power, pineal function, and the risk of breast cancer. Faseb J. 1992;6:853–860. [PubMed]

3. Straif K, Baan R, Grosse Y, et al. Carcinogenicity of shift-work, painting, and fire-fighting. The Lancet Oncology. 2007;8:1065–1066. [PubMed]

4. Band PR, Le ND, Fang R, et al. Cohort study of Air Canada pilots: mortality, cancer incidence, and leukemia risk. Am J Epidemiol. 1996;143:137–143. [PubMed]

5. Irvine D, Davies DM. British Airways flightdeck mortality study, 1950–1992. Aviat Space Environ Med. 1999;70:548–555. [PubMed]

6. Pukkala E, Aspholm R, Auvinen A, et al. Cancer incidence among 10,211 airline pilots: a Nordic study. Aviat Space Environ Med. 2003;74:699–706. [PubMed]

7. Kubo T, Ozasa K, Mikami K, et al. Prospective cohort study of the risk of prostate cancer among rotating-shift workers: findings from the Japan collaborative cohort study. Am J Epidemiol. 2006;164:549–555. [PubMed]

8. Conlon M, Lightfoot N, Kreiger N. Rotating shift work and risk of prostate cancer. Epidemiology. 2007;18:182–183. [PubMed]

9. Pukkala E, Ojamo M, Rudanko SL, Stevens RG, Verkasalo PK. Does incidence of breast cancer and prostate cancer decrease with increasing degree of visual impairment. Cancer Causes Control. 2006;17:573–576. [PubMed]

10. Feychting M, Osterlund B, Ahlbom A. Reduced cancer incidence among the blind. Epidemiology. 1998;9:490–494. [PubMed]

11. Kondratov RV, Gorbacheva VY, Antoch MP. The role of mammalian circadian proteins in normal physiology and genotoxic stress responses. Curr Top Dev Biol. 2007;78:173–216. [PubMed]

12. Fu L, Lee CC. The circadian clock: pacemaker and tumour suppressor. Nat Rev Cancer. 2003;3:350–361. [PubMed]

13. Zhu Y, Zheng T, Stevens RG, Zhang Y, Boyle P. Does "Clock" Matter in Prostate Cancer? Cancer Epidemiology, Biomarkers and Prevention. 2006;15:3–5. [PMC free article] [PubMed]

14. Stanford JL, Wicklund KG, McKnight B, Daling JR, Brawer MK. Vasectomy and risk of prostate cancer. Cancer Epidemiol Biomarkers Prev. 1999;8:881–886. [PubMed]

15. Agalliu I, Salinas CA, Hansten PD, Ostrander EA, Stanford JL. Statin use and risk of prostate cancer: results from a population-based epidemiologic study. Am J Epidemiol. 2008;168:250–260. [PMC free article] [PubMed]

16. de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, Altshuler D. Efficiency and power in genetic association studies.[see comment] Nature Genetics. 2005;37:1217–1223. [PubMed]

17. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]

18. Core Genotyping Facility, Office of Cancer Genomics, Division of Cancer Epidemiology and Genetics NCI. National Institutes of Health The Cancer Genetic Markers of Susceptibility Project. Bethesda (MD):

19. Tobler AR, Short S, Andersen MR, et al. The SNPlex genotyping system: a flexible and scalable platform for SNP genotyping. J Biomol Tech. 2005;16:398–406. [PMC free article] [PubMed]

20. De la Vega FM, Lazaruk KD, Rhodes MD, Wenz MH. Assessment of two flexible and compatible SNP genotyping platforms: TaqMan SNP Genotyping Assays and the SNPlex Genotyping System. Mutat Res. 2005;573:111–135. [PubMed]

21. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68:978–989. [PubMed]

22. Zhu Y, Brown HN, Zhang Y, Stevens RG, Zheng T. Period3 structural variation: a circadian biomarker associated with breast cancer in young women. Cancer Epidemiol Biomarkers Prev. 2005;14:268–270. [PubMed]

23. Zhu Y, Stevens RG, Leaderer D, et al. Non-synonymous polymorphisms in the circadian gene NPAS2 and breast cancer risk. Breast Cancer Res Treat. 2008;107:421–425. [PMC free article] [PubMed]

24. Chu LW, Zhu Y, Yu K, et al. Variants in circadian genes and prostate cancer risk: a population-based study in China. Prostate Cancer Prostatic Dis. 2008;11:342–348. [PubMed]

25. Storch KF, Lipan O, Leykin I, et al. Extensive and divergent circadian gene expression in liver and heart. Nature. 2002;417:78–83. [PubMed]

26. Fu L, Pelicano H, Liu J, Huang P, Lee C. The circadian gene Period2 plays an important role in tumor suppression and DNA damage response in vivo. Cell. 2002;111:41–50. [PubMed]

27. Matsuo T, Yamaguchi S, Mitsui S, Emi A, Shimoda F, Okamura H. Control Mechanism of the Circadian Clock for Timing of Cell Division in Vivo. Science. 2003;302:255–259. [PubMed]

28. Granda TG, Liu XH, Smaaland R, et al. Circadian regulation of cell cycle and apoptosis proteins in mouse bone marrow and tumor. Faseb J. 2005;19:304–306. [PubMed]

29. Lee E, Salic A, Kirschner MW. Physiological regulation of [beta]-catenin stability by Tcf3 and CK1epsilon. J Cell Biol. 2001;154:983–993. [PMC free article] [PubMed]

30. van de Wetering M, de Lau W, Clevers H. WNT signaling and lymphocyte development. Cell. 2002;109 Suppl:S13–S19. [PubMed]

31. Morin PJ. beta-catenin signaling and cancer. Bioessays. 1999;21:1021–1030. [PubMed]

32. Gery S, Komatsu N, Baldjyan L, Yu A, Koo D, Koeffler HP. The circadian gene per1 plays an important role in cell growth and DNA damage control in human cancer cells. Mol Cell. 2006;22:375–382. [PubMed]

33. Hoffman AE, Zheng T, Ba Y, Zhu Y. The circadian gene NPAS2, a putative tumor suppressor, is involved in DNA damage response. Mol Cancer Res. 2008;6:1461–1468. [PMC free article] [PubMed]

34. Brown CE, Warren S, Chute RN, Ryan KJ, Todd RB. Hormonally induced tumors of the reproductive system of parabiosed male rats. Cancer Res. 1979;39:3971–3976. [PubMed]

35. Gann PH, Hennekens CH, Ma J, Longcope C, Stampfer MJ. Prospective study of sex hormone levels and risk of prostate cancer. J Natl Cancer Inst. 1996;88:1118–1126. [PubMed]

36. Roddam AW, Allen NE, Appleby P, Key TJ. Endogenous sex hormones and prostate cancer: a collaborative analysis of 18 prospective studies. J Natl Cancer Inst. 2008;100:170–183. [PubMed]

37. Thompson IM, Goodman PJ, Tangen CM, et al. The influence of finasteride on the development of prostate cancer. N Engl J Med. 2003;349:215–224. [PubMed]

38. Soronen P, Laiti M, Torn S, et al. Sex steroid hormone metabolism and prostate cancer. J Steroid Biochem Mol Biol. 2004;92:281–286. [PubMed]

39. Francavilla A, Eagon PK, DiLeo A, et al. Circadian rhythm of hepatic cytosolic and nuclear estrogen and androgen receptors. Gastroenterology. 1986;91:182–188. [PMC free article] [PubMed]

40. Plymate SR, Tenover JS, Bremner WJ. Circadian variation in testosterone, sex hormone-binding globulin, and calculated non-sex hormone-binding globulin bound testosterone in healthy young and elderly men. J Androl. 1989;10:366–371. [PubMed]

41. Guignard MM, Pesquies PC, Serrurier BD, Merino DB, Reinberg AE. Circadian rhythms in plasma levels of cortisol, dehydroepiandrosterone, delta 4-androstenedione, testosterone and dihydrotestosterone of healthy young men. Acta Endocrinol (Copenh) 1980;94:536–545. [PubMed]

42. Chu LW, Zhu Y, Yu K, et al. Correlation between circadian gene variants and serum levels of sex steroids and insulin-like growth factor-I. Cancer Epidemiol Biomarkers Prev. 2008;17:3268–3273. [PubMed]

43. Zhu Y, Leaderer D, Guss C, et al. Ala394Thr polymorphism in the clock gene NPAS2: a circadian modifier for the risk of non-Hodgkin's lymphoma. Int J Cancer. 2007;120:432–435. [PMC free article] [PubMed]

44. Johansson C, Willeit M, Smedh C, et al. Circadian clock-related polymorphisms in seasonal affective disorder and their relevance to diurnal preference. Neuropsychopharmacology. 2003;28:734–739. [PubMed]

45. Lam RW, Levitan RD. Pathophysiology of seasonal affective disorder: a review. J Psychiatry Neurosci. 2000;25:469–480. [PMC free article] [PubMed]