Introduction
Shift work has been referred to as “a work activity scheduled outside standard daytime hours, where there may be a handover of duty from one individual or work group to another” (Health and Safety Executive 2006). Typically, shift work has been associated with industries that require 24-hour operation, such as essential public services or for practical purposes. In recent years there has been an increase in the number of shift workers in other industries (Health and Safety Executive 2006) with approximately 19% of the working population engaged in shift work as their main job (NHS Digital 2013) in the UK.
It has been argued that the introduction of shift work within wider industries does not consider the health and wellbeing costs to the individual shift worker (Folkard and Tucker 2003). Much of the research that investigates the health impacts of shift work has previously centred around circadian disruption particularly among night shift workers, which can result in disturbed sleep and excessive sleepiness during the work shift (Wright et al. 2013). However, adverse health behaviours have been found among shift workers which also put them at higher risk of disease (Dorrian et al. 2017; Loef et al. 2018; Suwazono et al. 2008) and shift work has been associated with higher risk of diseases, for example, type 2 diabetes, cardiovascular disease, stroke and cancer (Kecklund and Axelsson 2016), than non-shift workers.
Recent studies have investigated epigenetic changes in the form of circulating DNA methylation (DNAm) as an objective measure for evaluating the potential health impact of shift work (White et al. 2019). While these studies have typically focused on assessing individual sites in the genome (cytosine-phosphate-guanine “CpG” sites), recently DNAm scores derived from methylation levels at numerous CpG sites across the epigenome have been developed which can act as proxies for lifestyle exposures and may predict health outcomes (McCartney et al. 2018). Self-reported health behaviours (especially one-off measures) are subject to measurement error and bias (Prince et al. 2020), and using more objective measures of long-term exposure could identify more robust associations with shift work.
One group of DNAm scores aims to capture the epigenetic clock. DNA methylation age (DNAm Age) has been derived to provide an accurate estimate of biological age across a range of tissues, and at different life stages (Fransquet et al. 2019). Estimated DNAm Age exceeding true chronological age is known as “age acceleration” and studies suggest that DNAm age acceleration is associated with age-related health outcomes independent of chronological age (Marioni et al. 2015). As well as the clocks developed based on age, more recent “second generation” epigenetic clocks have been developed based on lifestyle factors and biomarkers which have been found to be highly predictive of both health and lifespan (Levine et al. 2018; Lu et al. 2019). Other DNAm scores which predict modifiable health, lifestyle and socio-economic factors include scores developed for alcohol consumption, smoking status, BMI and education (McCartney et al. 2018).
This study aimed to investigate associations of evening/night shift work participation and a series of blood based DNAm predictors of ageing, BMI, smoking, alcohol and education within the Generation Scotland (GS) and Understanding Society (UKHLS) studies.
Materials and Methods
Generation Scotland (GS)
The Generation Scotland: Scottish Family Health Study is a prospective cohort study comprising participants from the general population across five regions of Scotland. The recruitment protocol and cohort characteristics are described in detail elsewhere (Smith et al. 2013) and in the Supplementary Methods.
Blood DNAm was profiled using the Infinium MethylationEPIC BeadChip (Illumina Inc.) in two sample sets from GS: Set 1 and Set 2 (Supplementary Methods).
Methylation data were available for 777,193 CpGs measured in Set 1 (n = 2,578 unrelated individuals) and 773,860 CpGs measured in Set 2 (n = 4,450 unrelated individuals). The data sets had DNAm profiled at separate time points and quality control and normalisation was carried out separately. Details on the lifestyle variables and how they were derived can be found in the Supplementary Methods.
Understanding Society (UKHLS)
The UK Household Longitudinal Study (UKHLS) (also known as Understanding Society) is a longitudinal panel survey of 40 000 UK households from England, Scotland, Wales and Northern Ireland. The recruitment protocol and cohort characteristics are described in detail elsewhere (Chen et al. 2016; Hannum et al. 2013) and in the Supplementary Methods.
Blood DNAm was profiled using the same Infinium MethylationEPIC BeadChip (Illumina Inc.) as used by GS. Methylation data was available for 837,487 CpGs in 1,175 individuals (Supplementary Methods). Details on the lifestyle variables and how they were derived can be found in the Supplementary Methods.
Evening/Night Shift Work
In both Generation Scotland and UKHLS we defined evening/night shift work based on work reported outside normal daytime working hours (typically 9:00 AM-5:00 PM).
In Generation Scotland, participants were asked at baseline (when the samples for DNAm profiling were obtained) how many hours in a typical week they would work in the evening or overnight, from which a binary variable for “evening/night shift work” was generated based on those working ≥20 hours/week compared with those reporting to work 0 hours in the evening or overnight (i.e. excluding 1–19 hours of evening/night shift work).
In UKHLS, participants were asked to describe the times of day they usually worked at multiple time points, with those who reported working evenings or at night classed as “evenings/night shift workers.” Current evening/night shift work was defined based on reported shifts at the timepoint when the samples for DNAm profiling were obtained. More details of the evening/night shift work variables and how they were derived can be found in the Supplementary Methods.
DNA Methylation Scores
We derived four DNAm scores related to BMI, smoking, alcohol consumption and education which were based on CpG sites identified in a previous study (McCartney et al. 2018). Details of the DNAm scores are shown in Table 1. For each individual, DNAm scores were calculated as the sum of methylation values at each CpG multiplied by the effect sizes obtained in a previous study (McCartney et al. 2018).
We also derived six epigenetic biomarkers of ageing using previously published approaches (Chen et al. 2016; Hannum et al. 2013; Horvath 2013; Levine et al. 2018). (Supplementary Methods). In each case, age acceleration was defined as the residual obtained from regressing predicted age, as estimated by the epigenetic clock, on chronological age. This measure of age acceleration is independent of chronological age.
All methylation scores were standardized (mean = 0, standard deviation = 1) in both GS and UKHLS.
Statistical Analyses
We first assessed whether lifestyle factors were associated with evening/night shift work. We performed linear regression of each phenotype (as the outcome variable) in relation to evening/night shift work (as the exposure variable). This was with the exception for education, which was treated as the exposure variable given that education precedes evening/night shift work. In GS, two models were run: Model 1 with adjustment for age and sex, and Model 2 with additional adjustment for other self-reported phenotypes (e.g. for smoking, the model was also adjusted for BMI, alcohol, and education). In UKHLS, two models were run: Model 1 with adjustment for age, sex, blood processing day and batch and Model 2 with additional adjustment for other self-reported phenotypes.
We subsequently assessed whether the lifestyle factors could be proxied with DNAm scores within GS and UKHLS. We performed linear regression of the phenotypes (exposure variable) and lifestyle DNAm scores for alcohol, smoking, education, and BMI (outcome variable). For these models, a logistic regression was performed. In GS, two models were run: Model 1 with adjustment for age, sex and 20 methylation principal components (PCs) – the latter are used to identify and remove unwanted sources of variation (e.g. batch effects). Model 2 was additionally adjusted for other self-reported phenotypes (e.g. for smoking related DNAm scores, the model was also adjusted for BMI, alcohol, and education). In UKHLS, two models were run: Model 1 with adjustment for age, sex, blood processing day and batch; and Model 2 with additional adjustment for other self-reported phenotypes.
We finally assessed whether evening/night shift work was related to the lifestyle DNAm scores. We performed linear regression of each lifestyle methylation score (outcome) in relation to evening/night shift work (exposure). This again was with the exception for education, which was treated as the exposure. In GS, three models were run: Model 1 with adjustment for age, sex and 20 methylation PCs; Model 2 had additional adjustments for other self-reported phenotypes (e.g. for smoking related DNAm scores, model adjusted for BMI, alcohol, and education); and Model 3 with further adjustment for the corresponding phenotype (e.g. for smoking-related DNAm scores, the model was also adjusted for smoking). This was performed to examine whether the DNAm score was associated with shift work independently of the corresponding phenotype.
In UKHLS, three models were run: Model 1 with adjustment for age, sex, blood processing day and batch; Model 2 with additional adjustment of other self-reported phenotypes; and Model 3 with further adjustment for the corresponding phenotype.
We assessed whether evening/night shift work was related to any of the six epigenetic ageing measures. We performed linear regression of each epigenetic age acceleration (EEA) measure (as the outcome variable) in relation to evening/night shift work.
In GS, two models were run: Model 1 with adjustment for sex and 20 methylation PCs and Model 2 with additional adjustment for smoking, BMI, education, and alcohol. In UKHLS, two models were run: Model 1 with adjustment for sex, blood processing day and batch; and Model 2 with additional adjustment for smoking, BMI, education, and alcohol.
We conducted a series of meta-analyses of GS Sets 1 and 2 and UKHLS using an inverse-variance weighted fixed effects approach. For this, we used the binary measure of evening/night shift work derived in GS and the current evening/night shift work variable derived in UKHLS. The I2 statistic was used to assess heterogeneity across the studies.
We conducted meta-analyses of: i) the self-reported or conventionally measured (in the case of BMI) lifestyle factors in relation to evening/night shift work, ii) the lifestyle DNAm scores in relation to evening/night shift work and iii) the epigenetic age acceleration measures in relation to evening/night shift work.
Results
Baseline Characteristics
Summary characteristics of the participants from GS and UKHLS are presented in Table 2. 7,028 individuals with DNAm data were included from GS (n = 2,578 in Set 1 and n = 4,450 in Set 2) and 1,175 individuals from UKHLS. Set 1 and 2 of GS had a younger mean age (50.0 ± 12.5 and 51.4 ± 13.2) compared to UKHLS (58.0 ± 15.0). There was a comparable balance of men and women across the three datasets: 38.6% males in GS Set 1, 43.7% in GS Set 2 and 41.6% in UKHLS. BMI was also broadly comparable: 27.4 ± 5.5 kg/m2, 26.8 ± 5.0 kg/m2 and 28.1 ± 6.2 kg/m2, in GS Set 1, GS Set 2 and UKHLS, respectively. GS Set 1 had more current smokers whilst UKHLS had more former smokers and GS Set 2 had more never smokers. There was a higher proportion of daily drinkers in UKHLS (16.0%) compared with GS (12.4% in Set 1 and 13.2% in Set 2). There was also a higher proportion of less than monthly drinkers in UKHLS (25.2%) compared with GS (17.2% in Set 1 and 16.2% in Set 2). Years of full-time education was comparable between Set 1 and 2 of GS (13.6 ± 3.4 and 13.8 ± 3.4, respectively) whilst in UKHLS it was slightly lower (12.3 ± 5.1). In UKHLS, 1.6% of participants (n = 18) were currently working evening/night shifts while 8.8% (n = 103) had reported working evening/night shifts over the previous 11 years. In GS Set 1, 8.1% (n = 127) of participants reported working at evening/night for >20 hours per week at the time of sampling, compared with 7.9% (n = 193) in GS Set 2.
Lifestyle Factors
How Does Evening/Night Shift Work Relate to Conventionally Measured Lifestyle Factors?
In GS, evening/night shift work was inversely associated with alcohol frequency in both GS Set 1 and GS Set 2 (Model 1: Effect= − 0.056; 95% CI − 0.095, −0.017; p = 0.006 category change per hour for GS Set 1, and − 0.039; − 0.072, − 0.006; p = 0.019 for GS Set 2), although associations attenuated with adjustment for the other self-reported phenotypes in GS Set 2 (Model 2: − 0.024; − 0.057, 0.009; p = 0.151) (Table S1). Evening/night shift work was positively associated with smoking status in GS Set 1 (Model 1: 0.024; 0.006, 0.042; p = 0.009 category change per hour) but were more weakly associated in GS Set 2 (Model 2: 0.012; −0.006, 0.030; p = 0.216), and when adjusted for the other phenotypes (Table S1). There were positive associations between evening/night shift work and BMI in both GS Set 1 (Model 1: 0.15; 0.03, 0.28; p = 0.02 kg/m2 per hour) and GS Set 2 (Model 2: 0.23; 0.14, 0.32; p = 1 × 10−6), although these associations attenuated when adjusted for the other phenotypes. There was an inverse association between evening/night shift work and education in GS Set 2 (Model 1: −0.07; −0.13, −0.02; p = 0.01 hour per year) which was less apparent in GS Set 1 (Model 1: −0.05; −0.12, 0.03; p = 0.21), although there was stronger evidence of association of the binary evening/night shift work variable with education, based on in both datasets (Table S1).
In UKHLS, no strong associations were found between current, ever, or previous evening/night shift work and any of the lifestyle phenotypes, although effect estimates were typically in the same direction as in GS (Table S2).
In the GS-UKHLS combined meta-analysis there was evidence of a positive association between evening/night shift work and BMI (0.79; 0.02, 1.56; p = 0.04 kg/m2 difference between those who did and did not work evening/night shifts) and an inverse association between evening/night shift work and education (−0.18; −0.30, −0.07; p = 0.002 log odds per year of education) in the fully adjusted models (Model 2) (Figure 1). There was little evidence for an association between evening/night shift work and either alcohol intake or smoking status in a meta-analysis across the three datasets (0.00; −0.19, 0.20; p = 0.97 and 0.04; −0.19, 0.27; p = 0.73 category change between those who did and did not work evening/night shifts, respectively). There was little evidence for heterogeneity between the study estimates (I2 < 5%). Results of the minimally adjusted models (Model 1) are shown in Figure S1).
How are DNA Methylation Biomarkers Associated with Lifestyle Factors in GS and UKHLS?
There were positive associations between each DNAm score and its respective lifestyle phenotype in both GS and UKHLS (Table S3 and S4, respectively).
Is Evening/Night Shift Work Associated with DNA Methylation Biomarkers?
In GS, there were no clear associations between evening/night shift work and the alcohol DNAm score (Table S5). Evening/night hours were positively associated with the smoking DNAm score in both GS Set 1 and GS Set 2 (Model 1: 0.04 SD; 0.02, 0.06; p = 1.28 × 10−4 and 0.02 SD; 0.00, 0.04; p = 0.01), although associations attenuated with further covariate adjustment in Models 2 and 3. Number of evening/night hours was also positively associated with the BMI DNAm score in both GS Set 1 and GS Set 2 (Model 1: 0.03 SD; 0.01, 0.05; p = 0.01 and 0.03 SD; 0.01, 0.04; p = 0.04). Associations attenuated with further covariate adjustment in Models 2 and 3 (Table S5). Number of evening/night hours was inversely associated with the education DNAm score in GS Set 1 (Model 1: −0.20 hours; −0.33, −0.08; p = 0.001) and more weakly in GS Set 2 (Model 1: −0.09 hours; −0.19, 0.01; p = 0.08). With further covariate adjustment, the association persisted in GS Set 1 but was attenuated in GS Set 2 (Table S5).
None of the evening/night shift work measures were strongly associated with the lifestyle DNAm scores when UKHLS was assessed independently (Table S6). When estimates from GS were combined in a meta-analysis with UKHLS, there was little evidence of evening/night shift work associations with the alcohol and BMI DNAm scores in the fully adjusted models (Model 3) (Figure 2). However, there was some evidence of an inverse association between evening/night shift work and the education DNAm score in all three datasets (−0.24; −0.49, 0.01; p = 0.06 log odds per year) (Figure 2). There was also weak evidence for a positive association between evening/night shift work and the smoking DNAm score (0.06 SD, −0.03, 0.15; p = 0.18). There was little evidence for heterogeneity between the study estimates in the meta-analysis (I2 = 0%). Results of the minimally adjusted models (Models 1 and 2) are shown in Figures S2 and S3.
Epigenetic Ageing
There was evidence of an association between evening/night shift work and GrimAge in GS, but not for the other epigenetic clocks. For GrimAge, the number of evening/night hours was associated with higher age acceleration in both GS Set 1 (Model 1: 0.19 years, 0.11, 0.27; p = 1.18 × 10−5) and GS Set 2 (Model 1: 0.10 years, 0.04, 0.16; p = 9.45 × 10−4) (Table S7).
Associations with GrimAge acceleration remained, although were partially attenuated on adjustment for smoking, BMI, education and alcohol (Table S7). None of the evening/night shift work measures were strongly associated with epigenetic age acceleration in UKHLS (Table S8).
In the GS-UKHLS meta-analysis, evening/night shift work was associated with a 0.80 year (0.42, 1.18; p < 0.001) increase in GrimAge acceleration (Figure 3). There was also weak evidence of association with PhenoAge acceleration (0.46 years; 0.00, 0.92; p = 0.05). The other four epigenetic clocks showed limited evidence of association. There was low heterogeneity between the study estimates in the meta-analysis (I2 < 50%).
Discussion
We conducted analyses to investigate associations between evening/night shift work and both phenotypic and DNAm markers in two cohorts. When we assessed phenotypic traits, we found that evening/night shift work was associated with higher BMI and lower education. When assessing DNAm predictors of the same traits, there was similar evidence of association of evening/night shift work with BMI and education DNAm scores. While the association with the BMI score was attenuated after adjusting for the corresponding phenotype, evening/night shift work was nominally associated with education and smoking in fully adjusted models. Furthermore, two of the epigenetic age measures, GrimAge and PhenoAge, demonstrated higher age acceleration among evening/night shift workers.
The observational associations of evening/night shift work with lower education and higher BMI have been previously reported with comparable effect sizes (Daghlas et al. 2021; van Drongelen et al. 2011). While we did not find evidence of a phenotypic association with reported smoking, the association between evening/night shift work and the smoking methylation score is consistent with previous studies reporting that smoking behaviour is more common with shift work comparison to day work (Bae et al. 2017; Boini et al. 2022; Bushnell et al. 2010; Knutsson and Bøggild 2010; Knutsson and Nilsson 1998; Knutsson et al. 1988). One study also found that shift workers were also more likely to start smoking in comparison to their counterpart day workers (Cho et al. 2019; van Amelsvoort et al. 2006). We also found little evidence of an association between evening/night shift work and either self-reported or methylation measures of alcohol; however, previous research has found that shift workers were more likely to drink heavily (Dorrian et al. 2017).
There is growing evidence to suggest that DNAm-based measures are useful for health and lifestyle profiling (Hillary and Marioni 2021). Furthermore, several studies have shown that methylation predictors can provide a more accurate measurement of exposure than those based on self-report (Pośpiech et al. 2024; Zhang et al. 2016). For example, previous studies have shown that smoking methylation scores may provide a more accurate measure of true exposure compared with self-reported smoking (Lu et al. 2019; Zhang et al. 2016), possibly due to erroneous self-reporting, the broad categories for reporting exposure, or because DNAm is able to capture long-term biological changes as a result of smoking as well as secondary smoking. This is supported in part by our findings that the education and smoking scores were still weakly associated with night shift work even after adjusting for the corresponding self-reported exposures.
With respect to evening/night shift work, a number of studies have reported associations with individual CpG sites through epigenome-wide association analysis (EWAS) (Wackers et al. 2023), but few have looked at the relationship between evening/night shift work and DNAm profiles which capture markers of ageing and lifestyle exposure. Some small-scale studies have also reported epigenetic age acceleration in evening/night shift workers (n < 200) (Carugno et al. 2021; Wackers et al. 2023). White et al. (2019) investigated associations between epigenetic ageing in >2000 participants from the Sister Study (females only) and found that the length of time working night shifts was associated with increased PhenoAge, although GrimAge was not investigated.
Circadian oscillators have been found to contribute to epigenetic ageing (Koncevičius et al. 2024; Oh et al. 2018) and there is emerging evidence that DNAm age estimators relate to circadian rhythm (Horvath and Raj 2018). However, is should be noted that we did not find that evening/night shift work was associated with the Hannum and Horvath clocks, which were designed to estimate chronological age (Hannum et al. 2013; Horvath 2013). This absence of association is in accordance with findings from White et al. (2019) and suggests that intrinsic circadian processes are unlikely to underlie the associations observed (Levine et al. 2018; Lu et al. 2019).
Strengths and Limitations
One of the main strengths of the study is the use of GS, an epidemiological cohort study with a large sample size with DNAm data which has also captured information on evening/night shift work. Meta-analysing associations with those from the UKHLS dataset also improved power to detect associations between evening/night shift work, lifestyle factors and DNAm predictors and indicated consistency in associations across studies. Furthermore, we have the use of both self-reported and DNAm markers for the same exposure, so were able to directly compare.
There were some limitations to the study. The DNAm scores used in this study were developed in GS Set 1, which might lead to overfitting of the models. However, by using the second set of participants in GS and the independent UKHLS datasets, we hope that this should minimise any potential overfitting issues. There was also limited evidence for heterogeneity in the associations observed between GS set 1 and the other studies.
We specifically assessed evening/night shift work rather than other forms of shift work (such as morning or rotating shift work) in relation to DNA methylation. This is because of the previous evidence suggesting that work outside of normal daytime hours is likely to be particularly disruptive to biological processes and to have implications for adverse health (Kecklund and Axelsson 2016).
We acknowledge that the definition we have used to classify shift workers in our study differs from that used by the International Agency for Research on Cancer (IARC) Monographs to classify night shift work as “work occurring during the regular sleeping hours of the general population” (Ward et al. 2019). However, we were unable to separate evening from night shift work in the Generation Scotland study based on the question asked and so to best harmonise the data sources, we opted to combined those individuals working evening or at night in UKHLS as well. While this may have obscured some associations which might be more evident among pure night shift workers (e.g. those working between 11 pm and 7am), our definition captures a more diverse group of workers, which might lead to more generalizable findings across industries with varying shift schedules.
Another limitation was the lack of data on the intensity of evening/night shift work, e.g. whether the participant works 7pm-7 am three days a week or 7pm-10 pm every day which may have a differential biological impact (Stevens et al. 2011). Whilst rotating shift work (shifts which rotate or change according to a set schedule) has also been linked with circadian disruption, we did not specifically investigate this group of workers.
We are unable to make conclusions regarding causality of the associations observed. We cannot exclude reverse causation since the analysis in GS and that based on current evening/night shift work in UKHLS was assessed cross-sectionally. Given education preceeds shift work status, the association implies that those with less education engage in occupations that include night shift work. This is in line with a previous study which found an inverse association between a genetic (rather than epigenetic) risk score for higher education and shift work participation (Daghlas et al. 2021).
To support our findings, similar analysis should be performed in larger cohorts with DNAm data. Future studies should also look at evaluating whether these biomarkers could provide insights into the potential effects of night shift work on subsequent health outcomes, e.g. cardiometabolic diseases and cancer.
Conclusions
In over 8,000 participants from two cohort studies, evening/night shift work was associated with both phenotypic and DNA methylation-based measures of higher BMI and lower education. DNAm predictors of smoking and ageing were also related to evening/night shift work. Epigenetic measures may provide insights into the health and lifestyle profiles of night shift workers.