Article Text

Risk of leukaemia mortality from exposure to ionising radiation in US nuclear workers: a pooled case-control study
  1. Robert D Daniels,
  2. Stephen Bertke,
  3. Kathleen M Waters,
  4. Mary K Schubauer-Berigan
  1. National Institute for Occupational Safety and Health (NIOSH), Division of Surveillance, Hazard Evaluations, and Field Studies (DSHEFS), Industrywide Studies Branch (IWSB), Cincinnati, Ohio, USA
  1. Correspondence to Dr Robert D Daniels, National Institute for Occupational Safety and Health (NIOSH), Division of Surveillance, Hazard Evaluations, and Field Studies (DSHEFS), Industrywide Studies Branch (IWSB), 4676 Columbia Parkway, Mailstop R-13, Cincinnati, OH 45226, USA; rtd2{at}


Objective To follow-up on earlier studies of the leukaemogenicity of occupational ionising radiation exposure.

Methods We conducted a nested case-control analysis of leukaemia mortality in a pooled cohort of US nuclear workers followed through 2005. Each case was matched to four controls on attained age. Exposures were estimated from available records. General relative risk models were used to estimate the excess relative risk (ERR) of leukaemia, excluding chronic lymphocytic (CLL), acute myeloid leukaemia, chronic myeloid leukaemia and CLL while controlling for potential confounders. Preferred exposure lags and time-windows of risks were calculated using joint maximum likelihood. Dose-response was also examined using linear, linear-quadratic, categorical and restricted cubic spline models.

Results There were 369 leukaemia deaths in 105 245 US nuclear workers. The adjusted ERR for non-CLL leukaemia was 0.09 (95% CI −0.17 to 0.65) per 100 mGy. Elevated non-CLL risks were observed from exposures occurring 6–14 years prior to attained age of cases (ERR per 100 mGy=1.9; 95% CI <0 to 8.0). Lagged models indicated non-linearity of risk at very low (<10 mGy) and high (>100 mGy) doses, which contributed to the imprecision of results in linear models. Similar risk attenuation was not evident in time-windows-based models.

Conclusions Risk estimates were in reasonable agreement with previous estimates, with the temporality of non-CLL leukaemia risk as a dominant factor in dose-response analyses. Future research should focus on methods that improve evaluations of the dose-response, particularly in the low-dose range.

Statistics from

What this study adds

  • Questions remain concerning the translation of risks observed in Japanese atomic bomb survivors to working populations who are exposed to low doses of ionising radiation delivered at low rates of exposure.

  • With 369 leukaemia deaths identified in a pooled cohort of over 100 000 nuclear workers from six major US facilities, this study is among the largest leukaemia studies ever conducted with low-dose radiation. The large study size enabled analysis by leukaemia subtype and detailed examination of dose-response characteristics that have not been completed prior.

  • Positive but imprecise estimates of leukemia risk, which increased with dose, were observed for myeloid leukaemias but not for chronic lymphocytic leukaemia (CLL). These findings were consistent with previous studies that suggest a weak leukaemogenicity of ionising radiation for non-CLL subtypes.


