Согласно современным представлениям теломеры — это важнейшие концевые участки хромосом, которые выполняют большое количество функций, основная из которых — сохранение целостности генома. Отмечено, что длина теломер у каждого человека варьирует в широком диапазоне даже в пределах одной клетки и каждой хромосомной пары, что можно объяснить неравномерным укорочением теломер [1][2]. Такой феномен распределения длины теломер на отдельных плечах индивидуальных хромосом называется «теломерным профилем хромосомы» [3]. Теломерный профиль может являться потенциальным биомаркером для прогнозирования рисков развития некоторых заболеваний, а также для оценки воздействия внешних неблагоприятных факторов, в том числе ионизирующего излучения [4]. В связи с этим исследования, направленные на изучение влияния радиационного воздействия на длину теломер, являются актуальными. Облучение людей при атомных взрывах в Хиросиме и Нагасаки в дозах, условно разделенных на низкие (5–700 мЗв) и высокие (>700 мЗв), приводило к укорочению теломер, и эти изменения сохранялись даже в отдаленные сроки (через 50–68 лет) после воздействия в высоких дозах [5]. В некоторых исследованиях, напротив, у облученных людей установлено увеличение длины теломер, что связывали с повышенной экспрессией теломеразы [6].

При оценке относительной длины теломер у работников атомной отрасли, подвергавшихся внутреннему α- (0,05 Гр; 0,05–0,1 Гр; >0,1 Гр) и внешнему γ- (1 Гр; 0,1–1,5 Гр; >1,5 Гр) облучению, было установлено, что снижение этого показателя обнаруживалось лишь после облучения в низких дозах как при внешнем гамма-облучении (при дозах <1 Гр), так и при внутреннем альфа-облучении (при дозах <0,05–0,1 Гр). При более высоких дозах этот показатель не отличался от значений у необлученного контроля [7]. Исследование, проведенное с участием лиц, подвергавшихся хроническому облучению в широком диапазоне доз (0–4,7 Гр), выявило снижение показателя относительной длины теломер при хроническом облучении (0,6–4,7 Гр) относительно группы сравнения для отдельных хромосомных плеч мета- и акроцентрических хромосом. Не было отмечено монотонного снижения длины теломер на хромосомах одной клетки или в разных плечах одной хромосомы [1]. Это свидетельствует о том, что длина теломер является динамической структурой и определяется не просто хронологическим возрастом, а совместным влиянием множества факторов эндогенной и экзогенной природы.

В научном сообществе обсуждаются вопросы, касающиеся изучения связи длины теломер с радиочувствительностью клеток и организма [8]. Так, исследования показывают, что некоторые заболевания, для которых характерна клиническая повышенная радиочувствительность (атаксия-телеангиэктазия, синдром ломкой Х-хромосомы, анемия Фанкони), отличаются нарушениями в поддержании структуры теломер [9].

Изучение радиочувствительности фибробластов человека при облучении рентгеновским in vitro излучением показало, что клоны клеток, сохранившие жизнеспособность после воздействия облучением в дозе 4 Гр, характеризовались увеличенной длиной теломер по сравнению с необлученными клонами. Авторами была выдвинута гипотеза о том, что клетки с более длинными теломерами являются более радиорезистентными, а радиочувствительность коррелирует с потерей теломер, а не с их средней длиной [4].

В ряде исследований показано, что количество сверхкоротких теломер, а не их средняя длина является наиболее информативным прогностическим маркером развития некоторых заболеваний [10]. Например, установлена взаимосвязь количества сверхкоротких теломер с различными заболеваниями человека, такими как хроническая обструктивная болезнь легких (ХОБЛ) и идиопатический легочный фиброз [11].

Все существующие методы исследования теломер хромосом подразделяются на две группы: молекулярные и молекулярно-цитогенетические. Наиболее подходящим для индивидуальной оценки теломерных участков является молекулярно-цитогенетический метод количественной флуоресцентной гибридизация in situ (Q-FISH) в метафазных хромосомах. Q-FISH — метод, в ходе которого используют флуоресцентно меченные пептидные нуклеиновые кислоты, которые специфически гибридизируются с денатурированной ДНК теломер. С помощью метода Q-FISH можно исследовать теломерные участки в каждой хромосоме каждой клетки, так как в процессе анализа проводится кариотипирование хромосом [1].

