Next Article in Journal
What Are the Place and Modalities of Surgical Management for Pancreatic Neuroendocrine Neoplasms? A Narrative Review
Previous Article in Journal
CARMA3 Promotes Colorectal Cancer Cell Motility and Cancer Stemness via YAP-Mediated NF-κB Activation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Systematic Review

Long-Term Nightshift Work and Breast Cancer Risk: An Updated Systematic Review and Meta-Analysis with Special Attention to Menopausal Status and to Recent Nightshift Work

by
Christine Schwarz
1,
Ana María Pedraza-Flechas
2,
Roberto Pastor-Barriuso
3,4,
Virginia Lope
3,4,
Nerea Fernández de Larrea
3,4,
José Juan Jiménez-Moleón
5,
Marina Pollán
3,4 and
Beatriz Pérez-Gómez
1,3,4,*
1
Public Health and Preventive Medicine Teaching Unit, National School of Public Health, Instituto de Salud Carlos III, Avda. Monforte de Lemos 5, 28029 Madrid, Spain
2
Medicine and Health Sciences, Universidad del Rosario, Calle 12C No 6-25, Bogotá D.C. 111221, Colombia
3
Department of Epidemiology of Chronic Diseases, National Centre for Epidemiology, Instituto de Salud Carlos III, Avda. Monforte de Lemos 5, 28029 Madrid, Spain
4
Consortium for Biomedical Research in Epidemiology and Public Health (CIBER en Epidemiología y Salud Pública—CIBERESP), 28029 Madrid, Spain
5
Department of Preventive Medicine and Public Health, School of Medicine, University Granada, Avda. de la Investigación No 11, 18016 Granada, Spain
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(23), 5952; https://doi.org/10.3390/cancers13235952
Submission received: 6 October 2021 / Revised: 13 November 2021 / Accepted: 22 November 2021 / Published: 26 November 2021
(This article belongs to the Section Systematic Review or Meta-Analysis in Cancer Research)

Abstract

:

Simple Summary

Night shift work (NSW) may disturb circadian rhythms. This could influence the risk of breast cancer (BC), but research papers have reported conflicting results. We reviewed, summarized, and combined the results of those studies that measured the effect of long-term nightshift work (≥15 years of NSW) in BC with special attention to menopausal status and time since retirement age. Women with long-term NSW had 13% more risk of BC than women without NSW. Postmenopausal women showed no increased risk, while premenopausal women had a 27% higher risk. Women with a higher probability of recent long-term NSW (women under retirement age) had a 23% higher risk than women without NSW. We concluded that long-term NSW may increase BC risk, especially in women before menopause or shortly after NSW discontinuation.

Abstract

This systematic review discusses long-term NSW and female BC risk, with special attention to differences between pre- and postmenopausal BC, to test the association with recent NSW. The review follows PRISMA guidelines (Prospero registry: CRD42018102515). We searched PubMed, Embase, and WOS for case–control, nested case–control, and cohort studies addressing long-term NSW (≥15 years) as risk exposure and female BC as outcome until 31 December 2020. Risk of bias was evaluated with the Newcastle–Ottawa scale. Eighteen studies were finally included (eight cohorts; five nested case–control; five case–control). We performed meta-analyses on long-term NSW and BC risk; overall and by menopausal status; a subanalysis on recent long-term NSW, based on studies involving predominantly women below retirement age; and a dose–response meta-analysis on NSW duration. The pooled estimate for long-term NSW and BC was 1.13 (95%CI = 1.01–1.27; 18 studies, I2 = 56.8%, p = 0.002). BC risk increased 4.7% per 10 years of NSW (95%CI = 0.94–1.09; 16 studies, I2 = 33.4%, p = 0.008). The pooled estimate for premenopausal BC was 1.27 (95%CI = 0.96–1.68; six studies, I2 = 32.0%, p = 0.196) and for postmenopausal BC 1.05 (95%CI = 0.90–1.24,I2 = 52.4%; seven studies, p = 0.050). For recent long-term exposure, the pooled estimate was 1.23 (95%CI = 1.06–1.42; 15 studies; I2 = 48.4%, p = 0.018). Our results indicate that long-term NSW increases the risk for BC and that menopausal status and time since exposure might be relevant.

1. Introduction

Many women’s working lives are characterized by some degree of night work, either as permanent shifts or as shifts within rotating work schedules. The European Union defines “night time” as any period of seven or more hours- as defined by national law, which includes the period between midnight and 5.00 a.m. [1]. Two large surveys provide information about the prevalence of night work in western societies [2,3]. According to the “Sixth European Working Conditions Survey” (2015), 14% of working women carried out night work [3], while in the “American Time Use Survey” of the US Bureau of Labor Statistics (2017/2018), 3.9% of working women reported night work, 2.1% rotating shifts, and 2.5% irregular shifts [2].
Humans share with other organisms the fundamental properties of a biological circadian system [4]. Night work has been studied as a possible risk factor for several diseases because physiological functions such as temperature and glucose homeostasis, hormone secretion, and immunological functions follow endogenous circadian rhythms [5]. These rhythms are controlled by the central circadian pacemaker, located in the suprachiasmatic nucleus of the hypothalamus [4]. Cellular circadian clocks generate roughly 24 h oscillations in the function of thousands of genes, allowing the organism to anticipate regular changes, to optimize energy, and to temporally separate incompatible processes [6]. As light is the most important stimulus that may reset circadian rhythms daily, by direct neuronal connections between the retina and the pacemaker, mistimed exposure to artificial light, usual in nightshift work (NSW), can cause important changes in the phase and amplitude of circadian rhythms, with deleterious effects in health [7].
In 2019, the International Agency for Research on Cancer, confirming a previous evaluation, classified NSW as probably carcinogenic to humans (Group 2A on a scale with a range from 1 “carcinogenic to humans” to 3 “not classifiable as to its carcinogenicity to humans”). Its decision was based on limited evidence from studies in humans and strong evidence from experimental findings, which showed its effects in immunosuppression, chronic inflammation, and cell proliferation [8]. The tumors with positive associations with NSW were mainly breast, prostate, and colorectal cancer [8].
This issue has awakened enormous interest, given the global burden of breast cancer (BC) [9] and the considerable prevalence of NSW. Unfortunately, the results of individual epidemiological studies have so far been inconsistent. Additionally, a number of meta-analyses have come to different conclusions [10,11,12]. One factor that could be related to the inconsistency of results in observational studies is the duration of NSW. Most of the studies that evaluated different categories of duration observed an elevated BC risk for long-term NSW (≥15/≥20 years) [13,14,15,16,17,18,19], in contrast to the moderate [16,17,19,20] or absent excess risk observed for shorter exposures (<15 years) [12,14,15,18,21,22].
Another important issue that might help to understand the inconsistent results is the possible role of menopausal status in this presumptive association. Pre- and postmenopausal breast cancer have different temporal trends, can differ in their molecular profile, and even have opposite relationships with some of the known risk factors for this tumor, such as obesity (i.e., obese women have a lower risk of premenopausal cancer and higher risk of postmenopausal tumors than women of normal weight) [23,24,25].
Another possible reason for these discrepancies might be related to a critical point: the lag time between the last NSW exposure and the end of follow-up for the occurrence of BC in some cohort studies [26,27]. Some authors have hypothesized that NSW may act essentially as a promoter in BC; if that were the case, when NSW exposure ceases, the associated BC risk increase should wane [13,28], while the increased risk associated with NSW should be mostly visible in recently exposed women. This pattern of decreasing risk since the cessation of exposure has also been proposed for the suggested association between NSW and other tumors [29]. As far as we know, no meta-analysis has focused on this interesting issue, although a pooled analysis of case–control studies found a higher premenopausal BC risk in women with recent exposure than in those who had stopped night work more than two years before [30].
For this reason, we carried out an overall meta-analysis of observational studies on long-term NSW and BC, a dose–response meta-analysis, and a sub-meta-analysis on recent long-term NSW. Therefore, this meta-analysis allows exploring both factors, the duration of NSW, and recent exposure in BC risk assessment. Given the possible biological differences between pre- and postmenopausal BC associated with long-term NSW, we also carried out subgroup analyses by menopausal status.

2. Materials and Methods

2.1. Literature Search and Eligibility Criteria

This review was conducted according to the PRISMA guidelines [31]. The protocol is registered in PROSPERO [32] (CRD42018102515), where methodological details can be consulted. We updated the search for original articles in PubMed/MEDLINE, Embase, and the Web of Science, using the terms (“nightshift” OR “night shift” OR “shiftwork” OR “shift work” OR “nightwork”) combined with “Breast cancer” from inception until 31 December 2020.
Our final search strategy was complemented with a manual search based on the references cited in the papers selected for full-text review. Additionally, we also went through the reference lists of the published meta-analyses. We took into consideration results in English, Spanish, German, or French. For this review, we considered eligible all original case–control, nested case–control, and cohort studies that compared, as risk exposure, long-term NSW (defined as ≥15 years) to having never done NSW. These studies had to provide original data and to have female BC risk as the outcome.

2.2. Data Extraction

