Selection of reference groups in the Life Span Study of atomic bomb survivors

Comparison of All Solid Cancer Mortality and Incidence Dose-Response in the Life Span Study of Atomic Bomb Survivors, 1958-2009

Recent analysis of all solid cancer incidence (1958-2009) in the Life Span Study (LSS) revealed evidence of upward curvature in the radiation dose response among males but not females. Upward curvature in sex-averaged excess relative risk (ERR) for all solid cancer mortality (1950-2003) was also observed in the 0-2 Gy dose range. As reasons for non-linearity in the LSS are not completely understood, we conducted dose-response analyses for all solid cancer mortality and incidence applying similar methods [1958-2009 follow-up, DS02R1 doses, including subjects not-in-city (NIC) at the time of the bombing] and statistical models. Incident cancers were ascertained from Hiroshima and Nagasaki cancer registries, while cause of death was ascertained from death certificates throughout Japan. The study included 105,444 LSS subjects who were alive and not known to have cancer before January 1, 1958 (80,205 with dose estimates and 25,239 NIC subjects). Between 1958 and 2009, there were 3.1 million person-years (PY) and 22,538 solid cancers for incidence analysis and 3.8 million PY and 15,419 solid cancer deaths for mortality analysis. We fitted sex-specific ERR models adjusted for smoking to both types of data. Over the entire range of doses, solid cancer mortality dose-response exhibited a borderline significant upward curvature among males (P = 0.062) and significant upward curvature among females (P = 0.010); for solid cancer incidence, as before, we found a significant upward curvature among males (P = 0.001) but not among females (P = 0.624). The sex difference in magnitude of dose-response curvature was statistically significant for cancer incidence (P = 0.017) but not for cancer mortality (P = 0.781). The results of analyses in the 0-2 Gy range and restricted lower dose ranges generally supported inferences made about the sex-specific dose-response shape over the entire range of doses for each outcome. Patterns of sex-specific curvature by calendar period (1958-1987 vs. 1988-2009) and age at exposure (0-19 vs. 20-83) varied between mortality and incidence data, particularly among females, although for each outcome there was an indication of curvature among 0-19-year-old male survivors in both calendar periods and among 0-19-year-old female survivors in the recent period. Collectively, our findings indicate that the upward curvature in all solid cancer dose response in the LSS is neither specific to males nor to incidence data; its evidence appears to depend on the composition of sites comprising all solid cancer group and age at exposure or time. Further follow up and site-specific analyses of cancer mortality and incidence will be important to confirm the emerging trend in dose-response curvature among young survivors and unveil the contributing factors and sites.

The possible impact of passive smoke exposure on radiation-related risk estimates for lung cancer among women: the life span study of atomic bomb survivors

PURPOSE

Analyses of the Life Span Study cohort of atomic bomb survivors have shown a statistically significant sex difference in the excess risk of incident lung cancer due to radiation exposure, with the radiation-related excess relative risk per gray (ERR/Gy) for women approximately 4 times that for men, after accounting for active smoking. We sought to determine the extent to which this risk difference could be explained by adjustment for passive smoke exposure, which is a known risk factor for lung cancer that was not measured among Life Span Study participants, and which could be particularly influential among female never-smokers.

MATERIALS AND METHODS

The Life Span Study includes survivors of the atomic bombings of Hiroshima and Nagasaki and city residents who were not in either city at the time of the bombings, matched to survivors on city, sex, and age. First primary lung cancers were identified from population-based cancer registries between 1958 and 2009. Data on active smoking were obtained from mailed surveys and in-person questionnaires (1965-1991). We calculated passive smoke exposure for female never-smokers by attributing smoking pack-years at various intensities (5-50%) based on smoking patterns among men, stratified by city, birth year, radiation dose, and lung cancer status. Poisson regression models with additive and multiplicative interactions between radiation dose and smoking were used to estimate sex-specific radiation-related excess relative risks for lung cancer.

RESULTS