Настоящее исследование посвящено изучению длины теломер у облученных лиц, подвергшихся хроническому радиационному воздействию в широком диапазоне доз вследствие проживания на загрязненных радионуклидами территориях, в отдаленные сроки наблюдения (через 60–65 лет). Облучение людей было сочетанным: внутренним — за счет поступления и накопления в организме радионуклидов 89,90Sr и внешним — за счет γ-облучения. Критическим органом при воздействии излучения был красный костный мозг (ККМ). Дозы облучения на ККМ в основном были сформированы до 1960 г. [12].

Цель работы — изучение длины теломер методом Q-FISH у облученных лиц, а также выявление доли сверхкоротких и сверхдлинных теломерных участков хромосом.

Исследование проведено с участием 43 доноров (женщин) различных возрастных групп (21–28, 60–67, 71–83 года). В него включали лиц, у которых не было в анамнезе аутоиммунных и онкологических заболеваний, хронических воспалительных заболеваний в фазе обострения. Исключали лиц, принимающих цитостатики и антибиотики. Информация о состоянии здоровья облученных лиц была получена из базы данных «Человек», индивидуальные дозы облучения были рассчитаны по TRDS-2016 в биофизической лаборатории, информация о наличии онкопатологии в анамнезе обследуемых лиц была предоставлена эпидемиологической лабораторией [13].

На первом этапе было проведено исследование зависимости относительной длины теломер от возраста в контрольной группе лиц (характеристика доноров представлена в таблице 1). Были сформированы три возрастные группы доноров:

Доноры старшей возрастной группы проживали на загрязненных радионуклидами территориях, но их накопленная доза облучения ККМ не превышала 0,05 Гр за период жизни, в связи с чем далее (на II и III этапах исследования) этих доноров мы будем называть необлученными лицами или группой сравнения.

На II этапе в группе сравнения был определен референс средней длины теломер.

На III этапе, учитывая установленные референсные значения, изучали длину теломер у облученных лиц, в том числе в зависимости от возраста и дозы облучения ККМ.

Для оценки влияния ионизирующего излучения на длину теломер были сформированы две группы облученных доноров: подвергшиеся облучению в диапазоне доз 0,14–0,86 Гр (10 человек в возрасте 69–80 лет) составили группу облученных в средних дозах; группу облученных в высоких дозах составили 12 человек в возрасте 71–84 лет с дозами облучения ККМ 1,04–3,10 Гр (табл. 2).

Также для исключения влияния возрастного фактора был проведен анализ методом «случай — контроль». Для этого проводили отбор возможных контролей для каждого донора («случая») из группы облученных в средних дозах по критерию соответствия возраста на момент обследования. В качестве влияющего фактора рассматривалась доза облучения ККМ. Были сформированы 4 пары доноров близкого возраста. В каждой паре один из доноров имел установленную накопленную дозу облучения ККМ в диапазоне 0,61–0,86 Гр, второй донор из пары либо не подвергался аварийному облучению, либо имел дозу ≤0,05 Гр (табл. 3).

Для исследования «случай — контроль» в диапазоне высоких доз на ККМ было отобрано 5 пар доноров близкого возраста. Один из доноров в паре имел дозу свыше 1 Гр, второй донор либо не подвергался аварийному облучению, либо имел дозу ≤0,05 Гр (табл. 3).

Предметом исследования был ядерный хроматин стимулированных фитогемагглютинином (ФГА) Т-лимфоцитов периферической крови. Выбор Т-лимфоцитов обусловлен простотой получения исходного материала и достаточно высокой концентрацией клеток. Т-лимфоциты не делятся в периферической крови и находятся в G0-фазе клеточного цикла, являясь синхронизированной в организме популяцией клеток. Клетки-предшественники лимфоцитов подвергались облучению непосредственно в ККМ.

