Abstract
Background Shift work is associated with increased cardiometabolic disease risk, but whether this association is influenced by cardiometabolic risk factors driving selection into shift work is currently unclear. We addressed this question using Mendelian randomization (MR) in the UK Biobank.
Methods We created genetic risk scores (GRS) associating with nine cardiometabolic risk factors (including education, body mass index [BMI], smoking, and alcohol consumption), and tested associations of each GRS with self-reported current frequency of shift work and night shift work amongst employed UKB participants of European ancestry (n=190,573). We used summary-level MR sensitivity analyses and multivariable MR to probe robustness of the identified effects, and tested whether effects were mediated through sleep timing preference.
Results Genetically instrumented lower educational attainment and higher body mass index increased odds of reporting frequent shift work (odds ratio [OR] per 3.6 years [1-SD] decrease in educational attainment=2.40, 95% confidence interval [CI]=2.22-2.59, p=4.84 × 10−20; OR per 4.7kg/m2 [1-SD] increase in BMI=1.30, 95%CI=1.14-1.47, p=5.85 × 10−05). Results were unchanged in sensitivity analyses allowing for different assumptions regarding horizontal pleiotropy, and the effects of education and BMI were independent in multivariable MR. No causal effects were evident for the remaining factors, nor for any exposures on selection out of shift work. Sleep timing preference did not mediate any causal effects.
Conclusions Educational attainment and BMI may influence selection into shift work, which may have implications for epidemiologic associations of shift work with cardiometabolic disease.
Key messages
Although it has been hypothesized that cardiometabolic risk factors and diseases may influence selection into shift work, little evidence for such an effect is currently available.
Using Mendelian randomization, we assessed whether cardiometabolic risk factors and diseases influenced selection into or out of shift work in the UK Biobank.
Our results were consistent with a causal effect of both higher BMI and lower educational attainment on selection into current shift work, with stronger effects seen for shift work that is more frequent and includes more night shifts.
Using multivariable Mendelian randomization, we found that effects of higher BMI and lower education were independent. Sleep timing preference had a null effect on shift work selection and therefore did not mediate these effects.
Selection through education and BMI may bias the relationship of shift work with cardiometabolic disease. Social mechanisms underlying these effects warrant further investigation.
Introduction
Shift work is increasingly common and is a risk factor for cardiometabolic diseases1 including type 2 diabetes (T2D)2 and coronary artery disease (CAD)3. In these studies the baseline characteristics of shift workers often systematically differ from non-shift workers, including: higher rates of smoking2,4,5, lower rates of alcohol consumption2,4, lower educational attainment6, and greater body mass index (BMI)6–8. Such baseline differences may reflect either a selection effect9 of these factors on shift work, a causal effect of shift work on these factors, or confounding by shared influences. To date, few studies have examined evidence for a potential selection effect10,11.
Clarifying which factors influence participation in shift work is of public health interest. First, identifying factors that influence selection into shift work would offer basic insights into which populations are at highest risk for working in unfavorable conditions12. Second, shift work may mediate the effects of cardiometabolic risk factors, such as lower educational attainment13, on adverse health outcomes. Alternatively, factors influencing selection into shift work may confound observational relationships between shift work and adverse health outcomes. Estimating the magnitude of this selection effect could facilitate development of methodology to better account for this confounding and improve our understanding of the adverse effects of shift work on health.
Observational data have limitations for identifying selection factors, due to timing of exposure measurement, unmeasured confounding, and effects of shift work on the risk factor (i.e. reverse causation). One method to overcome these limitations is Mendelian randomization (MR)14. MR uses genetic variants as proxies for epidemiologic exposures to estimate causal effects. This approach is well-suited for identifying selection effects on shift work because genetic variants are precisely measured and less prone to bias through confounding or reverse causality15,16.
Previous MR analyses supported a causal effect of BMI on lower income and socioeconomic status12,17, but MR has not been applied to characterize the determinants of selection into shift work. We used MR to assess the causal effect of cardiometabolic risk factors and diseases on shiftwork selection in the UK Biobank18.
Methods
Population
The UK Biobank is a population-based cohort study that enrolled over 500,000 volunteers aged 40-69 from 2006-2010; participation rate: 5.5%18. Data collection included questionnaire and nurse interview information, anthropometric and physiological measurements, and genomic data as previously described18,19. All participants provided written informed consent, and data used in this study were de-identified.
Exposures
The exposures were single genetic variants or weighted genetic risk scores (GRS) derived from genome-wide association studies (GWAS) for cardiometabolic risk factors and diseases: alcohol consumption, smoking heaviness, body mass index (BMI), waist-to-hip ratio adjusted for BMI, type 2 diabetes, coronary artery disease, and educational attainment (Table 1)20,21,22,23,24,13,25,26,27. We focused on risk factors with publicly available GWAS summary statistics, and with putative causal effects on cardiometabolic outcomes, as supported by prior MR studies31,20,28–31. A single nucleotide polymorphism (SNP) in the alcohol dehydrogenase 1B (ADH1B) gene robustly associated with reduced alcohol consumption20 was used to proxy alcohol consumption. Consistent with previous analyses20, we coded this missense variant (rs1229984) under a dominant model (cases combining homozygotes and heterozygotes), with the effect allele oriented to reduced alcohol consumption. A single missense variant (rs16969968) in the Cholinergic Receptor Nicotinic Alpha 5 Subunit (CHRNA5) gene, associated at genome-wide significance with one additional cigarette smoked per day21, was used as the genetic instrument for smoking heaviness (limiting the sample to ever smokers32). For all other exposures, we used SNP associations from GWAS meta-analyses that did not include UKB data to generate multi-SNP GRS22,23,25,27. Figure 1 shows a causal diagram by which these exposures may influence shiftwork selection.
We utilized risk profiling in PLINK33 v1.9 to construct beta-weighted GRS using individual-level data in the UKB for participants of White British ancestry; weights were obtained from the respective GWAS. We regressed each GRS on the respective exposure to determine the strength of each genetic instrument, using an F-statistic >10 to indicate minimal influence of weak instrument bias34.
Outcomes
We derived the study outcomes using current shift work characteristics at baseline from currently employed UKB participants. Participants were asked to indicate shift work participation by answering the following question: ‘Does your work involve shift work?’ defined as ‘a work schedule that falls outside of the normal daytime working hours of 9am-5pm. This may involve working afternoons, evenings or nights or rotating through these kinds of shifts.’ Response options included: ‘never/rarely,’ ‘sometimes,’ ‘usually,’ ‘always,’ ‘prefer not to answer,’ and ‘do not know.’ ‘Never/rarely’ served as the reference group, ‘sometimes’ was an intermediate group, and ‘usually’ and ‘always’ were collapsed into ‘frequent shift work;’ totaling three groups in this analysis. ‘Prefer not to answer’ and ‘do not know’ were coded as missing.
Participants indicating that they currently worked shifts were further queried, ‘Does your work involve night shifts,’ defined as ‘a work schedule that involves working through the normal sleeping hours, for instance working through the hours from 12am to 6am.’ Response options were identical to those for the primary shift work question and coded analogously. In the night shift work analyses we compared employed individuals not participating in shift work (reference group) with individuals working ‘shift work, but no night shift work,’ ‘some night shift work,’ and ‘frequent night shift work;’ resulting in four groups in this analysis.
All UKB participants with email addresses (n∼330,000) were invited to complete an online follow-up questionnaire on lifetime employment information. This questionnaire queried information on each job worked over the participant lifetime (n=118,699). For each job, participants were asked ‘Did you ever work shifts (day and/or night shifts) for this job?’. To determine whether associations generalized beyond current shiftwork participation, we used this variable to derive an outcome of lifetime history of shiftwork, regardless of current employment status (n=83,133; 22,987 ever shift workers / 60,146 never shift workers). To identify exposures influencing selection out of shift work, we tested the outcome of history of quitting shift work. Cases were defined as participants reporting working a shiftwork job and subsequently working a non-shift work job (n=22,987; 14,584 cases of participants who ever quit shift work / 8,403 controls who worked shifts and never quit) as shown in Figure 2.
Individual-level analyses
We first tested the association of each GRS35 with categories of the shiftwork outcomes using multinomial or binomial logistic regression adjusted for age, sex, and the top ten principal components (PCs) of ancestry. We transformed the effects of binary exposures (GRS for CAD and T2DM) to reflect a doubling in the odds of the exposure on the odds of the outcome36 (i.e. the effect of lifelong increased risk for disease on shift work selection).
We then probed the robustness of results through a series of secondary and sensitivity analyses. In order to include individuals who previously worked a job with shifts, but would no longer be included in the “current shift work” analysis, we also tested associations with lifetime history of shiftwork participation. To identify exposures influencing selection out of shift work, we tested associations with history of leaving shiftwork. To assess the impact of population stratification37, we controlled for i) 40 PCs of ancestry, and ii) the 22 UK Biobank assessment centers. Prior work suggested sex differences in the causal effects of BMI on socioeconomic outcomes17, so we tested for GRS interactions with sex using a model with an interaction term. We used UKB job codes to stratify jobs as ‘skilled’ and ‘unskilled’, using the classification laid out in Howe et al12, and tested for interactions of job skill with the GRS using a model with an interaction term.
Finally, we expanded the alcohol consumption and smoking heaviness instruments to include SNPs identified in a recent large-scale meta-analysis38. We used GRS weights from summary statistics excluding UKB38. These instruments explain more variance in the exposures relative to the single SNP instruments, but may be at greater risk for pleiotropy due to the inclusion of multiple variants with uncharacterized function39.
Summary-level MR analyses and sensitivity analyses
Genetic associations may be explained by horizontal pleiotropy, which is when variants influence the outcome through paths independent of the exposure of interest. We therefore conducted summary-level MR analyses to further assess whether the genetic associations represented casual relationships. We first estimated causal effects using inverse-variance weighted (IVW) random-effects regression, which is robust to balanced pleiotropy40. As a global test for pleiotropy, we calculated Cochran’s Q for heterogeneity. We then used the following pleiotropy-robust sensitivity analyses: MR Egger41, weighted median42, and MR-PRESSO43 regressions (see Supplementary methods 1). We used the estimate of the MR Egger model intercept to further assess for balanced horizontal pleiotropy. Due to low power44, this method is generally more appropriate for ruling pleiotropy in rather than ruling it out.
These analyses were implemented using the TwoSampleMR R package v0.4.2245. For effects identified in univariable MR, we undertook summary-level multivariable MR (MVMR)45 simultaneously controlling for each exposure using the TwoSampleMR software45. We confirmed instrument strength using the Qstrength adapted for MVMR46. The Qstrength divided by the number of variants is analogous to the F statistic, with values >10 indicating minimal weak instrument bias. Prior MR analyses identified a nominal effect of greater BMI on earlier sleep timing preference so we tested for a causal effect of early sleep timing preference on selection into shift work using GWAS estimates that did not overlap with UKB47. We hypothesized that earlier sleep timing preference would be negatively associated with current shift work, reflecting a negative selection effect. We planned to test MVMR of the exposures and sleep timing preference only if the first-stage association of sleep timing preference with shift work demonstrated an effect.
Results
Individual level analyses
The median age of employed participants in the UKB was 53 years [IQR 47-59], and 49% were male. Shift workers were more likely to be male and report lower indicators of socioeconomic and health status (Table 2). Each GRS was strongly associated with the exposures (Supplementary Table 1).
Genetically predicted educational attainment increased odds of working some shift work (odds ratio per 1-standard deviation (SD) / 3.6y (OR) 1.32, 95% confidence interval (CI) 1.10-1.59, p=2.64×10−03), frequent shift work (OR=2.22, 95%CI=1.89-2.63, p<2 × 10−16), and frequent night shift work (OR 2.94, 95%CI=2.27-3.85, p< 2×10−16; Figure 4). (Figures 3-4). Genetically predicted BMI increased odds of working ‘some’ shift work (OR 1.21, 95% CI 1.08-1.35, p=9.20 × 10−04), ‘frequent’ shift work (OR 1.25, 95% CI 1.14-1.38, p=7.40 × 10−06), and frequent night shift work (OR 1.44, 95% CI 1.23-1.68, p=4.00 × 10−06; Figures 3-4).
Both educational attainment (OR 1.40, 95% CI 1.20-1.64, p=2.68 × 10−05) and BMI (OR 1.17, 95% CI 1.06-1.28, p=1.68 × 10−03) increased the odds of working at least one shiftwork job across the life course (Supplementary Table 2). No other exposures associated with shift work participation (Figures 3-4). No exposures were associated with history of leaving shift work (Supplementary Table 3).
Secondary and sensitivity analyses in individual level analyses
The BMI and education GRS associations with shift work did not interact with sex or job skill classification. Additional adjustment for population stratification through controlling for 40 principal components and UKB assessment center did not influence the BMI or education effects (Supplementary Table 4). Results were similar when using an expanded GRS for alcohol consumption (OR=1.07, 95%CI=0.88-1.32, p=0.50) and smoking heaviness (OR=1.03, 95%CI=1.00-1.05, p=0.02) on the outcome of frequent shift work. There was no effect of genetic propensity to early sleep timing preference on frequent shift work (OR 0.99, 95% CI 0.96-1.04, p=0.79) or frequent night shift work (OR 0.96, 95% CI 0.90-1.02, p=0.16).
Summary-level MR analyses
We conducted summary-level MR analyses to estimate instrumental causal effects robust to horizontal pleiotropy. IVW random-effects models showed consistent effects of higher BMI (OR 1.30, 95% CI 1.14-1.47, p=5.58×10−05) and lower educational attainment (OR 2.40, 95% CI 2.22-2.59, p=4.84×10−20) on shift work (Table 3). Four lines of evidence suggest that the observed effects were not explained by horizontal pleiotropy. First, there was minimal evidence of heterogeneity (Table 3). Second, similar MR effects were observed across the Egger and weighted median analyses (Table 3). Third, there was no evidence of unbalanced pleiotropy as estimated by the Egger intercepts for the effects of BMI (intercept=-0.004, p=0.37) and of education (intercept=0.015, p=0.07) on shift work, although this test is generally underpowered. Fourth, MR-PRESSO did not detect pleiotropic outlier SNPs.
We next turned to multivariable MR to assess the independence of exposure effects on shift work. The genetic instruments for educational attainment (Qstrength=2,407, Fadjusted=15.4) and BMI (Qstrength=4,393, Fadjusted=28) were conditionally strong in multivariable MR. The effects of BMI (OR=1.18, 95%CI=1.08-1.31, p=1.94×10−03) and of education (OR=2.56, 95%CI=2.12-3.13, p=1.48×10−23) on self-reported ‘frequent’ shift work were independent in MVMR.
Discussion
These analyses revealed a causal and independent effect of lower educational attainment and higher BMI on current and lifetime selection into shift work, particularly night shift work. These effects were not mediated by sleep timing preference. There was minimal evidence for selection through other tested exposures.
Our findings provide evidence for a causal effect of lower educational attainment and higher BMI, but not of other cardiometabolic risk factors, on selection into shift work within this population. The effect sizes were larger when shift work involved night shifts, and with increasing frequency of shift work, supporting a dose-dependent effect. Selection into shift work (including alcohol consumption, BMI, smoking, lipids, and glucose) by cardiometabolic risk factors has been evaluated in previous observational studies (N∼2,800) 10,11. These studies identified an association of smoking with future participation in shift work, but not of other cardiometabolic traits. This contrasts with our results, where we did not identify strong evidence for an effect of smoking heaviness on shift work selection. Given the inverse effect of educational attainment on shift work, and the known effect of educational attainment on reduced smoking30, it is possible that educational attainment confounded the previously observed relationship of smoking with shift work. Although shift workers typically report lower educational attainment6,48, prior work has not evaluated the role of educational attainment in shift work selection. Our findings suggest that this imbalance is driven by a causal effect of lower educational attainment on selection into shift work.
We observed a robust effect of higher BMI on shift work selection. This is in contrast to null associations in prior observational work10,11, There are several potential explanations for these differences. First, the observational estimates10,11 had wide confidence intervals, and could not exclude moderate effects of BMI on selection. Second, the influence of BMI on selection into shift work may accrue over time, such that effects are only observed later in life. Such a phenomenon would affect the present analysis, given the use of an older sample in contrast to the younger samples studied in prior studies of selection into shift work. However, this finding is consistent with the growing evidence for a causal effect of higher BMI on lower socioeconomic status12,17.
There are several potential pathways by which education and BMI may influence selection into shift work. Lower educational attainment is related to lower-skilled employment49, which is generally more typical of shift work50. It is less clear what mediates the effect of BMI on selection into shift work. One potential mechanism is weight stigma, which is a highly prevalent51, well-characterized driver of employment inequities52. Women experience a greater burden of this discrimination than men53, but we observed no interaction by sex, suggesting that men and women experience similar selection pressures into shift work on the basis of differences in adiposity and educational attainment. Although we hypothesized that differences in sleep timing preference due to variation in BMI47 or education may be a mediating mechanism, we did not observe an effect of earlier sleep timing preference on night shift odds. This is surprising given that night shift workers who tend to sleep earlier may experience the greatest burden of circadian misalignment54,55, and would therefore be more likely to select out of night shift work. A possible interpretation of this result is that selection effects of other social factors, such as socioeconomic status or educational attainment, preclude leaving a job on the basis of tolerability. We did not identify factors influencing selection out of shift work, but this sample was an order of magnitude smaller than the main sample, and we may have been underpowered to observe modest effects.
The present results have several potential implications for shift work research. The associations of shift work with cardiometabolic outcomes may be influenced by confounding by BMI and educational attainment. This does not invalidate the contribution of circadian misalignment to the excess disease risk conferred by shift work, but confounding may distort these estimates. Cross-sectional studies of shift work using obesity as an outcome may be at high risk for bias, as any observed associations will reflect a sum of forward and reverse causal effects. 56 Analyses restricted to shift workers may also be affected by collider bias57, which can induce false positive associations. Alternatively, shift work may mediate the effects of education and BMI on health outcomes. Recent MR mediation analyses showed that only half the variance of the effect of education on coronary artery disease risk is explained by conventional risk factors30, and shift work is an independent risk factor for cardiovascular disease3. Exposure to shift work may therefore explain part of the variance in this relationship. More broadly, MR may be a useful tool in occupational epidemiology to determine whether reverse causality may drive epidemiologic associations.
Our study has some limitations. We studied shift work in a middle-older aged population, and our results may therefore not generalize to younger populations. Given that UKB is a relatively healthy population58, our results may not generalize to sicker populations and may not be applicable to other populations or cohorts59 with different occupation structures or ethnicities. Moreover, the identified effects may be driven by selection bias in UKB given a 5% responder rate in the cohort58. Factors influencing selection into UKB may induce collider bias57 and false positive associations of exposures with selection into shift work. While we did not examine mediation through job-related factors, such as tolerability or fatigue, we were able to examine sleep timing preference and it did not appear to mediate the observed effects. Finally, despite consistent results in sensitivity analyses, the MR effect estimates may be biased by horizontal pleiotropy due to effects of genetic instruments on pathways influencing the outcome independently of the studied exposures. Triangulating results with analyses leveraging other forms of natural experiments, such as educational reform in the UK60, will strengthen confidence in causality.
In conclusion, higher BMI and lower educational attainment influenced selection into shift work and night shift work, while no other cardiometabolic risk factors were associated with selection into or out of shift work.
Data Availability
The following data availability statement is made by the UK Biobank (https://www.ukbiobank.ac.uk/scientists-3/): “UK Biobank is an open access resource. The Resource is open to bona fide scientists, undertaking health-related research that is in the public good. Approved scientists from the UK and overseas and from academia, government, charity and commercial companies can use the Resource.” Summary statistics used to generate the genetic instruments are all publicly available.
Funding
R01 DK105072 (Saxena, Dashti, Vetter), R21OH011052 (Schernhammer, Vetter), R01 DK107859 (Dashti, Saxena), R01DK102696 (Saxena), MGH Research Scholar Fund (Saxena), Diabetes UK 17/0005700 (Rutter), The University of Manchester (Research Infrastructure Fund). RCR and GDS are members of the MRC Integrative Epidemiology Unit at the University of Bristol funded by the Medical Research Council (MM_UU_00011/1). RCR is a de Pass Vice Chancellor’s Research Fellow at the University of Bristol.
Disclosures
C.V., during the conduct of the study, received research support from the NIH, was a scientific advisory board member of Circadian Light Therapy Inc., and served as a paid consultant to the US Department of Energy outside the submitted work.
Acknowledgements
This research has been conducted using the UK Biobank Resource (UK Biobank application number 6818). We would like to thank the participants and researchers from the UK Biobank who contributed or collected data.