94227592553Occup Environ MedOccup Environ MedOccupational and environmental medicine1351-07111470-792623000827456386110.1136/oemed-2012-100906HHSPA717197ArticleRisk of leukaemia mortality from exposure to ionising radiation in US nuclear workers: a pooled case-control studyDanielsRobert DBertkeStephenWatersKathleen MSchubauer-BeriganMary KNational Institute for Occupational Safety and Health (NIOSH), Division of Surveillance, Hazard Evaluations, and Field Studies (DSHEFS), Industrywide Studies Branch (IWSB), Cincinnati, Ohio, USACorrespondence 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@cdc.gov248201521920121201309920157014148Objective

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.

INTRODUCTION

Most information on the leukaemogenicity of ionising radiation stems from studies of radiation therapy patients and survivors of the atomic bomb blasts in 1945.12 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.37 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 al34 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.

METHODSStudy 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,34 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.

Exposure

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.1011 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.41113 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.48 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.4814 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 form

RR=(1+αj=13γjDj)ei=1nβiXi 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 cut-points 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

RESULTSDescriptive 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.

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.

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.

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).

DISCUSSION

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,34 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,2425 ionising radiation42628 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.

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.428 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.223233 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.48 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 benzene2434 and ionising radiation.3536 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.4041 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.4248 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.

CONCLUSIONS

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.

Supplementary Material

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.

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.

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.

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

