Оценка радиационных рисков заболеваемости и смертности по основным классам неонкологической заболеваемости среди участников ликвидации последствий аварии на Чернобыльской АЭС

К концу 1990-х годов исследования японской когорты жителей, переживших атомные бомбардировки городов Хиросима и Нагасаки в 1945 г. [1], показали, что основным отдалённым эффектом действия ионизирующей радиации в малых дозах на здоровье человека является увеличение частоты злокачественных новообразований (ЗНО) среди облучённых лиц. Однако попытки оценки радиационных эффектов для здоровья японского населения, касающиеся неонкологических (нераковых) болезней, начались ещё в 1964 г., когда было продемонстрировано увеличение смертности по причине туберкулёза среди мужчин, находившихся в радиусе до 2 км от эпицентра взрыва, по сравнению с мужчинами, находившимися в 2-3 км от эпицентра [2]. Исследование

отдалённых радиационных эффектов для нераковых заболеваний было также представлено в 1973 г. Belsky и соавт. [3], где статистически значимое влияние радиации было обнаружено для заболеваний щитовидной железы и глаза, а также для гипотензии в группе наиболее облучённых в молодом возрасте лиц.

Исследование 1992 г. Shimizu и соавт. [4] в японской когорте показало существование статистически значимой ассоциации между дозой облучения и показателем смертности от всех, за исключением ЗНО, причин. Впоследствии в японской когорте были обнаружены радиационные риски для болезней щитовидной железы, печени, сердечно-сосудистых и цереброваскулярных болезней [5].

Международная комиссия по радиологической защите (МКРЗ) в 2007 г. [6] отмечала «потенциальную важность наблюдений за частотой нераковых заболеваний» и констатировала недостаточность доказательств существования их радиационных рисков при дозах облучения менее 0,2 Гр.

В когорте российских участников ликвидации последствий аварии на Чернобыльской АЭС (ликвидаторов), зарегистрированных и наблюдающихся с 1986 г. в системе Национального радиационно-эпидемиологического регистра (НРЭР) [7], дозы облучения не превышали 1,5 Гр, а средняя доза в когорте составила около 0,1 Гр [8]. Численности российских когорт ликвидаторов для анализа радиационных рисков неонкологических заболеваний (86,5-89,5 тыс. человек) практически совпадают с численностью японской когорты, использовавшейся для оценки таких рисков среди лиц, переживших атомные бомбардировки 1945 г. (86,5 тыс. человек) [9].

Статистически значимый радиационный риск цереброваскулярных заболеваний (ЦВЗ), ишемической болезни сердца и эссенциальной гипертензии был обнаружен в российской когорте ликвидаторов уже в первые 20 лет после аварии на Чернобыльской АЭС [10, 11], т.е. до 2006 г. Ещё через 5 лет по наблюдениям НРЭР были получены оценки радиационный рисков смертности от болезней системы кровообращения (БСК) в целом и от ЦВЗ в частности [12, 13]. При этом коэффициенты избыточных относительных радиационных рисков ERR/Зв (от англ. Excess Relative Risk) для ЦВЗ, ишемической болезни сердца и БСК в российских облучённых когортах оценивались в диапазоне от 0,1 до 1,5 и оказались выше, чем в облучённых когортах других стран [14]. В период наблюдения 1986-2012 гг. статистически значимые радиационные риски наблюдались для ЦВЗ [15], сердечно-сосудистых заболеваний [16] и гипертензий [17].

Систематическое исследование радиационных рисков по всем основным группам зарегистрированных диагнозов неонкологических заболеваний ликвидаторов (закодированных в соответствии с Международной классификацией болезней (МКБ-10) [18, 19]) за период наблюдений с 1986 по 2019 гг. [20] показало наличие радиационных рисков заболеваемости для 9 трёхзначных рубрик диагнозов МКБ-10 [18, 19]): болезни крови (рубрики МКБ-10: D50-D89), эндокринной системы (E00-E90), психические расстройства (F00-F99), болезни нервной системы (G00-G99), БСК (I00-I99), органов дыхания (J00-J99), органов пищеварения (K00-K93), костно-мышечной системы (M00-M99) и мочеполовой системы (N00-N99). При этом радиационные риски смертности с высокой степенью достоверности (p<0,001) идентифицировались только для БСК (I00-I99), и на грани достоверности (p=0,011) - для болезней органов пищеварения (K00-K93).

Количественно оценки коэффициентов радиационных рисков на единицу дозы находились на уровне аналогичных оценок для солидных ЗНО [21], поэтому радиационные риски неонкологических заболеваний вносят существенный вклад в общий радиационный ущерб здоровью человека.