Забор крови производился в вакуумные пластиковые пробирки для взятия крови с гепарином лития.

Культуральную смесь готовили из следующих компонентов: 2 мл цельной гепаринизированной крови, 5 мл среды RPMI-1640 («Панэко», Россия), 0,5 мл отборной эмбриональной телячьей сыворотки (РАА Laboratories, Австрия), ФГА добавляли в итоговой концентрации 20 мг/мл («Панэко», Россия). Клетки культивировали в СО2-инкубаторе (SANYO MCO-18AIC, Япония) при 37,5 °С в течение 54 ч. За 3 ч до окончания культивирования вводили колхицин в итоговой концентрации 0,1 мг/мл («Панэко», Россия). Далее проводили гипотоническую обработку метафазных клеток, фиксацию метафазных пластинок (в соотношении 3:1 части этанола медицинского 95% и ледяной уксусной кислоты) и получение препаратов хромосом. Для флюоресцентного окрашивания методом Q-FISH стекла обрабатывали растворами РНК-азы (100 нг/мл) и пепсина, после чего проводили денатурацию ДНК зонда и препарата. Гибридизация проходила в соответствии с протоколом производителя, к которому прилагались оригинальные растворы.

Для оценки длины теломерных районов хромосом использовали зонды производства DAKO (Дания), для окраски центромерного района второй хромосомы, являющегося реперным сигналом, применяли центромерные зонды (Metasystems, Германия). Анализ флуоресцентно окрашенных препаратов проводили на флуоресцентном микроскопе AxioImager Z2 (Zeiss; Германия) с фильтрами DAPI и SpO (spectrum orange).

Флуоресцентный сигнал регистрировали и измеряли его интенсивность с помощью специального программного обеспечения для анализа изображения Isis. Анализировали 20–30 клеток на одного обследуемого. Результат получали отдельно для коротких (p-) и длинных (q-) плеч каждой из 46 хромосом. Длина теломер являлась относительной, поскольку выражалась в процентах отношения длины теломерного (Т) сигнала к длине центромерного (С) — (Т/C%). Подробно метод оценки длины теломерных районов описан в работе [14].

Для описания результатов использовали стандартные методы описательной статистики с вычислением медианы, ее 95% доверительного интервала (95% ДИ) и 5-го и 95-го процентилей (%). Доверительные интервалы рассчитывали процедурой бустрепа при n = 1000. Сравнение значений в группах осуществляли с помощью непараметрического критерия Манна – Уитни.

В качестве порога для сверхкоротких теломер было принято значение, ниже которого расположены 5% всех измерений длины теломерного участка в группе необлученных лиц. Аналогичным образом был рассчитан референсный показатель для сверхдлинных теломер: это все значения, превышающие 95% распределения измерений. Полученные референсные показатели были использованы для сравнительного анализа с данными относительной длины теломер у лиц, подвергшихся радиационному воздействию в разных диапазонах доз. При сравнении частоты встречаемости сверхкоротких и сверхдлинных теломер в исследовании «случай — контроль» и при сравнении с референсным значением использовали метод χ². Для определения связи длины теломер с возрастом и дозой был проведен многофакторный анализ. Для определения связи дозы облучения и количества сверхкоротких и сверхдлинных теломер у человека использовали коэффициент корреляции по Спирмену, а также регрессионный анализ. Статистический анализ выполнен с помощью программного обеспечения Past 4.01 (Oyvind Hammer; Великобритания), Statistica 10 (StatSoft, США), Sigmaplot (SYSTAT Software, USA).

В ходе работы получены данные измерений длины теломер для трех возрастных групп необлученных доноров (табл. 4).

В группе доноров 20–28 лет было выполнено 4331 измерение длины теломерных участков, минимальная длина теломер зафиксирована на уровне 0%, максимальная — 900,0%, а медианное значение длины теломер составляло 31,0%.

В то же время в группе доноров средней возрастной группы (60–67 лет) было проведено 22 311 измерений теломерных участков, значение (min–max) длины теломер составило 0–967,9%, а медианное значение — 13,0%.