During the study period, 2,446 first primary lung cancers were identified among 105,444 study participants. On average, male smokers started smoking 19.5 cigarettes per day at 21.5 years old. Partially attributing male smoking patterns to female never-smokers-to approximate passive smoke exposure-yielded lower radiation-related ERR/Gy estimates for women under a multiplicative radiation-smoking interaction model, leading to a lower female-to-male ratio of ERR/Gy estimates; however, this difference was evident only at very high passive smoke intensities. Under an additive radiation-smoking interaction model, the results were unchanged.

CONCLUSIONS

Our results are consistent with the possibility that failure to account for passive smoke might contribute, in small part, to the higher radiation risk estimates for lung cancer among women compared to men in the Life Span Study.

Lifetime risk of suicide among survivors of the atomic bombings of Japan

AIMS

The long-term physical health effects of the atomic bombings of Hiroshima and Nagasaki are well characterised, but the psychological effects remain unclear. Therefore, we sought to determine whether measures of exposure severity, as indirect measures of psychological trauma arising from exposure to the atomic bombings, are associated with suicide mortality among atomic bomb survivors.

METHODS

The Life Span Study is a prospective cohort study of 93 741 Japanese atomic bomb survivors who were located within 10 km of the hypocentre in Hiroshima or Nagasaki at the time of the bombings in 1945, and 26 579 residents of Hiroshima and Nagasaki who were not in either city at the time of the bombings, matched to survivors on city, sex and age. Measures of exposure severity included: proximity to the hypocentre, type of shielding between the survivor and the blast and self-reported occurrence of acute radiation and thermal injuries. Date of death was obtained from the Japanese National Family Registry system. Cause of death was obtained from death certificates. Adjusted hazard ratios (HRs) were estimated from Cox regression models overall and stratified by sex and age.

RESULTS

During the 60-year follow-up period (1950-2009), 1150 suicide deaths were recorded among 120 231 participants (23.6 per 100 000 person-years): 510 among 70 092 women (17.2 per 100 000 person-years) and 640 among 50 139 men (33.6 per 100 000 person-years). Overall, there was no association of proximity, type of shielding or the occurrence of acute injuries with suicide mortality. Among those <25 years of age at the time of the bombings, increased suicide risk was observed for survivors outside v. shielded inside any structure (HR: 1.24; 95% confidence interval (CI): 1.03, 1.48; interaction p = 0.054) and for those who reported flash burns (HR: 1.32; 95% CI: 1.00, 1.73; interaction p = 0.025). Sex-stratified analyses indicated that these associations were limited to men. Among women, closer proximity to the hypocentre was associated with a non-significant increase in suicide risk, with a positive association between proximity and suicide risk observed among women <15 years of age (HR: 1.09 per km; 95% CI: 1.00, 1.18; interaction p = 0.067).

CONCLUSIONS

Proximity to the hypocentre, shielding and acute injury presence do not generally appear to influence suicide mortality among atomic bomb survivors. However, heterogeneity may exist by age and sex, with younger survivors potentially more sensitive to psychological trauma. Coupled with other studies, our results suggest the importance of long-term monitoring of mental health among young populations exposed to catastrophic events or mass trauma.

Radiation Risks for the Incidence of Kidney, Bladder and Other Urinary Tract Cancers: 1958-2009

As part of the recent series of articles to create a comprehensive description of the radiation risks of solid cancer incidence after ionizing radiation exposure, based on the atomic bomb survivors' Life Span Study (LSS), this work focuses on the risks of urinary tract cancer (UTC) and kidney cancer. Analyses covered a 52-year period of follow-up, through 2009, among 105,444 eligible survivors who were alive and cancer free in 1958. This represents an additional 11 years of follow-up since the last comprehensive report, with a total of 3,079,502 person-years. We observed 790 UTC and 218 kidney cancer cases. Adjusted for smoking, there was a strong linear radiation dose response for UTC. The sex-averaged excess relative risk per 1 Gy (ERR/Gy) was 1.4 (95% confidence interval, CI: 0.82 to 2.1). Both males and females showed significantly increased ERRs/Gy with female point estimates at a factor of 3.4 (95% CI: 1.4 to 8.6) greater than male estimates. UTC radiation risks were largely unmodified by age at exposure or attained age. The attributable fraction of UTC to radiation exposure was approximately 18% while that attributed to smoking was 48%. Kidney cancer showed an increased ERR due to smoking (0.56 per 50 pack-years; 95% CI -0.007 to 1.6; P = 0.054), but we did not observe any strong associations of kidney cancer with radiation exposure, although sex-specific dose responses were found to be statistically different.