Представленное в данной статье исследование обновляет результаты оценок радиационных рисков неонкологических заболеваний ликвидаторов, полученных за период наблюдений НРЭР 1988-2019 гг. [20], с учётом новых данных, накопленных к концу 2024 г.

Материалы и методы

В Единой федеральной базе данных (ЕФБД) НРЭР диагнозы заболеваний и причин смерти кодируются в соответствии с МКБ-10 [18, 19]. В данной работе исследуются избыточные радиационные риски заболеваемости и смертности, сгруппированные по трёхзначным рубрикам диагнозов МКБ-10.

Начало наблюдения для каждого члена когорты определяется датой его въезда в чернобыльскую зону. Дата регистрации у наблюдаемого члена когорты первого диагноза из выбранной трёхзначной рубрики считается датой диагноза случая заболевания или смерти, относящейся к этой рубрике. Время нахождения ликвидатора под риском заболеть исследуемой болезнью или умереть от исследуемой болезни рассчитывается от даты начала наблюдения до наименьшей из следующих дат: даты регистрации случая (заболевания или смерти из выбранной рубрики болезней) или даты выбытия из-под наблюдения по причинам, не связанным с исследуемой рубрикой болезней.

Для оценки радиационного риска заболеваемости, относящейся к выбранной рубрике диагнозов, когорта формировалась из лиц, не имевших диагнозов из этой рубрики до начала наблюдения. Поэтому в случае исследования заболеваемости для разных рубрик болезней численность когорт различается. Для оценки радиационных рисков смертности с изменением исследуемых причин смерти численность когорты не меняется.

Для исследования радиационных рисков неонкологических заболеваний на основе ЕФБД НРЭР были сформированы ретроспективные когорты ликвидаторов мужского пола с зарегистрированными индивидуальными дозами внешнего гамма-облучения всего тела, въехавших в зону чернобыльских работ в период с 26.04.1986 г. по 31.12.1987 г. в возрасте от 18 до 70 лет и наблюдавшихся в период с 1986 по 2024 гг. Средняя доза облучения в исследованных когортах составила 0,132 Гр, средний возраст при въезде в чернобыльскую зону – 33,6 года, а средний достигнутый возраст – от 39,6 до 46,9 года.

В табл. 1 приведены основные характеристики когорт ликвидаторов, сформированных для анализа радиационных рисков заболеваемости. Наименьшую численность (86 532 человека) имеет когорта для исследования болезней костно-мышечной системы и соединительной ткани (МКБ-10: M00-M99), а наибольшую численность (89 536 человек) – когорты для исследования новообразований in situ (D00-D09) и новообразований неопределённого характера (D37-D48). Эти же когорты имеют минимальный (39,6 лет) и максимальный (46,9 лет) средний достигнутый возраст за время наблюдения.

Численность когорты для анализа радиационных рисков смертности от неонкологических болезней составила 89 540 человек, число человеко-лет под наблюдением – 2 445 031,5, средний возраст при въезде в чернобыльскую зону – 33,6 года, а средний достигнутый возраст за время наблюдения – 48,5 года.

Для исследования радиационных рисков был использован когортный метод анализа [22] и соответствующие численные методы расчёта, реализованные в специализированном программном обеспечении радиационно-эпидемиологического анализа EPICURE [23].

Для проведения необходимых расчётов с помощью пакета EPICURE [23] данные о ликвидаторах, вошедших в когорту, были стратифицированы следующим образом: 10 страт по возрасту ликвидаторов на момент въезда в чернобыльскую зону (интервалы в годах: 18-24, 25-29, 30-34, 35-39, 40-44, 45-49, 50-54, 55-59, 60-64, 65-70); 38 страт по календарному году наблюдения, по одному году с 1986 по 2024 гг.; 3 страты по месту регистрации ликвидаторов (78 субъектов РФ были разделены на 3 страты); а также 16 страт по дозе облучения со следующими верхними разделительными границами между стратами (не включая эту границу) в Гр: (0,005; 0,020; 0,030; 0,040; 0,050; 0,070; 0,090; 0,100; 0,110; 0,125; 0,150; 0,175; 0,200; 0,225; 0,250; 1,500).

Таблица 1

Основные характеристики когорт ликвидаторов, сформированных для анализа радиационных рисков заболеваемости

