The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (131I)–exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood 131I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose–dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

The U.S. National Cancer Institute (NCI), in cooperation with the Ministries of Health of Belarus and of Ukraine, is involved in epidemiological studies of thyroid diseases presumably related to the Chornobyl accident, which occurred in Ukraine on 26 April 1986. Within the framework of these studies, individual thyroid absorbed doses, as well as uncertainties, have been estimated for all members of the cohorts (13,215 Ukrainians and 11,918 Belarusians), who were selected from the large group of children aged 0 to 18 whose thyroids were monitored for gamma radiation within a few weeks after the accident. Information on the residence history and dietary habits of each cohort member was obtained during personal interviews. The methodology used to estimate the thyroid absorbed doses resulting from intakes of (131)I by the Ukrainian cohort subjects is described. The model of thyroid dose estimation is run in two modes: deterministic and stochastic. In the stochastic mode, the model is run 1,000 times for each subject using a Monte Carlo procedure. The geometric means of the individual thyroid absorbed doses obtained in the stochastic mode range from 0.0006 to 42 Gy. The arithmetic and geometric means of these individual thyroid absorbed doses over the entire cohort are 0.68 and 0.23 Gy, respectively. On average, the individual thyroid dose estimates obtained in the deterministic mode are about the same as the geometric mean doses obtained in the stochastic mode, while the arithmetic mean thyroid absorbed doses obtained in the stochastic mode are about 20% higher than those obtained in the deterministic mode. The distributions of the 1000 values of the individual thyroid absorbed dose estimates are found to be approximately lognormal, with geometric standard deviations ranging from 1.6 to 5.0 for most cohort subjects. For the time being, only the thyroid doses resulting from intakes of (131)I have been estimated for all subjects. Future work will include the estimation of the contributions to the thyroid doses resulting from external irradiation and from intakes of short-lived ((133)I and (132)Te) and long-lived ((134)Cs and (137)Cs) radionuclides, as well as efforts to reduce the uncertainties.

In collaboration with the Ukrainian Research Center for Radiation Medicine, the U.S. National Cancer Institute initiated a cohort study of children and adolescents exposed to Chornobyl fallout in Ukraine to better understand the long-term health effects of exposure to radioactive iodines. All 13,204 cohort members were subjected to at least one direct thyroid measurement between 30 April and 30 June 1986 and resided at the time of the accident in the northern part of Kyiv, Zhytomyr, or Chernihiv Oblasts, which were the most contaminated territories of Ukraine as a result of radioactive fallout from the Chornobyl accident. Thyroid doses for the cohort members, which had been estimated following the first round of interviews, were re-evaluated following the second round of interviews. The revised thyroid doses range from 0.35 mGy to 42 Gy, with 95 percent of the doses between 1 mGy and 4.2 Gy, an arithmetic mean of 0.65 Gy, and a geometric mean of 0.19 Gy. These means are 70% of the previous estimates, mainly because of the use of country-specific thyroid masses. Many of the individual thyroid dose estimates show substantial differences because of the use of an improved questionnaire for the second round of interviews. Limitations of the current set of thyroid dose estimates are discussed. For the epidemiologic study, the most notable improvement is a revised assessment of the uncertainties, as shared and unshared uncertainties in the parameter values were considered in the calculation of the 1,000 stochastic estimates of thyroid dose for each cohort member. This procedure makes it possible to perform a more realistic risk analysis.

On April 26, 1986, the worst nuclear reactor accident to date occurred at the Chornobyl (Chernobyl) power plant in Ukraine. Millions of people in Ukraine, Belarus and Russia were exposed to radioactive nuclides, especially (131)I. Since then, research has been conducted on various subgroups of the exposed population, and it has been demonstrated that the large increase in thyroid cancer is related to the (131)I exposure. However, because of study limitations, quantified risk estimates are limited, and there remains a need for additional information. We conducted an ecological study to investigate the relationship between (131)I thyroid dose and the diagnosis of thyroid cancer in three highly contaminated oblasts in Northern Ukraine. The study population is comprised of 301,907 persons who were between the ages of 1 and 18 at the time of the Chornobyl accident and were living in 1,293 rural settlements in the three study oblasts. Twenty-four percent of the study population had individual thyroid dose estimates and the other 76% had "individualized" estimates of thyroid dose based on direct thyroid measurements taken from a person of the same age and gender living in the same or nearby settlement. Cases include 232 thyroid cancers diagnosed from January 1990 through December 2001, and all were confirmed histologically. Dose-response analyses took into account differences in the rate of ultrasound examinations conducted in the three study oblasts. The estimated excess relative risk per gray was 8.0 (95% CI = 4.6-15) and the excess absolute risk per 10,000 person-year gray was estimated to be 1.5 (95% CI = 1.2-1.9). In broad terms, these estimates are compatible with results of other studies from the contaminated areas, as well as studies of external radiation exposure.

The 1986 accident at the Chernobyl nuclear power plant remains the most serious nuclear accident in history, and excess thyroid cancers, particularly among those exposed to releases of iodine-131 remain the best-documented sequelae. Failure to take dose-measurement error into account can lead to bias in assessments of dose-response slope. Although risks in the Ukrainian-US thyroid screening study have been previously evaluated, errors in dose assessments have not been addressed hitherto. Dose-response patterns were examined in a thyroid screening prevalence cohort of 13,127 persons aged <18 at the time of the accident who were resident in the most radioactively contaminated regions of Ukraine. We extended earlier analyses in this cohort by adjusting for dose error in the recently developed TD-10 dosimetry. Three methods of statistical correction, via two types of regression calibration, and Monte Carlo maximum-likelihood, were applied to the doses that can be derived from the ratio of thyroid activity to thyroid mass. The two components that make up this ratio have different types of error, Berkson error for thyroid mass and classical error for thyroid activity. The first regression-calibration method yielded estimates of excess odds ratio of 5.78 Gy−1 (95% CI 1.92, 27.04), about 7% higher than estimates unadjusted for dose error. The second regression-calibration method gave an excess odds ratio of 4.78 Gy−1 (95% CI 1.64, 19.69), about 11% lower than unadjusted analysis. The Monte Carlo maximum-likelihood method produced an excess odds ratio of 4.93 Gy−1 (95% CI 1.67, 19.90), about 8% lower than unadjusted analysis. There are borderline-significant (p = 0.101–0.112) indications of downward curvature in the dose response, allowing for which nearly doubled the low-dose linear coefficient. In conclusion, dose-error adjustment has comparatively modest effects on regression parameters, a consequence of the relatively small errors, of a mixture of Berkson and classical form, associated with thyroid dose assessment.