Radiation effects on atherosclerosis in atomic bomb survivors: a cross-sectional study using structural equation modeling

Past reports indicated that total-body irradiation at low to moderate doses could be responsible for cardiovascular disease risks, but the mechanism remains unclear. The purpose of this study was to investigate the association between radiation exposure and atherosclerosis, an underlying pathology of cardiovascular diseases, in the Japanese atomic bomb survivors. We performed a cross-sectional study measuring 14 clinical-physiological atherosclerosis indicators during clinical exams from 2010 to 2014 in 3274 participants of the Adult Health Study cohort. Multivariable analyses were performed by using a structural equation model with latent factors representing underlying atherosclerotic pathologies: (1) arterial stiffness, (2) calcification, and (3) plaque　as measured with indicators chosen a priori on the basis of clinical-physiological knowledge. Radiation was linearly associated with calcification (standardized coefficient per Gy 0.15, 95 % confidence interval: CI [0.070, 0.23]) and plaque (0.11, 95 % CI [0.029, 0.20]), small associations that were comparable to about 2 years of aging per Gy of radiation exposure, but not with arterial stiffness (0.036, 95 % CI [− 0.025, 0.095]). The model fitted better and had narrower confidence intervals than separate ordinary regression models explaining individual indicators independently. The associations were less evident when the dose range was restricted to a maximum of 2 or 1 Gy. By combining individual clinical-physiological indicators that are correlated because of common, underlying atherosclerotic pathologies, we found a small, but significant association of radiation with atherosclerosis.

Radiation risk of central nervous system tumors in the Life Span Study of atomic bomb survivors, 1958-2009

Radiation exposure is among the few factors known to be associated with risk of central nervous system (CNS) tumors. However, the patterns of radiation risk by histological type, sex or age are unclear. We evaluated radiation risks of first primary glioma, meningioma, schwannoma, and other or not otherwise specified (other/NOS) tumors in the Life Span Study cohort of atomic bomb survivors. Cases diagnosed between 1958 and 2009 were ascertained through population-based cancer registries in Hiroshima and Nagasaki. To estimate excess relative risk per Gy (ERR/Gy), we fit rate models using Poisson regression methods. There were 285 CNS tumors (67 gliomas, 107 meningiomas, 49 schwannomas, and 64 other/NOS tumors) among 105,444 individuals with radiation dose estimates to the brain contributing 3.1 million person-years of observation. Based on a simple linear model without effect modification, ERR/Gy was 1.67 (95% confidence interval, CI: 0.12 to 5.26) for glioma, 1.82 (95% CI: 0.51 to 4.30) for meningioma, 1.45 (95% CI: - 0.01 to 4.97) for schwannoma, and 1.40 (95% CI: 0.61 to 2.57) for all CNS tumors as a group. For each tumor type, the dose-response was consistent with linearity and appeared to be stronger among males than among females, particularly for meningioma (P = 0.045). There was also evidence that the ERR/Gy for schwannoma decreased with attained age (P = 0.002). More than 60 years after the bombings, radiation risks for CNS tumors continue to be elevated. Further follow-up is necessary to characterize the lifetime risks of specific CNS tumors following radiation exposure.

Misclassification of primary liver cancer in the Life Span Study of atomic bomb survivors