Диагнозы болезней (рубрики МКБ-10)	Численность когорты	Число человеко-лет под наблюдением	Средний достигнутый возраст, лет
Некоторые инфекционные и паразитарные болезни (A00-B99)	89 279	2 289 991,5	46,07
Новообразования in situ (D00-D09)	89 536	2 443 729,5	46,88
Доброкачественные новообразования (D10-D36)	89 448	2 306 302,5	46,14
Новообразования неопределённого или неизвестного характера (D37-D48)	89 536	2 436 058,5	46,84
Болезни крови, кроветворных органов и отдельные нарушения,	89 523	2 372 843,5	46,49
вовлекающие иммунный механизм (D50-D89)
Болезни эндокринной системы, расстройства питания и нарушения обмена	88 948	1 564 116,5	42,08
веществ (E00-E90)
Психические расстройства и расстройства поведения (F00-F99)	88 641	1 612 985,0	42,36
Болезни нервной системы (G00-G99)	88 319	1 230 431,0	40,25
Болезни глаза и его придаточного аппарата (H00-H59)	88 453	1 721 492,0	43,01
Катаракты (H25-H26)	89 500	2 285 422,0	46,03
Болезни уха и сосцевидного отростка (H60-H95)	89 299	2 171 189,5	45,41
Болезни системы кровообращения (I00-I99)	87 467	1 199 147,5	40,07
Болезни органов дыхания (J00-J99)	87 381	1 183 732,0	40,04
Болезни органов пищеварения (K00-K93)	86 594	1 258 191,0	40,50
Болезни кожи и подкожной клетчатки (L00-L99)	89 288	2 193 637,5	45,54
Болезни костно-мышечной системы и соединительной ткани (M00-M99)	86 532	1 109 064,0	39,60
Болезни мочеполовой системы (N00-N99)	89 075	1 909 915,5	44,00
Травмы, отравления и некоторые другие воздействия внешних причин	89 216	2 018 945,0	44,57
(S00-T98)
Внешние причины заболеваемости и смертности (V00-Y98)	89 530	2 435 346,5	46,84

В использованной когортной модели риска наблюдаемый показатель заболеваемости (риск) А представлен суммой фонового (в отсутствие исследуемых факторов риска) показателя и показателя, связанного с исследуемыми факторами риска:

А = Ао + Ar .(1)

В общем случае избыточный риск A r представляется двумя моделями: аддитивной и мультипликативной. В аддитивной модели наблюдаемый риск представляется в следующем виде:

А = Ао + EAR,(2)

где EAR – избыточный абсолютный риск (от англ. Excess Absolute Risk).

В мультипликативной модели – в виде:

А = Аох RR = Аох (1 + ERR),(3)

где безразмерная величина ERR есть избыточный относительный риск.

В данном исследовании аддитивные модели риска оказались статистически не значимыми, и при оценке зависимости заболеваемости и смертности ликвидаторов от дозы внешнего облучения использовалась Пуассоновская регрессия [20] с моделью избыточного относительного риска ERR :

A nc ( e,r,c,d) = A o_inc ( e,r,c ) х RR ( d ) = A o_inc ( e,r,c ) х [ 1 + ERR ( d)] ,                      (4)

A mrT ( e,r,c,d) = A o_mrt ( e,r,c) x RR ( d) = A o_mrt ( e,r,c ) х [ 1 + ERR ( d)] ,                   (5)

где e – категориальная переменная, обозначающая возраст при облучении, в годах; r – категориальная переменная, обозначающая место регистрации ликвидатора; с – категориальная переменная, обозначающая календарный год наблюдения; d – категориальная переменная, обозначающая страту (группу) по дозе облучения ликвидаторов; RR(d)=1 + ERR(d) - относительный риск в дозовой группе d; ERR(d) - избыточный относительный риск в дозовой группе d; AiNc(e,r,c,d) -наблюдаемая в когорте заболеваемость по группам ликвидаторов, определяемым значениями категориальных переменных (е,r,с,d); ЛмRт(e,r,c,d) - наблюдаемая в когорте смертность по группам ликвидаторов (е,r,с,d); Лo_lNc(e,r,c) - фоновая (в отсутствие облучения) заболеваемость в когорте по группам ликвидаторов (е,r,с); Ao_MRT(e,r,c) - фоновая (в отсутствие облучения) смертность в когорте по группам ликвидаторов (е, r, с).

В данном исследовании получены оценки параметров линейной беспороговой (ЛБП) зависимости ERR от дозы облучения D :

ERR(D) = fl* D ,                                                                  (6)

где в - оцениваемый параметр, показатель избыточного относительного риска на единицу дозы, ERR /Гр; D - средняя доза облучения в группе ликвидаторов, определяемой значениями категориальных переменных ( е , r , с , d ) , в Гр.