Two members of our research team (A.P., C.S.) performed the original abstract/title screening, the full-text review, and the data extraction. In case of doubts or disagreements, a third author’s opinion (B.P.G.) allowed reaching a consensus. In the presence of diverse options of risk estimators, we used preferentially those including NSW exposure between 12.00 p.m. and 3.00 a.m. In cohort studies with several reports, we used the latest update if there were new cases or if the analysis was improved in terms of methodology. With respect to exposure ascertainment, we prioritized self-reported individual data over job-exposure matrix estimations. According to our inclusion criteria, we selected studies providing data corresponding to a minimum of 15 years of NSW, without restricting the exposure to a minimum number of nights per month. Authors from two studies were contacted to request additional data about the definition of nonexposure (information received) and mean age at BC diagnosis (information finally not available).

2.3. Risk of Bias Assessment of Individual Studies (Newcastle–Ottawa Scale)

The risk of bias of each study was evaluated by two researchers (A.P., C.S.) using the Newcastle–Ottawa scale (NOS) [33], an eight-item scale developed to assess the quality of nonrandomized studies for meta-analyses with a maximum punctuation of nine stars (see Table S1). It assigns four stars to items related to selection and representativeness of the study groups. Another three stars correspond to the evaluation of the ascertainment of the exposure or the outcome of interest for case–control or cohort studies, respectively, and some of these items require operative definitions. In this case, for cohort studies, we assigned one star to those with at least 15 years of follow-up (15 years of NSW exposure does not determine the follow-up time, since exposure may be assessed retrospectively) and another to those with a difference ≤15% in the proportion of subjects lost to follow-up between exposed and unexposed. For case–control studies, we gave one star to those with differences ≤15% in nonresponse rates between cases and controls. Finally, the remaining two stars evaluate the comparability of the groups in terms of control for potential confounders; for this report, one star was given to studies that controlled for age and, additionally, for at least four of a predefined set of covariables (namely, family history of BC, age at first birth or number of children, body mass index in postmenopausal women, hormone replacement therapy, and menopausal status), and a second star was rendered to those studies that included at least five of these additional covariables: alcohol intake, oral contraception, age at menarche, age at menopause, breastfeeding, physical activity, history of benign breast disease, and education/socioeconomic status.

2.4. Statistical Methods

2.4.1. Meta-Analysis on Long-Term NSW and BC Risk

  • Global meta-analysis of long-term NSW (for ≥15 years) and BC
We compared BC risk for women exposed to NSW for ≥15 years to that of unexposed women. We estimated summary risk ratios (RRs) and 95% confidence intervals (CIs) as the measure of effect. Individual effect estimates were extracted to construct a database with the logarithm of the risk estimator reported in each study (RR, hazard ratio (HR), or odds ratio (OR)) and its standard error (SE). If SEs were not provided in the publication, we calculated them from the CI, using the formula: SE = (lnCIupper limit-lnCIlower limit)/(2 × 1.96). When the report provided different categories of exposure but not an overall estimator, we calculated pooled RRs and variances using the formulas: lnRR pooled = (lnRR1/Var1 + lnRR2/Var2 + lnRRn/Varn)/(1/Var1 + 1/Var2 + 1/Varn) and Varpooled = 1/(1/Var1 + 1/Var2 + 1/Varn). The proportion of variance due to heterogeneity was estimated by the I2 statistic [34]. Heterogeneity between studies was assumed, and a random-effects model was applied using the method of DerSimonian and Laird, with the estimate of heterogeneity taken from the inverse-variance fixed-effects model.
To assess publication bias, funnel plot asymmetry was visually inspected and statistically tested with Begg’s and Egger’s tests (NFL, CS). We also carried out two sensitivity analyses: (a) excluding those studies with uncertain exposure definition (see Figure S1) and (b) restricting the analysis to those studies with NOS ≥ 7 (see Figure S2).
Program Codes Specific STATA© version 15 codes (metan, metafunnel, metabias) from StataCorp. 2017. Stata Statistical Software: Release 15. College Station, TX: StataCorp LLC were used for the meta-analysis and RevMan 5.4 test for differences in the predefined subgroup analyses [35].
  • Dose–response meta-analysis
For studies reporting disaggregated results by duration of NSW, we further conducted a random-effects dose–response meta-analysis of log-relative risks (either adjusted log hazard ratios for cohort studies or adjusted log odds ratios for case–control studies) on mean durations of NSW for each exposure category [36]. The dose–response meta-analysis was fitted through generalized least squares, which allowed for within-study correlations among log-relative risks that used the same unexposed group as reference [37]. Linear trends were allowed to vary randomly across studies, whose variance was estimated by the method of moments [38]. We evaluated departures from the linear trend by including a quadratic term for mean duration, as well as distinct linear trends in cohort, nested case–control, and case–control studies by including interaction terms between mean duration and study design indicators. The heterogeneity of log-relative risks with respect to the pooled linear trend was contrasted with the Cochran chi-squared test and quantified with the I2 statistic, which described the proportion of total variation in log-relative risks due to heterogeneity [34]. The dose–response meta-analysis was performed using the glst command in Stata, version 14 from StataCorp. 2015. Stata Statistical Software: Release 14. College Station, TX: StataCorp LP, and graphics were produced in R, version 3 (R Foundation for Statistical Computing).
  • Subgroup meta-analysis by menopausal status
We stratified the overall analysis by menopausal status, which was preferentially based on self-reported information. If this information was not available, given cut-off points between age 50 and 55 were accepted to classify women either as pre- or as postmenopausal. Again, a sensitivity analysis was performed, restricting the meta-analyses to those studies with NOS ≥ 7.

2.4.2. Meta-Analysis on Recent Long-Term NSW and BC

As the information on the timing of long-term NSW exposure was not readily available in the selected studies, we used an indirect approach based on the analysis of BC before retirement. We classified as having recent NSW those participants in the subgroup of studies in which, either the mean age of the participants at study entry plus the maximum follow-up time (cohort studies) or the mean age at BC diagnosis (case–control studies) was below the country-specific retirement age in that year; this age was obtained from the official government websites of each country or of the Organization for Economic Co-operation and Development [39]. The meta-analysis was performed following the same methodology that has been previously described.

3. Results

3.1. Search Results

The initial search yielded 1524 records found in PubMed, Embase, and Web of Science. After the selection process, 16 articles with data from 18 research studies were finally included in the global meta-analysis (Figure 1).
It is noteworthy to indicate that the early publications of the Nurses Health cohort studies (NHS I and II) [20,40,41] were excluded from the main meta-analysis, due to the availability of updated data with improved analyses from these projects, and that five studies included in previous meta-analyses did not meet our inclusion criteria, as they did not investigate long-term exposures (≥15 years of NSW) [28,42,43,44,45]. Our manual search did not identify any additional articles.

3.2. Study Characteristics

Table 1 summarizes the main characteristics of the studies included in our global meta-analysis, and the characteristics of other studies that did not meet our inclusion criteria but that were included in previous reviews; this way, the reader may have a rapid overview of the differences among the published meta-analyses in this area. In regard to those studies fulfilling our criteria, eight research projects, published in six articles, had a cohort design [12,13,14,22,46,47], five were nested case–control studies [17,18,21,48,49], and five were case–control studies [15,16,19,50,51]. The numbers of participants and BC cases in the table correspond to the risk estimators used for the present global meta-analysis. For risk estimators, all cohort studies reported HRs, except for one that reported RRs [20]; four of the five nested case–control studies reported ORs, and the remaining one reported HR [48], while all case–control studies reported ORs. Regarding NSW exposure, it should be noted that most of the studies only reported exposure status at baseline and did not update this information; therefore, the real number of women who might have finally been exposed to long-term NSW could be higher than reported, had we taken into account the additional years of NSW during follow-up.
Table 2 provides further details on exposure characteristics for the studies selected in this review, as well as indicates which are included in each of the meta-analyses presented in this report. The table summarizes the definition of exposure and nonexposure in each study. It should be noted that the operational approach used in several cases may raise doubts about the accuracy of the classification [14,21,22,47]. Thus, in the two Nurses’ Health Studies [14], the definition for nonexposure was “never rotating NSW”, with unclear classification for those women with permanent NSW. In another study among Norwegian nurses [49], any woman working at infirmaries was classified as exposed to NSW, while women working in managerial, teaching, physiotherapy, or at out-patients’ departments formed the nonexposed group. Additionally, in some reports only those with current NSW and long-term stability at the same job (“job tenure”) [22] or with NSW in the last 10 years [47], were considered as exposed, classifying as nonexposed those women with former NSW exposure.
Concerning the timing of exposure ascertainment, in cohort studies, it was mostly performed only at baseline, with just a few studies collecting data of additional years of NSW during follow-up [14,47]. In nested case–control studies and case–control studies, the ascertainment was retrospective, taking place after follow-up in the former and at recruitment in the latter. Self-reporting was the most common source of information to define NSW [12,13,14,15,16,17,18,19,47,49,50,51], while three studies [21,22,48] applied job-exposure matrices (JEM) for this purpose. Interestingly, the Shanghai Women’s Health cohort study compared both approaches, concluding that the use of JEM overestimated the prevalence of NSW [46].
As regards the duration of exposure, in 10 studies, the risk estimates provided by the authors corresponded to NSW exposures of at least 20 years [12,13,14,15,17,22,47,48,50] and in one study to exposures greater than 17 years [46].

3.3. Results of the Risk of Bias Assessment in the Selected Studies