Primary liver cancer is difficult to diagnose accurately at death, due to metastases from nearby organs and to concomitant diseases, such as chronic hepatitis and cirrhosis. Trends in diagnostic accuracy could affect radiation risk estimates for incident liver cancer by altering background rates or by impacting risk modification by sex and age. We quantified the potential impact of death-certificate inaccuracies on radiation risk estimates for liver cancer in the Life Span Study of atomic bomb survivors. True-positive and false-negative rates were obtained from a previous study that compared death-certificate causes of death with those based on pathological review, from 1958 to 1987. We assumed various scenarios for misclassification rates after 1987. We obtained estimated true positives and estimated false negatives by stratified sampling from binomial distributions with probabilities given by the true-positive and false-negative rates, respectively. Poisson regression methods were applied to highly stratified person-year tables of corrected case counts and accrued person years. During the study period (1958-2009), there were 1,885 cases of liver cancer, which included 383 death-certificate-only (DCO) cases; 1,283 cases with chronic liver disease as the underlying cause of death; and 150 DCO cases of pancreatic cancer among 105,444 study participants. Across the range of scenarios considered, radiation risk estimates based on corrected case counts were attenuated, on average, by 13-30%. Our results indicated that radiation risk estimates for liver cancer were potentially sensitive to death-certificate inaccuracies. Additional data are needed to inform misclassification rates in recent years.

Assessing the Relative Biological Effectiveness of Neutrons across Organs of Varying Depth among the Atomic Bomb Survivors

When assessing radiation-related risk among the atomic bomb survivors, choices in modeling approach can have an important impact on the results, which are then used to inform radiation protection standards throughout the world. The atomic bombings of Hiroshima and Nagasaki produced a mixed-field radiation exposure from two sources: neutrons and gamma rays. Neutrons are more densely ionizing and cause greater biological damage per unit absorbed dose, resulting in greater relative biological effectiveness (RBE) than gamma rays. To account for this, a combined weighted dose is typically calculated as the sum of the gamma-ray dose and 10 times the neutron dose in the Radiation Effects Research Foundation's reports of mortality, solid cancer incidence and other outcomes. In addition, the colon, which is often chosen as the whole-body representative organ in these analyses, is relatively deep in the body and therefore its dose calculation involves heavy body shielding of neutrons and a low neutron/gamma-ray ratio. With added follow-up and recently updated doses, we used a data-driven approach to determine the best-fitting neutron RBE for a range of organs of varying depth. Aggregated person-year tables of solid cancer incidence (1958-2009) from the Life Span Study were created with separate neutron and gamma-ray DS02R1 doses for several organs including breast, brain, thyroid, bone marrow, lung, liver and colon. Typical excess relative risk models estimating the linear effect of radiation dose were fitted using a range of neutron weights (1-250) to calculate combined dose for each organ, and model deviances were compared to assess fit. Furthermore, models using separate terms for gamma-ray and neutron dose were also examined, wherein the ratio of the neutron/gamma-ray linear terms indicated the best estimate of the RBE. The best-fitting RBE value for the traditional weighted colon dose was 80 [95% confidence interval (CI): 20-190], while the RBEs for other organs using weighted doses ranged from 25 to 60, with the best-fitting weights and confidence interval widths both incrementally increasing with greater depth of organ. Models using separate neutron- and gamma-ray-dose terms gave similar results to weighted linear combinations, with a neutron/gamma-ray term ratio of 79.9 (95% CI: 18.8-192.3) for colon. These results indicated that the traditionally modeled RBE of 10 may underestimate the effect of neutrons across the full dose range, although these updated estimates still have fairly wide confidence bounds. Furthermore, the colon is among the deepest of organs and may not be the best choice as a single surrogate organ dose, as it may minimize the role of the neutrons. Future work with more refined organ doses could shed more light on RBE-related information available in the Life Span Study data.

Effect of Heterogeneity in Background Incidence on Inference about the Solid-Cancer Radiation Dose Response in Atomic Bomb Survivors

A recent analysis of solid cancer incidence in the Life Span Study of atomic bomb survivors (Hiroshima and Nagasaki, Japan) found evidence of a nonlinear, upwardly curving radiation dose response among males but not among females. Further analysis of this new and unexpected finding was necessary. We used two approaches to investigate this finding. In one approach, we excluded individual cancer sites or groups of sites from all solid cancers. In the other approach, we used joint analysis to allow for heterogeneity in background-rate parameters across groups of cancers with dissimilar trends in background rates. Exclusion of a few sites led to the disappearance of curvature among males in the remaining collection of solid cancers; some of these influential sites have unique features in their background age-specific incidence that are not captured by a background-rate model fit to all solid cancers combined. Exclusion of a few sites also led to an appearance of curvature among females. Misspecification of background rates can cause bias in inference about the shape of the dose response, so heterogeneity of background rates might explain at least part of the all solid cancer dose-response difference in curvature between males and females. We conclude that analysis based on all solid cancers as a single outcome is not the optimal method to assess radiation risk for solid cancer in the Life Span Study; joint analysis with suitable choices of cancer groups might be preferable by allowing for background-rate heterogeneity across sites while providing greater power to assess radiation risk than analyses of individual sites.