Идентификация моделей и оценки коэффициента ERR /Гр проводились методом максимального правдоподобия, где логарифмическая функция правдоподобия l имела следующий вид [20]: l[l₀(e,r,c\p] = VU(mrln(^                                       (7)

где N =10 x 38 x 3 x 16=18240 - число ячеек в таблице для сгруппированных данных; m i - число случаев заболеваний или случаев смерти в i -ой ячейке; PY i - суммарное время под риском в i -ой ячейке, в человеко-годах, остальные обозначения соответствуют обозначениям в формулах (4-6).

Для исследуемой группы болезней суммарное время под риском PY i для каждой i -ой ячейки данных, сгруппированных по переменным ( е , r , с , d ), определяется как сумма индивидуальных времён наблюдения за членами когорты в этой ячейке, а средняя доза - как средняя индивидуальная доза облучения ликвидаторов, попавших в ячейку. Поскольку дата выхода из-под наблюдения регистрируется с точностью до целого года, то год выхода из-под наблюдения для индивидуума учитывается в общей сумме времён под риском с весом 0,5.

Непараметрические оценки относительного риска (RR, от англ. Relative Risk) заболеваемости и смертности получены для трёх групп ликвидаторов, различающихся дозами облучения: для ликвидаторов с дозами 0,050-0,100 Гр, 0,100-0,150 Гр и 0,150-1,5 Гр. В качестве контрольной группы (для которой RR=1 по определению) была принята группа ликвидаторов с дозами менее 0,050 Гр. Модели относительного риска имели следующий вид:

^ iwc (e, r, c, п) = A_0JNC⁽e, r, c) * RR_n ,                                                   (8)

^■ MRT ⁽e,r, c,n) = A o_MRT ^(e,^r, ^cy * ^RR n ,                                                   ⁽⁹⁾

где A o_inc и A o_mrt - показатели контрольной группы ликвидаторов с дозами менее 0,030 Гр, которая принималась за необлучённую; RR n - относительный риск n -ой дозовой группы ликвидаторов по отношению к контрольной группе ( n =0 - для контрольной группы ликвидаторов с величиной RR o =1 по определению, n =1 - для дозовой группы 0,030-0,050 Гр, n =2 - для дозовой группы 0,050-0,100 Гр, n =3 - для дозовой группы 0,100-0,150 Гр, n =4 - для группы ликвидаторов с дозами 0,150 Гр и более).

Результаты и обсуждение

В табл. 2 для 33 исследованных групп диагнозов в 16 классах приведено число впервые зарегистрированных случаев заболеваний в когорте ликвидаторов в период 1986-2024 гг., оценки коэффициентов избыточного относительного риска на единицу дозы (ERR/Гр), 95% доверительные интервалы (ДИ) этих оценок и величина p, определяющая статистическую значимость найденной оценки ERR/Гр на основе статистики отношения правдоподобий идентифицированной ЛБП модели и модели «нулевой гипотезы» с отсутствием радиационного риска (ERR/Гр=0) [22]. Учитывая большое число зарегистрированных случаев заболеваний, статистическую значимость найденной оценки ERR/Гр можно определить величиной p<0,01.

Для 9 классов из 16 исследованных групп диагнозов статистические ассоциации избыточной заболеваемости с дозой облучения ликвидаторов являются статистически значимыми: 1) болезни крови, кроветворных органов и отдельные нарушения, вовлекающие иммунный механизм (D50-D89); 2) болезни эндокринной системы, расстройства питания и нарушения обмена веществ (E00-E90); 3) психические расстройства и расстройства поведения (F00-F99); 4) болезни нервной системы (G00-G99); 5) БСК (I00-I99); 6) болезни органов дыхания (J00-J99); 7) болезни органов пищеварения (K00-K93); 8) болезни костно-мышечной системы и соединительной ткани (M00-M99); 9) болезни мочеполовой системы (N00-N99).

Таблица 2 Оценки коэффициентов избыточного относительного риска (ERR/Гр) заболеваемости в когорте ликвидаторов для различных рубрик диагнозов болезней

Наибольшими коэффициентами радиационного риска характеризовались болезни крови, кроветворных органов и отдельные нарушения, вовлекающие иммунный механизм (D50-D89), а также БСК (I00-I99) с оценками ERR/Гр=0,52 и ERR/Гр=0,56 соответственно.