В группе доноров в возрасте 71–83 года было получено 9732 измерения, минимальное значение длины теломер 0%, максимальное — 141,2%, медиана — 5,8%.

Значения длин сверхкоротких теломер для доноров младшей возрастной группы находились в диапазоне 0–6,3%, для доноров 60–67 лет — в диапазоне 0–1,8%, для доноров 71–83 лет — в диапазоне 0–0,7%. При этом диапазоны значений сверхдлинных теломер у доноров 20–28 лет составили 221,2–900,0%, у доноров 60–67 лет — 105,9–967,9%, у доноров 71–83 лет — 25,6–141,2%.

Было установлено, что медианные значения длины теломер снижались с увеличением возраста доноров. Попарное сравнение значений в исследуемых контрольных группах лиц показало статистически значимые (p = 0,0001) различия. С помощью проведенного регрессионного анализа показано наличие достоверной (p < 0,001; R² = 0,07), но крайне слабой линейной зависимости между относительной длиной теломер и возрастом доноров:

y = 86,1 – 26,7 × x, (1)

где y — относительная длина теломер, х — возраст обследуемых лиц.

Для дальнейшего анализа длины теломер у облученных лиц за референсную длину сверхкоротких теломер принято значение 0,7%, что соответствовало нижней границе 90% диапазона длины теломер у необлученных доноров возрастом 71–83 года. А сверх­длинные теломеры оценивали по верхней границе 90% диапазона длины теломер — 25,6%

В таблице 5 представлены данные измерений длины теломер для облученных доноров. Группой сравнения являются необлученные доноры в возрасте 71–83 года.

Как видно, медианное значение длины теломер у всех облученных составило 10,3%, что статистически значимо выше, чем в группе сравнения (5,8%; p = 0,0001). Вместе с тем 90% интервал длины теломер находился в диапазоне 1,4–56,5%.

Учитывая референсное значение длины сверхкоротких теломер, в группе облученных доноров (0,14–3,10 Гр) сверхкороткими были 616 измерений теломер, что составило 1,6% от числа всех измерений. В референсном диапазоне для сверхдлинных теломер у облученных доноров было 7491 измерение, что составило 19,5%. Таким образом, у облученных лиц наблюдали меньшее количество сверхкоротких теломер (χ² = 397,03; p < 0,0001) и большее количество сверх­длинных теломер (χ² = 1184; p < 0,0001) относительно группы сравнения.

В группе облученных на уровне средних доз (0,14–0,86 Гр) медиана длины теломер составила 11,9% (min–max: 0–643,4%). У доноров, облученных в высоких дозах, медиана составила 10,6%, (min–max: 0–332,1%), что статистически значимо ниже, чем у облученных в средних дозах (p = 0,0001).

Относительно референсных значений частота сверхкоротких теломер у доноров, облученных в средних дозах, составила 1,3%, а при высоких дозах облучения — 1,4%. Частота сверхдлинных теломер у доноров, облученных в средних дозах, составила 24%, в высоких дозах — 19,5%, что статистически значимо отличалось от значений в группе сравнения, p = 0,0001. Таким образом, у облученных лиц с высокой накопленной дозой на ККМ частота сверхкоротких теломер такая же, как и в группе со средними дозами, но количество сверхдлинных теломер выше.

Для установления связи длины теломер с возрастом и дозой облучения ККМ был проведен двухфакторный анализ (общая линейная модель). Была получена статистически незначимая модель: длина теломер во всей выборке облученных лиц не зависела от возраста и дозы облучения ККМ.

В таблице 6 представлены данные попарного сравнения для анализа «случай — контроль», когда «случаи» представлены лицами со средними (0,14–0,86 Гр) и высокими (1,04–3,10 Гр) накопленными дозами на ККМ.

У всех исследуемых пар из группы «средние дозы облучения» было выявлено статистически значимое повышение длины теломер у облученных доноров относительно их контроля.