Most information on the leukaemogenicity of ionising radiation stems from studies of radiation therapy patients and survivors of the atomic bomb blasts in 1945.1 ,2 In particular, the Life-Span Study (LSS) of the Japanese atomic bomb survivors forms the foundation for current worker protection standards. The LSS has the advantages of large size (∼87 000), all ages at exposure, both sexes, long follow-up period and doses over a range of interest (0–4 Gy with 60% less than 0.1 Gy), although questions remain concerning the translation of risks to workers who are intermittently exposed to low doses delivered at low dose rates. A preferred alternative is the direct estimation of risk in working populations; however, most occupational studies lack sufficient statistical power for examining rare outcomes such as leukaemia. As such, recent studies have estimated leukaemia risk in nuclear workers using pooled, or meta-analyses, that offer improved precision.3–7 In particular, Daniels and Schubauer–Berigan conducted5 a meta-analysis using information from occupational and environmental studies to estimate the excess relative risk (ERR) of leukaemia, excluding the chronic lymphocytic subtype (CLL). The summary estimate from 10 studies (ERR at 100 mGy=0.19; (95% CI 0.07 to 0.32) was in reasonable agreement with the LSS. Nevertheless, a meta-analytic approach relies on the quality of data available in the published literature; therefore, these analyses may be vulnerable to biases in the source information. Moreover, studies have typically relied on models that assume risk is linear or linear-quadratic in dose for a heterogeneous group of diseases comprising non-CLL leukaemia. Few studies have examined departures from linearity in the dose-response or risk differences by leukaemia subtype.

This study extends follow-up of workers previously studied by Schubauer–Berigan et al3 ,4 and expands the cohort to include workers from the Idaho National Laboratory.8 Previous analyses reported estimates of the ERR of mortality per unit dose for non-CLL leukaemia that reasonably agreed with the LSS4; however, there was evidence of risk attenuation at doses above 100 mGy and temporal patterns in risk that varied with time since exposure.4 The purpose of the current study is to: (1) improve upon previous estimates through additional follow-up and cohort expansion; (2) report on the risk by leukaemia subtype and (3) further evaluate the shape and temporal behaviour of the dose-response.


Study cohort

This research complied with the requirements of the Federal Policy for Protection of Human Subjects (10CFR745 or, where applicable, 45CFR46), and was reviewed by the National Institute for Occupational Safety and Health (NIOSH) Human Subjects Review Board to ensure that the rights and welfare of study subjects were protected.

The study population consisted of workers who were employed at the Hanford Site (Hanford) near Richland, Washington; the Savannah River Site (SRS) near Aiken, South Carolina; the Oak Ridge National Laboratory (ORNL) in Oak Ridge, Tennessee; the Los Alamos National Laboratory (LANL), including the Zia Company employees, in Los Alamos, New Mexico; the Idaho National Laboratory (INL) near Idaho Falls, Idaho; and the Portsmouth Naval Shipyard (PNS) in Kittery, Maine. Eligible workers were monitored for radiation exposure and employed at one or more of the primary facilities for at least 30 days. Hanford, LANL and ORNL workers were further restricted to those first employed on or after 1 January 1951 to reduce exposure errors believed associated with early dosimetry practices at Manhattan Engineering District Facilities.4 Vital status ascertainment was extended through 31 December 2005, adding 9–15 years of follow-up to the previous study.

Cases were decedents with leukaemia (ICD-9: 204–208; ICD-10: C91–C95, excluding leukaemic reticuloendotheliosis (C91.4)), identified as the underlying cause of death. Distinctions were made by major subtypes, namely: CLL, acute lymphoid leukaemia (ALL), chronic myeloid leukaemia (CML) and acute myeloid leukaemia (AML). Non-CLL cases comprised all-acute, and subacute leukaemias, and myeloid or monocytic leukaemias, excluding CLL-indeterminate cases that were defined as cases lacking sufficient information on cell type and severity to determine non-CLL status (ie, ICD-9 codes: 204.9, 208.1, 208.9; ICD-10 codes: C91.9, C95.9, C95.1). Similarly, AML-indeterminate cases were defined as acute leukaemias of unspecified cell type (ie, ICD-9 code: 208.0; ICD-10 code: C95.0). Four controls were selected for each case using incidence density sampling.9 Controls were matched to their respective case on attained age and by primary study group, whereby, Group I risk sets were formed from the previously studied group of Hanford, LANL/Zia, ORNL, PNS and SRS workers,3 ,4 and Group II risk sets contained the INL workers.8 This group matching was done to maximise utility of information collected on cases and controls selected for the previous study.


The primary exposure was low linear energy transfer (low-LET) ionising radiation from all occupational sources, including work-related medical x-ray examinations (WRX) conducted during routine physical examinations by the employer.10 ,11 Workplace exposures to neutrons and internally deposited plutonium were also considered as potential cofactors in dose-response analyses. External ionising radiation exposure measurements for each study participant were abstracted from site dosimetry databases and dose registries maintained by the US Department of Energy (DOE) and US Nuclear Regulatory Commission (NRC). All doses were quantified in terms of absorbed dose to haematopoietic tissue using methods described elsewhere.4 ,11–13 There were five exposure variables describing ionising radiation dose (ie, gamma, tritium, WRX, neutron and plutonium). Low-LET exposure was calculated by summing doses from whole-body gamma, tritium and WRX sources.

Benzene exposures were assessed using a job exposure matrix to assign task-specific benzene ‘scores’ that were based on combinations of qualitative ratings of exposure intensity and frequency over time. Potentially exposed workers were identified using available employment information on job assignment (eg, job title, division and department) and work location in each study facility (eg, area, building and room). Tasks were identified and characterised using site-monitoring records and process information. Task-based scores were calculated for each study subject and period of employment, which were then accumulated over the observation period.

Socioeconomic status

The methods used to assign socioeconomic status (SES) are described elsewhere.4 ,8 In general, the job title at first hire for each worker was related to an occupational code listed in the 1980 US Bureau of Census. The first job title was chosen because it is thought that initial employment provides a better marker for some lifestyle factors (eg, smoking) that are typically established in adolescence and early adulthood.4 ,8 ,14 Each of these codes was then assigned to one of six SES categories, namely (from high to low SES): professional, intermediate, skilled non-manual, skilled manual, partly skilled and unskilled.

Statistical methods

All statistical analyses were conducted using SAS software.15 Primary modelling used general relative risk models that represent a linear trend in radiation dose.16 Relative risk models followed the general formEmbedded Image

where α is the radiation risk coefficient; βi is the coefficient for the ith covariate of 1, 2 …n covariates indexed by multiplicative factors, Xi; and γj is the radiation weighting factor for the jth dose variable of three linear dose covariates, D. For all models, three radiation sources were considered: (1) low-LET (D1) irradiation; (2) plutonium deposition (D2); and (3) neutron radiation (D3). Setting the low-LET radiation weighting factor to unity (γ1=1) allowed solutions to radiation weights for plutonium and neutron doses that can be compared with standard values. Profile likelihood-based (PL-based) two-sided 95% CIs were calculated for model parameter estimates.

Effect modification was evaluated by likelihood ratio test (LRT) at the 0.05 critical level of the cross-product term of the potential effect modifier and radiation dose. Factors evaluated were: sex, race (White non-Hispanic/all others), birth cohort (tertiles), period of first hire (tertiles), facility of longest employment duration and study group (Group I vs Group II). Potential confounders were incorporated into models when causing >15% relative change in the overall dose parameter estimate. To minimise the potential bias from selection of cutpoints for categories of birth cohort and hire date, continuous variables were constructed from terms calculated by restricted cubic spline (RCS) models, with three knots at the 10th, 50th and 90th percentile of birth date and hire date, respectively.17 Because birth date and hire date were correlated (Pearson coefficient=0.59, p<0.001), models were restricted to include no more than one of these variables at a time. Final selection between birth date and hire date was based on a comparison of observed effect and the statistical significance of the covariates.

Baseline models were constructed for non-CLL leukaemia and three of the four major subtypes: AML, CML and CLL. There were too few cases to adequately examine ALL risks separately. The methods proposed by Richardson et al18 were used to simultaneously estimate the preferred exposure lag by joint maximum likelihood. A preferred time-window of exposure was similarly estimated.

The linearity of the dose-response was examined by: (1) evaluating the fit of a quadratic term added to the model; (2) examining the risk by dose categories; (3) fitting RCS models and (4) fitting piecewise linear spline (PLS) models. Lagged categorical models used categories of 0−<1.0, 1−<10, 10−<50, 50−<100, and ≥100 mGy, with the lowest category as referent. Time-windows-based models used dose categories with cutpoints defined by quintiles of the dose distribution of the cases for each outcome. Lagged RCS models set four knots at the 5th, 35th, 65th and 95th percentiles of the distribution of low-LET dose, while windows-based RCS models used three knots defined by the 10th, 50th and 90th percentiles of dose.17 PLS models were restricted to non-CLL (lagged and windows-based); both two- and three-piece solutions were evaluated. Model fits were judged by Akaike Information Criterion (AIC).19 To test improvement in fit over a linear model, LRTs were conducted comparing the adjusted linear-quadratic, RCS and PLS models with the adjusted linear model.

The effects of high-LET radiations were examined in sensitivity analyses using model solutions to radiation weights and by varying weighting schemes for biological effectiveness. Three sets of weights for low-LET, neutron and plutonium were considered: 1, 10 and 20, which is consistent with the protection quantity ‘equivalent dose’20; 1, 11 and 3.6, which are median radiation effectiveness factors suggested by Kocher et al21; and 1, 10 and 1, based on the upper bound suggested by the United States Environmental Protection Agency (US EPA) in 1999 for α emitters deposited on bone surfaces or in bone.22 The current EPA-recommended value of 2 (95% CI 1 to 3) for incorporated α-emitters is within the range analysed.23


Descriptive and univariate

There were 94 517 workers from Hanford, LANL, ORNL, PNS and SRS (ie, Group I) and 63 560 INL workers (ie, Group II). Applying the inclusion criteria, and then combining the two cohorts, resulted in a base cohort of 105 245. There were 30 266 deaths (28%), which included 8721 cancer deaths and 369 leukaemia cases. Among leukaemias, 291 cases were identifiable to one of four subtypes. Non-CLL leukaemia (n=264) comprised AML (n=150), CML (n=52), ALL (n=18), unspecified acute leukaemias (n=36) and myeloid leukaemias (n=2), and other non-CLL subtypes (n=6). There were 74 CLL deaths. (See Web Appendix, tables 1–3 for additional information on the study cohort and case distribution.)

Dose distributions were right-skewed; the average, median and maximum values of low-LET dose lagged 2 years were 26.5 and 7.2 mGy, and 634 mGy, respectively. There were 293 (15.9%) cases and controls exposed to tritium, 237 (12.8%) with neutron dose, 130 (7.0%) with plutonium dose and 360 (19.5%) who were assigned benzene scores (see Web Appendix, table 4 for additional information on exposures).

Univariable analysis of potential confounders revealed similar results across non-CLL leukaemia and subtypes. The most pronounced effect was observed across categories of birth year and year of first hire, where less risk was consistently observed in persons born earlier and/or hired earlier. Decreased risk was observed in women and non-Hispanic Whites compared with men and other races/ethnicities, respectively; however, among these factors, only CML in other races/ethnicities compared with non-Hispanic Whites was statistically significant (HR=4.1; 95% CI 1.3 to 12.7). More information on univariate analyses is available in the Web Appendix, table 5.

Multivariable analyses

Of variables tested for interaction with low-LET radiation, only race approached statistical significance (eg, non-CLL: Χ2=3.07, df=1, p=0.080). As observed previously,4 restricting analyses to Group I showed evidence of risk heterogeneity across facilities (Χ2=18.07, df=4, p=0.003). Positive dose-responses were observed in PNS workers (ERR per 100 mGy=3.7; 95% CI 0.76 to 20.0) and SRS (ERR per 100 mGy=2.5; 95% CI 0.44 to 8.3), which influenced the overall estimate for the group (ERR per 100 mGy=0.12; 95% CI −0.16 to 0.66). The estimate of non-CLL leukaemia risk for INL was also positive (ERR per 100 mGy=0.03; 95% CI −0.23 to 0.81), although markedly less than Group I. There was no evidence of significant risk differences between Group I and Group II (Χ2=0.086, df=1, p=0.769).

Including sex in the model reduced risk estimates for non-CLL leukaemia by over 27%. Adjusting for race alone increased the estimate by nearly 40%. Including hire year by tertiles resulted in a 40% increase in the estimate compared with just under a 15% increase when controlling for hire year using RCS. Adjustment for benzene, SES or birth cohort resulted in less than 15% relative change. Similar results were obtained for CML, where final models were adjusted for race, sex and hire year. The model of AML risk was not influenced by sex or race, but was increased by 28% when including RCS terms for birth cohort, and 17% when including hire date. Given the strong association between birth cohort and AML compared with hire date and AML, subsequent models were controlled for birth date using spline terms. CLL risk estimates were influenced by sex (fourfold decrease), race (twofold increase), and birth cohort (10-fold increase using RCS).

Solving for radiation weights resulted in relatively large values for both high-LET components and a slight reduction in the risk per unit dose (table 1). Fixing weights in sensitivity analyses moderately increased ERR estimates of non-CLL leukaemia per unit dose. These patterns were consistent in lagged models; however, weight solutions were substantially reduced in windows-based models. For example, non-CLL risk from exposures 6–14 years before the attained age of the case had solutions for plutonium and neutron weights of 0 (fixed bound) and 16, respectively, and the relative change in the ERR estimate was about 3%. Given evidence of only moderate high-LET effects and small percentages of exposed workers, we excluded plutonium and neutron doses from final models.

Table 1

Sensitivity analysis of non-chronic lymphocytic leukaemia risk from high and low linear energy transfer radiations combined (lagged 2 years, fixed), adjusted for sex, race and hire year

In general, positive, but non-statistically significant estimates of excess radiation-associated risk were observed in all lagged models except for CLL (table 2). There was no evidence of CLL radiogenicity in any model tested. The overall estimate of non-CLL risk was 0.09 (95% CI −0.17 to 0.65) per 100 mGy (lagged 2 year), excluding CLL-indeterminate cases, and adjusting for sex, race and hire year. Including CLL-indeterminate cases in the model reduced the observed non-CLL risk (ERR per 100 mGy=0.03; 95% CI −1.5 to 0.38). Among leukaemia subtypes, the highest risk in lagged models was observed in CML (ERR per 100 mGy=0.29; 95% CI <0 to 1.8; 6.3 year lag). Temporal risk patterns were clearly evident in time-windows-based models, and the most pronounced effect was observed for AML in the 6–14-years window (ERR per 100 mGy=7.0; 95% CI 0.079 to 32.0), although CIs were wide. AML cases in the 6–14-years window had mean and maximum cumulative doses of 3.8 and 71.2 mSv, respectively, while control mean and maximum doses were 2.1 and 137.4 mSv, respectively. Results from alternate analyses of exposure lags and time since exposure using methods comparable with the previous study are reported in tables 6–8 of the Web Appendix.

Table 2

Linear dose*-response modelling results by preferred lag and exposure window

Dose-response shape analyses

The results of fitting linear, linear-quadratic, categorical and RCS dose-response models of low-LET dose and non-CLL leukaemia are shown in figure 1. Figures for AML and CML dose-response modelling are available as figures 1 and 2 in the Web Appendix. Panel A shows the results of lagged models while Panel B shows results using the best-fit exposure window. Model convergence was not achieved for some time-windows-based models (ie, AML and CML linear-quadratic forms and CML RCS, see Web Appendix). The best fitting non-CLL and AML-lagged models were obtained from RCS (dashed lines) based on AIC and the LRT (table 3). The shape of the dose-response for AML was similar, whereby the lagged RCS and categorical plots show risk attenuation at cumulative doses below 10 mGy, followed by a sharp increase to 100 mGy before decreasing again. CML dose-response under a 5-year lag was best-fit to a linear model; however, inspection of RCS and categorical results suggested some attenuation at low and high doses, albeit diminished relative to AML. By contrast, there is little evidence of a departure from linearity in low-dose regions of windows-based analyses of all outcomes, although non-linear models suggested some attenuation at higher doses. Alternate analyses examining the non-linearity of the dose-response at high doses are reported in the table 9 of the Web Appendix.

Table 3

Model fit characteristics by outcome, exposure (preferred lag or window), and model type

Figure 1

Mortality from leukaemia, excluding the chronic lymphocytic subtype, and cumulative absorbed dose to bone marrow. (A) Two-year lag; (B) 6–14-year window. Solid line: linear; long dashed line: linear-quadratic; short dashed line: restricted cubic spline (RCS); filled squares: category midpoints (with 95% CI shown by error bars).

Two knot solutions (8.3 and 46.9 mGy) in the non-CLL PLS-lagged model indicated the points of departure from linearity in the low- and high-dose ranges. The slope estimates (ERR per 100 mGy) in each segment were: −6.6 (95% CI −9.8 to −0.92) for doses ≤8.3 mGy; 2.0 (95% CI 0.30 to 4.2), for doses >8.3 and <46.9 mGy; and −0.16 (95% CI: <−0.22 to 0.39) for doses ≥46.9 mGy. By contrast, the best fit of time-windows-based data was achieved in models using two segments with the knot at 20.1 mGy; however, the two segments were not statistically significantly different (p=0.86). The first segment was positive and statistically significant (ERR per 100 mGy=6.4; 95% CI 0.31 to 17.0). The second segment was negative and lacked statistical significance (ERR per 100 mGy=−0.82; 95% CI −1.7 to 4.4).


With 369 leukaemia cases in a population of 105 245 US nuclear workers, the current study is among the largest examining the leukaemogenicity of protracted ionising radiation exposures. The large size improved power compared with the previous analyses,3 ,4 and enabled analyses by leukaemia subtype. Estimates of risk per unit dose were consistent with results from previous studies (table 4), as were temporal patterns observed in populations exposed to known leukaemogens, such as benzene,24 ,25 ionising radiation4 ,26–28 and smoking,29 whereby leukaemia risks lessen with time since exposure. Nevertheless, statistical power remains limited, leading to results which, in many cases, could be explained by the role of chance.

Table 4

Leukaemia ERR at 100 mGy exposure by subtype and by study

A non-linear dose-response was most evident in lagged models, whereby risk attenuation was observed at very low doses (<10 mGy) and at doses greater than 100 mGy. By contrast, strong departure from linearity was not evident in time-windows-based models, where risks appear greater than in lagged models. The elevated risk in windows-based models suggested temporal patterns in dose-response with time since exposure that were similar to, but more pronounced than, that observed in Japanese atomic bomb survivors26 and in other nuclear worker cohorts.4 ,28 AML risk appears particularly susceptible to temporal effects, where a 30-fold change in risk was observed in the windows-based model compared with a lagged model. Without accounting for these temporal effects in studies of protracted exposures with long follow-up since periods of highest exposure, dose-response analyses may provide distorted results, as evident in the lagged models.

We examined the combined effects from low- and high-LET radiations and determined that there was little evidence of high-LET effects in this study. Solutions to weights in time-windows-based models appeared in fair agreement with expectation.22 ,32 ,33 By contrast, weights in lagged models were larger than expected, although CIs were wide. Exposure misclassification is a possible cause of spurious results given the sparseness and poor quality of available plutonium and neutron exposure measurements. It is also conceivable that other characteristics unique to the few plutonium workers or persons with measurable neutron dose may have led to inflated weight estimates. Last, the cohort was assembled from nuclear facilities selected to minimise high-LET dose contributions4; therefore, only a small number of study participants were significantly exposed to plutonium or neutrons.

Risk attenuation in the low-dose range was not evident in previous studies of workers employed at these facilities.4 ,8 The addition of low-exposed cases through extended follow-up may be a source of the non-linear dose-response. Restricting the current study to non-CLL cases with follow-up through 1990 (n=112 cases) reduced the non-linear low-dose effect and the best fit was obtained from a linear model (ERR per 100 mGy=0.6; 95% CI −0.06 to 2.0). Diminishing leukaemia risk with continued follow-up has been observed in populations exposed to benzene24 ,34 and ionising radiation.35 ,36 Rinsky et al34 reported a decline in the leukaemia standardised mortality ratio from 3.37 (95% CI 1.54 to 6.41) to 2.56 (95% CI 1.43 to 4.22) after adding 15 years of follow-up to a cohort of white male benzene-exposed rubber hydrochloride workers. Silver et al24 also examined these workers using yearly follow-up information and found a fivefold decrease in the risk estimates from 1961 to 1996. In studies of ionising radiation, Boice et al35 reported decreased non-CLL leukaemia mortality after adding 9 years of follow-up to the study of Rocketdyne radiation workers,37 as did Krestinina et al36 after adding 6 years of follow-up to the Techa River cohort.38 Exposure profiles for these populations are temporally similar to the current study, where exposures occurred earlier in the follow-up period, and continued follow-up resulted in negligible additional dose. Given that leukaemia risk diminishes with time since exposure, it is less likely that new cases, who have not accrued additional dose in the extended follow-up period, are a consequence of ionising radiation exposure.

Risk attenuation at higher exposures is often observed in occupational studies, and may result from a number of causes, including a depletion of susceptible persons, biologic saturation, exposure misclassification or a healthy worker survivor bias.39 A healthy worker effect is less likely because leukaemia is not strongly related to lifestyle factors. Biologic saturation is also unlikely given that ‘high doses’ in this study are well below known thresholds for deterministic effects, and below doses that appear within the linear range in higher-dose studies. Exposure misclassification, although inevitable, is believed to be reduced by the availability of individual measurement data and the exclusion of early workers who were prone to highly uncertain dose estimates. Nevertheless, this source of uncertainty cannot be completely ruled out as cause for the observed attenuation. Some suggest that variations in risk in an exposed group may be attributable to small numbers of individuals who are predisposed to radiation-induced leukaemia.40 ,41 In chronically exposed populations, those most susceptible to the disease present early after relevant exposure, and persons who are resistant remain leukaemia-free and continue to accumulate dose, which may explain the observed risk attenuation at higher doses and the dilution of risk with increased follow-up.

Of course, results from low-dose studies are highly susceptible to distortion by confounding and exposure measurement errors. Benzene, a known leukaemogen and a confounder in the previous study, was not associated with leukaemia risk in this study. A possible explanation is the exclusion of early workers who were at greatest risk of benzene exposure in the previous study. SES was used as a surrogate for lifestyle factors; therefore, a potential for incomplete control of smoking effects cannot be ruled out, although smoking is only weakly related to leukaemia. Regarding exposure, extensive searches for exposure records were conducted to estimate career doses; nevertheless, some exposure misclassification from incomplete or erroneous exposure data is likely. For example, we found that cases with exposures of 1.0 mGy, or less, were typically professionals who were hired later, were older at hire, and worked fewer years prior to termination compared with cases in the higher-dose categories. These patterns suggested a potential for previous dose during employment elsewhere and the lack of this information in dose-response analyses may have led to biased estimates. To examine this potential, non-CLL cases in the lowest-dose category (n=74) were matched to exposure data maintained on claimants under the Energy Employee Occupational Illness Compensation Programme Act. There were eight claimants identified, and four had additional exposure information in claimant files. Had these data been available for our study, a change in exposure category would have resulted in every case.

Although we have offered several possible explanations, the cause for the non-linear dose-response remains unknown. In the absence of known biases, one could interpret the attenuated risk at low doses as evidence of a threshold or protective effect. There has been considerable debate regarding the validity of the linear-no-threshold (LNT) dose-response model because the carcinogenicity of low- and very low-dose radiation exposure is largely uncertain.42–48 Although threshold or protective effects at low doses cannot be ruled out by our findings, we did not observe attenuated risk in time-windows-based models that very likely represent the most relevant period of risk following exposure.


Positive, but imprecise estimates of excess relative risk were observed for all leukaemia outcomes except the CLL subtype. Non-CLL leukaemia risk was time-dependent, where peak excess risk per unit dose was observed at select times prior to disease onset. Marked improvements in estimates were observed in time-windows analyses, where linear dose-response models prevailed over other forms tested. Overall, the point estimates and temporal patterns observed in this study were in reasonable agreement with previous results; however, some dilution from extended follow-up without concomitant increases in dose was apparent.

Risk models using exposure lags exhibited attenuated risk in the low-dose (<10 mGy) and high-dose (>100 mGy) regions. The low-dose attenuation disappeared in models based on dose accrued in time-windows associated with maximum leukaemia risk prior to attained age of the case. Risk attenuation in the low-dose region is particularly difficult to explain given that similar results were not observed in previous studies; however, there is some evidence that extended follow-up and exposure misclassification in the low-dose range may be partially at fault. Low-dose effects, which are most important for risk assessment, are much more difficult to discern compared with the effects at high doses given limited statistical power and results that are highly susceptible to bias. Future research should focus on methods that improve evaluations of the dose-response, particularly in the low-dose range. In particular, better collection of career dose information may greatly reduce the potential for exposure misclassification. For example, future studies could be markedly improved by a consolidation of exposure information from NRC, DOE and military sources into a national dose registry.


The study also served to partially satisfy the requirements for doctoral dissertation research by the corresponding author, who is grateful to the dissertation committee: Dr S Pinney, Dr M Schubauer-Berigan, Dr C R Buncher, Dr R Hornung, and Dr H Spitz, for their encouragement and valuable input.


Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Files in this Data Supplement:


  • Contributors All coauthors had full access to the data used in the study, and take full responsibility for the integrity of the data and data analyses.

  • Funding This research was funded by the National Institute for Occupational Safety and Health (NIOSH) following competitive award as an intramural research project under the National Occupational Research Agenda (NORA), and made possible through the cooperation and support of the DOE and their employees and contractors.

  • Disclaimer The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Institute for Occupational Safety and Health.

  • Competing interests None.

  • Ethics approval The NIOSH Institutional Review Board (IRB).

  • Provenance and peer review Not commissioned; externally peer reviewed

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.