В классе болезней крови (D50-D89) коэффициент радиационного риска повышается с 0,52 до 0,86 при исключении из этого класса анемий, связанных с питанием (D50-D53), и гемолитических анемий (D55-D59). Здесь необходимо обратить особое внимание на эритробластопению и другие апластические анемии (D60-D61) с очень высоким коэффициентом ERR/Гр=6,85, который превышает коэффициент риска для лейкозов в этой же когорте ликвидаторов (для лейкозов, исключая хронический лимфолейкоз, ERR/Гр=4,41) [24].

В классе БСК высокими коэффициентами радиационных рисков, выше среднего по классу значения 0,56, выделяются ишемическая болезнь сердца (I20-I25) с оценкой ERR/Гр=0,67 и атеросклероз (I70) с оценкой ERR/Г р=0,71.

Стоит отметить, что некоторые инфекционные и паразитарные болезни (A00-B99) характеризуются отрицательной оценкой радиационного риска ERR/Гр=-0,29, находящейся на близком к статистической значимости уровне (p=0,014).

Несмотря на найденные здесь статистические связи заболеваемости с дозой облучения ликвидаторов, МКРЗ [25] считает, что для перечисленных выше болезней радиобиологические механизмы действия радиации в настоящее время не выяснены, хотя некоторые гипотезы являются предметом исследований.

Для таких состояний здоровья, как травмы, отравления и прочие внешние причины болезней (S00-T98, V00-Y98), какие-либо механизмы влияния облучения в малых дозах на их развитие предположить крайне сложно, что подтверждается в табл. 2 отсутствием статистических связей с дозой для этих диагнозов.

Для заболеваний глаза (H00-H59), включая катаракту (H25-H26), для болезней уха (H60-H95), для доброкачественных новообразований (D10-D36), а также для болезней кожи и подкожной клетчатки (L00-L99) радиационный риск не обнаружен, несмотря на большое число установленных диагнозов в этих рубриках за почти 40-летний период наблюдения (1986-2024 гг.).

В табл. 3 для 17 групп неонкологических диагнозов из 11 классов болезней и двух объединённых групп, которые являлись причинами смерти ликвидаторов, приведено число зарегистрированных случаев смерти в когорте ликвидаторов за период 1986-2024 гг., оценки коэффициентов избыточного относительного риска на единицу дозы (ERR/Гр), 95% ДИ этих оценок и величина p, определяющая статистическую значимость найденной оценки ERR/Гр [22]. Здесь, как и для радиационных рисков заболеваемости, статистическую значимость найденной оценки ERR/Г р можно определить величиной p<0,01.

Как видно из табл. 3, из 11 исследованных классов болезней, являвшихся причинами смертности ликвидаторов, только для двух классов, БСК (I00-I99) и болезней органов пищеварения (K00-K93), радиационные риски являлись статистически значимыми с оценками ERR/Г р=0,49 и 0,73 соответственно. В совокупности для всех причин смерти (A00-Z99) и для всех причин смерти, исключая ЗНО (A00-B99, D00-Z99), радиационные риски также были статистически значимыми (ERR/Гр=0,36 и 0,30 соответственно).

Таблица 3

Оценки коэффициентов избыточного относительного риска (ERR/Гр) смертности в когорте ликвидаторов для различных рубрик диагнозов причин смерти

Диагнозы болезней (рубрики МКБ-10)	Число случаев смерти	ERR/Гр (95% ДИ)	Величина p
Некоторые инфекционные и паразитарные болезни (A00-B99) Болезни крови, кроветворных органов и отдельные нарушения, вовлекающие иммунный механизм (D50-D89) Болезни эндокринной системы, расстройства питания и нарушения обмена веществ (E00-E90) Психические расстройства и расстройства поведения (F00-F99) Болезни нервной системы (G00-G99) Болезни системы кровообращения (I00-I99) Болезни, характеризующиеся повышенным кровяным давлением (I10-I15) Ишемическая болезнь сердца (I20-I25) Болезни органов дыхания (J00-J99) Болезни органов пищеварения (K00-K93) Болезни печени (K70-K76) Болезни поджелудочной железы (K85-86) Болезни мочеполовой системы (N00-N99) Травмы, отравления и некоторые другие воздействия внешних причин (S00-T98) Внешние причины смертности (V00-Y98) Все причины смерти, исключая ЗНО (A00-B99, D00-Z99) Все причины смерти (A00-Z99)	666 37 324 440 729 18 080 733 10 786 2 340 2 861 1 779 350 266 6 882 169 34 825 42 651	0,94 (-0,16; 2,42) 2,42 (н.н.; 19,01) 0,74 (-0,46.; 2,58) - 0,60 (н.н.; 0,66) - 0,03 (н.н.; 1,02) 0,49 (0,28; 0,71) 2,50 (1,09; 4,39) 0,42 (0,16; 0,70) 0,12 (-0,38; 0,70) 0,73 (0,20; 1,34) 0,98 (0,29; 1,82) 1,71 (0,04; 4,29) - 0,18 (н.н.; 1,79) -0,17 (0,45; 0,15) -0,71 (н.н.; 1,20) 0,30 (0,16; 0,45) 0,36 (0,23; 0,50)	0,103 0,340 0,268 0,306 > 0,5 < 0,001 < 0,001 0,001 > 0,5 0,005 0,004 0,044 > 0,5 0,293 0,384 < 0,001 < 0,001