В группе лиц, облученных в высоких дозах, у четырех пар из пяти было выявлено статистически значимое повышение длины теломер у облученных доноров относительно их контроля. В паре № 1 обнаружено статистически значимое снижение длины теломер у облученного донора относительно его контроля. Но при этом в данной паре «случай — контроль» частота сверхкоротких теломер не отличалась от контрольных значений (0,8% против 0,7%; χ² = 0,03; p = 0,86), а частота сверхдлинных теломер была значимо ниже (1,4% против 2,9% χ² = 8; p = 0,005) (табл. 7).

В таблице 7 представлены частоты сверхкоротких и сверхдлинных теломер для каждой пары «случай — контроль» (на основе референсных значений) среди лиц со средней и высокой дозой облучения на ККМ.

В работе получены статистически значимые данные об изменении частоты сверхкоротких теломер у трех из четырех пар (p << 0,0001), за исключением пары № 3 (p = 0,19). При этом у всех исследуемых пар частота встречаемости сверхкоротких теломер у облученных доноров статистически значимо ниже, чем в контроле. В то же время частота сверхдлинных теломер во всех парах облученных доноров в дозах (0,14–0,86 Гр) статистически значимо выше по сравнению с группой контроля (p << 0,0001).

Как видно, у лиц, облученных в диапазоне доз 1,04–3,10 Гр, частота сверхкоротких теломер статистически значимо ниже относительно контроля для четырех пар из пяти. В первой паре не отмечено статистически значимых различий в частоте сверхкоротких теломер. Также в этой же паре у облученного донора более низкое значение частоты сверхдлинных теломер по сравнению с контролем. У остальных исследуемых четырех пар наблюдали увеличение частоты сверхдлинных теломер у облученных лиц относительно контроля.

Далее для доноров, отобранных в исследование «случай — контроль» по методу ранговой корреляции Спирмена, оценили взаимосвязь между долями сверхкоротких теломер и дозой облучения. Результаты показали, что корреляционной зависимости сверхдлинных теломер от дозы на ККМ отмечено не было, но имеется тенденция уменьшения частоты сверхкоротких теломер с увеличением дозы облучения ККМ. Установлена статистически значимая обратная корреляционная зависимость умеренной степени выраженности между частотой сверхкоротких и сверхдлинных теломер у облученных доноров (rs = -0,654, p = 0,004).

Далее была подобрана наиболее подходящая регрессионная зависимость частоты сверхкоротких теломер от дозы на ККМ (рис.).

На рисунке показано, что с увеличением дозы облучения количество сверхкоротких теломер экспоненциально уменьшается.

y = 4,3 e–1,5D [ R² = 0,23, p = 0,0036], (2)

где y — частота сверхкоротких теломер, D — доза облучения ККМ, Гр.

Таким образом, было показано, что у облученных лиц снижена частота сверхкоротких теломер и увеличена частота сверхдлинных теломер относительно группы сравнения. Отмечено экспоненциальное снижение частоты сверхкоротких теломер с увеличением дозы облучения ККМ.

Настоящая статья продолжает проект по изучению влияния хронического радиационного воздействия на теломерные участки хромосом. В рамках многолетнего исследования была изучена длина теломерных участков отдельных хромосомных пар у хронически облученных лиц, оценивалась частота инверсий с захватом теломерного участка и частота потерь теломерных участков метафазных хромосом в культивированных Т-клетках периферической крови у облученных жителей Южного Урала [14]. Для этих целей был предложен и апробирован метод флуоресцентной окраски теломерных районов хромосом. Использование метода Q-FISH обосновано тем, что он дает возможность исследовать теломерные участки в каждой хромосоме клетки донора. Это позволяет получить представление о теломерном профиле.

В представленной работе впервые проведен анализ влияния радиационного и нерадиационного фактора (возраста) на частоту сверхкоротких и сверхдлинных теломер. Частоту сверхкоротких и сверхдлинных теломер у облученных лиц оценивали исходя из референсных значений длины теломерных районов, рассчитанных для группы сравнения схожего возрастного диапазона. Референсное значение длины теломер в группе сравнения (возраст 71–83 года) составило 0,7–25,6%. Теломеры, длина которых ниже 0,7%, считались сверхкороткими, а если выше 25,6% — сверх­длинными.

В группе облученных лиц показано статистически значимое повышение медианной длины теломер относительно группы сравнения. При этом как у облученных лиц, так и в группе сравнения наблюдали широкую вариабельность длины теломер, что согласуется с ранее полученными данными [1][15]. Эта вариабельность может свидетельствовать о циркулировании в периферической крови Т-клеток разного репродуктивного возраста.

Согласно результатам в общей выборке облученных лиц частота сверхкоротких теломер была ниже, а сверхдлинных — выше по сравнению с группой сравнения. Эти данные подтверждаются и оценкой с помощью метода «случай — контроль» с подбором сопоставимых по возрасту доноров. Более длинные теломеры у облученных лиц могут свидетельствовать как о молекулярно-генетических особенностях, связанных с повышенной теломеразной активностью [16] или эпигенетическими изменениями других ферментов, регулирующими репликацию и репарацию теломер [17][18], так и об отборе более радиорезистентных клеток и гибели радиочувствительных [4]. Также здесь мог сыграть роль фактор отбора доноров для цитогенетического исследования, когда в выборку попадают наиболее здоровые лица относительно генеральной совокупности, что является необходимым критерием отбора, поскольку исключаемые заболевания, лечебные процедуры и лекарственные препараты могут искажать цитогенетические данные [1].

Помимо этого, интересным является статистически значимое снижение медианной длины теломер в группе доноров с высокой дозой облучения ККМ относительно группы доноров, облученных в средних дозах. Данное различие в первую очередь было связано с тем, что в группе доноров с высокой дозой облучения ККМ снижена частота сверхдлинных теломер, а частота сверхкоротких находится на том же уровне, что у лиц со средними дозами. В то же время частота сверхдлинных теломер не зависела от дозы облучения ККМ, а частота сверхкоротких теломер экспоненциально снижалась с увеличением дозы. Представленные закономерности требуют дальнейшего более глубокого изучения. Если в последующих исследованиях подтвердится связь частоты сверхкоротких теломер с дозой облучения, то этот показатель может использоваться как индикатор радиационного воздействия или биомаркер радиочувствительности клеток [19][20].

В настоящем исследовании представлены данные о влиянии хронического радиационного воздействия на длину теломер у женщин, полученные с использованием метода Q-FISH. Этот метод позволяет определить относительную длину теломерных участков посредством расчета отношения интенсивности флюоресценции теломерного сигнала к интенсивности флюоресценции центромерного сигнала второй хромосомы, выраженную в процентах (Т/С%). Сверхкороткие и сверхдлинные теломеры у облученных лиц выявляли относительно референсного значения длины теломер, рассчитанного в группе сравнения.

Анализ полученных результатов показал статистически значимое повышение медианной длины теломер у облученных лиц относительно группы сравнения. Повышение длины теломер связано со снижением частоты сверхкоротких теломер и повышением частоты сверхдлинных теломер у облученных лиц. Исследования, выполненные методом «случай — контроль» с подбором сопоставимых по возрасту доноров, подтвердили, что данная закономерность сохраняется независимо от возраста, подчеркивая значимость радиационного фактора. Зависимость частоты сверхкоротких теломер от дозы облучения ККМ описывается нелинейной регрессией: наблюдается экспоненциальное снижение частоты сверхкоротких теломер с увеличением дозы. Этот феномен подчеркивает важность дальнейшего изучения роли сверхкоротких теломер как потенциального маркера радиационного воздействия.

Вклад авторов. Все авторы подтверждают соответствие своего авторства критериям ICMJE. Наибольший вклад распределен следующим образом: Я.В. Кривощапова — идея и план исследования, проведение лабораторных исследований, получение первичных данных, статистическая обработка результатов, написание текста статьи; Ю.Р. Ахмадуллина — статистическая обработка и интерпретация результатов, написание, редактирование и подготовка итоговой версии статьи.