The detailed risk of bias assessment for each selected study, according to the Newcastle–Ottawa scale (NOS), is available in the Table S1. In summary, NOS scores ranged from four [22] to eight stars [14,16,18]. Among cohort studies, the highest scores were obtained by NHS I and II [14], penalized only by their focus on a single occupation. The other cohort studies studied the general population, except for one, which included twins [13]. All but three cohort studies [22,28,46] achieved our highest quality criteria for control of confounders. Only NHS I and II reported the proportion of participants lost to follow-up for exposed and nonexposed groups and had a follow-up length over our threshold [14], but two other cohort studies presented an overall number of losses to follow-up below 5% [46,47]. Three of the nested case–control studies received high scores in the NOS (≥7) [17,18,49], and all case–control studies reached at least six points. The highest score among them corresponded to a Canadian study [16], only penalized for not blinding the interviewer, a common problem that often has a difficult logistic solution in this type of study; in fact, only one case–control study used interviewers blinded to the disease status [50]. The other common quality problem present in case–control-studies [15,19,50] laid on differences in the response rate between cases and controls.

3.4. Results of the Meta-Analyses

3.4.1. Long-Term NSW and BC Risk

  • Global meta-analysis on long-term NSW and BC risk
The meta-analysis of long-term NSW and BC risk (Figure 2) showed a global significant increased risk of BC, although there was moderate heterogeneity of effects between studies (pooled effect estimate = 1.13, 95%CI = 1.01–1.27, 18 studies, I2 = 56.8%, p = 0.002). When we took into account study designs, the meta-analysis yielded no association between NSW and BC in cohort studies (pooled effect estimate = 1.04, 95%CI = 0.91–1.18, I2 = 36.3%, p = 0.139); for nested case–control studies, it showed an almost significant risk increase, although with high heterogeneity (pooled effect estimate = 1.34, 95%CI = 0.96–1.88, I2 = 82.2%, p < 0.001); and, finally, in case–control studies, the pooled risk indicated an excess risk of BC among exposed women (pooled effect estimate = 1.22, 95%CI = 1.01–1.48, I2 = 0.0%, p = 0.613). We performed sensitivity analyses: (a) excluding the studies with uncertain exposure definition, which showed similar results (pooled effect estimate = 1.17, 95%CI = 0.99–1.39, I2 = 61.8%, p = 0.002) (Figure S1); and (b) including only studies that scored ≥7 stars at the NOS, which yielded the same overall risk estimate for long-term NSW and BC (pooled effect estimate = 1.13, 95%CI = 1.00–1.28, I2 = 48.2%, p = 0.031) (see Figure S2).
The funnel plot to evaluate a possible relation between the RR of BC and the study precision showed an imbalance (Figure 3), which was statistically significant (Egger’s test p = 0.03, Begg’s test p = 0.08). The asymmetry was not present in a funnel plot limited to the high-quality studies (NOS ≥ 7) (see Figure S3) (Egger’s test p = 0.33, Begg’s test p = 0.24).
  • Dose–response meta-analysis
The dose–response meta-analysis was based on 16 studies. Risk trends across categories of exposure duration varied substantially from one study to another, although most of them tended to show an increase in risk with rising duration of NSW (Figure 4). The pooled risk for BC increased by 4.7% for every 10 years of NSW duration (pooled effect estimate = 1.05, 95%CI = 0.94–1.09), with moderate heterogeneity across studies (I2 = 33.4%, p = 0.008), and no significant differences among cohort, nested case–control, and case–control studies (p for interaction between exposure duration and study design = 0.45). Adding a quadratic term to the model for exposure duration did not significantly improve model fit (p = 0.389).
  • Subgroup analysis by menopausal status
Menopausal status was generally self-reported, although in some studies, this information was not available. Overall, long-term NSW was associated with an almost significant increased risk for pre-/perimenopausal BC (pooled effect estimate = 1.27, 95%CI = 0.96–1.68, six studies, I2 = 32.0%, p = 0.196) (Figure 5), mostly due to the results of cohort studies, while no association was found for postmenopausal BC (pooled effect estimate = 1.05, 95%CI = 0.90–1.24, seven studies, I2 = 52.4%, p = 0.050) (Figure 6), although the test for subgroup differences was not significant (p = 0.27). However, in our sensitivity analysis, restricted to high-quality studies, the difference between both subgroups was more marked (pooled effect estimate premenopause = 1.57; 95%: 1.11–2.21; pooled effect estimate postmenopausee = 1.08; 95%: 0.94–1.24; test for subgroup differences: p = 0.05) (see Figures S4 and S5).

3.4.2. Recent Long-Term NSW and BC Risk

This last meta-analysis included data from 15 studies. Among them, data from the NHS1 after 10 years of follow-up could be included, as the biological age of menopause normally is below the retirement age of 65 years [20]. The updated HRs for women from the NHS2 after 10 years of follow-up [14] and the HRs for women Shanghai textile workers under 50 years [48] could also be included. Moreover, we included a Danish nested case–control study [17], which provided mean age at the beginning of the study, and the sum of the years of follow-up resulted in ages below the retirement age in Denmark at that time. The meta-analysis of recent long-term NSW and BC showed a significant pooled risk of BC of around 23% (95%CI = 1.06–1.42), with moderate heterogeneity between studies (I2 = 48.4%, p = 0.018) (Figure 7) (test of interaction between NSW and BC before and after retirement: p = 0.001). The results were similar in all study design groups (20% for cohort studies, 95%CI = 0.88–1.65, 29% for nested case–control studies, 95%CI = 0.96–1.73, and 22% for case–control studies, 95%CI = 1.01–1.48).

3.5. Comparison with Former Meta-Analyses

Table 3 shows the main summary effect estimates of previous meta-analyses, dose–response meta-analyses, and sub-meta-analyses for the association between long-term NSW and BC risk, globally and in pre- and postmenopausal women, to allow the reader to have a better perspective of the results that have been published in this matter.

4. Discussion

4.1. The Association between Long-Term NSW and BC Risk

The global meta-analysis showed that the risk of BC in women exposed to ≥15 years of NSW is 13% higher than that of those who had not worked nightshifts, although there was moderate heterogeneity among the studies. In the dose–response meta-analysis, the BC risk increased by 4.7% for every 10 years of exposure to NSW. The analysis by menopausal status suggested that this association was limited to premenopausal BC, while no association was observed for BC risk after menopause; this difference was quite more marked when we restricted the analysis to high-quality studies. We also found an excess of risk among those women exposed to recent long-term NSW, but due to the scarcity of studies, we did not evaluate differences by menopausal status in this subgroup. It should be noted that the use of aggregated data on age and retirement as a proxy for recent exposure, due to the lack of individual information, renders it very difficult to separate the possible effects of the time lag (between exposure cessation and outcome) from those of menopausal status and age.

4.2. Comparison with Former Meta-Analyses and Relevant Studies

Our global findings are similar to the results of a meta-analysis carried out by Jia et al. [10], who reported a 15% increased BC risk for long-term NSW. Additionally, He et al. [54] found a 6% increased BC risk for every 10 years of exposure in their dose–response meta-analysis, similar to our result. Considering only the results from cohort studies, our pooled estimate was comparable with the null results of previous meta-analyses [12,54], while another one detected an 8.8% risk increase after 20 years of NSW exposure [11]. For premenopausal BC and long-term NSW, we found a slightly lower increased risk than the one reported by Cordina-Duverger et al. [30] for an exposure of at least twenty years after pooling the results of case–control studies, though in that study, no elevated risk was observed for exposure durations between 10 and 20 years. As this review includes, to our knowledge, the first intent of carrying out a sub-meta-analysis of BC risk before retirement and long-term NSW, we could not compare the results related to recent NSW exposure.

4.3. Interpretation of the Results and Possible Mechanisms of the Association

Our analyses of long-term NSW suggest that this exposure was associated with BC before retirement or menopause. It is conceivable that, during regular NSW exposure, changes in sex hormone secretion [65] or other biological mechanisms might drive the growth of malignant tumor cells with a promoter effect; if that were the case, once exposure ceases, the BC risk might drop to normal levels. Similar effects had been proposed for the action of hormone replacement therapy [66] but were later partly questioned by a meta-analysis [67]. Alternatively, some cancer-relevant epigenetic, genetic, or microenvironmental changes caused or favored by long-term NSW might contribute to transforming breast epithelial cells into malignant cells before menopause, supported by the higher cell proliferation rates at this stage [68]. It is also plausible that the NSW-related risks persist during more years, but cumulative DNA damage with age, caused by a wide spectrum of reasons [69], could mask the deleterious effects of NSW in older women. In regard to long-term NSW, having been working in night shifts for at least 15 years makes it more likely to be exposed in windows of increased susceptibility, that is, in life periods during which harmful factors can more easily lead to cancer. One of these susceptibility periods is the time before a first birth [70]. The early life etiological model of BC [70] suggests that childbirth and breastfeeding cause differentiation of mammary stem cells so that they are consequently less prone to proliferate and cause BC. Thus, NSW exposure after the first childbirth is expected to be less harmful than before this life event. However, the few risk estimates for NSW exposure before the first childbirth found in the literature are contradictory [19,45,50]. Menopausal transition might also be a period of higher susceptibility [71], but unfortunately, the impact of NSW during this susceptibility window has not been studied yet.

4.4. Strengths and Limitations