Radiation and Risk of Liver, Biliary Tract, and Pancreatic Cancers among Atomic Bomb Survivors in Hiroshima and Nagasaki: 1958-2009

The Life Span Study (LSS) of atomic bomb survivors has consistently demonstrated significant excess radiation-related risks of liver cancer since the first cancer incidence report. Here, we present updated information on radiation risks of liver, biliary tract and pancreatic cancers based on 11 additional years of follow-up since the last report, from 1958 to 2009. The current analyses used improved individual radiation doses and accounted for the effects of alcohol consumption, smoking and body mass index. The study participants included 105,444 LSS participants with known individual radiation dose and no known history of cancer at the start of follow-up. Cases were the first primary incident cancers of the liver (including intrahepatic bile duct), biliary tract (gallbladder and other and unspecified parts of biliary tract) or pancreas identified through linkage with population-based cancer registries in Hiroshima and Nagasaki. Poisson regression methods were used to estimate excess relative risks (ERRs) and excess absolute risks (EARs) associated with DS02R1 doses for liver (liver and biliary tract cancers) or pancreas (pancreatic cancer). We identified 2,016 incident liver cancer cases during the follow-up period. Radiation dose was significantly associated with liver cancer risk (ERR per Gy: 0.53, 95% CI: 0.23 to 0.89; EAR per 10,000 person-year Gy: 5.32, 95% CI: 2.49 to 8.51). There was no evidence for curvature in the radiation dose response (P=0.344). ERRs by age-at-exposure categories were significantly increased among those who were exposed at 0-9, 10-19 and 20-29 years, but not significantly increased after age 30 years, although there was no statistical evidence of heterogeneity in these ERRs (P = 0.378). The radiation ERRs were not affected by adjustment for smoking, alcohol consumption or body mass index. As in previously reported studies, radiation dose was not associated with risk of biliary tract cancer (ERR per Gy: -0.02, 95% CI: -0.25 to 0.30). Radiation dose was associated with a nonsignificant increase in pancreatic cancer risk (ERR per Gy: 0.38, 95% CI: <0 to 0.83). The increased risk was statistically significant among women (ERR per Gy: 0.70, 95% CI: 0.12 to 1.45), but not among men.

Integrating the evidence from evidence factors in observational studies

A sensitivity analysis for an observational study assesses how much bias, due to nonrandom assignment of treatment, would be necessary to change the conclusions of an analysis that assumes treatment assignment was effectively random. The evidence for a treatment effect can be strengthened if two different analyses, which could be affected by different types of biases, are both somewhat insensitive to bias. The finding from the observational study is then said to be replicated. Evidence factors allow for two independent analyses to be constructed from the same dataset. When combining the evidence factors, the Type I error rate must be controlled to obtain valid inference. A powerful method is developed for controlling the familywise error rate for sensitivity analyses with evidence factors. It is shown that the Bahadur efficiency of sensitivity analysis for the combined evidence is greater than for either evidence factor alone. The proposed methods are illustrated through a study of the effect of radiation exposure on the risk of cancer. An R package, evidenceFactors, is available from CRAN to implement the methods of the paper.

Radiation Risks of Uterine Cancer in Atomic Bomb Survivors: 1958-2009

BACKGROUND

Ionizing radiation is known to be capable of causing cancer of many organs, but its relationship with uterine cancer has not been well characterized.

METHODS

We studied incidence of uterine cancer during 1958-2009 among 62 534 female atomic bomb survivors. Using Poisson regression analysis, we fitted excess relative risk (ERR) models to uterine cancer rates adjusted for several lifestyle and reproductive factors. Person-years at risk were also adjusted for the probability of prior hysterectomy, because it could affect the subsequent risk of uterine cancer. We assessed the modifying effect of age and other factors on the radiation risk. For analysis of the modifying effect of age at radiation exposure around menarche, we compared the radiation risk for several exposure-age categories as well as using parametric models.

RESULTS

There were 224 uterine corpus cancers and 982 cervical cancers. We found a significant association between radiation dose and risk of corpus cancer (ERR per Gray [ERR/Gy] = 0.73, 95% confidence interval [CI] = 0.03 to 1.87) but not for cervical cancer (ERR/Gy = 0.00, 95% CI = -0.22 to 0.31). For corpus cancer, we found statistically significant heterogeneity in ERR/Gy by age (P heterogeneity  = .001) with elevated risk for women exposed to radiation between ages 11 and 15 years (ERR/Gy = 4.10, 95% CI = 1.47 to 8.42) and no indication of a radiation effect for exposures before or after this exposure-age range.

CONCLUSIONS

The current data suggest that uterine corpus is especially sensitive to the carcinogenic effect of radiation exposure occurring during the mid-pubertal period preceding menarche. There is little evidence for a radiation effect on cervical cancer risk.

Population Density in Hiroshima and Nagasaki Before the Bombings in 1945: Its Measurement and Impact on Radiation Risk Estimates in the Life Span Study of Atomic Bomb Survivors

In the Life Span Study cohort of atomic bomb survivors, differences in urbanicity between high-dose and low-dose survivors could confound the association between radiation dose and adverse outcomes. We obtained data on the population distribution in Hiroshima and Nagasaki before the 1945 bombings and quantified the impact of adjustment for population density on radiation risk estimates for mortality (1950-2003) and incident solid cancer (1958-2009). Population density ranged from 4,671 to 14,378 people/km2 in the urban region of Hiroshima and 5,748 to 19,149 people/km2 in the urban region of Nagasaki. Radiation risk estimates for solid cancer mortality were attenuated by 5.1% after adjustment for population density, but those for all-cause mortality and incident solid cancer were unchanged. There was no overall association between population density and adverse outcomes, but there was evidence that the association between density and mortality differed according to age at exposure. Among survivors who were 10-14 years of age in 1945, there was a positive association between population density and risk of all-cause mortality (per 5,000-people/km2 increase, relative risk = 1.053, 95% confidence interval: 1.027, 1.079) and solid cancer mortality (per 5,000-people/km2 increase, relative risk = 1.069, 95% confidence interval: 1.025, 1.115). Our results suggest that radiation risk estimates from the Life Span Study are not sensitive to unmeasured confounding by urban-rural differences.

Incidence of Breast Cancer in the Life Span Study of Atomic Bomb Survivors: 1958-2009

The importance of reproductive history in breast tissue development and etiology of sporadic breast cancer in females is well established. However, there is limited evidence of factors, other than age, that modify risk of radiation-related breast cancer. In this study, we evaluated breast cancer incidence in the Life Span Study cohort of atomic bomb survivors, adding 11 years of follow-up and incorporating reproductive history data. We used Poisson regression models to describe radiation risks and modifying effects of age and reproductive factors. Among 62,534 females, we identified 1,470 breast cancers between 1958 and 2009. Of 397 new cases diagnosed since 1998, 75% were exposed before age 20. We found a strong linear dose response with excess relative risk (ERR) of 1.12 per Gy [95% confidence interval (CI): 0.73 to 1.59] for females at age 70 after exposure at age 30. The ERR decreased with increasing attained age ( P = 0.007) while excess absolute rate (EAR) increased with attained age up to age 70 ( P < 0.001). Age at menarche was a strong modifier of the radiation effect: for a given dose, both the ERR and EAR decreased with increasing age at menarche ( P = 0.007 and P < 0.001). Also, independently, age-at-exposure effects on ERR and EAR differed before and after menarche ( P = 0.043 and P = 0.015, respectively, relative to log-linear trends), with highest risks for exposures around menarche. Despite the small number of male breast cancers (n = 10), the data continue to suggest a dose response (ERR per Gy = 5.7; 95% CI: 0.3 to 30.8; P = 0.018). Persistently increased risk of female breast cancer after radiation exposure and its modification pattern suggests heightened breast sensitivity during puberty.