Committee to Assess Health Risks from Exposure to Low Levels of Ionizing Radiation, National Research CouncilHealth Risks from Exposure to Low Levels of Ionizing Radiation: BEIR VII—Phase 2Washington, DCNational Academies Press2006406United Nations Scientific Committee on the Effects of Atomic RadiationUNSCEAR 2006 Report: Effects of Ionizing Radiation, Annex A: Epidemiological Studies of Radiation and CancerNew York, NY2009310Schubauer-BeriganMKDanielsRDFlemingDAChronic lymphocytic leukaemia and radiation: findings among workers at five US nuclear facilities and a review of the recent literatureBr J Haematol200713979980817922878Schubauer-BeriganMKDanielsRDFlemingDARisk of chronic myeloid and acute leukemia mortality after exposure to ionizing radiation among workers at four U.S. nuclear weapons facilities and a nuclear naval shipyardRadiat Res20071672223217390730DanielsRDSchubauer-BeriganMKA meta-analysis of leukaemia risk from protracted exposure to low-dose gamma radiationOccup Environ Med2011684576420935290MuirheadCRO’HaganJAHaylockRGMortality and cancer incidence following occupational radiation exposure: third analysis of the national registry for radiation workersBr J Cancer20091002061219127272CardisEVrijheidMBlettnerMThe 15-country collaborative study of cancer risk among radiation workers in the nuclear industry: estimates of radiation-related cancer risksRadiat Res200716739641617388693Schubauer-BeriganMKMacievicGVUtterbackDFAn epidemiologic study of mortality and radiation-related risk of cancer among workers at the Idaho National Engineering and Environmental Laboratory, a US Department of Energy facilityCincinnati, OHNational Institute for Occupational Safety and Health2005220BeaumontJJSteenlandKMintonAA computer program for incidence density sampling of controls in case-control studies nested within occupational cohort studiesAm J Epidemiol1989129212192910063DanielsRDKubaleTLSpitzHBRadiation exposure from work-related medical X-rays at the Portsmouth Naval ShipyardAm J Ind Med2005472061615712259AndersonJLDanielsRDBone marrow dose estimates from work-related medical x-ray examinations given between 1943 and 1966 for personnel from five U.S. nuclear facilitiesHealth Phys2006905445316691102DanielsRDLodwickCJSchubauer-BeriganMKAssessment of plutonium exposures for an epidemiological study of US nuclear workersRadiat Prot Dosimetry2006118435516081492DanielsRDSchubauer-BeriganMKBias and uncertainty of penetrating photon dose measured by film dosemeters in an epidemiological study of US nuclear workersRadiat Prot Dosimetry20051132758915769802YiinJHSchubauer-BeriganMKSilverSRRisk of lung cancer and leukemia from exposure to ionizing radiation and potential confounders among workers at the Portsmouth Naval ShipyardRadiat Res20051636031315913392SAS Software, Version 9.2, Copyright (2002–2003) SAS and all other SAS Institute Inc. product or service names are registered trademarks or trademarks of SAS InstituteCary, NC, USASAS Institute IncLangholzBRichardsonDBFitting general relative risk models for survival time and matched case-control analysisAm J Epidemiol20101713778320044379HarrellFERegression modeling strategies with applications in linear models, logistic regression, and survival analysisNew York, NYSpringer-Verlag New York Inc2001RichardsonDBColeSRChuHLagging exposure information in cumulative exposure-response analysesAm J Epidemiol201117414162222047823AkaikeHInformation theory and an extension of the maximum likelihood principlePetrovBNCFSecond International Symposium on Information TheoryBudapest197326781International Commission on Radiological ProtectionConversion coefficients for use in radiological protection against external radiation. Adopted by the ICRP and ICRU in September 1995Ann ICRP1996261205KocherDCApostoaeiAIHoffmanFORadiation effectiveness factors for use in calculating probability of causation of radiogenic cancersHealth Phys20058933215951689PuskinJSNelsonCBPawelDJEstimating radiological cancer risks, addendum: uncertainty analysisWashington, DCU.S. Environmental Protection Agency199937Environmental Protection AgencyEPA Radiogenic cancer risk models and projections for the US populationWashington, DCU.S. Environmental Protection Agency2011175SilverSRRinskyRACooperSPEffect of follow-up time on risk estimates: a longitudinal examination of the relative risks of leukemia and multiple myeloma in a rubber hydrochloride cohortAm J Ind Med200242481912439871FinkelsteinMMLeukemia after exposure to benzene: temporal trends and implications for standardsAm J Ind Med2000381710861761RichardsonDSugiyamaHNishiNIonizing radiation and leukemia mortality among Japanese atomic bomb survivors, 1950–2000Radiat Res20091723688219708786RichardsonDBWingSLeukemia mortality among workers at the Savannah River SiteAm J Epidemiol200716610152217660455ShilnikovaNSPrestonDLRonECancer mortality risk among workers at the Mayak nuclear complexRadiat Res20031597879812751962RichardsonDBTerschurenCPohlabelnHTemporal patterns of association between cigarette smoking and leukemia riskCancer Causes Control200819435017906957VrijheidMCardisEAshmorePIonizing radiation and risk of chronic lymphocytic leukemia in the 15-country study of nuclear industry workersRadiat Res2008170661518959468CardisEGilbertESCarpenterLEffects of low doses and low dose rates of external ionizing radiation: cancer mortality among nuclear industry workers in three countriesRadiat Res1995142117327724726GroganHASinclairWKVoillequePGRisks of fatal cancer from inhalation of 239,240 plutonium by humans: a combined four-method approach with uncertainty evaluationHealth Phys2001804476111316075SpiersFWVaughanJThe toxicity of the bone seeking radionuclidesLeuk Res198913347502747266RinskyRAHornungRWSilverSRBenzene exposure and hematopoietic mortality: a long-term epidemiologic risk assessmentAm J Ind Med2002424748012439870BoiceJDCohenSSMummaMTUpdated mortality analysis of radiation workers at rocketdyne (atomics international), 1948–2008Radiat Res20111762445821381866KrestininaLPrestonDLDavisFGLeukemia incidence among people exposed to chronic radiation from the contaminated Techa River, 1953–2005Radiat Environ Biophys20104919520120012750BoiceJDCohenSSMummaMTMortality among radiation workers at Rocketdyne (atomics international), 1948–1999Radiat Res20061669811516808626KrestininaLYPrestonDLOstroumovaEVProtracted radiation exposure and cancer mortality in the Techa River CohortRadiat Res20051646021116238437StaynerLSteenlandKDosemeciMAttenuation of exposure-response curves in occupational cohort studies at high exposure levelsScand J Work Environ Health2003293172412934726NakamuraNA hypothesis: radiation-related leukemia is mainly attributable to the small number of people who carry pre-existing clonally expanded preleukemic cellsRadiat Res20051632586515733032HallEJBrennerDJWorgulBGenetic susceptibility to radiationAdv Space Res2005352495315934202BrennerDJSachsRKEstimating radiation-induced cancer risks at very low doses: rationale for using a linear no-threshold approachRadiat Environ Biophys200644253616470411TubianaMAurengoAAverbeckDThe debate on the use of linear no threshold for assessing the effects of low dosesJ Radiol Prot2006263172416926474ScottBRLow-dose radiation risk extrapolation fallacy associated with the linear-no-threshold modelHum Exp Toxicol200827163818480143BreckowJLinear-no-threshold is a radiation-protection standard rather than a mechanistic effect modelRadiat Environ Biophys2006442576016468065WallBFKendallGMEdwardsAAWhat are the risks from medical X-rays and other low dose radiation?Br J Radiol2006792859416585719CharlesMWLNT—an apparent rather than a real controversy?J Radiol Prot200626325916926475ChadwickKHLeenhoutsHPRadiation risk is linear with dose at low dosesBr J Radiol20057881015673520

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).

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

Radiation weights*
Linear ERR per 100 mGy (95% CI)Relative change in ERR (%)
PlutoniumNeutron
000.09 (−0.14 to 0.51)NA
171 (solved)406 (solved)0.07 (<0 to 0.49)−22
20100.15 (−0.06 to 0.57)67
1100.11 (−0.13 to 0.54)22
3.6110.13 (−0.12 to 0.56)33