The exhaustive literature search, including studies published up to December 2020, allowed us to provide updated estimates of risk associated with long-term NSW, as well as to perform an updated dose–response meta-analysis. In contrast to a previously published report [30], we included cohort, nested case–control, and case–control studies. For this review, we assessed the risk of bias of included studies and carried out a sensitivity analysis based on their quality scores. Additionally, we explored the consistency of our results after excluding those studies with uncertain exposure definitions or ascertainment. However, probably, the major strength of this meta-analysis resides in its deeper analysis of the risk for long-term NSW, thanks to the sub-meta-analyses carried out by menopausal status, a factor involved in clinical and epidemiological differences in BC, and by recent exposures, which allow exploring the hypothesis under debate of a reduction of risk after the cessation of NSW exposure. Noteworthy, the subgroup analysis by menopausal status was mainly based on within studies comparisons. Our approach on recent long-term NSW, even with its mentioned limitations, represents an original analysis trying to overcome the scarcity of information about time since NSW cessation. Although we are not able in this review to decipher the reasons for the association found between long-term NSW and breast cancer during working age, we think it is an interesting new finding that deserves further investigation.
Notwithstanding, this review also has limitations. The asymmetry of the funnel plot may suggest that our estimates can be affected by publication bias or by differences in the quality of smaller studies that may spuriously inflate their effects [72]. However, our sensitivity analysis, restricted to the studies of higher methodological quality (NOS ≥ 7), yielded the same overall risk as the global meta-analysis, while asymmetry, according to the funnel plot and the Egg test, was no longer significant. Another limitation comes from the heterogeneity observed in the meta-analysis. As previously mentioned, we explored several putative factors that might contribute to explaining it (e.g., study design and study quality, menopause, exposure definition); however, other design characteristics, such as NSW pattern (e.g., permanent versus rotating shift work or frequency of shifts) or NSW assessment (i.e., self-report versus JEM), may also contribute to it. It has been described that self-report can be used as a reliable method for NSW exposure [73], while, in contrast, JEM might result in misclassification bias towards the null [46]. As the three included studies [21,22,48] that used JEMs received high weights in the meta-analysis (in total 23%), it is possible that we might be underestimating the risk due to NSW. The distribution of BC molecular subtypes might also differ among studies, being an additional uncontrolled potential source of heterogeneity, although it is not yet known whether NSW may have a different role in the pathogenesis of BC per subtype. Cordina-Duverger recently described higher risks of NSW linked to tumors that were both HER+ and ER+ [30]. Finally, we could only carry out an indirect approach to the risk associated with recent long-term NSW. This approach assumes that NSW is current or recent if the woman—really, the mean age of the participants—has not reached retirement age. However, we do not have data to test this assumption, which surely includes a certain misclassification. Additionally, it is difficult to disentangle the effect of menopause, as both are age related.

5. Conclusions

Our meta-analysis supports the hypothesis that long-term NSW is a risk factor for BC. Even though the possible effect of the publication bias is small, we prudently classify the evidence as low. Among the factors that may modify this risk, which are possible sources of the observed heterogeneity, are menopause and time since last NSW exposure. Both factors point towards higher risks in younger women. In contrast, no association was found between long-term NSW and BC in postmenopausal women. Future investigators on NSW and this tumor should also take into consideration the timing of exposure and the menopausal status at BC diagnosis.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13235952/s1, Figure S1: Analysis excluding five studies with uncertain exposure definitions, Figure S2: Analysis of long-term NSW and BC restricted to those studies that reached at least seven points in the Newcastle Ottawa Scale (NOS ≥ 7), Figure S3: Funnel plot limited to studies of high quality (NOS ≥ 7), Figure S4: Analysis of long-term NSW and premenopausal BC limited to studies of high quality (NOS ≥ 7), Figure S5: Analysis of long-term NSW and postmenopausal BC limited to studies of high quality (NOS ≥ 7), Table S1: Risk of bias assessment by the Newcastle–Ottawa scale for each individual study.

Author Contributions

Conceptualization, B.P.-G., C.S. and A.M.P.-F.; methodology, B.P.-G., A.M.P.-F. and C.S. validation, B.P.-G., N.F.d.L., R.P.-B. and V.L; formal analysis, A.M.P.-F., R.P.-B. and C.S.; investigation, A.M.P.-F., R.P.-B. and C.S. resources, B.P.-G. and M.P., data curation, A.M.P.-F. and C.S.; writing—original draft preparation, C.S. and B.P.-G.; writing—review and editing R.P.-B., V.L., N.F.d.L., J.J.J.-M. and M.P.; visualization, A.M.P.-F., B.P.-G., R.P.-B. and C.S.; supervision, B.P.-G., M.P. and J.J.J.-M.; project administration, B.P.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Publicly available datasets were analyzed in this study. The analyzed data can be found in the reported publications (see references of Table 2).

Conflicts of Interest

The authors declare no conflict of interest. The research results do not represent any official position of the Instituto de Salud Carlos III.