Среди БСК (I00-I99) наибольший коэффициент радиационно обусловленной смертности (ERR/Гр=2,50), в 5 раз превышающий средний по классу, имели болезни, характеризующиеся повышенным кровяным давлением (I10-I15). В классе болезней органов пищеварения (K00-K93) коэффициенты радиационно обусловленной смертности выше среднего наблюдались для болезней печени (K70-K76) и болезней поджелудочной железы (K85-86) с оценками ERR/Гр=0,98 и ERR/Гр=1,71 соответственно.

Последняя строка табл. 3 показывает, что для всех причин смерти, включая ЗНО, радиационный риск увеличивается до значения ERR/Гр=0,36. Это связано с тем, что радиационный риск смертности от ЗНО в российской когорте ликвидаторов за период наблюдения 1992-2024 гг. является статистически значимым с величиной ERR/Гр=0,79 [21].

Учитывая, что оценки радиационных рисков БСК подтверждаются ранее проведёнными японскими исследованиями [26] и исследованиями многих других облучённых когорт по всему миру, включая когорты персонала атомной отрасли [14], существование радиационных рисков БСК (I00-I99) к настоящему времени является доказанным на эпидемиологическом уровне исследований. Кроме того, существование радиационных рисков болезней органов пищеварения (K00-K93) также подтверждается оценками радиационных рисков в японской когорте [26].

Оценки ERR/Гр для смертности от БСК (I00-I99) и от болезней органов пищеварения (K00-K93) в российской когорте ликвидаторов превышают соответствующие оценки в японской когорте [26] в несколько раз: 0,49 и 0,73 в российской когорте против ERR/Гр=0,11 для БСК и органов пищеварения в японской когорте. В отличие от японской когорты, в которой обнаружен высокий радиационный риск смертности от болезней органов дыхания (J00-J99) со значением коэффициента ERR/Гр=0,21 [26], в российской когорте ликвидаторов за весь период наблюдения 1986-2024 гг. такой риск не обнаружен (для оценки ERR/Гр p>0,5).

Несмотря на то, что в российской когорте ликвидаторов радиационный риск смертности от болезней крови (D50-D89) не выявлен из-за малого числа случаев смерти по этой причине (37 случаев), косвенным подтверждением оценки радиационного риска заболеваемости ликвидаторов болезнями крови (D50-D89) на уровне ERR/Г р=0,52 является высокий радиационный риск смертности от этих болезней в японской когорте (ERR/Г р=1,7) [26].

Следует отметить, что для смертности от инфекционных и паразитарных болезней (A00-B99) в японской когорте [26] получена статистически значимая отрицательная оценка ERR/Гр=-0,02, что косвенно подтверждается статистически значимой отрицательной оценкой ERR/Г р=-0,29 для заболеваемости этими болезнями в российской когорте ликвидаторов (табл. 2, первая строка).

Непараметрические оценки относительного риска RR для заболеваемости ликвидаторов БСК (I00-I99) были проведены по модели (8), а для смертности ликвидаторов от БСК - по модели (9), в четырёх дозовых интервалах: 0,030-0,050 Гр, 0,050-0,100 Гр, 0,100-0,150 Гр и более 0,150 Гр, относительно контрольной группы ликвидаторов с дозами облучения 0-0,030 Гр.

Статистически значимая оценка RR для заболеваемости БСК получена для ликвидаторов с дозами 0,100-0,150 Гр, а также 0,150 Гр и более: RR o,ioo-o,i5o =1,O4 при 95% ДИ (1,01; 1,08),

RR o,i5o+ =1,10 при 95% ДИ (1,08; 1,13). Статистически значимая оценка RR для смертности от БСК получена только для ликвидаторов с дозами 0,150 Гр и более: RR o,i5o+ =1,10 при 95% ДИ (1,05; 1,16).