Radiation weights used to construct weighted absorbed dose to the red bone marrow.

ERR, excess relative risk; NA, not applicable.

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

OutcomeNo. casesAdjusted forBest fit exposure lag ERR per 100 mGy (95% CI)
Best fit exposure window ERR per 100 mGy (95% CI)
Lag (years)Lag modelWindow (years)Window model
Non-CLL264Sex, race, hire year20.09 (−0.17 to 0.65)6 to 141.9 (<0 to 8.0)
AML150Birth cohort1.10.22 (−0.19 to 1.2)6 to 147.0 (0.079 to 32.0)
CML52Sex, race, hire year6.30.29 (<0 to 2.6)6.3 to 143.8 (−1.2 to 26.0)
CLL74Sex, race, birth cohort24−0.032 (−0.16 to 0.74)11 to 14−2.8 (NC to NC)

‘Dose’ was defined as absorbed dose to red bone marrow from low linear energy transfer radiation exposure.

Excluding leukaemias of indeterminate subtype.

AML, acute myeloid leukaemia; CML, chronic myeloid leukaemia; CLL, chronic lymphocytic leukaemia; ERR, excess relative risk; NC, not calculable.

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

Outcome Adjust forNon-CLL (n=264) race, sex and hire date
AML (n=150) birth date
CML (n=52) race, sex and hire date
Model Fit statistic*Lag (2 years)
Windows (6–14 years)
Lag (1.1 years)
Windows (6–14 years)
Lag (6.3 years)
Windows (6.3–14 years)
AICLRT p valueAICLRT p valueAICLRT p valueAICLRT p valueAICLRT p valueAICLRT p value
Linear850.9NA847.6NA480.9NA473.7NA168.0NA166.2NA
Linear-quadratic852.10.382846.90.104481.70.281NCNC170.00.962NCNC
RCS844.70.006845.60.046479.40.065473.50.136170.60.488NCNC
PLS846.20.013847.60.135NMNMNMNM
Categorical852.2NA849.8NA486.1NA477.7NA170.9NA169.6NA

LRT degrees of freedom (df): Linear-quadratic all=1; RCS lag=2, window=1; and PLS lag=4, window=2.

Best fitting PLS models were obtained using two knots (8.3 and 46.9 mGy) in the 2-year-lagged model and one knot (20.1 mGy) in the windows-based model.

AIC, Akaike Information Criterion; AML, acute myeloid leukaemia; CML, chronic myeloid leukaemia; CLL, chronic lymphocytic leukaemia; LRT, likelihood ratio test; NA, not applicable; NC, not calculable; NM, not modelled; PLS, piecewise linear spline; RCS, restricted cubic spline.

Leukaemia ERR at 100 mGy exposure by subtype and by study

StudyNon-CLLAMLCMLCLL
Current study0.09 (95% CI −0.17 to 0.65)0.22 (95% CI −0.19 to 1.2)0.29 (95% CI <0 to 2.6)−0.03 (95% CI −0.16 to 0.74)
264 decedents150 decedents52 decedents74 decedents
LCCS340.26 (95% CI <−0.10 to 1.03)NRNR−0.20 (95% CI <0 to 1.4)
184 decedents43 decedents
15-country study*7300.196 (90% CI −0.026 to 0.59)−0.41 (90% CI −0.76 to 0.17)1.0 (90% CI −0.086 to 4.02)−0.17 (95% CI −0.89 to 0.40)
196 decedents81 decedents45 decedents47 decedents
Three-country study310.22 (90% CI 0.013 to 0.57)0.34 (90% CI <0 to 1.49)1.1 (90% CI 0.29 to 3.09)−0.09 (90% CI <0 to 0.73)
119 decedents32 decedents28 decedents27 decedents
NRRW60.17 (90% CI 0.006 to 0.43)0.12 (90% CI −0.12 to 0.65)0.33 (90% CI 0.04 to 0.93)<−0.19 (90% CI <−0.19 to 0.12)
198 decedents102 decedents44 decedents69 decedents
LSS260.47 (90% CI 0.35 to 0.64)0.43 (90% CI 0.27 to 0.66)0.64 (90% CI 0.30 to 1.37)NR
310 decedents124 decedents58 decedents

AML and CLL results were calculated using a log-linear relative risk model of the form RR(d)=exp(βd) and ERR=RR−1.

Based on a linear ERR model without effect modification of the form ERR(d)=βd.

Also included seven CLL cases and 121 leukaemias not specified as AML or CML.

AML, acute myeloid leukaemia; CML, chronic myeloid leukaemia; CLL, chronic lymphocytic leukaemia; ERR, excess relative risk; LCCS, NIOSH Leukaemia Mortality Case control Study; LSS, Life-Span Study of Japanese Atomic Bomb Survivors; NC, not calculable; NR, not reported; NRRW, UK National Registry of Radiation Workers.

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.