References

  1. Directive 2003/88/EC of the European Parliament and of the Council of 4 November 2003 Concerning Certain Aspects of the Organisation of Working Time, EUR-LEX, Access to European Union Law; European Parlament and European Council. Available online: https://eur-lex.europa.eu/legal-content/en/TXT/?uri=CELEX:32003L0088 (accessed on 11 March 2020).
  2. Job Flexibilites and Work Schedules-2017-2018. Data from the American Time Use Survey. Available online: https://www.bls.gov/news.release/flex2.nr0.htm (accessed on 11 March 2020).
  3. Parent-Thirion, A.; Biletta, I.; Cabrita, J.; Llave Vargas, O.; Vermeylen, G.; Wilczynska, A.; Wilkens, M. 6th European Working Conditions Survey: Overview Report; EF; 2017 Update; Publications Office of the European Union: Luxembourg, 2017; ISBN 978-92-897-1596-6. [Google Scholar]
  4. Czeisler, C.A.; Gooley, J.J. Sleep and Circadian Rhythms in Humans. Cold Spring Harb. Symp. Quant. Biol. 2007, 72, 579–597. [Google Scholar] [CrossRef] [Green Version]
  5. Czeisler, C.A.; Duffy, J.F.; Shanahan, T.L.; Brown, E.N.; Mitchell, J.F.; Rimmer, D.W.; Ronda, J.M.; Silva, E.J.; Allan, J.S.; Emens, J.S.; et al. Stability, Precision, and Near-24-Hour Period of the Human Circadian Pacemaker. Science 1999, 284, 2177. [Google Scholar] [CrossRef] [PubMed]
  6. Panda, S. Circadian Physiology of Metabolism. Science 2016, 354, 1008–1015. [Google Scholar] [CrossRef] [Green Version]
  7. Duffy, J.F.; Czeisler, C.A. Effect of Light on Human Circadian Physiology. Sleep Med. Clin. 2009, 4, 165–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. International Agency for Research on Cancer (IARC). Night Shift Work. In Monographs on the Identification of Carcinogenic Hazards to Humans; IARC: Lyon, France, 2015; Volume 124, ISBN 9789283201625. [Google Scholar]
  9. Ginsburg, O.; Bray, F.; Coleman, M.P.; Vanderpuye, V.; Eniu, A.; Kotha, S.R.; Sarker, M.; Huong, T.T.; Allemani, C.; Dvaladze, A.; et al. The Global Burden of Women’s Cancers: A Grand Challenge in Global Health. Lancet 2017, 389, 847–860. [Google Scholar] [CrossRef]
  10. Jia, Y.; Lu, Y.; Wu, K.; Lin, Q.; Shen, W.; Zhu, M.; Huang, S.; Chen, J. Does Night Work Increase the Risk of Breast Cancer? A Systematic Review and Meta-Analysis of Epidemiological Studies. Cancer Epidemiol. 2013, 37, 197–206. [Google Scholar] [CrossRef] [PubMed]
  11. Lin, X.; Chen, W.; Wei, F.; Ying, M.; Wei, W.; Xie, X. Night-Shift Work Increases Morbidity of Breast Cancer and All-Cause Mortality: A Meta-Analysis of 16 Prospective Cohort Studies. Sleep Med. 2015, 16, 1381–1387. [Google Scholar] [CrossRef]
  12. Travis, R.C.; Balkwill, A.; Fensom, G.K.; Appleby, P.N.; Reeves, G.K.; Wang, X.-S.; Roddam, A.W.; Gathani, T.; Peto, R.; Green, J.; et al. Gathani Toral, Peto Richard, Green Jane, Key Timothy J., Beral Valerie. Night Shift Work and Breast Cancer Incidence: Three Prospective Studies and Meta-analysis of Published Studies. J. Natl. Cancer Inst. 2016, 108, djw169. [Google Scholar] [CrossRef]
  13. Åkerstedt, T.; Knutsson, A.; Narusyte, J.; Svedberg, P.; Kecklund, G.; Alexanderson, K. Night Work and Breast Cancer in Women: A Swedish Cohort Study. BMJ Open 2015, 5, e008127. [Google Scholar] [CrossRef]
  14. Wegrzyn, L.R.; Tamimi, R.M.; Rosner, B.A.; Brown, S.B.; Stevens, R.G.; Eliassen, A.H.; Laden, F.; Willett, W.C.; Hankinson, S.E.; Schernhammer, E.S. Rotating night shift work and risk of breast cancer in the nurses’ health studies. Am. J. Epidemiol. 2017, 186, 532–540. [Google Scholar] [CrossRef] [Green Version]
  15. Pesch, B.; Harth, V.; Rabstein, S.; Baisch, C.; Schiffermann, M.; Pallapies, D.; Bonberg, N.; Heinze, E.; Spickenheuer, A.; Justenhoven, C.; et al. Night Work and Breast Cancer—Results from the German GENICA Study. Scand. J. Work. Environ. Health 2010, 36, 134–141. [Google Scholar] [CrossRef] [Green Version]
  16. Grundy, A.; Richardson, H.; Burstyn, I.; Lohrisch, C.; SenGupta, S.K.; Lai, A.S.; Lee, D.; Spinelli, J.J.; Aronson, K.J. Increased Risk of Breast Cancer Associated with Long-Term Shift Work in Canada. Occup. Environ. Med. 2013, 70, 831–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Hansen, J.; Stevens, R.G. Case–Control Study of Shift-Work and Breast Cancer Risk in Danish Nurses: Impact of Shift Systems. Eur. J. Cancer 2012, 48, 1722–1729. [Google Scholar] [CrossRef] [PubMed]
  18. Hansen, J.; Lassen, C.F. Nested Case–Control Study of Night Shift Work and Breast Cancer Risk among Women in the Danish Military. Occup. Environ. Med. 2012, 69, 551–556. [Google Scholar] [CrossRef]
  19. Papantoniou, K.; Castaño-Vinyals, G.; Espinosa, A.; Aragonés, N.; Pérez-Gómez, B.; Ardanaz, E.; Altzibar, J.M.; Sanchez, V.M.; Gómez-Acebo, I.; Llorca, J.; et al. Breast Cancer Risk and Night Shift Work in a Case-Control Study in a Spanish Population. Eur. J. Epidemiol. 2016, 31, 867–878. [Google Scholar] [CrossRef] [Green Version]
  20. Schernhammer, E.S.; Laden, F.; Speizer, F.E.; Willett, W.C.; Hunter, D.J.; Kawachi, I.; Colditz, G.A. Rotating Night Shifts and Risk of Breast Cancer in Women Participating in the Nurses’ Health Study. J. Natl. Cancer Inst. 2001, 93, 1563–1568. [Google Scholar] [CrossRef] [Green Version]
  21. Lie, J.-A.S.; Roessink, J.; Kjærheim, K. Breast Cancer and Night Work among Norwegian Nurses. Cancer Causes Control 2006, 17, 39–44. [Google Scholar] [CrossRef]
  22. Koppes, L.L.J.; Geuskens, G.A.; Pronk, A.; Vermeulen, R.C.H.; de Vroome, E.M.M. Night Work and Breast Cancer Risk in a General Population Prospective Cohort Study in The Netherlands. Eur. J. Epidemiol. 2014, 29, 577–584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Word Cancer Resarch Fund International (WCRFI). American Institute for Cancer Research Breast Cancer. Available online: https://www.wcrf.org/dietandcancer/breast-cancer/ (accessed on 2 November 2021).
  24. Schoemaker, M.J.; Nichols, H.B.; Wright, L.B.; Brook, M.N.; Jones, M.E.; O’Brien, K.M.; Adami, H.-O.; Baglietto, L.; Bernstein, L.; Bertrand, K.A.; et al. Association of Body Mass Index and Age With Subsequent Breast Cancer Risk in Premenopausal Women. JAMA Oncol. 2018, 4, e181771. [Google Scholar] [CrossRef]
  25. Liao, S.; Hartmaier, R.J.; McGuire, K.P.; Puhalla, S.L.; Luthra, S.; Chandran, U.R.; Ma, T.; Bhargava, R.; Modugno, F.; Davidson, N.E.; et al. The Molecular Landscape of Premenopausal Breast Cancer. Breast Cancer Res. 2015, 17, 1–13. [Google Scholar] [CrossRef] [Green Version]
  26. Stevens, R.G. RE: Night Shift Work and Breast Cancer Incidence: Three Prospective Studies and Meta-Analysis of Published Studies. J. Natl. Cancer Inst. 2017, 109. [Google Scholar] [CrossRef] [Green Version]
  27. Schernhammer, E.S. RE: Night Shift Work and Breast Cancer Incidence: Three Prospective Studies and Meta-Analysis of Published Studies. JNCI J. Natl. Cancer Inst. 2017, 109. [Google Scholar] [CrossRef]
  28. Knutsson, A.; Alfredsson, L.; Karlsson, B.; Åkerstedt, T.; Fransson, E.I.; Westerholm, P.; Westerlund, H. Breast Cancer among Shift Workers: Results of the WOLF Longitudinal Cohort Study. Scand. J. Work. Environ. Health 2013, 39, 170–177. [Google Scholar] [CrossRef] [Green Version]
  29. Kogevinas, M.; Espinosa, A.; Papantoniou, K.; Aragonés, N.; Pérez-Gómez, B.; Burgos, J.; Gómez-Acebo, I.; Llorca, J.; Peiró, R.; Jimenez-Moleón, J.J.; et al. Prostate Cancer Risk Decreases Following Cessation of Night Shift Work. Int. J. Cancer 2019, 145, 2597–2599. [Google Scholar] [CrossRef] [PubMed]
  30. Cordina-Duverger, E.; Menegaux, F.; Popa, A.; Rabstein, S.; Harth, V.; Pesch, B.; Brüning, T.; Fritschi, L.; Glass, D.C.; Heyworth, J.S.; et al. Night Shift Work and Breast Cancer: A Pooled Analysis of Population-Based Case-Control Studies with Complete Work History. Eur. J. Epidemiol. 2018, 33, 369–379. [Google Scholar] [CrossRef]
  31. Page, M.J.; Moher, D.; Bossuyt, P.M.; Boutron, I.; Hoffmann, T.C.; Mulrow, C.D.; Shamseer, L.; Tetzlaff, J.M.; Akl, E.A.; Brennan, S.E.; et al. PRISMA 2020 Explanation and Elaboration: Updated Guidance and Exemplars for Reporting Systematic Reviews. BMJ 2021, 372, n160. [Google Scholar] [CrossRef] [PubMed]
  32. NIHR: National Institute for Health Research. PROSPERO International Prospective Register of Systematic Reviews. Available online: https://www.crd.york.ac.uk/prospero/ (accessed on 30 August 2020).
  33. Wells, G.A.; Shea, B.; O’Connell, D.; Paterson, J.; Welch, V.; Losos, M. Tugwell The Newcastle-Ottawa Scale (NOS) for Assessing the Quality of Nonrandomised Studies in Meta-Analyses. Available online: http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp (accessed on 28 November 2020).
  34. Higgins, J.P.T.; Thompson, S.G. Quantifying Heterogeneity in a Meta-Analysis. Stat. Med. 2002, 21, 1539–1558. [Google Scholar] [CrossRef] [PubMed]
  35. Review Manager(RevMan) [Computer Program]. The Cochrane Collaboration; Version 5.4; RevMan: New York, NY, USA, 2020. [Google Scholar]
  36. Shi, J.Q.; Copas, J.B. Meta-Analysis for Trend Estimation. Stat. Med. 2004, 23, 3–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Greenland, S.; Longnecker, M.P. Methods for Trend Estimation from Summarized Dose-Response Data, with Applications to Meta-Analysis. Am. J. Epidemiol. 1992, 135, 1301–1309. [Google Scholar] [CrossRef]
  38. DerSimonian, R.; Laird, N. Meta-Analysis in Clinical Trials Revisited. Contemp. Clin. Trials 2015, 45, 139–145. [Google Scholar] [CrossRef] [Green Version]
  39. Organisation for Economic Co-Operation and Development Pensions at a Glance-2013: Country Profiles. Available online: https://www.oecd.org/els/public-pensions/PAG2013-profile-country.pdf (accessed on 21 November 2020).
  40. Schernhammer, E.S.; Kroenke, C.H.; Laden, F.; Hankinson, S.E. Night Work and Risk of Breast Cancer. Epidemiology 2006, 17, 108–111. [Google Scholar] [CrossRef]
  41. Schernhammer, E. 0359 Nighshift Work and Breast Cancer Risk—Good News, Bad News? Occup. Environ. Med. 2014, 71, A121. [Google Scholar] [CrossRef] [Green Version]
  42. Tynes, T.; Hannevik, M.; Andersen, A.; Vistnes, A.I.; Haldorsen, T. Incidence of Breast Cancer in Norwegian Female Radio and Telegraph Operators. Cancer Causes Control 1996, 7, 197–204. [Google Scholar] [CrossRef]
  43. Davis, S.; Mirick, D.K.; Stevens, R.G. Night Shift Work, Light at Night, and Risk of Breast Cancer. J. Natl. Cancer Inst. 2001, 93, 1557–1562. [Google Scholar] [CrossRef] [Green Version]
  44. O’Leary, E.S. Shift Work, Light at Night, and Breast Cancer on Long Island, New York. Am. J. Epidemiol. 2006, 164, 358–366. [Google Scholar] [CrossRef] [Green Version]
  45. Menegaux, F.; Truong, T.; Anger, A.; Cordina-Duverger, E.; Lamkarkach, F.; Arveux, P.; Kerbrat, P.; Févotte, J.; Guénel, P. Night Work and Breast Cancer: A Population-Based Case-Control Study in France (the CECILE Study). Int. J. Cancer 2013, 132, 924–931. [Google Scholar] [CrossRef] [PubMed]
  46. Pronk, A.; Ji, B.T.; Shu, X.O.; Xue, S.; Yang, G.; Li, H.L.; Rothman, N.; Gao, Y.T.; Zheng, W.; Chow, W.H. Night-Shift Work and Breast Cancer Risk in a Cohort of Chinese Women. Am. J. Epidemiol. 2010, 171, 953–959. [Google Scholar] [CrossRef] [Green Version]
  47. Jones, M.E.; Schoemaker, M.J.; McFadden, E.C.; Wright, L.B.; Johns, L.E.; Swerdlow, A.J. Night Shift Work and Risk of Breast Cancer in Women: The Generations Study Cohort. Br. J. Cancer 2019, 121, 172–179. [Google Scholar] [CrossRef] [PubMed]
  48. Li, W.; Ray, R.M.; Thomas, D.B.; Davis, S.; Yost, M.; Breslow, N.; Gao, D.L.; Fitzgibbons, E.D.; Camp, J.E.; Wong, E.; et al. Shift Work and Breast Cancer among Women Textile Workers in Shanghai, China. Cancer Causes Control 2015, 26, 143–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Lie, J.-A.S.; Kjuus, H.; Zienolddiny, S.; Haugen, A.; Stevens, R.G.; Kjaerheim, K. Night Work and Breast Cancer Risk Among Norwegian Nurses: Assessment by Different Exposure Metrics. Am. J. Epidemiol. 2011, 173, 1272–1279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Fritschi, L.; Erren, T.C.; Glass, D.C.; Girschik, J.; Thomson, A.K.; Saunders, C.; Boyle, T.; El-Zaemey, S.; Rogers, P.; Peters, S.; et al. The Association between Different Night Shiftwork Factors and Breast Cancer: A Case–Control Study. Br. J. Cancer 2013, 109, 2472–2480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Tsc, L.A.; Wang, F.; Chan, W.C.; Wu, Ç.; Li, M.; Kwok, C.H.; Leung, S.L.; Yu, W.C.; Yu, I.T. 0076 Long-Term Nightshift Work and Breast Cancer Risk in Hong Kong Women: Results Update. Occup. Environ. Med. 2014, 71, A7–A8. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, F.; Yeung, K.L.; Chan, W.C.; Kwok, C.C.H.; Leung, S.L.; Wu, C.; Chan, E.Y.Y.; Yu, I.T.S.; Yang, X.R.; Tse, L.A. A Meta-Analysis on Dose-Response Relationship between Night Shift Work and the Risk of Breast Cancer. Ann. Oncol. 2013, 24, 2724–2732. [Google Scholar] [CrossRef] [PubMed]
  53. Ijaz, S.I.; Verbeek, J.H.; Seidler, A.; Lindbohm, M.-L.; Ojajärvi, A.; Orsini, N.; Costa, G.; Neuvonen, K. Night-Shift Work and Breast Cancer—A Systematic Review and Meta-Analysis. Scand. J. Work. Environ. Health 2013, 39, 431–447. [Google Scholar] [CrossRef] [PubMed]
  54. He, C.; Anand, S.T.; Ebell, M.H.; Vena, J.E.; Robb, S.W. Circadian Disrupting Exposures and Breast Cancer Risk: A Meta-Analysis. Int. Arch. Occup. Environ. Health 2015, 88, 533–547. [Google Scholar] [CrossRef]
  55. Tse, L.A. Project details. Research Grants Council of Hong Kong. Enquire Project Details by General Public. Available online: https://cerg1.ugc.edu.hk/cergprod/scrrm00542.jsp?proj_id=474811&old_proj_id=null&proj_title=&isname=&ioname=&institution=&subject=&pages=1&year=&theSubmit=474811 (accessed on 28 July 2020).
  56. Social Security Benefits Planner: Retirement. Retirement Age Calculator. Available online: https://www.ssa.gov/planners/retire/ageincrease.html (accessed on 28 November 2020).
  57. The State Council. The People’s Republic of China China Plans to Rise Age of Retirement. Available online: http://english.www.gov.cn/state_council/ministries/2016/03/01/content_281475299145295.htm (accessed on 10 November 2020).
  58. Atav, T.; Jongen, E.; Rabaté, S. Centraal Planbureau. In The Effects of the Increase in the Retirement Age in the Netherlands; Netherlands Bureau for Economic Policy Analysis (CPB): The Hague, The Netherlands, 2019; Available online: https://www.cpb.nl/sites/default/files/omnidownload/CPB-Discussion-Paper-408-The-effects-of-the-increase-in-the-retirement-age-in-the-Netherlands.pdf (accessed on 28 November 2020).
  59. The Government of Sweden. Ministry of Health and Social Affairs. Government Offices of Sweden Sustainable Pensions: Improved Basic Protection for Pensioners and a Gradual Increase in Retirement Age; Ministry of Health and Social Affairs: Stockholm, Sweden, 2019. Available online: https://www.government.se/articles/2019/06/1.---sustainable-pensions-improved-basic-protection-for-pensioners-and-a-gradual-increase-in-retirement-age/ (accessed on 28 November 2020).
  60. Check Your State Pension Age—GOV.UK. Available online: https://www.gov.uk/state-pension-age (accessed on 28 November 2020).
  61. Old Age Pensions. 2008. Available online: https://www.regjeringen.no/en/topics/pensions-and-welfare/innsikt/pensjoner/old-age-pensions/id449656/ (accessed on 28 November 2020).
  62. OECD; Ageing and Employment Policy Team Denmark. Key Policies to Promote Longer Working Lives. Country Note 2007 to 2017. Secretary General of the OECD. 2018. Available online: https://www.oecd.org/els/emp/Denmark%20Key%20policies_Final.pdf (accessed on 28 November 2020).
  63. Daniels, D.; Buckmaster, L.; Yeend, P. Parliament of Australia; Social Policy Section Social Security and Other Legislation Amendment (Pension Reform and Other 2009 Budget Measures) Bill; Issuing Body; Item ID: 154004-1001; Parliamentary Library: Canberra, Australia, 2009.
  64. GovHK Retirement Planning and Tips. Available online: https://www.gov.hk/en/residents/employment/mpf/employees/planandtips.htm (accessed on 28 November 2020).
  65. Papantoniou, K.; Pozo, O.J.; Espinosa, A.; Marcos, J.; Castano-Vinyals, G.; Basagana, X.; Pages, E.J.; Mirabent, J.; Martin, J.; Such Faro, P.; et al. Increased and Mistimed Sex Hormone Production in Night Shift Workers. Cancer Epidemiol. Biomark. Prev. 2015, 24, 854–863. [Google Scholar] [CrossRef] [Green Version]
  66. Beral, V. Million Women Study Collaborators Breast Cancer and Hormone-Replacement Therapy in the Million Women Study. Lancet 2003, 362, 419–427. [Google Scholar] [CrossRef]
  67. Collaborative Group on Hormonal Factors in Breast Cancer. Type and Timing of Menopausal Hormone Therapy and Breast Cancer Risk: Individual Participant Meta-Analysis of the Worldwide Epidemiological Evidence. Lancet 2019, 394, 1159–1168. [Google Scholar] [CrossRef]
  68. Khan, S.A.; Stickles, S. Cell Proliferation and Apoptosis in the Normal Breast Epithelium of Pre, Peri, and Postmenopausal Women. In Hormonal Carcinogenesis III; Springer: New York, NY, USA, 2001; pp. 418–423. [Google Scholar]
  69. Milholland, B.; Auton, A.; Suh, Y.; Vijg, J. Millholland Age-Related Somatic Mutations in the Cancer Genome. Oncotarget 2015, 6, 24627–24635. [Google Scholar] [CrossRef] [Green Version]
  70. Trichopoulos, D.; Lagiou, P.; Adami, H.-O. Towards an Integrated Model for Breast Cancer Etiology: The Crucial Role of the Number of Mammary Tissue-Specific Stem Cells. Breast Cancer Res. 2005, 7, 13–17. [Google Scholar] [CrossRef] [Green Version]
  71. Terry, M.B.; Michels, K.B.; Brody, J.G.; Byrne, C.; Chen, S.; Jerry, D.J.; Malecki, K.M.C.; Martin, M.B.; Miller, R.L.; Neuhausen, S.L.; et al. Trentham-Dietz A12; Breast Cancer and the Environment Research Program (BCERP). Environmental Exposures during Windows of Susceptibility for Breast Cancer: A Framework for Prevention Research. Breast Cancer Res. 2019, 21, 96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Page, M.J.; Higgins, J.P.T.; Clayton, G.; Sterne, J.A.C.; Hróbjartsson, A.; Savović, J. Empirical Evidence of Study Design Biases in Randomized Trials: Systematic Review of Meta-Epidemiological Studies. PLoS ONE 2016, 11, e0159267. [Google Scholar] [CrossRef] [PubMed]
  73. Härmä, M.; Koskinen, A.; Ropponen, A.; Puttonen, S.; Karhula, K.; Vahtera, J.; Kivimäki, M. Validity of Self-Reported Exposure to Shift Work. Occup. Environ. Med. 2017, 74, 228–230. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study flow diagram.
