A cross-sectional analysis of light at night, neighborhood sociodemographics and urinary 6-sulfatoxymelatonin concentrations: implications for the conduct of health studies
© Hurley et al.; licensee BioMed Central Ltd. 2013
Received: 17 July 2013
Accepted: 26 August 2013
Published: 30 August 2013
There is accumulating evidence that circadian disruption, mediated by alterations in melatonin levels, may play an etiologic role in a wide variety of diseases. The degree to which light-at-night (LAN) and other factors can alter melatonin levels is not well-documented. Our primary objective was to evaluate the degree to which estimates of outdoor environmental LAN predict 6-sulftoxymelatonin (aMT6s), the primary urinary metabolite of melatonin. We also evaluated other potential behavioral, sociodemographic, and anthropomorphic predictors of aMT6s.
Study participants consisted of 303 members of the California Teachers Study who provided a 24-hour urine specimen and completed a self-administered questionnaire in 2000. Urinary aMT6s was measured using the Bühlmann ELISA. Outdoor LAN levels were estimated from satellite imagery data obtained from the U.S. Defense Meteorological Satellite Program’s (DMSP) Operational Linescan System and assigned to study participants’ geocoded residential address. Information on other potential predictors of aMT6s was derived from self-administered surveys. Neighborhood socioeconomic status (SES) was based on U.S. Census block group data.
Lower aMT6s levels were significantly associated with older age, shorter nights, and residential locations in lower SES neighborhoods. Outdoor sources of LAN estimated using low-dynamic range DMSP data had insufficient variability across urban neighborhoods to evaluate. While high-dynamic range DMSP offered much better variability, it was not significantly associated with urinary aMT6s.
Future health studies should utilize the high-dynamic range DMSP data and should consider other potential sources of circadian disruption associated with living in lower SES neighborhoods.
KeywordsCircadian disruption Light at night Melatonin aMT6s Socioeconomic status Women
Light at night
California Teachers Study
There is growing evidence that environmental light pollution may play an etiologic role in a variety of diseases, including depression, cardiovascular disease, and cancer [1–10]. Such health effects are largely thought to be mediated via circadian disruption driven by alterations in melatonin production and secretion [11, 12]. Melatonin, an endogenous hormone produced by the pineal gland, is the primary hormonal modulator of circadian regulation in mammals and is strongly influenced by exposure to visible light, with levels peaking during the darkness of night . Predictors of melatonin levels in humans generally have not yet been well-characterized [11, 14–18]. While there is fairly strong evidence that indoor artificial light-at-night (LAN) is sufficient to suppress melatonin levels in humans, it is not known whether outdoor environmental light pollution is sufficient to exert a similar effect [1, 3, 19]. Fully understanding such predictors is a necessary prerequisite for implementing epidemiologic studies on this topic.
The primary objective of the current analyses was to evaluate the degree to which estimates of outdoor environmental light-at-night (LAN) predict levels of 6-sulfatoxymelatonin (aMT6s), the primary urinary metabolite of melatonin, among 303 California women. Secondarily, we evaluated whether selected anthropomorphic characteristics and behavioral and sociodemographic factors predict urinary concentrations of aMT6s.
Creatinine adjusted urinary aMT6s concentrations ranged from 0.58 ng/mg-creatinine to 102.06 ng/mg-creatinine with a mean of 20.66 ng/mg-creatinine, a standard deviation of 17.78 ng/mg-creatinine and an interquartile range of 17.74 ng/mg-creatinine.
Characteristics of study population (n = 303)
Age at urine collection, years (mean,standard deviation)
BMI at urine collection, kg/m2 (mean,standard deviation)
Height at CTS baseline survey, inches (mean,standard deviation)
Packyears of smoking (mean,standard deviation)
Pillyears of aspirin use (mean,standard deviation)
Coffee consumption, average grams/day (mean,standard deviation)
Coffee consumption, average times/day (mean,standard deviation)
Strenuous exercise, 3-years prior to baseline, hours/wk (mean,standard deviation)
Strenuous exercise, lifetime prior to baseline, hours/wk (mean,standard deviation)
Age at first full-term pregnancy, years (mean,standard deviation)
Parity, total number of live and still births (mean,standard deviation)
Asian & Pacific Islander
Menopausal Status, at baseline(n,%)
Premenopausal or Perimenopausal
Alcohol Consumption, year prior to baseline (n, %)
Tobacco Smoke Exposure, at baseline (n, %)
Former active smoker
Current active smoker
Oral contraceptive use, at baseline (n, %)
Antidepressant use, daily use for at least 2 months at baseline
Hormone therapy use, at baseline (n, %)
Peri/post-menopausal, never used
Peri/post-menopausal, prior use
Peri/post-menopasual, current use
Regular aspirin use, years at baseline
Regular aspirin use, days/week at baseline
Daily Use of Calcium Blockers, in 2000
Hours of sleep/night
Neighborhood Socioeconomic Status (mean, range)
Median Family Income ($)
Percent of population below poverty level
6.8 (0.0 – 37.6)
Percent of adults age 25+ with college degree
34.3 (1.5 – 79.0)
Percent of adults without a high school degree
14.6 (1.2 – 56.8)
Percent of adults employed in a professional occupation
Estimates of outdoor LAN among study participants (n = 303)
Outdoor LAN satellite values (annual average radiance)
Low-dynamic range, 2000 average (DN units)
43.4 (4 – 63)
High-dynamic range, 2006 average (scaled radiance units)
128.4 (3.8 – 465.2)
Length-of-night on day of urine collection (hours)
11 (9.2 – 14.4)
The results from our random forests analyses to identify the most important predictors of urinary aMT6s concentrations, based on the scaled variable importance values (mean divided by the standard deviation), identified the following 8 variables as potentially important predictors of aMT6s: hormone therapy use; oral contraceptive use; length-of-night; coffee consumption; strenuous exercise in the past 3 years; age; menopausal status; and neighborhood SES. Overall, however, none of these variables, individually or combined, explained a substantial proportion of the variability in measured aMT6s levels (data not shown).
Regression coefficients (β) for final model estimating log 2 creatinine-adjusted aMT6s
β, standard error
t statistic (2-sided p-value)
Ages > 55 years
Although the outdoor LAN variable was not selected through our stepwise model building approach, because it was of primary interest to us, we evaluated it separately in age-adjusted regression models and found it to have a very small inverse association with aMT6s that was not significant in either the full study population (β = −0.0028, se = 0.008, p = 0.73) or among those over age 55 (β = −0.0062, se = 0.012, p = 0.62).
Overall, results from our study failed to identify factors that substantially predict urinary measures of aMT6s. The only factors that were significantly related to aMT6s levels were age, length-of-night, and SES and together they explained very little of the variance in aMT6s levels. Nevertheless, our study does provide some important findings relevant to future investigations of health outcomes related to melatonin and/or LAN exposures.
There is increasing interest in the use of satellite imagery data to evaluate potential effects of outdoor environmental light pollution in wildlife and in human populations, including a number of studies that have investigated potential links to cancer incidence [20–22]. Our analyses demonstrated that the widely-available low-dynamic range annual satellite imagery data are insufficient for distinguishing areas with differing LAN values within suburban and urban areas. Since many health outcomes, including many cancers, are more common in urbanized areas [23–27], use of the low-dynamic range satellite data will likely be inadequate for investigating LAN exposures within such areas. While the high-dynamic range data are currently only available for 2006, efforts are underway to develop such high-dynamic range data for other years . Researchers interested in this topic would be prudent to pursue such data.
The other novel finding generated by our analyses was the significance of neighborhood SES. Our analysis showed that aMT6s concentrations tended to be low among women living in the lowest SES areas and increased with greater neighborhood SES up to a point after which the concentration leveled off and possibly declined slightly at the very highest levels of SES. Furthermore, this relationship, while statistically significant in the full study population, appeared to be of lesser importance in the women over age 55. Overall, very little is known about how melatonin varies by SES. Similar to our findings, a study of approximately 200 Canadian women reported that urinary measures of aMT6s were lower among the lowest educated women (i.e. high school diploma or less), highest in the middle category of education, with a slight decrease in the highest educated group . In contrast, Burgess and colleagues in their evaluation of predictors of salivary measures of melatonin found education was not a significant predictor .
Interpretation of our findings on SES can only be speculative. It is possible that the neighborhood measure of SES is serving as a proxy for a number of behavioral and/or anthropomorphic factors that influence melatonin levels. We did, however, incorporate variables in our model selection process for many of the most important of these factors including alcohol consumption, smoking, BMI, use of exogenous hormones and other medications, coffee consumption, and physical activity – none of which were independent and significant predictors of aMT6s. Another possibility is that low SES neighborhoods may have environmental stressors (e.g., noise, crime) other than light pollution that we did not measure in our study but could disrupt circadian rhythms and lead to lower levels of aMT6s. Given that a number of the health outcomes of high interest with respect to circadian disruption display strong SES-related risk gradients, research into other features of the built environment represents an important area for future inquiry.
In addition to SES, age and length-of-night were the only other significant predictors of aMT6s identified in our analysis. These findings are generally consistent with the limited body of literature on this topic. Two studies have reported seasonal differences in melatonin levels with higher levels associated with greater length-of-night [15, 17] and lower levels in the summer months , although two other studies reported no relationship with month  or season . Most studies have reported declines in melatonin levels with increasing age, albeit most have reported this relationship to be linear [15, 17, 18, 31–33]. Consistent with our results, two studies noted more dramatic age-related declines among older individuals [34, 35]; conversely, a few other studies suggested that the age-related declines in melatonin are greater, or are limited to, early adulthood and then level-off [18, 36–38]. Measurements of melatonin in older adults living under controlled light–dark conditions in one study reported no age-related declines in melatonin levels, leading the authors to suggest that changes in sleep behaviors might be responsible for the age-related declines reported in other observational studies . In our study, however, we observed no differences in sleep duration by age group (data not show). Our results, in the context of the somewhat conflicting literature on this, underscore the need to carefully consider the shape of the age-related risks when conducting health-related research on this topic. The fact that we observed a shift in the relationship between age and aMT6s around age 55 suggests this may be especially important when evaluating health risks in women, whose risks often change in their early- to mid-50s after menopause. This may be particularly relevant to breast cancer which exhibits a different constellation of risk relationships in pre- and post-menopausal women and for which LAN exposures have been postulated to be a potential risk factor.
Our evaluation of outdoor LAN, which showed a very modest inverse, but not statistically significant, relationship to aMT6s concentrations, was hindered by a number of obstacles. While we used the best satellite imagery data available to objectively estimate outdoor LAN values, because the high-dynamic range data was only available for 2006, the LAN estimates were not temporally congruent with the urine in which the aMT6s was measured (collected in 2000). An examination of the low-dynamic range data for all years spanning this time period (2000–2006), however, suggested relative stability in LAN values over this time for the study area. Furthermore, the Spearman rank correlation between the 2000 low-dynamic range data and the 2006 high-dynamic range data was 0.96.
Another reason for the lack of an association seen between outdoor LAN and urinary measures of aMT6s could be that while the satellite imagery data may be an adequate predictor of outdoor ambient light, it may not reflect light exposures experienced at night when participants are likely to spend the majority of their time indoors. Intervening factors, such as the use of blinds/curtains, time spent indoors, other sources of indoor light, etc. all are likely contributors to LAN exposures. The importance of this is underscored by findings from a recent study that compared estimates of LAN from satellite imagery data to calibrated photometric measurements obtained from personal monitors and reported that satellite imagery data did not correlate with personal photometric measurements . Thus, it is important that studies aimed at evaluating outdoor LAN exposures do so with full consideration of these other factors.
Finally, there are a number of limitations in the estimates of our melatonin levels that are worth noting. The urine specimens were collected over a 24-hour time period. This precludes our ability to examine the timing of peak melatonin concentration which is likely an important factor in circadian disruption . Furthermore, while there is a good deal of evidence that first morning urine captures peak night-time melatonin excretion, the degree to which 24-hour urine specimens reflect this is not well-documented [13, 31]. The aMT6s assay, however, is a well-validated biomarker, demonstrating good correlation with serum melatonin levels , and sufficient stability over time [13, 15, 31, 41, 42] to serve as a reasonable estimate of chronic levels.
There is growing evidence that circadian disruption, mediated by alterations in melatonin secretion driven by night-time light exposures may play an etiologic role in a large array of diseases. The use of satellite imagery data to estimate ambient measures of LAN offers an innovative and inexpensive source of data to test emerging hypotheses on this topic. The results from our analyses highlight some important cautionary limitations to the use of such data as well as point to some additional avenues of pursuit, including the elucidation of factors associated with neighborhood SES that may play a role in circadian disruption.
Study participants consisted of 303 members of the California Teachers Study (CTS) residing in the San Francisco Bay Area, who as part of a special substudy conducted in 2000, completed a self-administered questionnaire, an in-person interview, and supplied a 24-hour urine sample. The CTS is a prospective cohort study initiated in 1995, consisting of 133,479 female professional California school employees for whom extensive information has been collected through a baseline questionnaire and three subsequent mailed questionnaires . Members of the special substudy included a random sample of 528 CTS participants who were ≤ 85 years old when the cohort was established in 1995, and who resided in the substudy area (comprised of western Alameda, Santa Clara, San Mateo, Santa Cruz, Monterey, and northern San Benito counties). Forty-four (8%) of these women were not contacted because they had died, moved out of the substudy area, or could not be located. Of the 484 women invited to participate, 328 agreed, of whom 303 provided a 24-hour urine sample. Further details of the special substudy can be found elsewhere . The use of human subjects was reviewed and approved by the California Health and Human Services Agency, committee for the Protection of Human Subjects and the Cancer Prevention Institute of California’s Institutional Review Board.
Data and specimen collection
Participants were given urine collection kits and instructed to collect all urine produced over the next 24-hour period. Participants were asked to use cold packs to keep the sample chilled and, if willing, to refrigerate it during the 24-hour collection period. After 24-hours, staff collected and stored the samples at -20° Celsius until they could be delivered to our laboratory (within 2 weeks) where it was thawed, aliquoted into 10 mL vials, and frozen at -70° Celsius. Prior to aliquoting, total urine volume was measured and recorded.
Urinary concentration of aMT6s was measured by Pacific Toxicology Laboratory (Chatsworth, CA, USA) using Bühlmann aMT6s enzyme-linked immunosorbent assay (ELISA) kits purchased through ALPCO (Salem, NH, USA). The aMT6s ELISA is a competitive immunoassay using a capture antibody technique that has a lower detection limit of 0.8 ng/mL and an analytic high sensitivity limit of 40 ng/mL. The functional sensitivity of the assay was 1.5 ng/mL, while the intra-and inter-assay precision were 7.1% and 11.9%, respectively. Urinary creatinine levels were obtained from laboratory assays using a modified-rate Jaffe method conducted as part of a prior analysis by the University of Alabama, Birmingham in 2000. aMT6s levels were creatinine-adjusted by dividing the aMT6s concentration by the creatinine concentration.
Four samples were below the detection limit and five samples exceeded the analytic high sensitivity limit (3% of samples total). Using PROC MI with a TRANSFORM LOG(var) statement in SAS 9.3 we imputed the values below and above the limit of detection, setting range parameters of 0–0.8 ng/mL and 40–60 ng/mL respectively. The upper range of 60 ng/mL for those samples above the analytic high sensitivity limit was based upon upper levels observed in 24-hour samples from women in a previous study . We used the first imputation to assign values of aMT6s to these nine participants.
Outdoor environmental LAN data
Information on outdoor exposures to ambient LAN was derived from nighttime satellite imagery data obtained from the U.S. Defense Meteorological Satellite Program’s Operational Linescan System (DMSP-OLS) . This database contains annual composites, made after excluding the outer quarters of the satellite swath, light from the sun and moon, glare, clouds, and atmospheric lightning. Ephemeral events, such as fires, are also discarded. While these images capture only a fraction of the light originating from the earth’s surface, they accurately represent the relative levels of nighttime illumination at ground level . The imagery data is georectified to a 30 arc-second grid (equivalent to approximately one-square kilometer). Two kinds of radiance information were available: adaptive gain, low dynamic range data; and fixed-gain, high dynamic range data. The low dynamic range data are available for multiple years, but consist of unit-free (called DN, or “Digital Number” units) data with a dynamic range of six bits, i.e., integers ranging from 0 (background noise) to 63. The fixed gain, high dynamic range data, is available only for 2006 and is obtained by combining information from several fixed-gain sensors. It consists of floating point numbers (“scaled radiance” units) that, if needed, can be scaled by 1.51586 × 10-10 to obtain data that are watts/sr/cm2 (see the F16_2006_radiance_readme.txt file in reference  for caveats with respect to this scale factor). Such scaling was not performed for these analyses. Residences from study participants were assigned the average nighttime radiance value for the grid cell in which they were located using spatial analysis tools available in the raster, rgdal, and sp packages within R . Initially we assigned to residences both the low-dynamic range data for 2000 and the high-dynamic range data for 2006.
As an additional measure of LAN we generated a “length-of-night” variable, created by linking the date of urine collection with daily sunrise and sunset data for San Francisco for the year 2000 obtained from the Naval Observatory . From these data we calculated the hours of darkness on the day of urine collection.
Anthropomorphic and behavioral factors
Information on anthropomorphic and behavioral factors of interest was derived from the CTS surveys. Factors were considered as potential predictors for these analyses based on a review of the limited literature on this topic [15–18, 31–35, 37, 38, 42, 48–51] and included: age, race/ethnicity, body mass index (BMI), height, parity, number of live births, age at first full-term pregnancy, pack-years of smoking, alcohol consumption, oral contraceptive use, hormone therapy use, menopausal status, strenuous physical activity, average hours of sleep per night, age at menarche, coffee consumption, and the use of aspirin, calcium blockers and antidepressants.
Information on neighborhood socioeconomic characteristics and urbanization were derived from U.S. Census 2000 data. We used the ArcGIS, version 9.2 (ESRI, Redlands, CA, USA) to geocode the participants’ home addresses to both a latitude/longitude and to a Census 2000 block group. Data on socioeconomic status (SES) included: percentage of adults over age 25 years having completed a college degree or higher; percentage of adults without a high school degree; income (median family income), occupation (percentage of adults employed in managerial/professional occupations), and poverty (percentage of population below the poverty line). The degree of neighborhood urbanization was characterized as either urban/suburban or rural, based on a previously-developed algorithm using a combination of census block group characteristics, details of which can be found elsewhere .
Given the extensive information collected on potential predictors of melatonin levels, our initial approach was to use Breiman’s random-forest-based (RF) variable importance measures, as described by Lunetta  to identify the most important predictors of log2 transformed, creatinine adjusted aMT6s from the large number of factors for which we had information. (The aMT6s levels were highly skewed to the right, necessitating the log2 transformation.) This approach allowed us to compose a list of candidate variables, ordered by their relative importance, taking into account both multiple testing and potential interactions. The top ranking variables were then considered in linear models to further quantify predictors of aMT6s concentrations. Final model selection via R’s stepAIC function was chosen in a backwards elimination process using the Akaike information criterion (AIC) to quantify the model fit.
We first focused on determining the most appropriate model specification for age. A plot, obtained by fitting a Generalized Additive Model (GAM) with an intercept and a spline term for age, demonstrated a decrease in aMT6s with increasing age, with a possible downward curvature at age 55. An analysis of variance to compare null, linear, and quadratic linear models for the regression of the creatinine adjusted aMT6s concentrations on age, however, indicated that the best fitting model was one with age modeled as a simple linear term. Because restricting our analyses to women over age 55 resulted in a much stronger and significant linear effect for age than observed in the full study population, we conducted our subsequent analyses on the full study sample, as well as among women over the age of 55. To address the high degree of collinearity between the SES variables, we conducted a principal components analysis (PCA) to create a composite measure of SES based on the five individual components as described above. The loadings of the first principal component were then used as a composite measure in the random forests and regression models.
We express our appreciation to all the participants in the California Teachers Study and to the researchers, analysts and staff who have contributed so much for the success of this research project. We also thank Christine Duffy for coordinating the urine collection and storage, Minhthu Le for administrative support, and the CTS Steering Committee members who are responsible for the formation and maintenance of the cohort within which this study was conducted but who did not directly contribute to the current paper: Leslie Bernstein, Hoda Anton-Culver, Jessica Clague, Christina A. Clarke, Dennis Deapen, Pam Horn-Ross, James V. Lacey Jr, Yani Lu, Huiyan Ma, Susan L. Neuhausen, Hannah Park, Rich Pinder, Fredrick Schumacher, Sophia S. Wang, and Argyrios Ziogas. This research and the development of this manuscript was supported by funding from the California Breast Cancer Research Program (Grant #16IB-0071) and National Cancer Institute grants R01 CA77398 and K05 CA136967.
- Chepesiuk R: Missing the dark: health effects of light pollution. Environ Health Perspect 2009, 117:A20-A27.PubMedView Article
- Hede K: Cancer and the circadian clock: Has the time finally come? J Natl Cancer Inst 2009, 101:550–553.PubMedView Article
- Kantermann T, Roenneberg T: Is light-at-night a health risk factor or a health risk predictor? Chronobiol Int 2009, 26:1069–1074.PubMed
- Korkmaz A: Epigenetic actions of melatonin. J Pineal Res 2009, 46:117–118.PubMedView Article
- Roenneberg T, Lucas RJ: Light, endocrine systems, and cancer–a view from circadian biologists. Neuro Endocrinol Lett 2002,23(Suppl 2):82–83.PubMed
- Sahar S, Sassone-Corsi P: Metabolism and cancer: the circadian clock connection. Nat Rev Cancer 2009, 9:886–896.PubMedView Article
- Stevens RG: Circadian disruption and breast cancer: from melatonin to clock genes. Epidemiology 2005, 16:254–258.PubMedView Article
- Stevens RG: Working against our endogenous circadian clock: breast cancer and electric lighting in the modern world. Mutat Res 2009, 679:6–8.View Article
- Stevens RG: Electric light causes cancer? Surely you’re joking, Mr. Stevens. Mutat Res 2009, 682:1–6.PubMedView Article
- Stevens RG: Light-at-night, circadian disruption and breast cancer: assessment of existing evidence. Int J Epidemiol 2009, 38:963–970.PubMedView Article
- Dopfel RP, Schulmeister K, Schernhammer ES: Nutritional and lifestyle correlates of the cancer-protective hormone melatonin. Cancer Detect Prev 2007, 31:140–148.PubMedView Article
- Schernhammer ES, Schulmeister K: Melatonin and cancer risk: does light at night compromise physiologic cancer protection by lowering serum melatonin levels? Br J Cancer 2004, 90:941–943.PubMedView Article
- Mirick DK, Davis S: Melatonin as a biomarker of circadian dysregulation. Cancer Epidemiol Biomarkers Prev 2008, 17:3306–3313.PubMedView Article
- Arendt J: Melatonin: characteristics, concerns, and prospects. J Biol Rhythms 2005, 20:291–303.PubMedView Article
- Davis S, Kaune WT, Mirick DK, Chen C, Stevens RG: Residential magnetic fields, light-at-night, and nocturnal urinary 6-sulfatoxymelatonin concentration in women. Am J Epidemiol 2001, 154:591–600.PubMedView Article
- Grundy A, Sanchez M, Richardson H, Tranmer J, Borugian M, Graham CH, Aronson KJ: Light intensity exposure, sleep duration, physical activity, and biomarkers of melatonin among rotating shift nurses. Chronobiol Int 2009, 26:1443–1461.PubMedView Article
- Knight JA, Thompson S, Raboud JM, Hoffman BR: Light and exercise and melatonin production in women. Am J Epidemiol 2005, 162:1114–1122.PubMedView Article
- Schernhammer ES, Kroenke CH, Dowsett M, Folkerd E, Hankinson SE: Urinary 6-sulfatoxymelatonin levels and their correlations with lifestyle factors and steroid hormone levels. J Pineal Res 2006, 40:116–124.PubMedView Article
- Graham C, Cook MR, Gerkovich MM, Sastre A: Examination of the melatonin hypothesis in women exposed at night to EMF or bright light. Environ Health Perspect 2001, 109:501–507.PubMedView Article
- Kloog I, Haim A, Stevens RG, Barchana M, Portnov BA: Light at night co-distributes with incident breast but not lung cancer in the female population of Israel. Chronobiol Int 2008, 25:65–81.PubMedView Article
- Kloog I, Haim A, Stevens RG, Portnov BA: Global co-distribution of Light At Night (Lan) and cancers of prostate, colon, and lung in men. Chronobiol Int 2009, 26:108–125.PubMedView Article
- Kloog I, Stevens RG, Haim A, Portnov BA: Nighttime light level co-distributes with breast cancer incidence worldwide. Cancer Causes Control 2010, 21:2059–2068.PubMedView Article
- Doll R: Urban and rural factors in the aetiology of cancer. Int J Cancer 1991, 47:803–810.PubMedView Article
- Haenszel W, Marcus SC, Zimmerer EG: Cancer morbidity in urban and rural Iowa. Public Health Monogr 1956, 37:1–85.PubMed
- Parkin DM, Bray F, Ferlay J, Pisani P: Global Cancer Statistics, 2002. Ca Cancer J Clin 2005, 55:74–108.PubMedView Article
- Reynolds P, Hurley S, Goldberg DE, Anton-Culver H, Bernstein L, Deapen D, Horn-Ross PL, Peel D, Pinder R, Ross RK, West D, Wright WE, Ziogas A: Regional variations in breast cancer among California teachers. Epidemiology 2004, 15:746–754.PubMedView Article
- Reynolds P, Hurley SE, Quach AT, Rosen H, Von Behren J, Hertz A, Smith D: Regional variations in breast cancer incidence among California women, 1988–1997. Cancer Causes Control 2005, 16:139–150.PubMedView Article
- Ziskin DB, Hsu FC, Ghosh T, Elvidge C: Methods used for the 2006 radiance lights. Proceedings Of The Asia Pacific Advanced Network Meeting 2010 2010.
- Levallois P, Dumont M, Touitou Y, Gingras S, Masse B, Gauvin D, Kroger E, Bourdages M, Douville P: Effects of electric and magnetic fields from high-power lines on female urinary excretion of 6-Sulfatoxymelatonin. Am J Epidemiol 2001, 154:601–609.PubMedView Article
- Burgess HJ, Fogg LF: Individual differences in the amount and timing of salivary melatonin secretion. Plos One 2008, 3:E3055.PubMedView Article
- Travis RC, Allen DS, Fentiman IS, Key TJ: Melatonin and breast cancer: a prospective study. J Natl Cancer Inst 2004, 96:475–482.PubMedView Article
- Mahlberg R, Tilmann A, Salewski L, Kunz D: Normative data on the daily profile of urinary 6-sulfatoxymelatonin in healthy subjects between the ages of 20 and 84. Psychoneuroendocrinology 2006, 31:634–641.PubMedView Article
- Wetterberg L, Bratlid T, Von Knorring L, Eberhard G, Yuwiler A: A multinational study of the relationships between nighttime urinary melatonin production, age, gender, body size, and latitude. Eur Arch Psychiatry Clin Neurosci 1999, 249:256–262.PubMedView Article
- Bojkowski CJ, Arendt J: Factors influencing urinary 6-sulphatoxymelatonin, a major melatonin metabolite, in normal human subjects. Clin Endocrinol (Oxf) 1990, 33:435–444.View Article
- Youngstedt SD, Kripke DF, Elliott JA, Assmus JD: No association of 6-sulfatoxymelatonin with in-bed 60-Hz magnetic field exposure or illumination level among older adults. Environ Res 2002, 89:201–209.PubMedView Article
- Kennaway DJ, Lushington K, Dawson D, Lack L, Van Den Heuvel C, Rogers N: Urinary 6-sulfatoxymelatonin excretion and aging: new results and a critical review of the literature. J Pineal Res 1999, 27:210–220.PubMedView Article
- Zhao ZY, Xie Y, Fu YR, Bogdan A, Touitou Y: Aging and the circadian rhythm of melatonin: a cross-sectional study of chinese subjects 30–110 yr of age. Chronobiol Int 2002, 19:1171–1182.PubMedView Article
- Zhou JN, Liu RY, Van Heerikhuize J, Hofman MA, Swaab DF: Alterations in the circadian rhythm of salivary melatonin begin during middle-age. J Pineal Res 2003, 34:11–16.PubMedView Article
- Czeisler CA, Duffy JF, Shanahan TL, Brown EN, Mitchell JF, Rimmer DW, Ronda JM, Silva EJ, Allan JS, Emens JS, Dijk DJ, Kronauer RE: Stability, precision, and near-24-hour period of the human circadian pacemaker. Science 1999, 284:2177–2181.PubMedView Article
- Rea MS, Brons JA, Figueiro MG: Measurements Of Light At Night (Lan) for a sample of female school teachers. Chronobiol Int 2011, 28:673–680.PubMedView Article
- Davis S, Mirick DK, Stevens RG: Night shift work, light at night, and risk of breast cancer. J Natl Cancer Inst 2001, 93:1557–1562.PubMedView Article
- Schernhammer ES, Rosner B, Willett WC, Laden F, Colditz GA, Hankinson SE: Epidemiology of urinary melatonin in women and its relation to other hormones and night work. Cancer Epidemiol Biomarkers Prev 2004, 13:936–943.PubMed
- Bernstein L, Allen M, Anton-Culver H, Deapen D, Horn-Ross PL, Peel D, Pinder R, Reynolds P, Sullivan-Halley J, West D, Wright W, Ziogas A, Ross RK: High breast cancer incidence rates among California teachers: results from the California teachers study (United States). Cancer Causes Control 2002, 13:625–635.PubMedView Article
- Gunier RB, Reynolds P, Hurley SE, Yerabati S, Hertz A, Strickland P, Horn-Ross PL: Estimating exposure to polycyclic aromatic hydrocarbons: a comparison of survey, biological monitoring, and geographic information system-based methods. Cancer Epidemiol Biomarkers Prev 2006, 15:1376–1381.PubMedView Article
- Image and Data Processing By Noaa’s National Geophysical Data Center: DMSP data collected by the Us Air Force weather agency. http://ngdc.noaa.gov/eog/
- R Development Core Team: R: A Language And Environment For Statistical Computing. Vienna, Austria: R Foundation For Statistical Computing; 2009.
- Sun Or Moon Rise/Set Table: 2000. http://aa.usno.navy.mil/data/docs/RS_OneYear.php
- Hartman TJ, Mahabir S, Baer DJ, Stevens RG, Albert PS, Dorgan JF, Kesner JS, Meadows JW, Shields R, Taylor PR: Moderate alcohol consumption and 24-hour urinary levels of melatonin in postmenopausal women. J Clin Endocrinol Metab 2012, 97:E65-E68.PubMedView Article
- Mcpherson M, Janssen I, Grundy A, Tranmer J, Richardson H, Aronson KJ: Physical activity, sedentary behavior, and melatonin among rotating shift nurses. J Occup Environ Med 2011, 53:716–721.PubMedView Article
- Wu AH, Stanczyk FZ, Wang R, Koh WP, Yuan JM, Yu MC: Sleep duration, spot urinary 6-sulfatoxymelatonin levels and risk of breast cancer among chinese women in Singapore. Int J Cancer 2013, 132:891–896.PubMedView Article
- Youngstedt SD, Kripke DF, Elliott JA: Circadian phase-delaying effects of bright light alone and combined with exercise in humans. Am J Physiol Regul Integr Comp Physiol 2002, 282:R259-R266.PubMed
- Lunetta KL, Hayward LB, Segal J, Van Eerdewegh P: Screening large-scale association study data: exploiting interactions using random forests. BMC Genet 2004, 5:32.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.