Линейная зависимость RR от дозы облучения в рамках ЛБП модели, выражающейся формулами (4-5) и (6), и непараметрические оценки RR, рассчитанные по моделям (8) и (9), представлены на рис. 1 и 2. Линия, состоящая из длинных штрихов, соответствует среднему значению RR(d)=1+ERR(d); линии, состоящие из коротких штрихов, соответствуют нижней и верхней границе 95% ДИ для оценки RR, определяемой по оценкам ERR/Гр из табл. 2-3. Точки и вертикальные отрезки соответствуют оценкам RR и их 95% ДИ, рассчитанным по моделям (8) и (9).

Кроме того, вариабельность точечных непараметрических оценок RR в более узких дозовых интервалах представлена на рис. 1 и 2 небольшими окружностями: оценки RR относительно контрольной группы ликвидаторов с дозами менее 0,005 Гр получены для 15 последовательных непересекающихся дозовых интервалов с нижними границами 0,005; 0,020; 0,030; 0,040; 0,050; 0,070; 0,090; 0,100; 0,110; 0,125; 0,150; 0,175; 0,200; 0,225 и 0,250 Гр.

• 1,3    ...................................................................................................,..................................................................................................,..................................................................................................I..................................................................................................,..................................................................................................,..................................................................................................,

  ERR/Гр = 0,56

  1,2

  

  1,1

  1,0 «“

  ф

  о

  0,9

  
  

О

0,8  1.................................................................................................................................................................................................................................................................................................................................................................................................................................................................................

0,00             0,05             0,10             0,15             0,20             0,25             0,30

Доза, Гр

Рис. 1. Относительный радиационный риск (RR) заболеваемости ликвидаторов болезнями системы кровообращения (I00-I99) в рамках ЛБП модели (4, 6) и непараметрические оценки RR по дозовым интервалам в рамках модели (8).

1,3

0,8

0,00             0,05             0,10             0,15             0,20             0,25             0,30

Доза, Гр

Рис. 2. Относительный радиационный риск (RR) смертности ликвидаторов от болезней системы кровообращения (I00-I99) в рамках ЛБП модели (5, 6) и непараметрические оценки RR по дозовым интервалам в рамках модели (9).

Оценки RR для БСК, полученные в рамках ЛБП модели (4, 5), показывают очень хорошее согласие с непараметрическими оценками RR (8, 9) в рамках 95% ДИ этих оценок для доз облучения ликвидаторов 0,100-0,150 Гр и свыше 0,150 Гр в случае заболеваемости (рис. 1), а в случае смертности – для доз облучения ликвидаторов свыше 0,150 Гр (рис. 2).

На рис. 3 и 4 оценки коэффициентов избыточного относительного риска на единицу дозы (ERR/Гр) для заболеваемости БСК (рис. 3) и для смертности от БСК (рис. 4), а также 95% ДИ этих оценок представлены в ретроспекции для последовательно увеличивающихся периодов наблю- дения когорты, начиная с 1986 г.

ERR/Гр             НГ 95% ДИ ERR/Гр

ВГ 95% ДИ ERR/Гр

Рис. 3. Оценки коэффициентов избыточного относительного риска на единицу дозы (ERR/Гр) для заболеваемости БСК в когорте ликвидаторов, полученные за последовательно увеличивающиеся календарные периоды наблюдения когорты, начиная с 1986 г.

ERR/Гр         НГ 95% ДИ ERR/Гр          ВГ 95% ДИ ERR/Гр

Рис. 4. Оценки коэффициентов избыточного относительного риска на единицу дозы (ERR/Гр) для смертности от БСК в когорте ликвидаторов, полученные за последовательно увеличивающиеся календарные периоды наблюдения когорты, начиная с 1986 г.

Как видно из рис. 3 и 4, по мере накопления данных (с увеличением периода наблюдения) 95% ДИ оценки коэффициентов ERR/Гр сужаются, а нижняя граница этих интервалов устойчиво превышает нулевой уровень уже с 1987 г. для ERR заболеваемости и с 1992 г. – для ERR смертности. Следует отметить, что на протяжении всего периода накопления данных радиационные риски БСК обнаруживались устойчиво уже с 1987 г. (оценки коэффициентов ERR/Гр устойчиво превышали нулевой уровень).

Полученные в данном исследовании оценки коэффициентов радиационных рисков для неонкологических заболеваний ликвидаторов являются устойчивыми также и по отношению к неопределённостям в оценках доз облучения [27].

С другой стороны, проводившиеся ранее исследования влияния мультиморбидности состояния ликвидаторов на радиационные риски БСК [28], болезней органов пищеварения [29] и болезней органов дыхания [30] показали, что связанная с этим фактором вариабельность оценок коэффициентов ERR/Гр для неонкологических заболеваний может достигать нескольких раз, в отличие от практически отсутствующей (при учёте мультиморбидности) вариабельности оценок радиационных рисков ЗНО [31].