Figure 1. Study flow diagram.
Cancers 13 05952 g001
Figure 2. Forest plot of the meta-analysis on long-term nightshift work (≥15 years of NSW) and BC risk.
Figure 2. Forest plot of the meta-analysis on long-term nightshift work (≥15 years of NSW) and BC risk.
Cancers 13 05952 g002
Figure 3. Funnel plot of reports included in the meta-analysis of long-term NSW and BC risk. Dashed lines represent pseudo 95% confidence limits.
Figure 3. Funnel plot of reports included in the meta-analysis of long-term NSW and BC risk. Dashed lines represent pseudo 95% confidence limits.
Cancers 13 05952 g003
Figure 4. Dose–response meta-analysis with BC risk trends across categories of NSW duration. (For this analysis, the exposure categories of 10–20 years and ≥20 years of the EPIC-Oxford study were combined). Black line and shadowed area represent overall pooled RR and 95%CI.
Figure 4. Dose–response meta-analysis with BC risk trends across categories of NSW duration. (For this analysis, the exposure categories of 10–20 years and ≥20 years of the EPIC-Oxford study were combined). Black line and shadowed area represent overall pooled RR and 95%CI.
Cancers 13 05952 g004
Figure 5. Forest plot of the meta-analysis on long-term NSW (≥15 years) and pre-/perimenopausal BC.
Figure 5. Forest plot of the meta-analysis on long-term NSW (≥15 years) and pre-/perimenopausal BC.
Cancers 13 05952 g005
Figure 6. Forest plot of the meta-analysis on long-term NSW (≥15 years) and postmenopausal BC.
Figure 6. Forest plot of the meta-analysis on long-term NSW (≥15 years) and postmenopausal BC.
Cancers 13 05952 g006
Figure 7. Forest plot of the meta-analysis on recent long-term NSW (≥15 years) and BC risk.
Figure 7. Forest plot of the meta-analysis on recent long-term NSW (≥15 years) and BC risk.
Cancers 13 05952 g007
Table 1. Main characteristics of studies included in this or in previous meta-analyses and pooled analyses on long-term NSW exposure and BC.
Table 1. Main characteristics of studies included in this or in previous meta-analyses and pooled analyses on long-term NSW exposure and BC.
(a) Cohort Studies
Individual Studies: Main CharacteristicsPublished Meta-Analyses or Pooled Analyses on Long-Term NSW
Wang
(2013)
[52]
Jia
(2013)
[10]
Ijaz
(2013)
[53]
He
(2015)
[54]
Lin
(2015)
[11]
Travis
(2016)
[12]
Cordina
Duverger
(2018) [30]
Present
Meta-Analysis
Present Dose–Response Meta-Analysis
Author, Study, Year of PublicationCountryBase
Population
Ascert. of Breast Cancer≥15 yrs NSWNo NSWTotal Cases5Follow-Up
(yrs)
Lost to Follow-UpQD-R
per
5 yrs
≥15
yrs
D-R
per
5 yrs
D-R
per
10 yrs
10–20
yrs
>20
yrs
>20
yrs
>30
yrs
10–20
yrs
>20
yrs
≥15 yrsD-R
per
10 yrs
N 1Ca 2N 3Ca 4
Schernhammer NHS I (2001) [20]USAOccup (nurses)SR/HOSP max 10<10%7XXXXXXXX--------
Wegrzyn NHS I (2017) [14]USAOccup (nurses)SR/HOSP580442731,74623822809max 24<10%8--------------------XX
Schernhammer NHS I and II (2014) [41]USAOccup (nurses)SR/HOSP max 22<10%8------ --X------------
Schernhammer NHSII (2006) [40]USAOccup (nurses)SR/HOSP max 12<10%7XXXXXXX----------
Wegrzyn NHS II (2017) [14]USAOccup (nurses)SR/HOSP1623543,529950985max 24<10%8--------------------XX
Pronk (2010) SWHS based on JEM [46] 6CHNGeneral popREG av 9.03%6XX--XX--------------
Ponk (2010) SWHS based on SR [46]CHNWorking popREG57201951,238276295av 4.43%7----X------XX----XX
Knutsson WOLF (2013) [28]SEGeneral popREG av 12.4n/a6--------X--------------
Koppes Dutch Labor Force (2014) [22]NLWorking popREGn/an/a255,9002312n/aav 6.9n/a4--------XXX------X--
Åkerstedt Salt Study (2015) [13]SETwinsREG305189674354372av 8.7n/a6------------X------XX
Travis (2016) EPIC—Oxford [12]UKGeneral popREG461119,289153154av 3.1n/a7------------X------XX
Travis Million Women Study (2016) [12]UKGeneral popREG964789450,23241364225av 2.6n/a7------------XX----XX
Jones Generations Study (2019) [47]UKGeneral popREGn/a6084,88818451905median 9.54%7--------------------XX
(b)Nested Case-Control and Case-Control Studies
Individual Studies: Main CharacteristicsPublished Meta-Analyses or Pooled Analyses on Long-Term NSW
Wang
(2013)
[52]
Jia
(2013)
[10]
Ijaz
(2013)
[53]
He
(2015)
[54]
Lin
(2015)
[11]
Travis
(2016)
[12]
Cordina
Duverger
(2018) [30]
Present
Meta-Anal.
Present Dose–Response Meta-Anal.
CountryBase
Population
Control
Selection
Case
Ascert
RatioLT-NSW Exposed (n)Participation (%)QD-R
per
5 yrs
15
yrs
D-R
per
5 yrs
D-R
per
10 yrs
10–20
yrs
>20
yrs
>20
yrs
>30
yrs
10–20
yrs
>20
yrs
15
yrs
D-R
per
10 yrs
ControlsCasesControlsCases
NESTED CASE–CONTROL STUDIES
Tynes (1996) [42]NOROccup (radio/tele)n/aREG4–7:1 n/an/a5----X------------------
Lie (2006) [21]NOROccup (nurses)IDSREG4:1417125n/an/a5X--X--------------XX
Lie (2011) [49]NOROccup (nurses)IDSREG1.5:123117965747XXXX------------XX
Hansen (2012a) [17]DKOccup (nurses)IDSREG4:11243991927X--XX------------XX
Hansen (2012b) [18]DKOccup (military)IDSREG4:1291261678XXXX------------XX
Li (2015) [48]CHNOccup (textile)RNDREG2.8:1n/a576n/an/a5------------XX----XX
CASE–CONTROL STUDIES
Davis (2001) [43]USAGeneral popRNDREG1:1 75785X--XX----------------
O’ Leary (2006) LIBCSP/EBCLIS [44]USADiverse sourcesRNDHOSP1:1 83875X--XX----------------
Pesch GENICA (2010) [15]DEGeneral popRNDHOSP1:151267886XXXX--------XXXX
Menegaux (2013) [45]FRGeneral popRNDHOSP1:1 76796----XX--------XX----
Fritschi (2013) [50]AUSGeneral popRNDREG1.5:1538441587------X--------XXXX
Grundy (2013) [16]CABC screeningRNDMixed 71:15365V: 54,
K: 49
V: 57 K:598------X--------XXXX
Tsc (2014 [51,55]CHNHospitalizedCONSECHOSP1:1n/an/a93 (55)91 (55)6--------------------X--
Papantoniou (2015) [19]SPGeneral popRNDHOSP1:1979152727----------------XXXX
Ascert (ascertainment), yrs (years), Q (quality index of the Newcastle–Ottawa scale; score range: 0–9), D-R (dose–response meta-analysis), Occup (occupational), SR (self-reported), SR/HOSP (self-reported/hospital), max (maximum), SWHS (Shanghai Women’s Health Study), pop (population), REG (registry), av (average), n/a (not available), Ratio (ratio controls–cases), LT-NSW (long-term nightshift work), radio/tele (radio/telegraph operators), IDS (incidence density sampling), RND (randomly), V (controls and cases from Vancouver), K (controls and cases from Kensington), CONSEC (consecutive) 1 Number of women exposed to ≥15 yrs of NSW. 2 BC cases among women exposed to ≥15 yrs of NSW. 3 Number of women not exposed to NSW. 4 Number of BC cases among women not exposed to NSW. 5 Total number of cases included in the present meta-analysis. 6 This study used a job-exposure matrix and self-reporting to ascertain NSW exposure. 7 Registry + BC screening.
Table 2. Exposure characteristics and age of participants in relation to retirement age in each study country, in the included publications.
Table 2. Exposure characteristics and age of participants in relation to retirement age in each study country, in the included publications.
(a) Cohort Studies
Author, Year of Publication, StudyExposure DefinitionMin.
Freq. of
NSW/Month
Nonexposure
Definition
Exposure AscertainmentNSW
Duration
(yrs)
Age at Start of Follow-UpFollow-Up
Period
Mean Age at Start of Follow-Up + Full Follow-Up Period (yrs)Retirement
Age (yrs)
Estimated
yrs since
Retirement at the End of Follow-Up
Inclusion in our (sub)-Meta-Analyses
Direction and MomentSourceMeanRangeAll BCPremen
BC
Postmen
BC
“Recent Long-Term NSW”
Schernhammer (2001), NHS I [20]Yrs of rotating night shifts with at least
3 nights/month
3/mNever rotating NSW 1Retro. at start of follow-upSR≥1557.142–671988–199867.165 [56]2.1--XXX 2
Wegrzyn (2017), NHS I [14]See Schernhammer (2001)3/mNever rotating NSW 1Retro. at start of follow-upSR≥1557.142–671988–201281.165 [56]16.1X------
Wegrzyn (2017), NHS II update [14]See Schernhammer (2001)3/mNever rotating NSW 1Retro. at start of follow-up and 5 updatesSR≥203925–421989–201363.066 [56]Still working, but no NSWXX--X 3
Pronk (2010) SWHS [46]Starting work after 10 p.m.3/mNever night shiftsRetro. and two yrsafter start of follow-upSR>1752.540–702000–200761.550–55 [57]6.5–11.5X------
Koppes (2014), Dutch Labor Force [22]Occ/regular current work at night, between 0 and 6 a.m. combined with yrs of job tenureNoneNo current NSW≥20 yrs same jobBefore start of follow-upJEM≥2038 *15–641996–200951.965 [58]Still workingX----X
Åkerstedt (2015), Salt twin study [13]Yrs working at nights at least now and thenNoneNot worked nightsRetro. at start of follow-upSR21–4551.641–601998–201063.665 [59]Still workingX----X
Travis (2016), EPIC-Oxford [12]Worked regularly at night, on NSW or on call at night1/mNever night shiftsRetro. at start of follow-upSR≥2057.8n/a2010–201360.962 [60]Still workingX--XX
Travis (2016), MWS [12]Ever regularly worked at night or on night shifts at any time between 00 and 06 h3/mNever night shiftsRetro. at start of follow-upSR≥2068.6n/a2011–201371.260 [60]11.2X--X--
Jones (2019), Generations Study [47]Any job that regularly involved work in the late evening or night (between 10 p.m. and 7 a.m.)NoneNot being an NSWer within the last ten yrsAt start of follow-up and 1 update 6 yrs after recruitmentSR≥2045 **35–552003–201860.560–68 [60]Still workingX----X
(b) Nested Case-Control and Case-Control Studies
NESTED CASE CONTROL STUDIES
Author,
Year of Publication
Exposure DefinitionMin. Freq of
NSW/month
Non-Exposure
Definition
Source of Expo. AscertNSW
Duration
(yrs)
Age at Start of Follow-UpMean Age at dx of BC (yrs) [Range]Mean Age at Start of Follow-Up + Full Follow-Up Period (yrs)Retirement
Age (yrs)
Estimated
yrs since
Retirement at BC dx
Inclusion in our (sub) Meta-Analyses
MeanRangeAll BCPremen
BC
Postmen
BC
“Recent Long-Term NSW”
Lie (2006) [21]Work at infirmaries (based on Nurse Registry and census, only considering time after grad)NoneManagerial, teaching, physiotherapy, outpatient department worksite, other than infirmariesJEM15–2939.7 *27–8554 * (27–85)61.767 [61]Still workingXXXX
Lie (2011) [49]Rotating and permanent night work3/mNever night workSR≥15n/a20–7054.5 (35–74)n/a67 [61]Still workingX----X
Hansen (2012a) [17]Working ≥1 yr during hours between 7 p.m. and 9 a.m. not including overtimeNoneNever (<1 year) “after midnight shifts”SR≥2054 *31–69n/a5665 [62]Still working ***X----X
Hansen (2012b) [18]Working ≥1 yr during hours between 5 p.m. and 9 a.m., not including overtimeNoneNo yrs (<1 yr) of NSWSR≥15n/a22–75n/an/a n/aX------
Li (2015) [48]Jobs involving rotating NSW (22.00–06.00 h) according to factory processesNoneNo rotating NSW 4JEM20–2748.9 *30–6653.459.950 [57]3.4XXXX 5
CASE CONTROL STUDIES
Author
Year of publication, study
Exposure DefinitionMin.
freq. of NSW/Month
Nonexposure DefinitionSource of expo. AscertNSW Duration(yrs) Mean Age at dx of BC (yrs) [Range]Year of
BC- dx
Retirement
Age (yrs)
Estimated yrs since Retirement at BC DiagnosisInclusion in our (sub) Meta-Analyses
All BCPremen
BC
Postmen
BC
“Recent Long-Term NSW”
Pesch (2010) [15]Ever having worked in NSW for ≥1 year and working the full time period between 0.00 and 5.00NoneEmployed but never in shift work (day shifts only)SR≥20 54 (42–62)2000–200465 [63]Still workingX----X
Fritschi (2013) [50]Worked any number of hours between 0.00 and 5 a.m. (graveyard shift)NoneNever graveyard shiftSR≥20 57 * (18–80)2009–201164.5 [63]Still workingX----X
Grundy (2013) [16]Jobs that started/ended between 11 p.m. and 7 a.m.NoneNo yrs in jobs with start or end between 11 p.m. and 7 a.m.SR≥15 57 (<80)2005–201065 [39]Still workingXXXX
Tsc (2014) [51,55]Nightshift at least once per month for ≥1 year1/mPermanent day work 6SR≥15 55 (40–69)2011–201365 [64]Still workingXXXX
Papantoniou (2015)
MCC-Spain [19]
Partly/entirely working between 0 and 6 a.m.3/mNever night workSR≥15 56 (23–85)2008–201365 [39]Still workingX----X
Min (minimum), freq (frequency), yrs (years), Retirement age (legal retirement age in the study country at the moment of the study), Premen (premenopausal), Postmen (postmenopausal), m (month), retro (retrospective), SR (self-reported), NSWer (nightshift worker), Occ (occasionally), JEM (job-exposure matrix), n/a (not available), dx (diagnosis), yrs (years), expo ascert. (exposure ascertainment), dx (diagnosis), grad (graduation). 1 May include permanent NSW. 2 Only the RR for premenopausal was included. 3 The first ten years of follow-up were included, as no NSW was assumed since 2009 by study authors. 4 No jobs in the factory involved permanent NSW. 5 The BC risk for women < 50 years was included. 6 Kindly reported by e-mail from the author, * Calculated based on given categories of age or birth year; ** Median age; *** Mean age at BC diagnosis was missing in this study, but the mean age at start of follow-up + years of full follow-up was under the retirement age.
Table 3. Summary results of the meta-analyses, dose-response meta-analyses and sub-meta-analyses for the association between long-term NSW and BC risk.
Table 3. Summary results of the meta-analyses, dose-response meta-analyses and sub-meta-analyses for the association between long-term NSW and BC risk.
Meta-AnalysisWang
(2013) [52]
Jia
(2013) [10]
Ijaz
(2013) [53]
He
(2015) [54]
Lin
(2015) [11]
Travis
(2016) [12]
Cordina-Duverger
(2018) [30]
Present
Study
Inclusion
Criteria
(Study Types)
Cohort,
nested
case–control, and
case–control
Cohort,
nested
case–control
and
case–control
Cohort,
nested
case–control, and
case–control
Cohort,
nested
case–control
and
case–control
CohortCohort
and
nested
case–control
Case–controlCohort,
nested
case–control, and
case–control
Duration of exposure
10–20 years 1.07
(1.01–1.14)
0.98
(0.78–1.22)
≥15 years 1.15
(1.03–1.29)
1.13
(1.01–1.27)
≥20 years 1.09
(1.01–1.17)
1.01
(0.93–1.10)
1.10
(0.87–1.39)
≥30 years 1.00
(0.87–1.14)
Dose–response meta-analysis
per 5 years1.03
(1.01–1.05)
1.05
(1.01–1.10)
per 10 years 1.06
(0.98–1.15)
1.05
(0.94–1.09)
Sub-meta-analyses
Premenopausal women
10–20 years 1.05
(0.74–1.47)
≥15 years 1.27
(0.96–1.68)
≥20 years 1.34
(0.85–2.13)
Postmenopausal women
10–20 years 0.92
(0.68–1.23)
≥15 years 1.05
(0.90–1.24)
≥20 years 1.04
(0.80–1.36)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schwarz, C.; Pedraza-Flechas, A.M.; Pastor-Barriuso, R.; Lope, V.; de Larrea, N.F.; Jiménez-Moleón, J.J.; Pollán, M.; Pérez-Gómez, B. Long-Term Nightshift Work and Breast Cancer Risk: An Updated Systematic Review and Meta-Analysis with Special Attention to Menopausal Status and to Recent Nightshift Work. Cancers 2021, 13, 5952. https://doi.org/10.3390/cancers13235952

AMA Style

Schwarz C, Pedraza-Flechas AM, Pastor-Barriuso R, Lope V, de Larrea NF, Jiménez-Moleón JJ, Pollán M, Pérez-Gómez B. Long-Term Nightshift Work and Breast Cancer Risk: An Updated Systematic Review and Meta-Analysis with Special Attention to Menopausal Status and to Recent Nightshift Work. Cancers. 2021; 13(23):5952. https://doi.org/10.3390/cancers13235952

Chicago/Turabian Style

Schwarz, Christine, Ana María Pedraza-Flechas, Roberto Pastor-Barriuso, Virginia Lope, Nerea Fernández de Larrea, José Juan Jiménez-Moleón, Marina Pollán, and Beatriz Pérez-Gómez. 2021. "Long-Term Nightshift Work and Breast Cancer Risk: An Updated Systematic Review and Meta-Analysis with Special Attention to Menopausal Status and to Recent Nightshift Work" Cancers 13, no. 23: 5952. https://doi.org/10.3390/cancers13235952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop