Breast cancer risk and protracted low-to-moderate dose occupational radiation exposure in the US Radiologic Technologists Cohort, 1983-2008.

BACKGROUND: Although high-dose ionising radiation is associated with increased breast cancer risks, the association with protracted low-dose-rate exposures remains unclear. The US Radiologic Technologist study provides an opportunity to examine the association between low-to-moderate dose radiation and breast cancer incidence and mortality.
METHODS: One thousand nine hundred and twenty-two self-reported first primary cancers were diagnosed during 1983-2005 among 66 915 female technologists, and 586 breast cancer deaths occurred during 1983-2008 among 83 538 female cohort members. Occupational breast dose estimates were based on work histories, historical data, and, after the mid-1970s, individual film badge measurements. Excess relative risks were estimated using Poisson regression with birth cohort stratification and adjustment for menopause, reproductive history, and other risk factors.
RESULTS: Higher doses were associated with increased breast cancer incidence, with an excess relative risk at 100 mGy of 0.07 (95% confidence interval (CI): -0.005 to 0.19). Associations were strongest for technologists born before 1930 (excess relative risk at 100 mGy=0.16; 95% CI: 0.03-0.39) with similar patterns for mortality among technologists born before 1930.
CONCLUSIONS: Occupational radiation to the breast was positively associated with breast cancer risk. The risk was more pronounced for women born before 1930 who began working before 1950 when mean annual doses (37 mGy) were considerably higher than in later years (1.3 mGy). However, because of the uncertainties and possible systematic errors in the occupational dose estimates before 1960, these findings should be treated with caution.

The risk of breast cancer has been associated with moderate and high ionising radiation doses in atomic bomb survivors and radiotherapy patients (United Nations Scientific Committee on the Effects of Atomic Radiation, 2008). However, data on breast cancer risk among women exposed to protracted low-dose radiation exposure are limited. Exposure–response trends for breast cancer have been reported for women who underwent multiple diagnostic x-rays during childhood and adolescence to assess spinal curvature from scoliosis (Ronckers ) and for tuberculosis patients with repeated fluoroscopic examinations (Howe and McLaughlin, 1996; Little and Boice, 1999). A cohort study of Russians (Eidemuller ; Ostroumova ; Davis ) with low-dose-rate environmental exposures found indications of a dose–response association. Most nuclear worker studies are comprised solely of men; however, a study in the United Kingdom that included women found no evidence of an association between radiation dose and breast cancer rates (Muirhead ). Most other occupational studies, including those of medical radiation workers (Linet ) and flight crews (Zeeb ), were based on work history factors rather than dose estimates.

Using self-reported work history data in the US Radiologic Technologists (USRT) study, we previously reported that technologists who worked before 1950 experienced elevated breast cancer incidence (Doody ) and mortality (Mohan ) compared with those who first worked after 1959. Hormonal, reproductive, lifestyle, anthropometric, and other known risk factors account for ∼40% of breast cancer cases (Madigan ; Sprague ), and our earlier analyses of breast cancer incidence in the USRT Cohort (Mohan ; Doody ) and the WeCare study of second breast cancers (Langholz ; Bernstein ) are the only studies that examined individual radiation exposure and breast cancer risks adjusting for a broad range of known lifestyle and reproductive risk factors.

The current paper, based on data from the USRT study, describes the first comprehensive dose–response analysis for breast cancer risk among medical radiation workers using individualised annual dose estimates that were based on badge dose measurements available in more recent years covering ∼30% of individual worker total doses, on average, together with information on work history and practices. For each worker, we calculated mean annual breast doses from 1000 dose realisations that reflect the uncertainty in individual annual doses from various sources (Simon ). Dose estimates have been corroborated by an analysis of chromosome translocation frequencies in a subset of USRT cohort members that found a radiation dose–response similar to that reported in other radiation-exposed populations (Little ). Using risk factor information obtained from three detailed questionnaires administered at roughly 10-year intervals, we investigated and adjusted for confounding by a wide range of known or suspected breast cancer risk factors.

Materials and methods

Target population and data sources

Details of the study population and methods are given in Boice and Doody and can be found online at http://www.radtechstudy.nci.nih.gov (Radiation Epidemiology Branch, 2015). Briefly, in the mid-1980s, the US National Cancer Institute, in collaboration with the University of Minnesota and the American Registry of Radiologic Technologists (ARRT), began a cohort study of 146 022 radiologic technologists, including 106 953 women, who were certified for at least 2 years between 1926 and 1982 (Sigurdson ). Vital status was ascertained from ARRT recertification records and linkage with the Social Security Administration. Cause of death was ascertained from the National Death Index through 31 December 2008.

Work history, behaviour, and cancer incidence data were obtained from three surveys: the first administered between 1983 and 1989 to 98 233 female technologists (71% response); the second administered to 94 508 women between 1994 and 1998 (74% response); and the third administered between 2003 and 2005 to 78 506 women (74% response) who completed at least one previous questionnaire.

Study population and follow-up

The target population for both incidence and mortality follow-up included 83 358 female technologists who responded to at least one questionnaire. Incidence analyses were limited to 66 915 out of the 83 358 technologists who completed at least two questionnaires and reported no cancer other than non-melanoma skin cancer before completing their first questionnaire (2366 women were excluded). Follow-up extended from the first (baseline) questionnaire answered until the earliest of a self-reported diagnosis of a cancer other than non-melanoma skin cancer or the last questionnaire answered. Mortality analyses were based on all 83 538 women who completed at least one questionnaire with follow-up from that response date until the date of death or 31 December 2008.

Breast cancer outcomes and medical validation

Incident cases were self-reported first primary malignancies occurring between the first and last questionnaires answered. Medical records were sought for all cases and obtained for 1611 of 1934 (83%), of which all but 12 (99%) were confirmed. We excluded the latter. Because medical records were obtained for a high percentage of self-reported cases, nearly all of which were confirmed, the 323 self-reported breast cancers for which medical records could not be obtained were treated as cases. Thus, the study included 1922 incident cases. Mortality analyses were based on 586 breast cancer deaths, 80% of which had been reported as an incident cases on at least one questionnaire.

Dose estimation

Dose reconstruction was undertaken to estimate individual annual occupational absorbed doses to specific organs (Simon ). Input data included: 684 465 annual badge dose measurements for 83 719 female cohort members between 1960 and 1997 from the nation's largest commercial dosimetry provider, military radiation dose registries, civilian employers, and a major US hospital; detailed information on work procedures and protection practices obtained from the cohort surveys; and literature-derived historical data on badge doses, x-ray imaging technology, radiation protection standards, protection practices, radiation energies and filtration; and other factors affecting exposure (Simon , 2006). The system provides 1000 realisations of breast dose for each cohort member in each year worked. The individual annual dose estimates used for these analyses were the arithmetic means of the 1000 realisations. Tables 1 and 2 present summary information on the doses. Further details on the dosimetry system are given in (Simon ) and Supplementary Data A.

Table 1

Summary information on dose, years worked, and year of birth by period worked, US Radiologic Technologists study, 1983–2008

Period worked	Total years worked	Mean years worked in period	Mean annual breast dose (mGy)	Badge dose available (%)a	Mean birth year
Pre-1930	390	3	84.6	0	1905
1930–1939	4016	4.3	82.2	0	1912
1940–1949	18 469	4.4	26.5	0	1921
1950–1959	66 033	4	7.4	0	1933
1960–1969	176 351	4.7	2.8	1	1941
1970–1979	415 873	5.8	1.5	18	1949
1980–1989	497 654	7.5	0.6	61	1949
1990+	356 814	6.7	0.3	73	1950
Total	1 535 600	6.2	1.9	42	1947

Abbreviation: mGy=milliGray.

Percentage of annual dose estimates for which a badge dose reading was available.

Table 2

Summary information on dose, years worked, and year of birth by year first worked, US Radiologic Technologists study, 1983–2008

Year first worked	Women	Mean total years worked	Mean total cumulative breast dose (mGy)	Mean birth year
Pre-1930	131	28.1	1168	1905
1930–1939	814	21.0	556	1913
1940–1949	3339	20.1	180	1923
1950–1959	12 100	18.2	50	1935
1960–1969	25 620	19.9	27	1945
1970–1979	40 379	17.4	14	1954
1980+	1155	11.4	6	1951
Total	83 538	18.4	37	1947

Abbreviation: mGy=milliGray.

Statistical methods

Our a priori approach for the dose–response analyses was to focus on breast cancer incidence using internal comparisons of excess relative risk (ERR) stratified by birth cohort and year first worked based on findings from our earlier work-history analyses. We also assessed excess absolute risks. As a secondary component of our analysis, we examined dose–response ERR for breast cancer mortality. In sensitivity analyses, we evaluated dose–response ERR for breast cancer incidence using external comparisons with US incidence rates from SEER population-based cancer registries.

Poisson regression was used to calculate ERRs and 95% confidence intervals (CIs) for breast cancer incidence or mortality, in a person-year table stratified by age, 5-year lagged cumulative breast dose, and other factors (for more details see Supplement Data B). Rates were described using ERR models of the form , in which represents the background rates depending on age (a), birth cohort (b) and other factors (z0). The ERR depends on dose (d) and, as necessary, various effect modifiers (z1). The primary risk model is described in detail in Supplementary Data D.

We evaluated known or suspected breast cancer risk factors to identify those which improved model fit or had an appreciable effect on inference about the radiation dose–response. Factors included in the final models were attained age, birth cohort, number of live births, menopausal status, body mass index (BMI) at baseline, family breast cancer history, alcohol consumption, and use of hormone replacement therapy. Other potential risk factors, including race, marital status, smoking, and other reproductive factors were not included in the primary analyses as their inclusion neither improved model fit nor had any appreciable effect on inference about the radiation dose–response association. With the exception of menopause age and status, missing data for risk factors were represented by missing-value indicator variables. A single imputation method was used to assign menopause age to women who were premenopausal at the time of last reported menopause status (Supplementary Data C). Standardised incidence/mortality ratio (SIR, SMR) analyses based on SEER incidence and mortality rates (Surveillance Epidemiology and End Results (SEER) Program, 2015) were used for some supplementary analyses.

Maximum-likelihood estimates were computed using AMFIT (Preston ). Two-sided hypothesis tests and CI's were based on likelihood ratio tests (Cox and Hinckley, 1974).

Results

Incidence

Among 66 915 women with 1 089 502 years of follow-up and 1922 reported first primary breast cancers, incidence increased with age with a complex interaction between age and menopausal status. Supplementary Data D contains information on the fitted rate model and includes a plot illustrating how rates varied with age, birth cohort, and menopausal age. Incidence did not vary substantially by race or marital status, but was higher in those with more education, current and former smokers, greater reported alcohol consumption, higher BMI at baseline, and a family history of breast cancer (see Supplementary Table D1 for details). Risks were higher in nulliparous women, declined with increasing number of children, and were slightly increased among those using hormone replacement therapy (see Supplementary Table D2 for details). There was no association with age at menarche. Age-specific breast cancer rates have decreased markedly with increasing birth year in the USRT (Supplementary Table D1).

The magnitude of the estimated linear ERR dose–response slope depended on the nature of the birth cohort adjustment. Without birth cohort adjustment, there was a positive dose–response association (P<0.001) with an estimated ERR at 100 mGy of 0.12 (95% CI: 0.03–0.25). In contrast, using a simple log-linear adjustment for birth cohort, the estimated ERR at 100 mGy was 0.05 (95% CI: -0.02 to 0.16; P=0.4). Adjustment for birth cohort had little effect on other parameters in the risk model (Supplementary Table D1). In view of the unusual nature of the simple birth cohort trend, our primary analyses were stratified on birth cohort. We also examined the extent to which the breast cancer risk varied with time-dependent length of employment as a radiological technologist. These analyses (details not shown) provided no indication of a significant trend in risk with employment duration and the estimated magnitude of the trend was small, ∼1% change in risk per decade increase in employment length.

The upper portion of Table 3 and Figure 1A present results of a stratified analysis with separate ERR estimates for each of four birth cohorts (before 1930, 1930–1939, 1940–1949, and 1950+). Supplementary Table D2 contains information on the birth cohort-specific distributions of various risk factors. A linear dose–response association was observed for cohort members born before 1930 with an estimated ERR at 100 mGy of 0.16 (95% CI: 0.03–0.39). There was no indication of a positive dose–response association in the other birth cohorts and weak evidence of heterogeneity (P=0.08). When the ERR was assumed constant over strata, the estimate (ERR at 100 mGy=0.07; 95% CI: −0.10 to 0.19, P=0.09) was slightly larger than that in an unstratified model in which baseline rates were assumed to vary log-linearly with birth cohort. Table 4 presents ERR estimates within year first worked categories (before 1950, 1950–1959, 1960–1969, and 1970+). There was a dose–response association in the earliest category (ERR at 100 mGy=0.11) with no indication of increased risks for the other three groups, but only weak evidence of heterogeneity over first year worked categories (P=0.10). Estimated excess absolute rates exhibited similar variation by birth cohort (results not shown).

Table 3

ERR estimates for breast cancer incidence and mortality by birth cohort, US Radiologic Technologists study, 1983–2008

Year of birth	ERR at 100 mGya	95% CI	P-value	Fitted excess cases	Observed cases	Mean dose (mGy)
Incidence
Before 1930	0.16	0.03, 0.39	0.01	44.2	216	167
1930–1939	−0.1	<−0.3, 0.11	0.41	−21	406	49
1940–1949	−0.14	<−0.5, 0.22	0.4	−25.8	745	25
1950+	−0.25	<−0.5, 0.31	0.27	−27.9	555	14
		Phet=0.08
Common ERR	0.07	−0.005, 0.19	0.09	54.3	1922	30
Mortality
Before 1930	0.26	0.08, 0.63	<0.001	43.7	128	217
1930–1939	0.28	−0.08, 0.89	0.26	19.5	159	51
1940–1949	1.61	0.20, 4.50	0.01	52.6	190	25
1950+	0.61	<−0.5, 4.70	>0.5	8.6	109	14
		Phet=0.35
Common ERR	0.31	0.11, 0.67	<0.001	88.5	586	37

Abbreviations: BMI=body mass index; CI=confidence interval; ERR=excess relative risk; Phet=P-value of test for hetergeneity over birth cohorts.

Adjusted for attained age, birth cohort category, menopause status, family history of breast cancer, baseline BMI and use of hormone replacement therapy.

Figure 1

Birth cohort-specific fitted breast cancer incidence (A) and mortality (B) dose–response curves (solid lines) with categorical excess relative risk estimates (ERRs) (closed circles) and 95% confidence limits. The dashed lines are the fitted linear dose–response curves assuming that the ERR is the same in all birth cohorts.

Table 4

ERRs and 95% CIs for breast cancer incidence by year first worked, US Radiologic Technologists study, 1983–2008

Year first worked	ERR at 100 mGya	95% CI	P-value	Fitted excess cases	Observed cases	Mean dose (mGy)
Before 1950	0.11	0.03, 0.39	0.03	31	163	213
1950–1959	−0.04	<−0.3, 0.11	>0.5	−8.6	450	49
1960–1969	−0.19	<−0.5, 0.22	0.18	−25.8	684	26
1970+	−0.33	<−0.5, 0.33	0.22	−27.9	625	14
		Phet=0.10
Common ERR	0.07	−0.005, 0.19	0.07	54.3	1922	30

Abbreviations: BMI=body mass index; CI=confidence interval; ERR=excess relative risk; mGy=milliGray.

Adjusted for attained age, birth cohort category, menopause status, family history of breast cancer, baseline BMI and use of hormone replacement therapy.

The SEER-based SIR estimate was 1.08 (95% CI: 1.04–1.14). Birth cohort-specific ERR estimates in an SIR model were similar to those from parametric models (details not shown).

Mortality

The mortality analyses were based on 83 538 women in the USRT Cohort with 1 722 303 years of mortality follow-up who completed at least one study questionnaire. There were 586 breast cancer deaths, of which 119 were in women who self-reported breast cancer and were included in the cohort of 66 915 women in the incidence analyses who completed at least two study questionnaires and 349 among women who reported a breast cancer diagnosis before the date of completion of the first survey. Only 20% (118) of the breast cancer deaths occurred among women who had never reported breast cancer. While breast cancer was the most common cause of death among the 1922 breast cancer cases in the incidence analyses, most of the women in this group (87%) were still alive at the end of the mortality follow-up.

The associations between non-radiation lifestyle and reproductive factors and breast cancer mortality were generally similar to those seen for incidence. Age-specific breast cancer mortality rates for cohort members decreased over time (see Supplement Data E for details).

Without birth cohort adjustment, there was a positive dose–response association (P<0.001) with an ERR at 100 mGy of 0.71 (95% CI: 0.31–1.43) and no evidence of non-linearity (P>0.5). As shown in the lower portion of Table 3, in a stratified analysis in which the ERR was assumed to be the same in all birth cohorts, the estimated ERR at 100 mGy was 0.31 (P<0.001). Although there was no evidence of significant heterogeneity in the ERRs over birth cohort (P=0.35), as with incidence, the experience of the early workers is the primary determinant of the risk (see Table 3 and Figure 1B). As with incidence, when birth cohort associations were described using a log-linear birth cohort trend, evidence for a radiation association was greatly reduced: the estimated ERR at 100 mGy was essentially zero (95% CI: <−0.03 to 0.10; P>0.5). An SMR analysis indicated that mortality was slightly lower than expected (SMR=0.9; 95% CI 0.82–0.97) with birth cohort-specific ERR estimates similar to those based on internal comparisons (details not shown).

Discussion

Analyses of breast cancer incidence and mortality rates in the USRT Cohort stratified on birth cohort suggest that low-to-moderate-dose protracted occupational radiation exposure is associated with increased breast cancer risk. Statistically significant trends were seen for incidence and mortality among women born before 1930. The findings, especially for incidence, are less clear for more recent birth cohorts where lifetime cumulative doses are much lower.

There were no indications of a dose–response association for either incidence or mortality when the variation in age-specific rates with birth cohort was modelled using continuous log-linear or log-quadratic functions of birth year. However, using such simple parametric birth cohort trend models may underestimate the radiation dose–response association because of intrinsic confounding between estimated cumulative breast dose and birth year (r=−0.58) and the marked variation in age-specific rates with birth cohort. The effect of this confounding is exacerbated by the large uncertainties in individual dose estimates for the early (higher dose) years when badge dose readings were unavailable. Concern about intrinsic confounding is supported by the contrast between the decline in age-specific incidence rates with birth cohort in the USRT study and the increase with some leveling off in recent years seen in US cancer registry data (Devesa ; Edwards ) and in age period cohort analyses of Danish (Anderson ) and Norwegian (Weedon-Fekjaer ) population data.

Age-specific breast cancer mortality rates have declined over time in the USRT Cohort, a pattern broadly similar to that seen in age-standardised US breast cancer mortality rates (Tarone ; Jatoi ). However, the marked effect of adjusting for birth cohort using smooth models on the radiation dose–response further suggests confounding between dose and birth cohort.

Despite this confounding, both the incidence and mortality data indicated that radiation is associated with increased breast cancer risks in the earliest workers with little evidence of increased risk in more recent workers who received lower doses. The findings from the mortality analyses are consistent with results of the work-history-based analyses that showed increased breast cancer mortality for those who started working before 1950. As there was some indication of heterogeneity in the radiation ERRs with year first worked or birth cohort, the birth cohort-stratified analyses with a single ERR trend estimate appear to provide the most useful characterisation of the breast cancer dose–response association.

The USRT has a number of strengths. It is a large nationwide cohort with individualised dose estimates, long-term follow-up, detailed information on potential confounders, and a high confirmation rate for the self-reported incident cases. The individualised bone marrow dose estimates were used in a dose–response analysis of chromosome translocation frequency (Little ) that found radiation dose–response trends similar to those seen in other populations.

This study has a number of limitations. For technologists who began working in more recent decades, power to detect associations was limited owing to low cumulative occupational radiation doses. Incidences may have been under-reported during follow-up as this information was captured in follow-up surveys administered every 8 to 10 years. Underascertainment of breast cancers should not have led to biased risk estimates unless cohort members with higher doses were more likely to report a case or if dose was more strongly associated with aggressive and rapidly fatal breast tumour types that were less likely to have been captured via self-report. However, Freedman found that only ∼3% of breast cancers were not reported by survey respondents, suggesting that underascertainment is unlikely to be an important source of bias. As dose estimates and risk factor data were available only for technologists who survived to at least the mid-1980s or, for the incidence analyses, the mid-1990s, older and higher dose technologists are under-represented in these analyses. While this reduces the power of the study, it seems unlikely that it would bias the risk estimates for the study period.

Despite the carefully designed dosimetry system, uncertainties, and, more importantly, the potential for systematic errors in the dose estimates before 1960, and for most cohort members, the mid-1970s, for the early years when badge dose measurements were lacking and information on dose distributions is scarce are the most important limitations of this analysis. Only ∼30% of individual worker total doses were based on actual badge readings. The correlation between birth cohort and total cumulative dose, together with the uncertainty in the dose estimates (particularly for the early years), make it difficult to disentangle birth cohort and dose–response associations. Lifestyle differences between early workers and the general population may account for some of the unusual birth cohort trends, but the persistence of these differences after adjustment and the similarity of the risk factor associations to those observed in other studies (Colditz ), suggests that lifestyle differences are unlikely to account for the patterns observed. Table 5 contrasts the breast cancer incidence and mortality ERR estimates for the USRT population with those from several other cohort studies. Although the USRT mortality risk estimate is roughly two times that in atomic bomb survivor studies (Preston ; Ozasa ), the incidence ERR for the current analysis is similar to the risks for cohorts in the pooled analysis of atomic bomb survivors and patients exposed to radiation from medical sources (Preston ). The difference in the ERRs for Japanese and western women appeared to be largely due to marked differences in the baseline rates in Japan from those in western countries. The latest incidence ERR estimate in the cohort of Russian women who received chronic low-dose-rate environmental exposures (Davis ) is only slightly higher than that seen in this study, whereas the breast cancer mortality ERR in the Mayak plutonium production workers (Sokolnikov ) was somewhat lower than our USRT estimate.

Table 5

ERR estimates for breast cancer incidence and mortality in the USRT study (1983–2008) and previous studies

	Incidence			Mortality
	Cases	ERR at 100 mGy	95% CI	Deaths	ERR at 100 mGy	95% CI
USRT study	1922	0.07	−0.005, 0.19	586	0.31	0.11, 0.67
Atomic bomb survivorsa	1073	0.19b	0.12, 0.28	324	0.16	0.09, 0.24
Pooled analysisc	245	0.07	0.04, 0.12
Techa River Cohortd	109	0.19	−0.06, 0.61
Mayak Worker Cohorte				107	0.02	−0.01, 0.06

Abbreviations: CI=confidence interval; ERR=excess relative risk; mGy=milliGray; USRT=US Radiologic Technologists.

From Preston (incidence) and Ozasa (mortality).

Scaled to attained age 50 years to correspond to the average attained age at diagnosis in the current USRT data.

The estimate shown here was based on the results for the three US cohorts used in that analysis (the two Massachusetts tuberculosis cohorts and a US thymic irradiation cohort) for the pooled ERR model shown in Table 9 of Preston . The estimate corresponds to the ERR at attained age 50 years.

From Davis .

From Sokolnikov .

The USRT is the first large cohort of medical radiation workers with dose estimates based on a comprehensive, individualised dose reconstruction, extended follow-up, large numbers of female technologists and breast cancer cases, and extensive information on breast cancer risk factors and potential confounders. Despite the difficulty of disentangling variations in breast cancer rates with dose and birth cohort and the shortcomings of the early-worker dose estimates, a radiation dose–response association was observed for both incidence and mortality in those born before 1930. Occupational radiation dose to the breast was positively associated with breast cancer risk. The risk was more pronounced for women born before 1930 who began working before 1950 when mean annual doses (37 mGy) were considerably higher than in later years (1.3 mGy). However, because of the uncertainties and possible systematic errors in the occupational dose estimates before 1960 these findings should be treated with caution.