Поэтому практическое использование оценок коэффициентов радиационных рисков неонкологических заболеваний в задачах радиационной защиты требует более дифференцированных подходов к защите различных групп облучённых лиц, чем просто их разделение по полу и возрасту, как это сейчас принято при учёте стохастических эффектов в виде радиационно обусловленных ЗНО.

Заключение

Полученные результаты исследований радиационных рисков неонкологических болезней в российской когорте ликвидаторов аварии на Чернобыльской АЭС за период наблюдений 1986-2024 гг. позволяют сделать следующие выводы.

В исследованной российской когорте ликвидаторов наблюдались статистически значимые коэффициенты радиационных рисков (ERR/Гр) заболеваемости для следующих 9 из 16 исследованных классов диагнозов МКБ-10 [18, 19]: 1) болезни крови, кроветворных органов и отдельные нарушения, вовлекающие иммунный механизм (D50-D89); 2) болезни эндокринной системы, расстройства питания и нарушения обмена веществ (E00-E90); 3) психические расстройства и расстройства поведения (F00-F99); 4) болезни нервной системы (G00-G99); 5) болезни системы кровообращения (I00-I99); 6) болезни органов дыхания (J00-J99); 7) болезни органов пищеварения (K00-K93); 8) болезни костно-мышечной системы и соединительной ткани (M00-M99); 9) болезни мочеполовой системы (N00-N99). Наибольшими коэффициентами радиационных рисков характеризовались болезни крови (D50-D89) с оценкой ERR/Г р=0,52 и болезни системы кровообращения (I00-I99) с оценкой ERR/Гр=0,56.

В классе болезней крови (D50-D89) необходимо обратить особое внимание на эритробластопению и другие апластические анемии (D60-D61) с очень высоким коэффициентом радиационно обусловленной заболеваемости, ERR/Гр=6,85, который в 1,5 раза выше коэффициента риска для лейкозов в этой же когорте ликвидаторов.

Наиболее важными заболеваниями при дальнейших расчётах возможного радиационного ущерба ликвидаторов являются болезни системы кровообращения (I00-I99) и болезни органов пищеварения (K00-K93), поскольку для этих болезней идентифицированы радиационные риски смертности с коэффициентами ERR/Гр=0,49 при 95% ДИ (0,28; 0,71) и ERR/Гр=0,73 при 95% ДИ (0,20; 1,34) соответственно. Эти оценки в 4 и 7 раз превышают оценки, полученные ранее в японской когорте лиц, переживших атомные бомбардировки в 1945 г., и по величине соответствуют оценкам ERR/Гр для солидных ЗНО, рекомендованных МКРЗ для обоснования норм и правил радиационной защиты. Среди болезней системы кровообращения (I00-I99) наибольший коэффициент радиационно обусловленной смертности, в 5 раз превышающий средний по классу, имели болезни, характеризующиеся повышенным кровяным давлением (I10-I15). В классе болезней органов пищеварения (K00-K93) наибольшие коэффициенты радиационно обусловленной смертности наблюдались для болезней печени (K70-K76) и болезней поджелудочной железы (K85-86).

Для смертности от болезней системы кровообращения, являющихся основной причиной смерти ликвидаторов из всех исследованных в данной работе, непараметрические оценки относительных радиационных рисков (RR) в дозовых интервалах статистически значимы при дозах облучения свыше 0,150 Гр. При дозах облучения менее 0,150 Гр оценённая ЛБП модель радиационных рисков смертности с коэффициентом ERR/Гр=0,49 является консервативной, т.е. приводит к завышенным оценкам радиационных рисков.

Оценки коэффициентов радиационных рисков неонкологических заболеваний, в первую очередь - болезней системы кровообращения (I00-I99) и болезней органов пищеварения (K00-K93), полученные в данном исследовании, должны быть использованы для оценки вклада неонкологических болезней в общий радиационный ущерб здоровью человека.

Принимая во внимание проведённые ранее исследования, которые показали значительное влияние мультиморбидности на вариабельность оценок коэффициентов радиационных рисков неонкологических заболеваний, практическое использование последних в задачах радиационной защиты требует более дифференцированных подходов к защите различных групп облучённых лиц, чем просто разделение по полу и возрасту, как это принято при учёте стохастических эффектов в виде радиационно обусловленных ЗНО. Влияние мультиморбидности на оценки радиационных рисков неонкологических заболеваний требует дальнейшего исследования.