На правах рукописи
Корытина Гульназ Фаритовна
Mолекулярно-генетический анализ наследственной предрасположенности
к хроническим заболеваниям органов дыхания
03.02.07 - генетика
Автореферат
диссертации на соискание ученой степени
доктора биологических наука
Уфа - 2012
Работа выполнена в Учреждении Российской Академии наук
Институте биохимии и генетики Уфимского научного центра РАН
Научный консультант | доктор медицинских наук, профессор Викторова Татьяна Викторовна |
Официальные оппоненты: | доктор биологических наук, профессор Иващенко Татьяна Эдуардовна, ГУ НИИ акушерства и гинекологии им. Д.О. Отта РАМН |
доктор медицинских наук, профессор Полоников Алексей Валерьевич, Курский государственный медицинский университет | |
доктор биологических наук, доцент Хидиятова Ирина Михайловна, Институт биохимии и генетики Уфимского научного центра РАН | |
Ведущая организация: | ГУ Медико-генетический Научный Центр РАМН, г. Москва |
Защита состоится л____марта_2012 г. в л_____ часов
на заседании Объединенного Диссертационного совета ДМ 002.133.01 при
Учреждении Российской Академии наук Институте биохимии и генетики
Уфимского научного центра РАН по адресу:
450054, г. Уфа, пр. Октября, 71. ИБГ УН - РАН
С диссертацией можно ознакомиться в библиотеке Уфимского научного
центра РАН (Уфа, пр. Октября, 71),
с авторефератом - на сайтах ВАК РФ и ИБГ УН - РАН
ibg.anrb.ru/dissov.html
e-mail: molgen@anrb.ru
Автореферат разослан л___________________2012 г.
Ученый секретарь
диссертационного совета Бикбулатова С.М.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность проблемы. Перспективы в изучении природы широко распространенных заболеваний связаны с познанием генома человека, с выявлением его особенностей, создающих предпосылки для формирования патологического фенотипа (Пузырев, 2009; Баранов, 2009). В развитии предрасположенности человека к многофакторной патологии важную роль играют как генетическая конституция индивидуума, так и действие факторов внешней среды. При этом соотношение средовых и генетических факторов сильно варьирует при различных заболеваниях (Пузырев, 2009; Баранов, 2009). В структуре общей заболеваемости населения России болезни органов дыхания традиционно занимают лидирующую позицию - на их долю приходится до 27,8% всей зарегистрированной и 42,7% выявленной впервые в жизни патологии (Щепин, 2001). Среди хронических респираторных заболеваний у взрослых значительная часть приходится на хроническую обструктивную болезнь легких (ХОБЛ), у детей и подростков - на хронические бронхолегочные заболевания (ХБЛЗ) (Чучалин, 2009). Профессиональный бронхит формируется у 4,5-24,6% лиц, работающих во вредных и неблагоприятных производственных условиях (Лещенко с соавт., 2004). К факторам, которые способствуют развитию заболеваний органов дыхания, относят курение табака, профессиональную вредность (воздействие аэроирритантов, пыль, химических веществ), атмосферные поллютанты (диоксид серы, азота и углерода, свинец, угарный газ, кремний, асбест, уголь), частые инфекции респираторного тракта, а также социально-экономический статус пациента (GOLD, 2010). Наследственная предрасположенность является важным внутренним фактором риска развития хронических заболеваний органов дыхания, с которыми тесно связаны особенности иммунологической реактивности, роста и развития легких (GOLD, 2010). Многие аспекты развития хронических заболеваний органов дыхания до сих пор не ясны, но известны ключевые звенья, патогенеза. Это, прежде всего, хроническое воспаление, запускающее весь каскад патогенетических реакций, приводящих к нарушению баланса протеолитических ферментов и их ингибиторов, и развитию окислительного стресса. Генетические механизмы формирования ХОБЛ в последние годы стали объектом широкомасштабных исследований во всем мире. В результате проведенного полногеномного анализа сцепления при ХОБЛ, показана высокая степень сцепления локусов на хромосоме 19q - TFGB1 (transforming growth factor, beta 1, трансформирующий фактор роста бета1) и хромосоме 2q - SERPINE2 (serpin peptidase inhibitor, clade E, серпиновый ингибитор пептидаз, тип E) (Silverman, 2002; Palmer, 2003). Из опубликованных к настоящему моменту данных полногеномных исследований генетических ассоциаций (GWAS), выявлено несколько локусов, связанных с развитием ХОБЛ, на 4-й и 15-й хромосомах в области 15q25.1 - CHRNA3, CHRNA5, (cholinergic receptor nicotinic alpha 3, холинергические никотиновые рецепторы альфа 3 и 5), IREB2 (iron-responsive element binding protein 2, белок железо-связывающего элемента 2), PSMA4 (proteasome subunit alpha type 4, субъединица альфа 4-го типа протеосомного комплекса), 4q31.21 - HHIP (hedgehog interacting protein, взаимодействующий протеин hedgehog сигнального пути), 4q22.1 - FAM13A (family with sequence similarity 13 member A, член Аа семейства со схожей последовательностью 13) и новый маркер rs7937 на хромосоме 19q13 (Pillai, 2009; Wilk, 2009; Repapi, 2009; Cho, 2010; Wan, 2010; Kong, 2010, 2011; Siedlinski, 2011; Cho, 2011). До Эры GWAS генетическая составляющая ХОБЛ достаточно активно изучалась с использованием анализа генов-кандидатов (Wang et al., 2009; Sadeghnejad et al., 2009; Haq et. al., 2010; Silverman et al., 2011). Данный метод до сих пор наиболее широко иcпользуется при проведении исследований по генетике ХОБЛ, поскольку позволяет сосредоточиться на одном или нескольких функционально значимых аллельных вариантах гена, кодирующих соответственно варианты белка, различающиеся по структуре и функциям, некоторые из которых могут быть вовлечены в развитие патогенетических изменений. Несмотря на значительное число исследований, выполненных с использованием разных подходов, молекулярно-генетические основы ХОБЛ до сих пор остаются во многом не ясными. На сегодняшний день перспективным является изучение не только генетических и молекулярных, но и эпигенетических факторов (метилирования ключевых генов, анализ активности гистонов, роли микро-РНК) (Silverman et al., 2011). Работа в различных направлениях позволит приблизиться к пониманию сложных процессов, приводящих к развитию болезни. Выбранные для настоящего исследования хронические заболевания органов дыхания - хроническая обструктивная болезнь лёгких, хронические бронхолегочные заболевания у детей, хронический профессиональный бронхит - занимают лидирующие позиции по распространенности, тяжести течения, серьезности осложнений и являются основными причинами высоких темпов роста заболеваемости, инвалидизации и смертности населения, как в Российской Федерации, так и во многих других странах, что свидетельствует об актуальности изучения генетических основ данных заболеваний и выяснения основных причин их возникновения.
Цель работы: на основе комплексного анализа полиморфных локусов генов-кандидатов и факторов окружающей среды охарактеризовать особенности формирования наследственной предрасположенности к развитию хронической обструктивной болезни легких, хроническим бронхолегочным заболеваниям у детей и профессиональному хроническому бронхиту
Задачи исследования:
1. Провести сравнительный анализ полиморфных вариантов генов ферментов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза в трех этнических группах жителей Республики Башкортостан (русские, татары, башкиры).
2. Изучить ассоциации полиморфных вариантов генов-кандидатов с предрасположенностью к хронической обструктивной болезни легких, хроническим бронхолегочным заболеваниям у детей и профессиональному хроническому бронхиту.
3. Выявить генетические локусы, ассоциированные с развитием хронических заболеваний органов дыхания, с учетом этнической принадлежности и полового диморфизма.
4. Исследовать взаимосвязь полиморфных вариантов генов-кандидатов с клиническими формами, тяжестью течения и возрастом манифестации хронической обструктивной болезни легких с учетом влияния факторов окружающей среды и гендерных различий.
5. Проанализировать характер взаимодействий между генами-кандидатами и средовыми факторами (статус и интенсивность курения, стаж работы во вредных производствах), определить их роль в этиопатогенезе хронических заболеваний органов дыхания.
6. На основе изучения моделей взаимодействия генов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза оценить значимость ключевых генов в формировании наследственной предрасположенности к хроническим заболеваниям органов дыхания.
Научная новизна: Впервые проведен комплексный анализ формирования наследственной предрасположенности к развитию хронических заболеваний органов дыхания, спровоцированных курением, профессиональной деятельностью и респираторными инфекциями, на основе анализа функционально значимых полиморфных локусов генов-кандидатов, изучения взаимодействия генетических и средовых факторов, выявления этнической гетерогенности полученных ассоциаций. Впервые получены данные о частотах аллелей и генотипов полиморфных локусов генов CYP1A2, CYP1B1, CYP2A6, CYP2E1, CYP2F1, CYP2J2, CYP2S1, NQO1, EPHX1, UGT2B7, GPX1, MMP1, 2, 3, 9, 12; TIMP2, 3; SERPINA1, SERPINA3, VDBP в этнических группах русских, татар и башкир, проживающих в Республике Башкортостан. Впервые установлена ассоциация полиморфных вариантов генов NQO1, MMP3, CYP2S1, CYP2F1, TIMP3, UGT2B7, MMP9 с развитием хронических заболеваний органов дыхания. Определены значимые ген-средовые взаимодействия с курением для локусов ADAM33 (12418A>G), TIMP3 (-1296T>C), VDBP (1296T>G), CYP1A1 (2454A>G), CYP2F1 (c.14_15insC), GPX1 (599C>T), SOD3 (691C>G), UGT2B7 (2146C>), EPHX1 (415A>G) и со стажем работы для локусов IL1RN (VNTR), VDBP (1307C>A) и CYP1A1 (3798T>C) при формировании хронических заболеваний органов дыхания. Выявлена ассоциация полиморфных локусов UGT2B7 (2146C>T), NAT2 (590G>A), CYP1A1 (2454A>G) с возрастом манифестации, LTA (252A>G), IL6 (-174G>C), CYP2E1 (-1053C>T), MMP9 (-1562C>T) с тяжестью течения, UGT2B7 (2146C>T), ADAM33 (13491 C>G), IL1RN (VNTR), VDBP (1296T>G) с клиническими вариантами хронической обструктивной болезни легких.
Научно-практическая значимость работы: Полученные результаты вносят вклад в понимание структуры наследственной предрасположенности к развитию хронических заболеваний органов дыхания и являются базой для познания механизмов их патогенеза, профилактики, диагностики и лечения. Результаты исследования широкого спектра генов-кандидатов с хроническими заболеваниями органов дыхания, анализ их взаимодействия с факторами среды, особенности ассоциации генетических маркеров в дифференцированных по полу и этнической принадлежности группах могут быть использованы в профилактической медицине для формирования групп повышенного риска. Данные диссертационной работы могут быть использованы в учебном процессе на биологических и медицинских факультетах ВУЗов, а также на курсах последипломного образования.
Основные положения, выносимые на защиту:
1. Полиморфизм генов ферментов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза в этнических группах (русских, татар, башкир), проживающих в Республике Башкортостан, характеризуется широким аллельным разнообразием. Межэтнические различия по распределению частот аллелей и генотипов полиморфных локусов генов CYP1A1 (3798T>C), CYP1A1 (2454A>G), CYP1A2 (-2467delT), CYP2A6 (del), CYP2E1 (-1053C>T), SOD1 (c.239+34A>C), IL1RN (VNTR), IL10 (-627C>A), VDBP (1296T>G).
2. С развитием хронической обструктивной болезни легких ассоциированы полиморфные варианты генов NQO1 (465C>T), CYP2A6 (del), CYP1A2 (-163C>A), CYP2S1 (13255A>G), EPHX1 (337T>C), GSTP1 (313A>G), SOD3 (691C>G), CAT (-262C>T), IL6 (-174G>C), VDBP (1296T>G и 1307C>A), MMP3 (-11715А>6А), ADAM33 (12418A>G). Развитие хронических бронхолегочных заболеваний у детей ассоциировано с полиморфными вариантами генов CYP2F1 (c.14_15insC), CYP1A1 (2454A>G, 3798>C), NQO1 (465C>T), GSTT1 (del), GSTP1 (313A>G), CAT (-262C>T), IL6 (-174G>C), IL1RN (VNTR), VDBP (1296T>G и 1307C>A), MMP1 (-1607 G>GG), MMP12 (-82A>G), MMP3 (-11715А>6А), MMP9 (2660A>G). Ассоциация с профессиональным хроническим бронхитом полиморфных вариантов генов VDBP (1296T>G и 1307C>A), UGT2B7 (2146C>T), CYP1A2 (-2467delT), GSTP1 (313A>G), NQO1 (465C>T), CAT (-262C>T), IL8 (-251T>A), ADAM33 (13491C>G), MMP9 (2660A>G), MMP1 (-519A>G). Спектр выявленных маркеров наследственной предрасположенности при каждой патологии существенно варьирует в зависимости от этнической принадлежности и полового диморфизма.
3. Полиморфизм генов UGT2B7 (2146C>T), NAT2 (857G>A), GSTM1 (del), CYP1A1 (2454A>G) вносит вклад в вариабельность возраста манифестации хронической обструктивной болезни легких; ассоциированы с тяжестью течения заболевания полиморфные варианты генов LTA (252A>G), IL6 (-174G>C), CYP2E1 (-1053C>T) MMP9 (-1562C>T); развитие эмфиземы легких ассоциировано с локусами UGT2B7 (2146C>T), ADAM33 (13491C>G), IL1RN (VNTR), VDBP (1296T>G).
4. Структура наследственной предрасположенности при формировании хронических заболеваний органов дыхания различается у индивидов подверженных курению и некурящих. При хронической обструктивной болезни легких и профессиональном хроническом бронхите наблюдаются значимые взаимодействия средовых факторов (статуса и интенсивности курения, стажа работы во вредных производствах) с полиморфными вариантами генов ADAM33 (12418A>G), TIMP3 (-1296T>C), VDBP (1296T>G), CYP1A1 (2454A>G), CYP2F1 (c.14_15insC), GPX1 (599C>T), SOD3 (691C>G), UGT2B7 (2146C>T), EPHX1 (415A>G), IL1RN (VNTR).
5. Характер межгенных взаимодействий исследованных полиморфных локусов генов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза варьирует в зависимости от конкретного заболевания, этнической принадлежности, полового диморфизма и статуса курения.
Апробация работы: Основные материалы диссертационной работы доложены и обсуждены на 14 - 21 Национальных конгрессах по болезням органов дыхания, (Москва, 2003, 2005, 2009, 2010; Санкт-Петербург, 2004, 2006; Казань, 2007; Екатеринбург, 2008; Уфа, 2011); на ежегодных международных конгрессах Европейского респираторного общества (Вена, 2003, 2009; Глазго, 2004; Копенгаген, 2005; Мюнхен, 2006; Стокгольм, 2007; Берлин, 2008; Барселона, 2010; Амстердам, 2011); на конференциях Всероссийского общества генетиков и селекционеров (Москва, 2004, 2009); на V и VI и Съездах Российского общества медицинских генетиков (Уфа, 2005; Ростов-на-Дону, 2010); на IX и Х Всероссийских Форумах с международным участием им. акад. В.И. Иоффе Дни иммунологии в Санкт-Петербурге (Санкт-Петербург, 2005, 2006); на XIX международном конгрессе аллергологов и клинических иммунологов (Мюнхен, 2005); на ежегодных конференциях Европейского общества по генетике человека (Копенгаген, 2005; Амстердам, 2006, 2011; Ницца, 2007; Барселона, 2008; Вена, 2009; Гётенбург, 2010); на международном конгрессе Геном человека Human Genome Meeting (Киото, 2005); на Всероссийском научном симпозиуме с международным участием Цитокины. Стволовая клетка. Иммунитет, (Новосибирск, 2005); на II Санкт-Петербургском международном экологическом форуме (Санкт-Петербург, 2008); на 3-м конгрессе Европейского общества педиатров Europaediatrics (Стамбул, 2008).
Публикации: По теме диссертации опубликовано 78 работ, из них 32 статьи в журналах из Перечня ВАК, в том числе, 15 в журналах из международной базы цитирования (PubMed), 1 статья в зарубежной печати. В виде тезисов и статей в материалах международных и российских научно-практических конференций опубликовано 45 работ.
Структура и объем диссертации: Диссертация изложена на 411 страницах машинописного текста. Данные проиллюстрированы таблицами (85), рисунками (51) и приложениями (2). Список литературы включает 577 источников. Благодарности. Автор выражает глубокую благодарность за неоценимую помощь в проведении совместных исследований и обсуждении результатов работы, директору ИБГ УН - РАН, д.б.н., проф. Вахитову В.А., зав. Отделом геномики ИБГ УН - РАН д.б.н., проф. Хуснутдиновой Э.К., главному пульмонологу РБ, зав. кафедрой пропедевтики внутренних болезней БГМУ д.м.н., проф. Загидуллину Ш.З. (БГМУ), зав. пульмонологическим отделением ГКБ №21 д.м.н, проф. Фархутдинову У.Р., к.м.н. Байнак О.В. (БГМУ), директору Уфимского научно-исследовательский института медицины труда и экологии человека, г. Уфы д.м.н., проф. Бакирову А.Б. и врачуЦпрофпатологу отделения профессиональной аллергологии и иммунореабилитации к.м.н., Мингазовой С.Р., зав. кафедрой детских болезней БГМУ д.м.н., проф. Эткиной Э.И., к.м.н. Бабенковой Л.И. (БГМУ), врачу-педиатру пульмонологического отделения РДКБ Нуриахметовой А.Ж.; аспиранту Department of Genetic Epidemiology, George-August University of Goettingen (Германия) Викторовой Е.В. Сотрудникам и коллегам лаборатории экологической генетики человека ИБГ УН - РАН к.б.н. Ахмадишиной Л.З., к.б.н. Кочетовой О.В., к.б.н. Янбаевой Д.Г., к.б.н. Данилко К.В., к.б.н. Целоусовой О.С., сотрудникам лаборатории молекулярной генетики человека ИБГ УН - РАН к.б.н. Хусаиновой Р.И., к.м.н. Карунас А.С. и всем коллегам лаборатории физиологической генетики и молекулярной генетики человека ИБГ УН - РАН.
Материалы и методы исследования
Объект исследования. Для оценки частот полиморфных вариантов генов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза в этнических группах жителей Республики Башкортостан была сформирована выборка из 794 человек, включающая индивидов разной этнической принадлежности: русские (N=376), татары (N=301) и башкиры (N=147). Для проведения анализа ассоциации генов-кандидатов с хроническими заболеваниями органов дыхания были сформированы три группы, подобранные по принципу случай-контроль, с учетом возраста, пола, этнической принадлежности и влияния факторов среды. Характеристики изученных групп приведены в таблицах 1, 2, 3. Диагноз ХОБЛ устанавливался лечащими врачами городских клинических больниц г. Уфы №21, 13, 22, 18 в соответствии с международной классификацией болезней 10-го пересмотра (МКБ-10) и с учетом рекомендаций GOLD (Global initiative for chronic obstructive lung disease, 2010). Больные ХОБЛ ранее не имели профессиональной вредности и не подвергались действию комплекса вредных производственных факторов. Критериями исключения из группы являлись наличие симптомов аллергических заболеваний, бронхиальной астмы, онкологических заболеваний, специфических инфекционных заболеваний органов дыхания (туберкулез). Группу контроля составили практически здоровые индивиды - жители Республики Башкортостан без патологии дыхательной системы в анамнезе и без профессионального контакта с вредными химическими веществами. Для выявления ассоциации полиморфных вариантов генов-кандидатов с тяжестью течения, возрастом манифестации и развитием эмфиземы легких внутри выборки больных выделяли следующие подгруппы: ХОБЛ с эмфиземой легких (N=202) и ХОБЛ без эмфиземы легких (N=189); группа больных с тяжелой ХОБЛ (N=296), включающая больных с IV стадией заболевания и осложнениями в виде пневмосклероза и/или эмфиземы легких, и ХОБЛ средней тяжести (N=95), которая включала больных со II и III стадиями ХОБЛ без осложнений; больные ХОБЛ с ранней манифестацией заболевания (до 40 лет, средний возраст манифестации составил 29.9210.68 лет) и группа с поздней манифестацией ХОБЛ (после 40 лет, средний возраст манифестации - 56.018.79 лет).
Таблица 1
Характеристика больных ХОБЛ и группа сравнения
Больные ХОБЛ n (%) | Контроль n (%) | |
мужчин женщин | 312 (79.79) 79 (20.21) | 417 (81.13) 97 (18.87) |
возраст (лет) Mm | 61. 3 12.7 | 55.8 7.22 |
татары русские | 183 (46.80) 208 (53.2) | 279 (54.28) 235 (45.72) |
индекс курения для курильщиков (пачка-лет) Mm | 35.521.41 | 13.198.087 |
статус курения: курильщики некурящие | 267 (68.29) 124 (31.71) | 348 (67.70) 166 (32.3) |
Всего | 391 | 514 |
Исследовано 257 детей с хроническими бронхолегочными заболеваниями (ХБЛЗ) в возрасте от 2 до 18 лет (табл. 2). Диагноз был верифицирован на основании комплексного клинико-инструментального исследования лечащими врачами Республиканской детской клинической больницы и клиники БГМУ г. Уфы. Более, чем у половины детей (65% обследованных) установлены необратимые морфологические изменения в виде деформации бронхов и пневмосклероза в одном или нескольких сегментах, сопровождающиеся рецидивами воспаления в бронхах и легочной ткани, что трактовалось как затяжная (рецидивирующая) пневмония (J18 по МКБ-10).
Таблица 2
Характеристика детей ХБЛЗ и группы сравнения
ХБЛЗ n (%) | Контроль n (%) | |
мальчики девочки | 123 (47.86) 134 (52.14) | 153 (45.67) 182 (54.33) |
возраст (лет) Mm | 12.44 3.96 | 15.513.44 |
русские татары | 94 (36.58) 163 (63.42) | 104 (31.05) 231 (68.95) |
всего | 257 | 335 |
У 35% обследованных с диагнозом хронический бронхит (J41 по МКБ-10) отмечались распространенные поражения бронхов, сопровождающиеся продуктивным кашлем (более 3-х месяцев в год), постоянными разнокалиберными хрипами в легких (в течение нескольких месяцев) при наличии 2-3 обострений заболевания в год на протяжении не менее 2-х лет подряд. Критерии исключения из исследования: дети с бронхиальной астмой, аллергическими заболеваниями и муковисцидозом. В качестве группы сравнения обследованы практически здоровые дети и подростки в возрасте от 4 до 18 лет. Все индивиды контрольной группы от 1 до 3-х раз в год переносили острые респираторные заболевания, не приводившие к развитию затяжного воспалительного процесса и госпитализации по поводу бронхолегочных заболеваний.
Диагностика профессионального хронического бронхита (ПХБ) (табл. 3) проводилась сотрудниками ФГУН УфНИИ Медицины труда и экологии человека Роспотребнадзора г. Уфы
Таблица 3
Характеристика больных с профессиональным хроническим бронхитом
и группы сравнения
ПХБ n (%) | Здоровые рабочие n (%) | |
мужчин женщин | 87 (71.31) 35 (28.69) | 158 (95.18) 8 (4.82) |
возраст (лет) Mm | 55.699.46 | 46.016.99 |
русские татары | 45 (36.89) 77 (63.11) | 84 (50.60) 82 (49.40) |
стаж работы во вредных условиях труда, лет | 21.708.27 | 16.576.77 |
превышение ПДК, раз | 7.5213.17 | н/д |
индекс курения у курильщиков, (PY) Mm | 19.9712.32 | 18.6212.55 |
статус курения: курильщики некурящие | 37 (30.33) 85 (69.67) | 104 (62.65) 62 (37.35) |
Всего | 122 | 166 |
в соответствии с нормативами, изложенными в сборнике нормативно-методической документации по профпатологии и гигиене труда от 1998 и 2005 гг. Исследовался характер действующего этиологического фактора и выполняемой работы, особенности клинической формы заболевания, конкретных санитарно-гигиенических условий производственной среды и трудового процесса, учитывался стаж работы во вредных условиях труда. По профессиональному составу в группу больных ПХБ вошли рабочие следующих специальностей: электрогазосварщики (34.43%); рабочие горнообогатительных комбинатов (42.64%); рабочие других редких профессий (22.93%). Комплекс воздействующих вредных производственных факторов включал в себя сварочный аэрозоль, токсические вещества, силикатную, металлическую, известковую и стеклянную пыль. В качестве группы сравнения были обследованы практически здоровые высокостажированные рабочие со сходными условиями труда: электрогазосварщики, машинисты буровых и насосных установок, проходчики, крепильщики, взрывники, горные мастера. Проведенное исследование было одобрено Комитетом по этике Учреждения Российской Академии Наук Института биохимии и генетики Уфимского научного центра РАН. От всех участников исследования было получено информированное добровольное согласие на использование биологического материала в планируемых исследованиях.
Методы исследования. ДНК выделяли из лейкоцитов периферической крови методом фенольно-хлороформной экстракции (Mathew, 1984). Нами были изучены функциональные полиморфные локусы генов, продукты которых потенциально могут быть вовлечены в патогенез заболеваний органов дыхания. Гены и локусы для генотипирования выбирались с учетом распространения аллельных вариантов, функционального характера полиморфного локуса, ассоциированного с изменениями уровня экспрессии гена, активности, количества соответствующего белкового продукта (табл. 4). Полиморфные локусы генов анализировали методом полимеразной цепной реакции (ПЦР) или рестрикционного анализа.
Статистическая обработка результатов. Для количественных показателей вычислялись средние величины и их стандартные ошибки (Mm). Сравнение групп проводили с помощью непараметрического U-теста Манна-Уитни. При сравнении частот качественных признаков использовался критерий χ2 Пирсона. Статистическая обработка проводилась в программе Biostat (Primer of Biostatistics version 4.03). Оценка частоты редкого аллеля (minor allele frequency, MAF), соответствие распределения частот генотипов равновесию ХардиЦВайнберга (тест χ2), анализ ассоциации с использованием базового аллельного теста и расчета
Таблица 4
Гены-кандидаты хронических заболеваний органов дыхания
№ | Ген | полиморфный локус | HGVS Names* | Ch | RefSNP | 1 | 2 | 3 |
Гены ферментов биотрансформации ксенобиотиков | ||||||||
1 | CYP1A1 | 2454A>G, I462V | c.1384A>G p.Ile462Val | 15 | rs1048943 | |||
2 | CYP1A1 | 3798T>C | g.1759A>G | 15 | rs4646903 | |||
3 | CYP1A2 | -163C>A | g.32035C>A | 15 | rs762551 | |||
4 | CYP1A2 | -2467delT | g.29731delT | 15 | rs35694136 | |||
5 | CYP1B1 | 4326C>G, L432V | - | 2 | rs1056836 | |||
6 | CYP2A6 | large deletion | g.10817_11469del653 | 19 | rs71790353 | - | - | |
7 | CYP2E1 | -1053C>T | g.3979C>T | 10 | rs2031920 | |||
8 | CYP2F1 | c.14_15insC | c.14_15insC | 19 | rs11399890 | |||
9 | CYP2J2 | -76G>T | g.4930G>T | 1 | rs890293 | |||
10 | CYP2S1 | 13106C>T, P466L | c.1397C>T p.Pro466Leu | 19 | rs34971233 | |||
11 | CYP2S1 | 13255A>G, 3' UTR | c.*31A>G | 19 | rs338583 | |||
12 | GSTM1 | large deletion | - | 1 | ||||
13 | GSTT1 | large deletion | c.528+1_529-1del204 | 22 | rs71748309 | |||
14 | GSTP1 | 313A>G , I105V | c.313A>G p.Ile105Val | 11 | rs1695 | |||
15 | GSTP1 | 341C>T, A114V | c. 341C>T p.Ala114Val | 11 | rs1138272 | |||
16 | NQO1 | 609C>T, P187S | c.445C>T p.Pro187Ser | 16 | rs1800566 | |||
17 | NQO1 | 465C>T, R139W | - | 16 | rs1131341 | |||
18 | EPHX1 | 337T>C, Y113H | - | 1 | rs1051740 | |||
19 | EPHX1 | 415A>G, H139R | c.416A>G p.His139Arg | 1 | rs2234922 | |||
20 | NAT2 | 481C >Т, L161L | c.481C>T p.Leu161= | 8 | rs1799929 | - | - | |
21 | NAT2 | 590G>A, R197Q | c.590G>A p.Arg197Gln | 8 | rs1799930 | - | - | |
22 | NAT2 | 857G>A, G286E | c.857G>A p.Gly286Glu | 8 | rs1799931 | - | - | |
23 | UGT2B7 | 2146C>T , H268Y | c.802T>C p.Tyr268His | 4 | rs7439366 | |||
Гены ферментов антиоксидантной защиты | ||||||||
24 | SOD1 | c.239+34A>C | c.239+34A>C | 21 | rs2234694 | - | ||
25 | SOD3 | 691C>G, R231G | c.691C>G p.Arg231Gly | 4 | rs1799895 | - | ||
26 | CAT | -262C>T | g.4760C>T | 11 | rs1001179 | |||
27 | CAT | 1167 C>T, D389D | c.1167C>T p.Asp389= | 11 | rs769217 | |||
28 | GPX1 | 599C>T, P197L | c.599C>T p.Pro200Leu | 3 | rs1050450 |
Продолжение таблицы на следующей странице
Продолжение таблицы 4
№ | Ген | полиморфный локус | HGVS Names | Ch | RefSNP | 1 | 2 | 3 |
Гены ферментов системы протеолиза | ||||||||
29 | MMP1 | -1607 G>GG | g.3471delG | 11 | rs1799750 | |||
30 | MMP1 | -519 A>G | g.4757A>G | 11 | rs494379 | |||
31 | MMP12 | -82 A>G | g.102745791A>G | 11 | rs2276109 | |||
32 | MMP3 | -1171 5А>6А | g.3395_3396insA | 11 | rs35068180 | |||
33 | MMP9 | -1562C>T | g.3430C>T | 20 | rs3918242 | |||
34 | MMP9 | 2660A>G, R279Q | c.836A>G p.Gln279Arg | 20 | rs17576 | |||
35 | MMP2 | -735C>T | g.4297C>T | 16 | rs2285053 | |||
36 | ADAM33 | 13491C>G | c.*449G>C | 20 | rs2787094 | |||
37 | ADAM33 | 12418A>G, M738T | c.2213A>G p.Met738Thr | 20 | rs2280091 | - | ||
Гены ферментов сывороточных и тканевых игибиторов протеаз | ||||||||
38 | TIMP2 | -418G>C | g.76921889C>G | 17 | rs8179090 | |||
39 | TIMP3 | -1296T>C | g.4892T>C | 22 | rs9619311 | |||
40 | SERPINA1 | 1237G>A | c.*224G>A | 14 | rs2073333 | |||
41 | SERPINA1 | 2313A>T, E264V | c.863A>T p.Glu288Val | 14 | rs17580 | |||
42 | SERPINA1 | 4628A>G, E342K | c.1096G>A p.Glu366Lys | 14 | rs28929474 | |||
43 | SERPINA3 | 25G>A, A15T | c.25G>A, p.Ala9Thr | 14 | rs4934 | |||
Гены медиаторов воспаления | ||||||||
44 | TNFA | -308G>A | g.4682G>A | 6 | rs1800629 | |||
45 | LTA | 252A>G | g.5438A>G | 6 | rs909253 | |||
46 | IL1B | -511C>T | g.4490T>C | 2 | rs16944 | - | ||
47 | IL1B | 3539C>T | - | 2 | rs1143634 | |||
48 | IL1RN | VNTR интрон 2 | - | 2 | rs71941886 | |||
49 | IL6 | -174G>C | g.4880C>G | 7 | rs1800795 | - | ||
50 | IL8 | -251T>A | g.74606024A>T | 4 | rs4073 | |||
51 | IL10 | -627C>A | g.4433A>C | 1 | rs1800872 | |||
52 | TLR4 | 896A>G, D299G | c.896A>G p.Asp299Gly | 9 | rs4986790 | - | ||
53 | VDBP | 1296T>G, D432E | c.1296T>G p.Asp432Glu | 4 | rs7041 | |||
54 | VDBP | 1307C>A, T436K | c.1307C>A p.Thr455Lys | 4 | rs4588 | |||
55 | TGFb | -509C>T | g.3473A>G | 19 | rs4803457 | - | - |
Примечание: Ch Цхромосома, HGVS Names - номенклатура Human Genome variation Society (2011 года), RefSNP - Reference SNP (National Center for Biotechnology Information (NCBI)), 1 - ХОБЛ, 2 - ХБЛЗ у детей, 3 - ПХБ
показателя отношения шансов OR для редкого аллеля каждого анализируемого локуса, статистическую значимость различий между группами по частотам аллелей и генотипов (тест χ2 на гомогенность выборок и значение P-value для теста) проводили в программе PLINK v. 1.07 ( Учитывая этническую гетерогенность выборок, были использованы тесты Кохрана-Мантеля-Хензеля и Бреслоу-Дэя, а также тест на гомогенность отношения шансов (OR), предусмотренные в пакете программ PLINK v. 1.07, которые позволяли оценить ассоциацию аллелей полиморфных локусов в стратифицированной выборке.
огистическая регрессия использовалась для выявления ассоциации полиморфных локусов и гаплотипов сцепленных локусов с учетом количественных (возраст, возраст манифестации заболевания, индекс курения, стаж работы) и бинарных признаков (пол, этническая принадлежность, статус курения), вводимых в уравнение регрессии в качестве независимых переменных. Статистически значимыми считали различия при P<0.05. Для минимизации статистической ошибки первого типа вводили поправку на множественность сравнений (поправка Бонферрони) и получали новое значение Рcor. Проверку гипотез о значимости независимых факторов проводили на основе коэффициента t-статистики и уровня значимости (P-value) для коэффициента t. Экспоненту отдельного коэффициента регрессии (beta) интерпретировали как отношение шансов (OR) для логистической модели с расчетом 95% доверительного интервала. Проверку гипотезы о существенности построенной модели в целом с учетом всех переменных проводили на основании теста отношения правдоподобия и его значимости Padj (уровень значимости с учетом всех переменных, включенных в уравнение регрессии). Для выбора лучшей модели использовали информационный критерий Акайке (AIC). Для каждого локуса из статистически значимых выбирались модели с наименьшим значением AIC. Логистическая регрессия использовалась также для оценки взаимодействия полиморфного локуса и фактора окружающей среды (статус и индекс курения, стаж работы во вредных производствах). Все расчеты проводили в программах PLINK v. 1.07 и SNPStats ( Частоты гаплотипов и величину стандартизированного коэффициента неравновесия по сцеплению D', а также различия в частотах гаплотипов между группами рассчитывали в программе Haploview 4.2 (
medical-and-population-genetics/haploview/haploview), анализ ассоциации гаплотипов сцепленных локусов в аддитивной модели проводили в программе SNPStats. Для исследования межгенных взаимодействий использовали метод MDR (Multifactor Dimensionality Reduction) (Moore J.H. et al., 2002; а также его модификацию GMDR (Lou X et al., 2007;
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
1. Анализ распространенности полиморфных вариантов генов-кандидатов в этнических группах Республики Башкортостан.
Принимая во внимание генетически неоднородную этническую структуру населения Республики Башкортостан, нами проведен анализ частот генотипов, аллелей и гаплотипов полиморфных локусов генов ферментов биотрансформации ксенобиотиков, антиоксидантной защиты, медиаторов воспаления, протеолитических ферментов и ингибиторов протеолиза в трёх этнических группах (русских, татар и башкир). Полученные данные сравнивались с различными популяциями мира. Полиморфные локусы изученных генов в исследуемых этнических группах характеризовались широким аллельным разнообразием. Статистически значимые различия в распределении частот аллелей между тремя изученными группами были выявлены по локусам CYP1A1 (3798T>C) (P=0.007), CYP1A1 (2454A>G) (P=0.0001), CYP1A2 (-2467delT) (P=0.017), CYP2A6 (del) (P=0.044), CYP2E1 (-1053C>T) (P=0.013), SOD1 (c.239+34A>C) (P=0.032), IL1RN (VNTR) (P=0.005), IL10 (-627C>A) (P=0.003), VDBP (1296T>G) (P=0.0001), связанные с увеличением частоты редких аллелей в этнической группе башкир. Башкиры существенно отличались по распределению частот полиморфных вариантов генов CYP1A1 (2454A>G), CYP1A1 (3798>C), CYP2A6 (del), CYP2E1 (-1053C>T), GSTP1 (313A>G), SOD1 (c.239+34A>C), MMP3 (-1117 5A>6A), IL10 (-627C>A), VDBP (1296T>G), IL1RN (VNTR) от этнической группы русских, тогда как между русскими и татарами выявлено сходство по большинству изученных локусов. По частотам аллелей исследованных полиморфных локусов русские и татары не отличались от других европеоидных популяций (Hapmap, тогда как c представителями монголоидной расы были выявлены значимые различия. С другой стороны, для этнической группы башкир различия в распределении частот генотипов и аллелей оказались статистически значимыми как с популяциями европеоидов и негроидов, так и с монголоидными популяциями (Hapmap,
2. Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием хронической обструктивной болезни легких
Результаты выявленных статистически значимых ассоциаций полиморфных вариантов генов-кандидатов с учетом таких факторов, как этническая принадлежность, возраст, статус и индекс курения, вводимых в уравнение регрессии в качестве независимых переменных резюмированы в табл. 5. Определена ассоциация гетерозиготного генотипа локуса NQO1 (465C>T) (23.27% против 10.58% в контроле, Padj=0.000001, Pcor=0.000002, OR=2.53) и гаплотипа C-T гена NQO1 по локусам 609C>T и 465C>T (Padj=0.00014 в целом, OR=2.77 95%CI 1.75-4.40) с развитием ХОБЛ. НАД(Ф)Н-хинон оксидоредуктаза, кодируемая геном NQO1, предотвращает образование свободных радикалов семихинона и активных форм кислорода (АФК), а также активирует нитрозамины и гетероциклические амины. Замена 465C>T приводит к изменению активности фермента (Park et al., 2011). Гомозиготный по частому аллелю генотип 6A6A локуса MMP3 (-1171 5А>6А) является маркером риска заболевания (Padj=0.0047, OR=2.42), в то время как гетерозиготный генотип маркирует устойчивость (OR=0.41). Гомозиготы по редкому аллелю 5A5A в нашей выборке выявлены не были. MMP3 является одним из ключевых ферментов семейства MMP, так как активирует и индуцирует синтез других ММР, таких как ММР1 и ММР9, он также взаимодействует с цитокинами, факторами роста, разлагает протеогликаны, ламинины и коллаген IV типа, образующие внеклеточный матрикс легочной ткани (Santus et al., 2009). Полиморфный вариант MMP3 (-1171 5А>6А) представляет инсерционно-делеционный полиморфизм, который изменяет сайт связывания транскрипционного фактора NFB, что влияет на экспрессию ММР3 (Juran et al., 2011). Ассоциация с частым генотипом может объясняться возможным сцеплением с другим функциональным локусом гена MMP3. Частота делеции гена CYP2A6 в группе больных была снижена (6.14% против 17.14% в контроле; Padj=0.00001, OR=0.27), маркером развития ХОБЛ является частый генотип без делеции (OR=3.16). Главной функцией фермента CYP2A6 является катализ процесса окисления никотина - основного соединения, который вызывает и поддерживает табачную зависимость. Изученная нами протяжённая делеция гена блокирует процесс экспрессии гена, что продлевает период полураспада никотина от двух до одиннадцати часов (Nelson, 2009). Вероятно, у носителей делеционного аллеля в гомозиготном состоянии удлиняется время метаболизма никотина, и это объясняет протективную роль данного генетического маркера в отношении развития целого комплекса заболеваний, напрямую связанных с курением, таких как ХОБЛ, рак легкого и т.д. (Gyamfi et al., 2005). Генотип GG локуса CYP2S1 (13255A>G) является маркером устойчивости к развитию ХОБЛ (Padj=0.0006, Pcor=0.0012, OR=0.26).
Таблица 5
окус | Генотип риска | с ХОБЛ в целом (N=827) | с ХОБЛ в группе русских (N=406) | с ХОБЛ в группе татар (N=421) | |||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
NQO1 (465C>T) | CT | 0.000001 | 2.53 (1.72-3.86) | 0.00001 | 3.74 (1.99-7.02) | 0.021 | 1.88 (1.09-3.25) |
CYP2A6 (del) | n | 0.00001 | 3.16 (1.89-5.33) | 0.0035 | 2.55(1.30-4.98) | 0.00001 | 5.48 (2.59-11.60) |
CYP2S1 (13255A>G) | GG | 0.0006 | 0.26 (0.11-0.59) | 0.013 | 0.28 (0.10-0.80) | 0.0064 | 0.23 (0.07-0.74) |
MMP3 (-1171 5А>6А) | 6A6A | 0.0047 | 2.42 (1.42-4.13) | 0.0069 | 2.86 (1.45-5.37) | 0.023 | 2.37 (1.09-5.12) |
EPHX1 (337T>C) | TC | 0.0022 | 1.72 (1.21-2.44) | 0.021 | 1.62 (1.07-2.45) | 0.012 | 1.84 (1.14-2.97) |
CAT (-262C>T) | СT | 0.016 | 0.40 (0.19-0.87) | 0.033 | 0.39 (0.16-0.95) | 0.041 | 0.39 (0.15-1.01) |
CAT (-262C>T) | TT | 0.011 | 1.61 (1.12-2.33) | - | - | - | - |
GSTP1 (313A>G) | AA | 0.0053 | 1.51 (1.13 - 1.99) | - | - | 0.018 | 1.69 (1.08-2.46) |
GSTP1( 313A>G) | AG | 0.0043 | 0.65 (0.49-0.88) | 0.04 | 0.65 (0.42-0.98) | 0.043 | 0.65 (0.43-0.99) |
CYP1A1 (3798T>C) | TC | 0.028 | 0.63 (0.42-0.96) | 0.0083 | 0.41 (0.21-0.81) | - | - |
IL6 (-174G>C) | GG | 0.0042 | 1.38 (1.12-1.84) | - | - | 0.018 | 1.41 (1.12-2.11) |
IL6 (-174G>C) | CC | 0.032 | 0.58 (0.35-0.96) | 0.0012 | 0.29 (0.13-0.64) | ||
CYP1A2 (-163A>C) | AA | 0.0062 | 1.47 (1.12-1.95) | - | - | 0.042 | 0.62 (0.39-0.98) |
ADAM33 (12418A>G) | AA | 0.018 | 1.49 (1.07-2.09) | - | - | 0.032 | 1.62 (1.11-2.48) |
SOD3 (691C>G) | аддитивная CC(0) CG(1) GG (2) | 0.021 | 1.81 (1.08-3.02) | - | - | 0.03 | 2.00 (1.05-3.80) |
CYP1A2 (-2467delT) | del/del | - | - | - | - | 0.014 | 2.42 (1.03-5.72) |
VDBP (1307C>A) | аддитивная CC(0) CA(1) AA (2) | 0.015 | 1.33 (1.06-1.68) | - | - | 0.061 | - |
Ассоциация полиморфных вариантов генов-кандидатов с развитием ХОБЛ
Примечание: Здесь и далее в таблицах и тексте: Padj - значимость для теста отношения правдоподобия лог-регрессионной модели с учетом возраста, этнической принадлежности, статуса и индекса курения. ORadj - отношение шансов с учетом всех факторов. *аддитивная модель на дозу редкого аллеля - увеличение дозы редкого аллеля в ряду: гомозигота по частому аллелю (0) - гетерозигота (1) - гомозигота по редкому аллелю (2). N - количество индивидов (больные - контроль), включенных в анализ ассоциации.
Выявлено значимое увеличение доли гетерозигот по локусу EPHX1 (337T>C) (45.01% против 35.7%, в контроле; Padj=0.0022, Pcor=0.0044, OR=1.72) и частоты гаплотипа C-A гена EPHX1 по локусам 337T>C и 415A>G среди больных ХОБЛ (Padj=0.00056 в целом, OR=2.16 95%CI 1.48-3.14). Показана ассоциация гомозиготного по редкому аллелю генотипа TT полиморфного локуса CAT (-262C>T) (Padj=0.011, Pcor=0.022, OR=1.61) и гаплотипа C-T гена CAT по локусам -262C>T и 1167C>T (P=0.032 в целом, OR=1.43 95%CI 1.08 - 1.90) с ХОБЛ. Гомозиготный по частому аллелю генотип AA локуса GSTP1 (313A>G) связан с развитием заболевания (Padj=0.0053, Pcor=0.0106, OR=1.51), тогда как гетерозиготный генотип маркирует устойчивость к формированию ХОБЛ (Padj=0.0043, Pcor=0.009, OR=0.62 95%CI 0.44-0.89). Частота генотипа AA полиморфного локуса CYP1A2 (-163A>C) в группе больных значимо выше, чем в контроле (50.9% против 41.3%, P=0.0062, Pcor=0.0124, OR=1.47), гетерозиготный генотип локуса CYP1A2 (-163A>C) является маркером устойчивости к развитию ХОБЛ (P=0.0044, Pcor=0.0088, OR=0.67 95%CI 0.51-0.88). Гомозиготы по редкому аллелю GG локуса SOD3 (691C>G) были выявлены только в группе больных с частотой 4.21% (Padj=0.0008, Pcor=0.0016); статистически значимая ассоциация локуса SOD3 (691C>G) показана в аддитивной модели (Padj=0.021, OR=1.81). Выявлено повышение доли гетерозигот по локусу ADAM33 (12418A>G) среди здоровых индивидов (29.76% против 21.99% в группе больных, Padj=0.022, Pcor=0.044, OR=0.61 95%CI 0.40-0.93); риск развития ХОБЛ возрастает у гомозигот по частому аллелю AA (P=0.018, Pcor=0.036, OR=1.49), а гаплотип G-G гена ADAM33 по локусам 13491 C>G и 12418 A>G маркирует устойчивость к развитию ХОБЛ (Padj=0.0071 в целом, OR=0.40 95%CI 0.23-0.68). Значимые ассоциации с ХОБЛ были выявлены для локуса VDBP (1307C>A) в аддитивной модели (P=0.015, OR=1.33) и гаплотипа C-T (GC*F) гена VDBP по локусам 1307C>A и 1296T>G (P=0.00036 в целом, OR=0.71 95%CI 0.56 - 0.89). Маркером риска развития ХОБЛ является гомозиготный по частому аллелю генотип GG локуса IL6 (-174G>C) (Padj=0.0042, Pcor=0.0084, OR=1.38).
Анализ ассоциации полиморфных вариантов генов-кандидатов с ХОБЛ
с учетом этнической принадлежности
Тесты Кохрана-Мантеля-Хензеля и Бреслоу-Дэя, а также тест на гомогенность отношения шансов в стратифицированных по этнической принадлежности выборках (русских и татар) показали гетерогенность выявленных ассоциаций по некоторым локусам. Анализ ассоциации в дифференцированных по этнической принадлежности группах проводился с учетом статуса и индекса курения, возраста и пола обследованных. В табл. 5 представлены полученные данные. Определены ключевые генетические маркеры, ассоциации с которыми были подтверждены как в этнической группе русских, так и татар NQO1 (465C>T), MMP3 (-1171 5А>6А), CYP2A6 (del), CYP2S1 (13255A>G), EPHX1 (337T>C), CAT (-262C>T), GSTP1 (313A>G). Необходимо отметить, что для локусов NQO1 (465C>T) и MMP3 (-1171 5А>6А) уровень значимости был ниже в этнической группе русских, тогда как для татар более информативными были локусы CYP2A6 (del) и CYP2S1 (13255A>G). Ассоциация локусов CYP1A1 (3798T>C) и LTA (252A>G) с развитием ХОБЛ выявлена только в этнической группе русских. В выборке татар ассоциация с развитием заболевания показана для локусов SOD3 (691C>G), IL6 (-174G>C), ADAM33 (12418A>G), CYP1A2 (-163C>A), CYP1A2 (-2467delT) и гаплотипа C-T гена VDBP по локусам 1307C>A и 1296T>G (P=0.017 в целом, OR=0.67 95%CI (0.49 - 0.92).
Анализ ассоциации полиморфных вариантов генов-кандидатов с ХОБЛ
с учетом пола обследованных
Большая часть больных ХОБЛ в нашем исследовании представлена мужчинами (80%), однако выявление специфических генетических маркеров заболевания, связанных с полом, является важным этапом молекулярно- генетического анализа наследственной предрасположенности к заболеванию. Известно, что ХОБЛ, при прочих равных условиях, более тяжело протекает у женщин (Чучалин, 2009), что указывает на генетическую модификацию предрасположенности к заболеванию у мужчин и женщин. Анализ ассоциации полиморфных локусов генов-кандидатов с ХОБЛ в дифференцированных по полу выборках проводился с учетом возраста, этнической принадлежности, статуса и индекса курения. Результаты представлены в табл. 6. Выявлены генетические маркеры, показавшие ассоциацию с ХОБЛ как в выборке мужчин, так и в выборке женщин - это полиморфные варианты локусов NQO1 (465C>T), MMP3 (-1171 5А>6А), CYP2A6 (del), CYP2S1 (13255A>G), EPHX1 (337T>C), CYP1A2 (-163C>A). С развитием ХОБЛ у мужчин ассоциируют полиморфные варианты локусов GSTP1 (313A>G), VDBP (1307C>A), ADAM33 (12418A>G), CAT (-262C>), CAT (1167C>T). Ассоциации с заболеванием у женщин установлены для локусов IL6
(-174G>C), CYP1A1 (3798T>C), SERPINA3 (25G>A). Для исключения возможного влияния на выявленные ассоциации этнической гетерогенности нами был проведен анализ ассоциации полиморфных локусов генов-кандидатов с развитием ХОБЛ у мужчин и женщин в стратифицированных по этнической принадлежности группах. Соотношение русских и татар среди мужчин (323 (49.01%) : 336 (50.99%)), среди женщин (83 (49.40%) : 85 (50.60%)). Ассоциация локуса NQO1 (465C>T) с низким уровнем значимости (P=0.0007) как у мужчин, так и у женщин, выявлена только у русских. Локус CYP2A6 (del) ассоциировал с ХОБЛ у мужчин обеих этнических групп и у женщин татарской этнической принадлежности.
Таблица 6
Ассоциация полиморфных вариантов локусов генов-кандидатов с развитием ХОБЛ у мужчин и женщин
окус | Генотип риска | с ХОБЛ в группе мужчин (N=659) | с ХОБЛ в группе женщин (N=168) | ||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
NQO1( 465C>T) | CT | 0.0004 | 2.19 (1.41-3.41) | 0.0005 | 5.01 (1.79-14.02) |
CYP2A6 ( del) | n | 0.0001 | 3.57 (2.03-6.26) | 0.0004 | 4.36 (1.63-11.46) |
CYP2S1 (13255A>G) | GG | 0.0001 | 0.22 (0.09-0.51) | 0.011 | 0.17 (0.03-0.80) |
MMP3 (-1171 5А>6А) | 6A6A | 0.041 | 1.87 (1.01-3.45) | 0.0029 | 5.17 (1.64-16.3) |
EPHX1 (337T>C) | TC | 0.022 | 1.61 (1.07-2.43) | 0.022 | 2.09 (1.10-3.95) |
CYP1A2 (-163A>C) | AC | 0.03 | 0.71 (0.52-0.97) | 0.026 | 1.85 (1.05-3.25) |
CAT (-262C>T) | CT | 0.0005 | 2.15 (1.39-3.33) | - | - |
GSTP1(313A>G) | AA | 0.0014 | 1.72 (1.26-2.40) | - | - |
GSTP1 (313A>G) | AG | 0.0004 | 0.47 (0.30-0.72) | - | - |
CAT (1167C>T) | аддитивная CC(0) CT(1) TT(2) | 0.018 | 1.39 (1.06-1.83) | - | - |
VDBP (1307C>A) | аддитивная CC(0) CA (1)AA (2) | 0.012 | 1.39 (1.07-1.79) | - | - |
ADAM33 (12418A>G) | AA | 0.016 | 1.57 (1.08-2.28) | - | - |
CYP1A1 (3798T>C) | TT | - | - | 0.0053 | 2.82 (1.29-6.18) |
CYP1A1 (3798T>C) | TC | - | - | 0.0045 | 0.31 (0.14-0.72) |
SERPINA3 (25G>A) | аддитивная GG(0)AG(1)AA(2) | - | - | 0.016 | 0.58 (0.37-0.91) |
IL6 (-174G>C) | GG | - | - | 0.0091 | 2.05 (1.12-3.87) |
Примечание: Padj - значимость для теста отношения правдоподобия лог-регрессионной модели с учетом возраста, этнической принадлежности, статуса и индекса курения. ORadj - отношение шансов с учетом всех факторов. N - количество индивидов (больные - контроль), включенных в анализ ассоциации.
окус MMP3 (-1171 5А>6А) ассоциировал с ХОБЛ у женщин, как русских, так и татар, у мужчин значимых ассоциаций в дифференцированных по этнической принадлежности группах выявлено не было. Ассоциации с полиморфными вариантами гена GSTP1 (313A>G) были выявлены у мужчин русских и татар. Локус CYP2S1 (13255A>G) ассоциировал с развитием ХОБЛ у мужчин в этнической группе татар. Полиморфный вариант локуса CYP1A1 (3798T>C) ассоциировал с развитием ХОБЛ у русских женщин. Локус CAT (-262C>T) подтвердили свою значимость в только группе мужчин татар. Локус IL6 (-174G>C) значимо ассоциировал с ХОБЛ только в группе женщин татар. Таким образом, особенности ассоциации полиморфных локусов генов-кандидатов с развитием ХОБЛ у мужчин и женщин, дифференцированных по этнической принадлежности, в целом отражает тот спектр генетических маркеров, который был определен отдельно в группах, подразделенных по этнической принадлежности и полу. Выявленные гендерные различия могут быть связаны с особенностями гормональной регуляции и спецификой воспалительных реакций у мужчин и женщин. Так, только у женщин риск развития ХОБЛ связан с двумя генами цитохромов P450 (CYP1A2 и CYP1A1), кодирующими ферменты метаболизма ксенобиотиков, а также эндогенной активации и деградации половых гормонов (Nelson, 2009) и геном IL6 (-174G>C), кодирующим ключевой цитокин острой фазы воспаления, оказывающего влияние на транскрипционный уровень С-реактивного белка (СРБ). Известно, что СРБ, в свою очередь, ассоциирует со снижением функций легких у здоровых лиц и у некурящих больных ХОБЛ (Shaaban et al., 2006; He et al., 2009). Для мужчин специфичной является ассоциация с генами, кодирующими ферменты антиоксидантной защиты (GSTP1, CAT, SOD3). Свободнорадикальное окисление играет первостепенную роль в патогенезе ХОБЛ. Большая часть мужчин с ХОБЛ являлись курильщиками с длительным стажем курения, поэтому, возможно, полученные ассоциации связаны с проявлением эффекта генов в сочетании с фактором курения.
Анализ ассоциации полиморфных вариантов генов-кандидатов с ХОБЛ с учетом курения, выявление ген-средовых взаимодействий
Комплексная оценка влияния генетических и средовых факторов на риск возникновения хронических заболеваний органов дыхания позволяет не только установить ключевые ген-средовые взаимодействия, формирующие основу предрасположенности к болезни, но и понять механизмы, посредством которых факторы внешней среды способны спровоцировать патологические изменения (Hunter, 2005; Полоников, 2006). Наибольший интерес для нашего исследования представляли факторы, влияние которых на организм непосредственно связано с поражением дыхательной системы. Это, прежде всего, статус и индекс курения, а также стаж работы во вредных производствах, характеризующий интенсивность воздействия пылевых и токсических веществ на органы дыхания. Нами были изучены ген-средовые взаимодействия статуса и индекса курения с полиморфными вариантами генов-кандидатов при формировании ХОБЛ. Статистически значимые взаимодействия со статусом курения определены для локусов ADAM33 (12418A>G) (Pinteract=0.026), TIMP3 (-1296T>C) (Pinteract=0.044), VDBP (1296T>G) (Pinteract=0.05104), CYP1A1 (2454A>G) (Pinteract=0.047), CYP2F1 (c.14_15insC) (Pinteract=0.049), GPX1 599C>T (Pinteract=0.041). Взаимодействия с индексом курения установлены для локуса SOD3 (691C>G) (Pinteract=0.05402). Анализ взаимодействий генотип-среда проводился также путем сопоставления величин отношения шансов, рассчитанных для полиморфных вариантов генов-кандидатов, в группах больных и контроля, сформированных по критерию наличия или отсутствия влияния средового фактора - курильщиков и некурящих. Результаты анализа представлены в табл. 7.
Таблица 7
Ассоциация полиморфных вариантов генов-кандидатов с развитием ХОБЛ в группах, дифференцированных по статусу курения
окус | генотип риска | с ХОБЛ в группе курильщиков (N=560) | с ХОБЛ в группе некурящих (N=267) | ||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
NQO1(465C>T) | CT | 0.00001 | 2.71 (1.62-4.53) | 0.0081 | 2.60 (1.26-5.38) |
CYP2A6 (del) | n | 0.0001 | 2.82 (1.55-5.13) | 0.00001 | 6.24 (3.01-10.75) |
CYP2S11(3255A>G) | GG | 0.0001 | 0.14 (0.05-0.44) | 0.018 | 0.34 (0.11-0.99) |
EPHX1 (337T>C) | TC | 0.02 | 1.51 (1.06-2.14) | 0.052 | - |
CAT (-262C>T) | TT | 0.0057 | 0.33 (0.14-0.76) | - | - |
CAT (-262C>T) | CT | 0.015 | 1.57 (1.09-2.25) | - | - |
GSTP1 (313A>G) | AA | 0.0039 | 1.66 (1.17-2.35) | - | - |
GSTP1 (313A>G) | AG | 0.0085 | 0.62 (0.43-0.89) | - | - |
VDBP (1307C>A) | аддитивная CC(0) CA(1) AA(2) | 0.02 | 1.39 (1.05-1.82) | - | - |
ILRN (VNTR) | 2/2 | 0.028 | 1.94 (1.06-3.56) | - | - |
MMP3 (-1171 5А>6А) | 6A6A | - | - | 0.0023 | 3.97 (1.53-10.2) |
IL6 (-174G>C) | GG | - | - | 0.042 | 1.74 (1.12-2.52) |
CYP1A2 (-163A>C) | AA | - | - | 0.0098 | 1.90 (1.16-3.11) |
CYP1A2 (-163A>C) | AC | - | - | 0.013 | 0.56 (0.34-0.94) |
ADAM33 (12418A>G) | AA | - | - | 0.0023 | 2.41 (1.34-4.32) |
ADAM33 (12418A>G) | AG | - | - | 0.0021 | 0.37 (0.19-0.71) |
TIMP3 (-1296T>C) | TT | - | - | 0.011 | 1.59 (1.10-2.53) |
CYP1A1 (3798T>C) | TT | - | - | 0.033 | 1.83 (1.04-3.22) |
CYP1A1 (3798T>C) | TC | - | - | 0.018 | 0.50 (0.28-0.90) |
Примечание: Padj значимость для теста отношения правдоподобия лог-регрессионной модели с учетом возраста, пола и этнической принадлежности. ORadj - отношение шансов с учетом всех факторов. N - количество индивидов (больные - контроль), включенных в анализ ассоциации.
Только три локуса ассоциировали с развитием ХОБЛ как у курильщиков, так и у некурящих индивидов - NQO1 (465C>T), CYP2A6 (del), CYP2S1 (13255A>G). Выявлены специфические генетические маркеры развития ХОБЛ в группах, дифференцированных по статусу курения. У курильщиков - это локусы EPHX1 (337T>C), CAT (-262C>T), GSTP1 (313A>G), VDBP (1307C>A), IL1RN (VNTR). У некурящих спектр выявленных ассоциаций представлен локусами MMP3 (-1171 5А>6А), IL6 (-174G>C), CYP1A2 (-163C>A), ADAM33 (12418A>G), TIMP3 (-1296T>C), CYP1A1 (3798T>C). Для курильщиков формирование хронической обструктивной болезни легких детерминировано, прежде всего, генами биотрансформации ксенобиотиков и антиоксидантой защиты; для некурящих - генами медиаторов воспаления и системы протеолиза-антипротеолиза. Таким образом, проявление эффектов полиморфных вариантов отдельных генов при развитии ХОБЛ существенно зависит от характера и интенсивности влияния факторов среды.
Анализ ассоциации полиморфных вариантов генов-кандидатов с возрастом манифестации ХОБЛ
Для выявления ассоциации полиморфных вариантов изученных генов с возрастом манифестации сравнивали две группы больных ХОБЛ, с ранней и поздней манифестацией заболевания, как между собой, так и отдельно каждую группу - с контрольной. Группы больных не различались по соотношению мужчин и женщин, по статусу и индексу курения, по тяжести заболевания. В группе ХОБЛ с ранней манифестацией частота генотипа СС локуса UGT2B7 (2146C>T) была значимо выше, чем в группе с поздней манифестацией (44.10% против 30.40%, P=0.0005, OR=2.83 95%CI 1.73-4.61). Среди больных с манифестацией заболевания после 40 лет частота гомозиготного генотипа по редкому аллелю AA локуса NAT2 (590G>A) составляет 11.21% против 3.37% у больных с ранней манифестацией ХОБЛ (P=0.018, OR=0.27 95%CI 0.06-0.99), а доля гетерозигот по локусу CYP1A1 (2454A>G) достигала 10% против 2.90% (P=0.0097, OR=0.27 95%CI 0.06-0.96). Сравнительный анализ с контролем показал, что для больных с ранней манифестацией значимые ассоциации были получены с полиморфными вариантами генов NQO1 (465C>T) генотип CT (P=0.0018, OR=3.31 95%CI 1.47-7.44), CYP2A6 (del) генотип без делеции (P=0.001, OR=3.06 95%CI 1.77-5.30), EPHX1 (337T>C) в аддитивной модели (P=0.0018, OR=1.84 95%CI 1.26-2.68), SOD3 (691C>G) в аддитивной модели (P=0.015, OR=2.06 95%CI 1.17-3.62), VDBP (1307C>A) генотип CA (P=0.0062, OR=1.85 95%CI 1.19-2.88), IL1RN (VNTR) генотип 2/4 (P=0.0024, OR=0.48 95%CI 0.29-0.78). Выявленные ассоциации полиморфных локусов изученных генов-кандидатов с возрастом манифестации ХОБЛ свидетельствуют о том, что данные генетические маркеры могут выступать в качестве модификаторов начала проявления первых симптомов ХОБЛ.
Анализ ассоциации полиморфных вариантов генов-кандидатов
с тяжестью ХОБЛ
С целью выявления генетических маркеров, ассоциированных с тяжестью заболевания, сравнивали больных с тяжелым течением ХОБЛ и больных со средне тяжелым течением заболевания. В нашем исследовании у 75.7% больных диагностировано тяжелое течение. Различий между группами больных по среднему возрасту, возрасту манифестации заболевания, соотношению мужчин и женщин, статусу и индексу курения выявлено не было. Ассоциация с тяжелым течением ХОБЛ получена для генотипа AA локуса LTA (252A>G) (Padj=0.037, OR=1.63 95%CI 1.12-2.34). Гетерозиготный генотип локуса LTA (252A>G) чаще встречался в группе больных со средней тяжестью ХОБЛ (Padj=0.0073, OR=0.53 95%CI 0.33-0.84). Маркерами тяжелого течения ХОБЛ являются генотип GG локуса IL6 (-174G>C) (Padj=0.015, OR=1.82 95%CI 1.11-2.96) и генотип СС локуса CYP2E1 (-1053C>T) (Padj=0.015, OR=2.90 95%CI 1.20-6.97). Ассоциация с тяжелой ХОБЛ была выявлена для локуса MMP9 (-1562C>T) (P=0.023 OR=1.88 в аддитивной модели), так как частота гомозигот по редкому аллелю (2% против 0%) и гетерозигот (23.6% против 15.8%) выше у больных с более тяжелым течением ХОБЛ. Таким образом, в результате анализа ассоциации всех изученных локусов генов-кандидатов выявлены генетические маркеры, которые связаны с более тяжелым течением ХОБЛ. Два локуса (IL6 и LTA) кодируют цитокины, играющие ключевую роль в патогенезе ХОБЛ и развитии хронического воспаления (Fishman et al., 1998). MMP9 участвует в процессах ремоделирования и деградации легочной ткани при ХОБЛ. Аллель T локуса -1562C>T гена MMP9, с которым выявлена ассоциация, снижает способность промотора связываться с репрессором транскрипции, таким образом, повышая экспрессию гена, что может играть важную роль в быстром развитии тяжелых осложнений заболевания.
Генетические маркеры развития эмфиземы легких при ХОБЛ
Согласно современным представлениям эмфизема легких (ЭЛ) является обязательным компонентом ХОБЛ. Однако степень выраженности эмфиземы различается у разных больных. Сравнительный анализ групп больных ХОБЛ с ЭЛ и без ЭЛ по таким показателям, как возраст больного, возраст манифестации ХОБЛ, пол, статус и индекс курения выявил достоверные различия между группами. Так, среди больных с ЭЛ было больше мужчин (P=0.0001) и курильщиков (P=0.0001). Индекс курения в этой группе был выше (P=0.0001). Средний возраст манифестации ХОБЛ у больных с ЭЛ также был больше (P=0.024). В связи с этим анализ ассоциации проводился с учетом статуса и интенсивности курения, возраста манифестации ХОБЛ и гендерных различий больных. Ассоциация с развитием ЭЛ была выявлена для полиморфных вариантов локусов UGT2B7 (2146C>T) (Padj=0.0072, OR=1.50 95%CI 1.11-2.04 в аддитивной модели), ADAM33 (13491 C>G) (Padj=0.013, OR=1.74 95%CI 1.43-2.79 для генотипа GG); IL1RN (VNTR) (Padj=0.042, OR=1.59 95%CI 1.11-2.57 для гетерозиготного генотипа); VDBP (1296T>G) (Padj=0.013, OR=0.58 95%CI 0.38-0.89) для генотипа TG). Таким образом, формирование ЭЛ детерминировано комплексом внешних факторов (курением и его интенсивностью), возрастам манифестации, полом и генетической конституцией индивида.
3. Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием хронических бронхолегочных заболеваний у детей
Результаты выявленных статистически значимых ассоциаций полиморфных вариантов генов-кандидатов с развитием ХБЛЗ у детей с учетом этнической принадлежности и пола резюмированы в табл. 8. Ассоциация с развитием ХБЛЗ с наименьшим уровнем значимости выявлена для локуса CYP2F1 (c.14_15insC). Маркером риска является гомозиготный по частому аллелю генотип без инсерции (Padj=0.000001, Pcor=0.000002, OR=3.16). Локусы CYP1A1 (2454A>G) и CYP1A1 (3798T>C) ассоциируют с ХБЛЗ в аддитивной модели (Padj=0.001, Pcor=0.002, OR=2.31 и Padj=0.0003, Pcor=0.0006, OR=1.72, соответственно), частоты гомозигот по редкому аллелю и гетерозигот по данным локусам в группе больных значимо выше; гаплотипы A-С и G-Cапо локусам 2454A>G и 3798T>C гена CYP1A1 являются маркерами риска заболевания (Padj=0.00022 в целом, OR=1.62 95%CI 1.16-2.26 и OR=2.28 95%CI 1.23-4.21). Среди больных детей доля индивидов с делецией гена GSTT1 составляет 35% против 21.3% в группе здоровых (Padj=0.0003, OR=1.98). Ассоциация с развитием ХБЛЗ у детей выявлена c гомозиготным по частому аллелю генотипом AA локуса GSTP1 (313A>G) (Padj=0.0046, Pcor=0.0092, OR=1.65), тогда как гетерозиготный генотип является маркером устойчивости (Padj=0.0012, Pcor=0.0024, OR=0.55). Локус NQO1 (465C>T) ассоциирует с развитием ХБЛЗ в аддитивной модели (Padj=0.006, OR=1.89), доля гетерозигот и гомозигот по редкому аллелю T среди больных выше, чем в контроле. Гомозиготный по частому аллелю генотип СС локуса CAT (-262C>T) чаще встречается среди больных детей (71.6% против 57.7% в контроле; Padj=0.0009, Pcor=0.0018, OR=1.84), а гаплотип T-C по локусам -262C>T и 1167C>T гена CAT является маркером устойчивости к заболеванию (Padj=0.0054 в целом, OR=0.65 95%CI 0.46-0.91). Маркером риска развития ХБЛЗ является гомозиготный по частому аллелю генотип AA локуса MMP12 (-82A>G) (Padj=0.007, OR=1.83). Ассоциация с ХБЛЗ получена с генотипом 6A6A локуса MMP3 (-1171 5А>6А) (Padj=0.0013, OR=2.75).
Таблица 8
Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием ХБЛЗ у детей
окус | генотип риска | c ХБЛЗ в целом (N=558) | с ХБЛЗ у русских (N=188) | с ХБЛЗ у татар (N=370) | |||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
CYP2F1 (c.14_15insC) | n/n | 0.000001 | 3.16 (2.10-4.77) | 0.0039 | 2.71 (1.35-5.46) | 0.000001 | 3.41 (2.04-5.67) |
CAT (-262C>T) | CC | 0.0009 | 1.84 (1.22-2.65) | 0.0063 | 2.41 (1.27-4.58) | 0.031 | 1.61 (1.04-2.49) |
IL6 (-174G>C) | аддитивная GG(0) GC(1) CC(2) | 0.0009 | 0.64 (0.49-0.84) | 0.035 | 0.62 (0.40-0.97) | 0.012 | 0.66 (0.47-0.91 |
IL6 (-174G>C) | GG | 0.0051 | 1.67 (1.16-2.39) | 0.019 | 2.06 (1.12-3.79) | - | - |
GSTP1 (313A>G) | AA | 0.0046 | 1.65 (1.17-2.34) | 0.0001 | 3.47 (1.85-6.54) | - | - |
GSTP1 (313A>G) | AG | 0.0012 | 0.55 (0.38-0.79) | 0.00001 | 0.25 (0.13-0.49) | - | - |
MMP3 (-1171 5А>6А) | 6A6A | 0.0013 | 2.75 (1.43-5.31) | 0.0015 | 5.98 (1.60-21.45) | - | - |
TIMP3 (-1296T>C) | TT | 0.033 | 1.48(1.11-2.01) | 0.02 | 2.05 (1.11-3.76) | - | - |
SERPINA1 (1237G>A) | GA | - | - | 0.018 | 3.28 (1.15-9.39) | - | - |
GSTT1 (del) | del | 0.0003 | 1.98 (1.36-2.89) | - | - | 0.0011 | 2.18 (1.36-3.48) |
MMP9 (2660A>G) | AG | 0.0004 | 0.52 (0.36-0.75) | - | - | 0.0003 | 0.44 (0.28-0.69) |
MMP9 (2660A>G) | AA | 0.017 | 1.54 (1.08-2.21) | - | - | 0.038 | 1.59 (1.02-2.47) |
CYP1A1 (3798T>C) | аддитивная TT(0) TC(1) TT(2) | 0.0003 | 1.72 (1.28-2.33) | - | - | 0.0004 | 1.92 (1.33-2.76) |
CYP1A1 (2454A>G) | аддитивная AA(0) AG(1) GG(2) | 0.001 | 2.31 (1.38-3.88) | - | - | 0.0002 | 3.14 (1.62-6.09) |
NQO1 (465C>T) | аддитивная CC(0) CT(1) TT(2) | 0.006 | 1.89 (1.19-3.02) | - | - | 0.036 | 1.80 (1.03-3.16) |
MMP12 (-82A>G) | AA | 0.007 | 1.83 (1.16-2.89) | - | - | 0.024 | 1.92 (1.07-3.46) |
VDBP (1296T>G) | GG | 0.0096 | 0.56 (0.35-0.87) | - | - | 0.02 | 1.64 (1.08-2.50) |
ILRN (VNTR) | 2/2 | 0.018 | 0.37 (0.16-0.89) | - | - | 0.007 | 0.22 (0.06-0.78) |
VDBP (1307C>A) | CA | - | - | - | - | 0.0092 | 1.77 (1.15-2.74) |
MMP1 (-1607 G>GG) | аддитивная GG/GG(0) GG/G(1) G/GC(2) | - | - | - | - | 0.0027 | 1.53 (1.16-2.03) |
Примечание: Padj - значимость для теста отношения правдоподобия лог-регрессионной модели с учетом пола и этнической принадлежности, ORadj - отношение шансов с учетом всех факторов
Риск развития заболевания связан с генотипом AA локуса MMP9 (2660A>G) (Padj=0.017, Pcor=0.034, OR=1.54), тогда как гетерозиготный генотип является протективным в отношении развития ХБЛЗ (Padj=0.0004, Pcor=0.0008, OR=0.52). Локус IL6 (-174G>C) ассоциирует с ХБЛЗ в аддитивной модели (Padj=0.0009, OR=0.64), поскольку доля гомозигот по редкому аллелю СС и гетерозигот выше среди здоровых детей; маркером риска ХБЛЗ у детей является гомозиготный по частому аллелю генотип GG (Padj=0.0051, Pcor=0.01, OR=1.67). Частота генотипа TT локуса TIMP3 (-1296T>C) в группе больных детей была выше (63.42% против 54.44% в контроле; Padj=0.033, Pcor=0.066, OR=1.48). Генотип 2/2 VNTR-локуса гена IL1RN и гомозиготный по редкому аллелю генотип GG локуса VDBP (1296T>G) чаще встречаются среди здоровых детей (Padj=0.018, Pcor=0.036, OR=0.37 и Padj=0.0096, Pcor=0.0019, OR=0.56, соответственно).
Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием ХБЛЗ у детей с учетом этнической принадлежности
Тесты Кохрана-Мантеля-Хензеля и Бреслоу-Дэя, а также тест на гомогенность отношения шансов в стратифицированных по этнической принадлежности выборках выявили гетерогенность ассоциации по некоторым локусам. Результаты анализа в дифференцированных по этнической принадлежности группах представлены в табл. 8. По трем локусам ассоциации подтверждены как в этнической группе русских, так и татар - CYP2F1 (c.14_15insC), CAT (-262C>T), IL6 (-174G>C). В этнической группе русских ассоциации были выявлены для локусов GSTP1 (313A>G), MMP3 (-1171 5А>6А), TIMP3 (-1296T>C), SERPINA1 (1237G>A).. В этнической группе татар ассоциация с ХБЛЗ у детей показана для локусов CYP1A1 (2454A>G), CYP1A1 (3798T>C), GSTT (del), NQO1 (465C>T), MMP12 (-82A>G), MMP9 (2660A>G), IL1RN (VNTR), VDBP (1296T>G), VDBP (1307C>A), MMP1 (-1607G>GG). У русских самые значимые ассоциации (P<0.001) с развитием ХБЛЗ у детей получены с полиморфными вариантами локуса GSTP1 (313A>G), а у татар с полиморфными вариантами генов CYP2F1 (c.14_15insC), CYP1A1 (2454A>G), CYP1A1 (3798T>C) и MMP9 (2660A>G).
Анализ ассоциации полиморфных вариантов генов-кандидатов с ХБЛЗ у детей с учетом пола обследованных
Установлено, что как для мальчиков, так и для девочек ассоциация с развитием ХБЛЗ сохраняется для локусов CYP2F1 (c.14_15insC), CYP1A1 (2454A>G), GSTP1 (313A>G), MMP9 (2660A>G), IL6 (-174G>C) (табл. 9). Выявлены генетические маркеры, связанные с развитием заболевания у мальчиков - это локусы GSTT1 (del), NQO1 (465C>T), MMP12 (-82A>G), TIMP3(-1296T>C), SOD3 (691C>G). Для девочек ассоциации были выявлены с полиморфными вариантами локусов CYP1A1 (3798T>C), CAT (-262C>T), MMP3 (-1171 5А>6А), IL1RN (VNTR), VDBP (1296T>G), CYP2S1 (13106C>T). Для исключения возможного влияния на выявленные ассоциации этнической гетерогенности нами был проведен анализ в группах мальчиков и девочек, стратифицированных по этнической принадлежности.
Таблица 9
Ассоциация полиморфных вариантов локусов генов-кандидатов с развитием ХБЛЗ у детей в зависимости от пола обследованных
окус | Генотип риска | с ХБЛЗ у мальчиков (N=264) | с ХБЛЗ у девочек (N=294) | ||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
CYP2F1 (c.14_15insC) | n/n | 0.00001 | 3.95 (2.16-7.23) | 0.0006 | 2.59 (1.47-4.55) |
MMP9 (2660A>G) | AG | 0.0086 | 0.49 (0.29-0.84) | 0.014 | 0.54 (0.33-0.89) |
GSTP1 (313A>G) | AG | 0.011 | 0.52 (0.31-0.87) | 0.04 | 0.59 (0.36-0.99) |
CYP1A1 (2454A>G) | аддитивная AA(0) AG(1) GG(2) | 0.014 | 2.55 (1.17-5.57) | 0.025 | 2.11 (1.06-4.19) |
IL6 (-174G>C) | аддитивная GG(0) GC(1) CC(2) | 0.018 | 0.64 (0.43-0.93) | 0.018 | 0.65 (0.45-0.93) |
CYP1A1 (2454A>G) | AG | 0.0028 | 3.42 (1.46-8.04) | - | - |
GSTT1 (de) | del | 0.001 | 2.45 (1.43-4.21) | - | - |
NQO1 (465C>T) | аддитивная CC(0) CT(1) TT(2) | 0.036 | 1.93 (1.02-3.64) | - | - |
MMP12 (-82A>G) | AA | 0.04 | 1.96 (1.02-3.78) | - | - |
SOD3 (691C>G) | аддитивная CC(0) CG(1) GG(2) | 0.007 | 2.63 (1.26-5.48) | - | - |
TIMP3 (-1296T>C) | TT | 0.011 | 2.95 (1.56-3.31) | - | - |
CYP1A1 (3798T>C) | аддитивная TT(0) TC(1) TT(2) | - | - | 0.002 | 1.93 (1.26-2.96) |
CAT (-262C>T) | CC | - | - | 0.0001 | 2.78 (1.67-4.62) |
GSTP1 (313A>G) | AA | - | - | 0.019 | 1.79 (1.09-2.93) |
MMP3(-1171 5А>6А) | 6A6A | - | - | 0.006 | 3.16 (1.31-7.62) |
VDBP (1296T>G) | TG | - | - | 0.021 | 1.74 (1.09-2.80) |
VDBP (1296T>G) | GG | - | - | 0.0012 | 0.33 (0.16-0.67) |
CYP2S1 (13106C>T) | аддитивная CC(0) CT (1) TT (2) | - | - | 0.0046 | 2.32 (1.27-4.24) |
Ассоциации с полиморфными вариантами гена GSTP1 (313A>G) были значимы только у русских как мальчиков, так и девочек. Делеция гена GSTT1 ассоциировала с развитием ХБЛЗ у мальчиков - татар по этнической принадлежности с уровнем значимости (P=0.0004). Ассоциации полиморфных локусов CYP1A1 (3798T>C), CYP1A1 (2454A>G), IL6 (-174G>C), MMP9 (2660A>G) также были значимы только у татар. Для локуса CYP2S1 (13106C>T) ассоциация была подтверждена только у русских девочек (P=0.014). Локус CAT (-262C>T) был связан с развитием ХБЛЗ у девочек обеих этнических групп, но уровень значимости (P=0.0001) был выше у русских. В результате проведённого анализа выявлен широкий спектр генетических маркеров развития ХБЛЗ у детей. Учитывая, что в данной группе больных отсутствует направленное и длительное влияние фактора среды, можно предположить, что риск развития заболеваний органов дыхания у детей обусловлен, в первую очередь, генетической компонентой, представленной генами всех изученных систем. Наиболее значимой является ассоциация с полиморфным вариантом гена CYP2F1 (c.14_15insC). Ген CYP2F1 (c.14_15insC) экспрессируется только в легочной ткани, функция его, помимо метаболизма токсических соединений и эндогенных молекул (холестерола, стероидных гормонов), до сих пор не изучена (Tornell et al., 2007). Важнейшими маркерами риска развития ХБЛЗ являются также полиморфные локусы генов, кодирующих протеолитические ферменты и ингибиторы протеолиза. В детском возрасте идет активный процесс роста и развития легких, а влияние средовых факторов связано, в основном, с частыми инфекциями и неспецифичным воздействием неблагоприятных экологических условий. Вполне вероятно, что у детей с определенной генетической конституцией существуют нарушения регуляции воспалительного ответа, что влияет на физиологические процессы роста легких, приводит к структурной перестройке бронхиального дерева и патологическому повреждению ткани легкого с последующим формированием эмфиземы и пневмосклероза.
4. Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием профессионального хронического бронхита
Анализ частот генотипов полиморфных локусов генов-кандидатов в группе здоровых высокостажированных рабочих выявил отклонения от равновесия Харди-Вайнберга, связанные с увеличением доли гомозигот по частому аллелю по локусам CAT (1167C>T) (25.30% против 35.95%; P=0.0003), MMP1 (-519A>G) (19.88% против 26.83%; P=0.002), CYP1A2 (-2467delT) (17.47% против 26.00%; P=0.00015), ADAM33 (13491C>G) (35.54% против 49.91%; P=0.00019), MMP9 (2660A>G) (23.49% против 29.23%; P=0.015), TIMP3 (-1296T>C) (34.34% против 45.99%; P=0.0012). Это может быть обусловлено отбором против гетерозигот и гомозигот по редким аллелям. Возможно, в связи с постоянным, направленным, происходящим в течение длительного времени воздействием определенных факторов производственной среды, среди работающих происходит отбор против аллелей, ассоциированных с высоким риском заболевания. Сравнительный анализ этой группы с популяционным контролем выявил значимые различия, в основном за счет увеличения доли гомозигот по частым аллелям по локусам генов CYP1A2 (-2467delT) (75.9% против 58% в контроле), CYP2E1 (-1053C>T) (94.6% против 88.8%), CYP2S1 (13106C>T) (94.6% против 81.2%), MMP1 (-1607 G>GG) (63.2% против 38.7%), MMP9 (2660A>G) (70.5% против 40.2%), VDBP (1296 T>G) (63.9% против 33.7%). Группа высокостажированных здоровых рабочих представляет интерес в плане выявления генетических маркеров устойчивости к воздействию вредных факторов производственной среды.
Результаты статистически значимых ассоциаций полиморфных вариантов генов-кандидатов с развитием ПХБ с учетом таких факторов, как этническая принадлежность, стаж работы во вредных производствах, статус и индекс курения представлены в табл. 10. В группе больных ПХБ выявлено одновременное увеличение доли редких гомозигот GG и гетерозигот полиморфного локуса VDBP 1296T>G (16.4% против 7.2% и 43.4% против 28.9%, соответственно), поэтому наиболее информативной моделью, описывающей риск развития заболевания у рабочих, является аддитивная модель (Padj=0.00005, OR=2.06); гаплотип С-G (GC*1S) по локусам 1307C>A и 1296T>G гена VDBP является маркером риска развития ПХБ (Padj=0.00066 в целом, OR=2.54 95%CI 1.60-4.04). Показана ассоциация гетерозиготного генотипа локуса MMP1 (519A>G) (Padj=0.0001, Pcor=0.0002, OR=2.52) с развитием профессионального бронхита. Генотип CG локуса ADAM33 (13491C>G) является маркером развития ПХБ у рабочих (Padj=0.0004, Pcor=0.0008, OR=2.52). Выявлена ассоциация гомозиготного по частому аллелю генотипа СС локуса NQO1 (465C>T) (Padj=0.001, Pcor=0.002, OR=4.11) и гаплотипа С-T по локусам 609C>T и 465C>T гена NQO1 (Padj=0.00087 в целом, OR=0.31 95%CI 0.12 - 0.80) с ПХБ. Полиморфный вариант СС локуса UGT2B7 (2146C>T) был связан с развитием заболевания у рабочих (Padj=0.0021, Pcor=0.0042, OR=2.34). Доля гетерозигот по локусу CYP1A2 (-2467delT) среди больных ПХБ была выше (30.33% против 17.47% в контроле; Padj=0.008, Pcor=0.016, OR=2.17). Гетерозиготный генотип локуса CYP2F1 (c.14_15insC) чаще выявляется у здоровых рабочих (Padj=0.0048, Pcor=0.01, OR=0.44). Маркерами устойчивости к развитию заболевания также являются гетерозиготный генотип локуса GSTP1 (313A>G) (Padj=0.0086, Pcor=0.017, OR=0.50) и гаплотип С-T гена CAT по локусам -262C>T и 1167C>T (Padj=0.0074 в целом, OR=0.54 95%CI 0.33 - 0.88).
Таблица 10
Ассоциация полиморфных вариантов генов-кандидатов с развитием ПХБ
окус | модель | c ПХБ в целом (N=288) | c ПХБ у руcских (N=129) | c ПХБ у татар (N=151) | |||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
MMP1 (-519 A>G) | AG | 0.0001 | 2.52 (1.36-4.69) | 0.040 | 2.28 (1.04-5.00) | 0.0092 | 2.62 (1.25-5.51) |
NQO1 (465C>T) | CC | 0.001 | 4.11 (1.91-8.85) | 0.0064 | 4.60 (1.10-16.63) | 0.0037 | 3.81 (1.44-10.11) |
CYP1A2 (-2467delT) | de l/ n | 0.008 | 2.17 (1.20-3.91) | - | - | - | - |
CAT (1167C>T) | TT | 0.015 | 0.27 (0.09-0.83) | - | - | - | - |
MMP9 (2660A>G) | AG | 0.015 | 0.27 (0.09-0.83) | - | - | - | - |
GSTP1 (313A>G) | AG | 0.0086 | 0.50 (0.29-0.84) | 0.0053 | 0.32 (0.14-0.74) | - | - |
CYP2F1 (c.14_15insC) | w / ins | 0.0048 | 0.44 (0.24-0.79) | 0.0007 | 0.20 (0.07-0.57) | - | - |
CYP2F1 (c.14_15insC) | w /w | - | - | 0.0011 | 4.07 (1.63-10.16) | - | - |
ADAM33 (13491C>G) | CG | 0.0004 | 2.52 (1.40-4.52) | 0.023 | 2.34 (1.12-4.90) | - | - |
TIMP3 (-1296T>C) | аддитивная TT (0) TC (1) CC (2) | - | - | 0.01 | 1.83 (1.15-2.93) | - | - |
IL8 (-251T>A) | аддитивная TT (0) TA (1) AA (2) | 0.014 | 1.51 (1.02-2.23) | - | - | - | - |
IL8 (-251T>A) | AA | - | - | 0.016 | 2.70 (1.18-6.19) | ||
VDBP (1296T>G) | аддитивная TT (0) TG (1) GG (2) | 0.00005 | 2.06 (1.36-3.13) | - | - | 0.0001 | 3.22 (1.86-5.59) |
UGT2B7 (2146C>T) | CC | 0.0021 | 2.33 (1.35-4.04) | - | - | 0.016 | 2.44 (1.16-5.16) |
Примечание: Здесь и далее в тексте: Padj - значимость для теста отношения правдоподобия лог-регрессионной модели с учетом этнической принадлежности, стажа работы, статуса и индекса курения. ORadj - отношение шансов с учетом всех факторов. N - количество индивидов (больные - контроль), включенных в анализ ассоциации.
Анализ ассоциации полиморфных вариантов генов-кандидатов
с развитием ПХБ с учетом этнической принадлежности
Тесты Кохрана-Мантеля-Хензеля, Бреслоу-Дэя, тест на гомогенность отношения шансов в стратифицированных выборках показали, что показатели отношения шансов не отличаются в выборках русских и татар для локусов UGT2B7 (2146C>T), MMP1 (-519A>G) и NQO1 (465C>T). Наряду с этим для локуса VDBP (1296T>G) показаны значимые различия по ассоциации в группе русских и татар (PBD=0.04602). Дальнейший анализ в дифференцированных по этнической принадлежности группах (см. табл. 10) подтвердил ассоциацию с развитием ПХБ следующих локусов в этнических группах русских и татар: NQO1 (465C>T) и MMP1 (-519A>G). Только в группе русских значимые ассоциации показаны для локусов GSTP1 (313A>G), гаплотипа С-T гена CAT по локусам -262C>T и 1167C>T (P=0.048 в целом, OR=0.49 95%CI 0.26 0.94), ADAM33 (13491C>G), CYP2F1 (c.14_15insC), TIMP3 (-1296T>C). В этнической группе татар самые значимые ассоциации выявлены с локусом VDBP (1296T>G) и гаплотипом С-G (GC*1S) по локусам 1307C>A и 1296T>G гена VDBP (P=0.00001 в целом, OR=4.86 95% CI 2.46-9.62). Маркерами риска развития ПХБ в данной этнической группе также являются полиморфные варианты локусов IL8 (-251T>A), UGT2B7 (2146C>T).
Анализ ассоциации генов-кандидатов с развитием ПХБ с учетом курения и стажа работы, выявление ген-средовых взаимодействий
Нами изучены взаимодействия статуса и индекса курения с полиморфными вариантами генов-кандидатов при развитии ПХБ. Статистически значимые взаимодействия со статусом курения определены для локусов UGT2B7 (2146C>T) (Pinteract=0.015), EPHX1 (415A>G) (Pinteract=0.04), GPX1 (599C>T) (Pinteract=0.037); с индексом курения для локуса CYP2F1 (c.14_15insC) (Pinteract=0.05). Анализ взаимодействий генотип-среда проводился также путем сопоставления величин отношения шансов, рассчитанных для полиморфных локусов генов-кандидатов в группах, сформированных по критерию наличия или отсутствия влияния средового фактора - у курильщиков и некурящих. Результаты анализа представлены в табл. 11. Ассоциация с развитием ПХБ как в группе курильщиков, так и некурящих выявлена для локусов VDBP (1296T>G), MMP1 (-519A>G), ADAM33 (13491C>G), GSTP1 (313A>G), UGT2B7 (2146C>T). Только у курильщиков ассоциация показана с полиморфными вариантами локусов CYP1A2 (-2467delT), CYP2F1 (c.14_15insC). Для некурящих ассоциации с развитием заболевания были выявлены с полиморфными вариантами генов MMP9 (2660A>G), NQO1 (465C>T), TIMP3 (-1296T>C). Среди внешнесредовых факторов, оказывающих влияние на формирование профессиональной патологии, большое значение имеет стаж работы во вредных производствах. Стаж работы является своего рода показателем длительности воздействия определенного фактора на организм работающего и несомненно, играет важную роль в формирование ПХБ. Нами была проведена оценка влияние этого фактора на полученные ассоциации с генетическими локусами, а также выявлены ген-средовые взаимодействия. Статистически значимые взаимодействия со стажем работы показаны для локусов IL1RN (VNTR) (Pinteraction=0.02656), VDBP (1307C>A) (Pinteraction=0.02258) и CYP1A1 (3798T>C) (Pinteraction=0.01657).
Таблица 11
Анализ ассоциации полиморфных вариантов генов-кандидатов с развитием ПХБ в группах, дифференцированных по статусу курения
окус | модель | с ПХБ у курильщиков (N=141) | с ПХБ у некурящих (N=147) | ||
Padj | ORadj (CI95%) | Padj | ORadj (CI95%) | ||
VDBP (1296T>G) | аддитивная TT (0) TG (1) GG (2) | 0.028 | 1.88 (1.07-3.31) | 0.0013 | 2.23 (1.33-3.73) |
MMP1 (-519 A>G) | AG | 0.025 | 2.54 (1.13-5.70) | 0.017 | 2.53 (1.15-5.57) |
ADAM33 (13491C>G) | CG | 0.014 | 2.62 (1.20-5.71) | 0.026 | 2.22 (1.09-4.52) |
GSTP1 313A>G | AG | 0.049 | 0.45 (0.20-1.02) | 0.05 | 0.51 (0.26-1.00) |
UGT2B7 (2146C>T) | аддитивная CC (0) CT (1) TT (2) | 0.03 | 0.53 (0.29-0.95) | - | - |
UGT2B7 (2146C>T) | CC | - | - | 0.0029 | 3.62 (1.42-9.02) |
CYP1A2 (-2467delT) | Del / n | 0.031 | 2.64 (1.11-6.30) | - | - |
CYP2F1 (c.14_15insC) | w / ins | 0.035 | 0.41 (0.17-0.98) | - | - |
MMP9 (2660A>G) | AG | - | - | 0.0067 | 3.27 (1.31-8.17) |
TIMP3 (-1296T>C) | СС | - | - | 0.041 | 2.29 (1.01-5.20) |
NQO1( 465C>T) | CC | - | - | 0.0002 | 7.04 (2.22-22.33) |
Полученные данные показывают, что при реализации ПХБ большую роль играют такие внешнесредовые факторы, как стаж работы, уровень ПДК, тип воздействия (пыль или токсическое вещество). Профессиональный бронхит при воздействии промышленных аэрозолей возникает преимущественно у лиц с определенной генетической конституцией и является следствием взаимодействия генетических и средовых факторов. Вклад генетической компоненты в развитие профессионального заболевания, в отличие от ХОБЛ и ХБЛЗ у детей, ограничен небольшим количеством полиморфных маркеров. Наиболее значимые ассоциации получены нами с полиморфными вариантами гена VDBP. Показано существование отбора против гетерозигот и гомозигот по редким аллелям в группе здоровых высокостажированных рабочих. В работах Спицына В.А. (2008) был отмечен эффект поляризации больных и здоровых рабочих по генетическим маркерам, когда в группе здоровых высокостажированных рабочих наблюдается снижение частоты гомозигот по редким аллелям. В нашей работе в группе рабочих, подвергающихся воздействию промышленных аэрозолей, но не имеющих ПХБ, выявлено увеличение доли гетерозиготных генотипов по локусам GSTP1 (313A>G), CYP2F1 (c.14_15insC), MMP9 (2660A>G) и частоты гаплотипа С-T гена CAT по локусам -262C>T и 1167C>T, которые могут рассматриваться как маркеры устойчивости организма к воздействию промышленных факторов.
5. Анализ межгенных взаимодействий при формировании наследственной предрасположенности к хроническим заболеваниям органов дыхания
Важным этапом ассоциативных исследований является выявление взаимодействий между генами с целью определения таких комбинаций генных локусов, которые имеют наибольшую патогенетическую значимость для развития болезни. В результате анализа межгенных взаимодействий полиморфных вариантов генов биотрансформации ксенобиотиков, антиоксидантной защиты, протеолиза, антипротеолиза и медиаторов воспаления определены модели, детерминирующие риск развития хронической обструктивной болезни легких, бронхолегочных заболеваний у детей и профессионального бронхита с учетом пола, этнической принадлежности индивидов и статуса курения. В табл. 12 представлены наиболее значимые 4-x локусные модели межгенных взаимодействий, выявленные при алгоритме выборочного поиска для каждой группы больных. На рис. 1 приведены результаты кластерного анализа межгенных взаимодействий при различных хронических заболеваниях органов дыхания.
Были выявлены три статистически значимые 4-х локусные модели, взаимодействие полиморфных вариантов которых лежит в основе предрасположенности к развитию ХОБЛ. Одна из них CYP2A6 (del), NQO1 (465C>T), ADAM33 (12418A>G), IL1B (3539C>T) состоит их двух независимых кластеров. Первый включает гены CYP2A6 (del) и IL1B (3539C>T), для которых характерно тесное взаимодействие и взаимное усиление эффектов (синергизм). Второй кластер - это гены NQO1 (465C>T) и ADAM33 (12418A>G), между которыми также выявлено сильное взаимодействие, но их влияние на развитие ХОБЛ дублируется. Анализ данной модели в таблицах сопряженности, представляющей собой комбинации всех возможных вариантов 4-х локусной модели, выявил несколько протективных и рисковых комбинаций. В дифференцированных по этнической принадлежности группах выявлено несколько моделей, наиболее значимые из которых представлены в табл. 12. Для русских значимой является модель, включающая следующие полиморфные локусы MMP1(-519 A>G), CYP1A2 (-163C>A), NQO1 (465C>T), ADAM33 (12418A>G).
Таблица 12
Модели межгенных взаимодействий при хронических
заболеваниях органов дыхания
Группа | Tr. Bal. Acc | Test. Bal. Acc | Sign Test (P) | Se | Sp | CVC | Pre. |
ХОБЛ в целом | CYP2A6 (del), NQO1 (465C>T), ADAM33 (12418A>G), IL1B (3539C>T) | ||||||
0.7412 | 0.6254 | 0.001 | 0.5627 | 0.9091 | 10/10 | 0.8609 | |
ХОБЛ русские | MMP1(-519 A>G), CYP1A2 (-163C>A), NQO1 (465C>T), ADAM33 (12418A>G) | ||||||
0.8402 | 0.8196 | 0.001 | 0.7981 | 0.8824 | 10/10 | 0.8715 | |
ХОБЛ татары | GSTP1 (313A>G), CAT (-262C>T), ADAM33 (12418A>G), IL6 (-174G>C) | ||||||
0.7732 | 0.6695 | 0.001 | 0.7760 | 0.7632 | 10/10 | 0.7662 | |
ПХБ в целом | NQO1 (465C>T), UGT2B7 (802C>T), MMP1 (-1607G>GG), ADAM33 (13491C>G) | ||||||
0.7159 | 0.6158 | 0.010 | 0.7623 | 0.6621 | 10/10 | 0.6929 | |
ПХБ русские | NQO1 (465C>T), UGT2B7 (802C>T), MMP1 (-1607 G>GG), MMP1 (-519 A>G) | ||||||
0.8461 | 0.7116 | 0.001 | 0.9211 | 0.7703 | 10/10 | 0.8046 | |
ПХБ татары | MMP1(-1607 G>GG), CYP1B1 (4326C>G), SERPINA3 (25G>A), VDBP (1296T>G) | ||||||
0.8836 | 0.8154 | 0.001 | 0.8289 | 0.9324 | 10/10 | 0.9246 | |
ХБЛЗ в целом | CYP2F1( c.14_15insC), NQO1 (465C>T), MMP12 (-82 A>G), VDBP (1296T>G) | ||||||
0.6878 | 0.615 | 0.001 | 0.8093 | 0.56 | 10/10 | 0.8128 | |
ХБЛЗ русские | GSTT (del), GSTP1 (313A>G), CAT (-262C>T), MMP3 (-1171 5А>6А) | ||||||
0.7636 | 0.695 | 0.001 | 0.7471 | 0.7742 | 10/10 | 0.9028 | |
ХБЛЗ татары | CYP2F1 (c.14_15insC), NQO1 (465C>T), MMP12 (-82 A>G), VDBP (1296T>G) | ||||||
0.7308 | 0.6972 | 0.001 | 0.8188 | 0.6377 | 10/10 | 0.8299 |
Примечание Tr.Bal.Acc. - тренировочная сбалансированная точность, Test.Bal.Acc. - тестируемая сбалансированная точность, Sign Test (P) - тест на значимость, Se. - чувствительность, Sp. - специфичность, CVС - повторяемость результата, Pre. (Precision) - точность модели.
ХОБЛ
ПХБ
ХБЛЗ у детей
Рис. 1. Дендрограммы межгенных взаимодействий при хронических заболеваниях органов дыхания. Здесь и далее на дендрограммах (рис. 2, 3) короткие линии указывают на сильное взаимодействие генных локусов; длинные - на слабую связь; красным и оранжевым цветом указывается синергизм, т.е. взаимное усиление эффектов между локусами; синим и зеленым - дублирование эффектов между локусами; коричневый цвет указывает на независимость эффектов отдельных локусов.
Кластерный анализ позволил выявить три кластера, при этом тесное взаимодействие показано для локусов CYP1A2 (-163C>A) и MMP1 (-519A>G), что, возможно, объясняется взаимным усилением экспрессии на уровне сигнального каскада арилгидрокарбонового рецептора (AhR-pathway) (Villano et al., 2006; Ishida et al., 2009). Локусы NQO1 (465C>T) и ADAM33 (12418A>G) представляют собой гены с независимым эффектом. В этнической группе татар развитие ХОБЛ определено взаимодействием локусов GSTP1 (313A>G), CAT (-262C>T), ADAM33 (12418A>G), IL6 (-174G>C). Модель представлена двумя независимыми кластерами: GSTP1 (313A>G) и IL6 (-174G>C) с дублирующими эффектами, и кластера взаимно усиливающих друг друга генов CAT (-262C>T) и ADAM33 (12418A>G). Анализ межгенных взаимодействий в группах, дифференцированных по статусу курения, позволил выделить специфические для курильщиков гены, определяющие развитие ХОБЛ. Одна из значимых комбинаций представлена генами GSTP1 (313A>G), MMP1 (-1607G>GG), IL8 (-251T>A), VDBP (1296T>G). Сбалансированная точность (Bal.Acc.) данной модели составила 0.87, чувствительность (Se) - 0.78, специфичность (Sp) - 0.94, воспроизводимость модели (CV Consistency) - 10/10, значимость P=0.001. Кластерный анализ вывил сильное взаимодействие между генами GSTP1 (313A>G) и NQO1 (465C>T), и генами IL8 (-251T>A) и VDBP (1296T>G), составляющими два отдельных кластера с геном MMP1 (-1607G>GG) (рис. 2). Такая схема взаимодействия отражает патогенез ХОБЛ - воспалительный компонент и взаимная активация через метаболические пути биотрансформации продуктов сигаретного дыма глутатион S-трансфераз, НАД(Ф)Н-хинон оксидоредуктазы 1 (NQO1) и матриксных металлопротеаз.
ХОБЛ курильщики
ХОБЛ некурящие
Рис. 2. Дендрограмма межгенных взаимодействий при формировании ХОБЛ у курильщиков и некурящих.
В группе некурящих определена модель, включающая следующие полиморфные локусы NQO1 (465C>T), ADAM33 (1241A>G), IL1B (3539C>T), VDBP (1296T>G). Сбалансированная точность модели составила 0.80, чувствительность (Se) - 0.83, специфичность (Sp) - 0.76, воспроизводимость модели (CV Consistency) - 10/10, значимость P=0.01. Кластерный анализ показывает, что эффект гена NQO1 (465C>T) не зависит от других генов, входящих в модель, тогда как три других гена образуют блок слабо взаимодействующих генов с дублирующими эффектами.
Анализ структуры межгенных взаимодействий при формировании ПХБ позволил выявить следующую модель NQO1 (465C>T), UGT2B7 (802C>T), MMP1 (-1607G>GG), ADAM33 (13491C>G). Основные статистические параметры модели приведены в табл. 12. Показано более тесное взаимодействие между локусами UGT2B7 (802C>T), NQO1 (465C>T) и MMP1 (-1607G>GG), эффекты которых синонимичны, тогда как вклад локуса ADAM33 (13491C>G) в развитие болезни не зависит от других ДНК-локусов (см. рис. 1). В рамках данной модели было выявлено несколько комбинаций генотипов повышенного и пониженного риска. Анализ межгенных взаимодействий в группах, дифференцированных по этнической принадлежности, выявил значимые комбинации ДНК-локусов, определяющих генетическую структуру предрасположенности к ПХБ у русских и татар (см. табл. 12). У русских выявлена следующая модель: NQO1 (465C>T), UGT2B7 (2146C>T), MMP1 (-1607G>GG), MMP1 (-519A>G). Структура взаимодействий между генами, входящими в модель, не отличается от таковой в общей группе. Показана тесная связь между генами NQO1 (465C>T) и MMP1(-1607G>GG). У татар лучшей моделью межгенных взаимодействий с наименьшей ошибкой предсказания является модель MMP1 (-1607G>GG), CYP1B1 (4326C>G), SERPINA3 (25G>A), VDBP (1296T>G). Кластерный анализ позволил установить тесное взаимодействие между локусами MMP1 (-1607G>GG) и VDBP (1296T>G), входящими в один кластер с локусом CYP1B1 4326C>G, тогда как влияние локуса SERPINA3 (25G>A) на риск ПХБ не зависит от других ДНК-локусов. Анализ межгенных взаимодействий в группах, дифференцированных по статусу курения, выявил специфические для курильщиков гены, определяющие развитие ПХБ в сочетании с действием факторов производственной среды. Самой лучшей моделью с наименьшей ошибкой предсказания была комбинация CYP1A2 (-2467del), MMP1 (-1607G>GG), ADAM33 (13491C>G), IL8 (-251T>A). Сбалансированная точность (Bal.Acc.) модели составила 0.81 чувствительность (Se) - 0.86, специфичность (Sp) - 0.76, воспроизводимость модели (CV Consistency) - 10/10, значимость P=0.01. Кластерный анализ выявил два независимых кластера генных локусов, определяющих развитие ПХБ у курильщиков (рис. 3). В первый кластер входят гены ADAM33 (13491C>G) и IL8 (-251T>A), суммирующий эффект взаимодействия которых намного сильнее вклада каждого отдельного локуса. Второй кластер представлен генами CYP1A2 (-2467delT) и MMP1 (-1607G>GG), сильный синергизм которых, возможно, связан с индукцией арилгидрокарбонового рецептора Ah-R. В модели присутствуют ДНК-локусы, которые были определены при анализе межгенных взаимодействий у курильщиков с ХОБЛ, что указывает на сходный механизм развития болезни у курильщиков и существование взаимодействий фактора среды (сигаретный дым) с генетическими факторами. У здоровых курильщиков и пациентов с ХОБЛ определяется повышенный уровень IL8 в бронхоальвеолярной жидкости, что позволяет предположить, что этот цитокин имеет большое значение в патогенезе хронических заболеваний легких (Rutgers et al., 2000). Уровень экспрессии ММР1 значительно возрастает при различных патологических процессах (Pearce et al., 2007).
ПХБ курильщики
ПХБ некурящие
Рис. 3. Дендрограмма межгенных взаимодействий при формировании ПХБ
у курильщиков и некурящих.
Кроме того, IL8 - хемоаттрактант нейтрофилов, являясь мишенью для ММР, активирует секрецию ММР нейтрофилами (Chakrabarti, Patel, 2005). Формирование синергизма, возможно, происходит уже на уровне взаимодействий транскрипционных факторов с участками промоторов генов цитокинов и металлопротеаз. Для некурящих развитие ПХБ связано с комбинацией ДНК-локусов, представленной тремя кластерами: TIMP3 (-1296T>C), LTA (252A>G), IL10 (-627C>A), VDBP (1292T>G). Самые сильные взаимодействия с эффектом синергизма определены между генами TIMP3 (-1296T>C) и IL10 (-627C>A), действие которых усиливается локусом LTA (252A>G), тогда как локус VDBP (1292T>G) является независимым фактором риска. У некурящих индивидов, для которых действие факторов среды связано, прежде всего, с воздействием токсических аэрозолей и промышленной пыли в рабочей зоне, развитие поражений органов дыхания детерминировано генами медиаторов воспаления и тканевым ингибитором протеаз.
Анализ межгенных взаимодействий при развитии ХБЛЗ у детей выявил 4-х локусную модель с наименьшей ошибкой предсказания, включающую локусы CYP2F1 (c.14_15insC), NQO1 (465C>T), MMP12 (-82A>G), VDBP (1296T>G). Все характеристики модели приведены в табл. 12. Кластерный анализ определил три кластера локусов с дублирующими эффектами (см. рис. 1). Показано тесное взаимодействие между локусами CYP2F1 (c.14_15insC) и MMP12 (-82A>G). Анализ в дифференцированных по этнической принадлежности выборках вывил модель межгенных взаимодействий, связанную с развитием ХБЛЗ у русских - GSTT (del), GSTP1 (313A>G), CAT (-262C>T), MMP3 (-1171 5А>6А). Все характеристики модели приведены в табл. 12. Кластерный анализ выявил сильное взаимодействие между генами GSTP1 (313A>G) и CAT (-262C>T), кодирующими ферменты 2-ой фазы биотрансформации ксенобиотиков и антиоксидантной защиты. С этими локусами слабо взаимодействует ген GSTT (del), при этом эффекты данных ДНК-локусов дублируются. Ген MMP3 (-1171 5А>6А) является независимым фактором риска ХБЛЗ у детей. Выявлено несколько моделей межгенных взаимодействий, связанных с развитием ХБЛЗ в этнической группе татар. Mодель CYP2F1 (c.14_15insC), NQO1 (465C>T), MMP12 (-82A>G), VDBP (1296T>G) характеризовалась минимальной ошибкой предсказания (0.30), высокой точностью (0.83), остальные оценки данной модели представлены в табл.12. Эта же модель является лучшей в общей группе, что связано с большой долей татар (63%) в изучаемой группе. Анализ межгенных взаимодействий при развитии ХБЛЗ у детей в группах, дифференцированных по полу, показал, что для мальчиков значимой моделью риска развития ХБЛЗ является CYP1A1 (3798T>C), CYP2F1 (c.14_15insC), GSTP1 (313A>G), VDBP (1296T>G). Сбалансированная точность (Bal.Acc.) модели составила 0.86, чувствительность (Se) - 0.77, специфичность (Sp) - 0.93, воспроизводимость модели (CV Consistency) - 10/10, значимость P=0.001. Кластерный анализ выявил два слабо взаимодействующих между собой кластера ДНК-локусов, определяющих развитие ХБЛЗ. В первый кластер вошли два локуса GSTP1 (313A>G) и VDBP (1296T>G), эффекты которых взаимно усиливаются. Второй кластер представлен двумя взаимодействующими генами CYP1A1 (3798T>C) и CYP2F1 (c.14_15insC) с дублирующими эффектами. Выявлены рисковые и протективные комбинации генотипов полиморфных локусов. Для девочек выявлена 4-х локусная модель, связанная с развитием ХБЛЗ - MMP1 (-1607G>GG), MMP9 (2660A>G), VDBP (1296T>G), VDBP (1307C>A). Сбалансированная точность (Bal.Acc.) модели составила 0.74, чувствительность (Se) - 0.62, специфичность (Sp) - 0.86, воспроизводимость модели (CV Consistency) - 10/10, значимость P=0.05. Кластерный анализ выявил три независимых кластера генных локусов. Первый представлен локусами MMP1 (-1607G>GG) и VDBP (1307C>A), для которых наблюдается явление синергизма. С данным кластером более слабо взаимодействует локус VDBP (1296T>G), что, по-видимому, можно объяснить сцеплением с локусом VDBP (1307C>A). Эффекты гена MMP9 (2660A>G) не зависят от других генов.
Структура межгенных взаимодействий при формировании хронической обструктивной болезни легких, бронхолегочных заболеваний у детей и профессионального бронхита имеет особенности в различных этнических группах, что подтверждает результаты анализа ассоциации отдельных локусов с заболеваниями. Известно, что при становлении этносов на протяжении значительного времени происходит их адаптация к условиям проживания в конкретном регионе и на уровне генома формируются своеобразные межгенные и ген-средовые взаимоотношения (Животовский, Хуснутдинова, 2007), что и отражается при проведении комплексных ассоциативных исследований многофакторных заболеваний в различных этнических группах. Анализ позволил установить сложный характер взаимодействий между генами ферментов биотрансформации ксенобиотиков, антиоксидантной защиты, протеолитических ферментов, ингибиторов протеолиза и медиаторов воспаления при развитии различных заболеваний органов дыхания. Была проведена оценка вклада каждого из исследуемых генов и выявлены особенности взаимодействия исследованных полиморфных локусов при развитии заболеваний. Несмотря на присутствие некоторых общих генов в моделях при различных заболеваниях органов дыхания, структура и характер взаимосвязей между локусами существенно варьировали в зависимости от конкретного заболевания, статуса курения, пола и этнической принадлежности.
Выводы
1. Сравнительный анализ трех этнических групп (русских, татар, башкир) жителей Республики Башкортостан выявил сходство этнических групп татар и русских по большинству изученных локусов. Башкиры значимо отличались по распределению частот аллелей полиморфных локусов генов CYP1A1 (2454A>G), CYP1A1 (3798>C), CYP2A6 (del), CYP2E1 (-1053C>T), GSTP1 (313A>G), SOD1 (c.239+34A>C), MMP3 (-1117 5A>6A), IL10 (-627C>A), VDBP (1296T>G), IL1RN (VNTR) от этнической группы русских.
2. Установлена ассоциация с развитием хронической обструктивной болезни легких в этнических группах русских и татар полиморфных локусов генов NQO1 (465C>T), CYP2A6 (del), CYP2S1 (13255A>G), MMP3 (-1171 5А>6А), EPHX1 (337T>C), GSTP1 (313A>G).. Выявлены полиморфные варианты, ассоциированные с развитием заболевания только в одной этнической группе: у русских - CYP1A1 (3798T>C), LTA (252A>G), CAT (-262C>T); у татар - SOD3 (691C>G), IL6 (-174G>C), ADAM33 (12418A>G), CYP1A2 (-163C>A), CYP1A2 (-2467delT), CAT (-262C>T), VDBP (1296T>G и 1307C>A). Показано, что структура наследственной предрасположенности к развитию хронической обструктивной болезни легких отличается у мужчин и женщин.
3. При хронической обструктивной болезни легких тяжесть течения заболевания модифицируют полиморфные варианты генов LTA (252A>G), IL6 (-174G>C), CYP2E1 (-1053C>T), MMP9 (-1562C>T). Развитие эмфиземы легких связано с полиморфными локусами генов UGT2B7 (2146C>T), ADAM33 (13491C>G), ILRN (VNTR), VDBP (1296T>G). Полиморфные варианты генов UGT2B7 (2146C>T), NAT2 (857G>A), GSTM1 (del), CYP1A1 (2454A>G) вносят вклад в вариабельность возраста манифестации заболевания.
4. Определены маркеры развития хронических бронхолегочных заболеваний у детей, включающие локусы, ассоциации с которыми были получены в этнических группах русских и татар - CYP2F1 (c.14_15insC), IL6 (-174G>C), CAT (-262C>T). Выявлены этноспецифические маркеры риска развития заболевания: у русских установлена ассоциация с полиморфными локусами генов GSTP1 (313A>G), MMP3 (-1171 5А>6А), TIMP3 (-1296T>C), SERPINA1 (1237G>A); у татар с полиморфными локусами генов GSTT1 (del), MMP9 (2660A>G), CYP1A1 (2454A>G и 3798T>C), NQO1 (465C>T), MMP12 (-82A>G), VDBP (1296T>G и 1307C>A), IL1RN (VNTR), MMP1 (-1607 G>GG). Показаны особенности ассоциации полиморфных локусов с учетом гендерных различий.
5. Развитие профессионального хронического бронхита ассоциировано с полиморфными локусами генов VDBP (1296T>G и 1307C>A), UGT2B7 (2146C>T), CYP2F1 (c.14_15insC), CYP1A2 (-2467delT), ADAM33 (13491C>G), IL8 (-251T>A), NQO1 (465C>T), MMP1 (-519A>G), TIMP3 (-1296T>C). Выявлено, что у индивидов русской этнической принадлежности риск развития профессионального заболевания ассоциирован с полиморфными вариантами генов GSTP1 (313A>G), CYP2F1 (c.14_15insC), ADAM33 (13491C>G), TIMP3 (-1296T>C); у татар с полиморфными вариантами генов IL8 (-251T>A), VDBP (1296T>G и 1307C>A).
6. Установлено существование отбора против гетерозигот и гомозигот по редким аллелям некоторых полиморфных локусов изученных генов в группе высокостажированных здоровых рабочих: CAT (1167C>T), CYP1A2 (-2467delT), CYP2E1 (-1053C>T), CYP2S1 (13106C>T), ADAM33 (13491C>G), MMP1 (-519A>G), MMP1 (-1607 G>GG), MMP9 (2660A>G), TIMP3 (-1296T>C), VDBP (1296 T>G). Генетическими маркерами устойчивости организма к воздействию производственных факторов являются генотип AG гена GSTP1 (313A>G), генотип n/ins гена CYP2F1 (c.14_15insC), генотип AG гена MMP9 (2660A>G) и гаплотип
С-T гена CAT по локусам -262C>T и 1167C>T.
7. Определены значимые ген-средовые взаимодействия при развитии хронической обструктивной болезни легких и профессионального хронического бронхита с курением для полиморфных локусов генов ADAM33 (12418A>G), TIMP3 (-1296T>C), VDBP (1296T>G), CYP1A1 (2454A>G), CYP2F1 (c.14_15insC), GPX1 (599C>T), SOD3 (691C>G), UGT2B7 (2146C>T), EPHX1 (415A>G); со стажем работы для локусов IL1RN (VNTR), VDBP (1307C>A), CYP1A1 (3798T>C). Спектр маркеров, связанных с формированием хронической обструктивной болезни легких и профессионального бронхита различается в группах подверженных курению и некурящих индивидов.
8. Выявлены наиболее информативные модели межгенных взаимодействий детерминирующих развитие хронической обструктивной болезни легких (CYP2A6 (del), NQO1 (465C>T), ADAM33 (12418A>G), IL1B (3539C>T)), профессионального бронхита (NQO1 (465C>T), UGT2B7 (802C>T), MMP1 (-1607G>GG), ADAM33 (13491C>G)) и бронхолегочных заболеваний у детей (CYP2F1 (c.14_15insC),NQO1 (465C>T), MMP12 (-82A>G), VDBP (1296T>G)). Этническая принадлежность является важным фактором, определяющим специфику взаимодействий между генами при конкретной патологии. Характер взаимодействий ДНК-локусов при формировании хронической обструктивной болезни легких и профессионального бронхита варьирует у курильщиков и некурящих индивидов.
9. В структуру наследственной предрасположенности к бронхолегочным заболеваниям у детей наибольший вклад вносят гены протеолитических ферментов и ингибиторов протеолиза, тогда как при развитии хронической обструктивной болезни легких и профессионального бронхита у взрослых важное значение имеют гены биотрансформации ксенобиотиков и антиоксидантной защиты, эффекты которых реализуются при взаимодействии с факторами окружающей среды. Вклад генетической компоненты снижается при одновременном возрастании роли факторов внешней среды в следующей последовательности: бронхолегочные заболевания у детей - хроническая обструктивная болезнь легких - профессиональный хронический бронхит.
СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ
1. Викторова Т.В., Корытина Г.Ф., Янбаева Д.Г., Макарова О.В., Якупова Э.В., Хуснутдинова Э.К.. Анализ полиморфизма гена N-ацетилтрансферазы 2 в популяциях волго-уральского региона и у больных хроническими обструктивнымим болезнями легких // Генетика. 2003. Т. 39. № 6. С. 855-857.
2. Макарова О.В., Викторова Т.В., Янбаева Д.Г., Корытина Г.Ф., Якупова Э.В., Каримова Л.К.. Полиморфизм генов метаболизма ксенобиотиков у рабочих нефтехимических производств // Генетика. 2003. Т. 39. №. 9. С. 1268-1274.
3. Викторова Т.В., Корытина Г.Ф., Янбаева Д.Г. Взаимодействие генетических и внешнесредовых факторов в процессе развития хронических обструктивных болезней легких. Медицинская генетика. 2003. №2. с. 50-59.
4. Корытина Г.Ф., Янбаева Д.Г., Викторова Т.В. Роль полиморфных вариантов генов цитохромов Р450 (CYP1A1, CYP2E1) и микросомальной эпоксидгидролазы (mEPHX) в патогенезе муковисцидоза и хронических заболеваний дыхательной системы // Молекулярная Биология. 2003. Т. 37. № 5. С. 784-792.
5. Корытина Г.Ф., Янбаева Д. Г., Викторова Т.В. Полиморфизм генов глутатион S-трансфераз у больных муковисцидозом и хроническими заболеваниями дыхательной системы. Генетика. 2004. Т. 40. № 3. С. 401-408.
6. Янбаева.Д.Г, Корытина Г.Ф., Викторова Т.В. Аллельные варианты генов суперсемейства TNF как маркеры тяжести течения хронической обструктивной болезнилегких и бронхоэктактической болезни // Генетика. 2004. Т. 40. №3. С. 409-418.
7. Корытина Г.Ф., Янбаева Д.Г., Бабенкова Л.И., Викторова Т.В. Ассоциация полиморфных вариантов в генах биотрансформации ксенобиотиков с тяжестью легочной патологии у больных муковисцидозом // Медицинская генетика. 2003. Т. 2. № 5. С. 227-232.
8. Янбаева Д.Г., Корытина Г.Ф., Викторова Т.В. Комплексный анализ генов ферментов протеолиза-антипротеолиза (PI, AACT, MMP1) у больных хроническими обструктивными заболеваниями органов дыхания // Молекулярная биология. 2004. Т. 38. № 6. C. 973-979.
9. Корытина Г.Ф., Бабенкова Л. И., Янбаева Д. Г., Эткина Э.И. ,Викторова Т.В. Полиморфизм генов провоспалительных цитокинов (TNFA, LTA, IL1B и IL1RN) у больных муковисцидозом и хроническими бронхолегочными заболеваними // Медицинская генетика. 2004. Т. 3. №7. С. 333-339.
10. Янбаева Д.Г., Байнак О.В., Корытина Г.Ф., Загидуллин Ш.З., Викторова Т.В.. Полиморфные варианты генов провоспалительных цитокинов как маркеры предрасположенности к хронической обструктивной болезни легких // Пульмонология. 2004. № 5. С. 17-22.
11. Korytina G.F., Yanbaeva D.G, Babenkova L.I., Etkina E.I., Victorova T.V. Genetic polymorphisms in the cytochromes P450 (1A1, 2E1), microsomal epoxide hydrolase and glutathione S-transferase M1, T1 and P1 genes and their relationship to severe chronic lung disease in children // Journal of Molecular Medicine. 2005. V. 83. № 9. P. 700-710.
12. Ахмадишина Л.З., Корытина Г.Ф., Мингазова С.Р., Янбаева Д. Г., Бакиров А.Б., Викторова Т.В. Роль полиморфизма генов CYP1A1, EPHX1, GSTM1, GSTT1 и GSTP1 в развитии хронических бронхитов профессионального генеза // Экологическая генетика. 2005. Т. 3. № 1. С. 11-17.
13. Янбаева Д.Г., Корытина Г.Ф., Загидуллин Ш.З., Викторова Т.В. Полиморфизм генов ферментов I и II фазы биотрансформации ксенобиотиков у больных хронической обструктивной болезнью легких // Молекулярная медицина. 2005. № 2. С. 58-64.
14. Ахмадишина Л.З., Янбаева Д.Г., Корытина Г.Ф., Загидуллин Ш.З., Викторова Т.В. Делеционный полиморфизм гена цитохрома P450 (CYP2A6) и его возможная роль в развитии хронической обструктивной болезни легких // Медицинская генетика. 2005. № 12. С. 578-582.
15. Корытина Г.Ф., Ахмадишина Л.З., Янбаева Д.Г., Викторова Т.В. Генотипы гена витамин D связывающего белка (dbp) у больных хронической обструктивной болезнью легких и здоровых жителей Республики Башкортостан // Молекулярная биология. 2006. Т. 40. № 2. С. 231-238.
16. Данилко К.В., Корытина Г. Ф., Ахмадишина Л. З., Янбаева Д. Г., Загидуллин Ш.З., Викторова Т. В. Полиморфизм генов цитокинов (IL1B, IL1RN, TNFA, LTA IL6, IL8, IL10) и их связь с развитием хронической обструктивной болезни легких // Молекулярная биология. 2007. Т. 41. № 1. С. 26-36.
17. Корытина Г.Ф., Ахмадишина Л.З., Янбаева Д. Г., Викторова Т. В. Полиморфизм промоторных областей генов матриксных металлопротеаз (MMP1, MMP9 и MMP12) у больных хронической обструктивной болезни легких и здоровых жителей Республики Башкортостан // Генетика. 2008. Т.44. № 2. С. 242-249.
18. Корытина Г.Ф., Ахмадишина Л.З., Янбаева Д.Г., Байнак О.В., Загидуллин Ш.З., Викторова Т. Ассоциация полиморфных вариантов генов ферментов матриксных металлопротеаз и антипротеаз с развитием и тяжестью течения хронической обструктивной болезни легки // Пульмонология. 2008. № 1 С. 33-39.
19. Кочетова О.В., Корытина Г.Ф., Ахмадишина Л.З., Исхакова Г.М., Викторова Т.В. Анализ полиморфизма гена цитохрома Р450 1А1 (CYP1A1) в этнических группах Республики Башкортостан // Генетика. 2008. Т. 44. № 12 С. 1454-1460.
20. Корытина Г. Ф., Ахмадишина Л. З., Кочетова О.В., Загидуллин Ш.З.,
Викторова Т. В. Ассоциация полиморфных маркеров генов семейства цитохрома P450 1 (CyP1a1 и CYP1a2) с развитием хронической обструктивной болезни легких в Республике Башкортостан // Молекулярная биология. 2008. Т. 42. №. 1. С. 32-41.
21. Ахмадишина Л. З., Корытина Г. Ф., Кочетова О.В., Мингазова С.Р., Бакиров А.Б., Викторова Т. В. Полиморфизм генов цитохромов Р 450 (CYP1A1, CYP1A2, CYP2E1) и риск развитии профессионального хронического бронхита // Медицинская генетика. 2007. Т. 6. № 6 (61) С. 32-38.
22. Корытина Г.Ф., Ахмадишина Л.З., Бабенкова Л.И., Викторова Т.В.Ассоциация полиморфных маркеров генов семейства цитохрома P450 и ферментов антиоксидантной защиты с развитием хронических заболеваний респираторной системы у детей // Медицинская генетика. 2007. Т. 6. № 8. С.42-48.
23. Корытина Г.Ф., Ахмадишина Л.З., Целоусова О.В., Загидуллин Ш.З., Викторова Т.В. Анализ полиморфных вариантов генов ферментов антиоксидантной защиты и их связь с развитием хронической обструктивной болезни легких у жителей Республики Башкортостан // Генетика. 2009. Т. 45. № 5. С. 850-858.
24. Корытина Г.Ф., Целоусова О.С., Ахмадишина Л.З., Бабенкова Л.И., Викторова Т.В. Анализ ассоциации полиморфных локусов генов медиаторов воспаления (IL1В, TNFA, LTA, IL8, IL6, IL1RN, IL10, TGFb, TLR4, DBP) с развитием хронических заболеваний респираторной системы у детей // Медицинская генетика. 2008. T. 7. № 2 (68). С.17-26.
25. Ахмадишина Л. З., Корытина Г. Ф., Кочетова О.В., Мингазова С.Р., Бакиров А.Б., Викторова Т. В. Профессиональный хронический бронхит: роль полиморфных вариантов генов ферментов-антиоксидантов в предрасположенности к заболеванию // Пульмонология. 2008. № 2. С. 68-73.
26. Кочетова О.В., Корытина Г.Ф., Ахмадишина Л.З., Каримова Л.К., Викторова Т.В. Применение новых технологий ДНК диагностики в медицине труда // Гигиена и санитария. 2008. № 5. С.43-47.
27. Кочетова О.В., Корытина Г.Ф., Ахмадишина Л.З., Каримова Л.К., Викторова Т.В. Роль полиморфных маркеров генов антиоксидантов (NQO1, GSTP1, GPX1) в формировании адаптивного ответа работающих к вредным условиям труда // Медицина труда и промышленная экология. 2010. №3. С. 17-24.
28. Целоусова О.С., Корытина Г.Ф., Ахмадишина Л.З., Бабенкова Л.И., Викторова Т.В. Роль генов матриксных металлопротеаз (MMP1, MMP9 и MMP12) в формировании предраcположенности к хроническим заболеваниям респираторной системы у детей // Медицинская генетика. 2009. №2. С.11-17.
29. Целоусова О.С., Ахмадишина Л. З., Корытина Г.Ф., Викторова Т.В. Ассоциативное исследование полиморфных вариантов генов систем протеолиза-антипротеолиза и медиаторов воспаления с развитием профессионального хронического бронхита // Медицинская генетика. 2009. Т. 8. № 11. (89). С. 28-37.
30. Целоусова О.С., Кочетова О.В., Ахмадишина Л.З., Корытина Г.Ф., Викторова Т.В. Полиморфные варианты генов цитохромов Р450 (CYP1A1, CYP2E1, CYP2D6) в развитии предрасположенности к профессиональному поражению печени // Бюллетень ВСН - СО РАМН. 2009. №1. (65) С.136-140.
31. Ахмадишина Л.З., Целоусова О.С., Корытина Г.Ф., Викторова Т.В. Роль полиморфных вариантов генов глутатион-S-трансфераз в развитии предрасположенности к профессиональному хроническому бронхиту // Бюллетень ВСН - СО РАМН. 2009. № 1. (65) С. 88- 92.
32. Корытина Г.Ф., Ахмадишина Л.З., Викторова Т.В. Частоты полиморфных вариантов генов CYP1B1 и CYP2F1 в трёх этнических группах жителей Республики Башкортостан и у больных хронической обструктивной болезнью легких // Молекулярная биология. 2010. Т. 44, №. 1. С. 33Ц41.
33. Ахмадишина Л. З., Корытина Г. Ф., Викторова Т. В. Полиморфные маркеры генов CYP1B1 (4326C>G), CYP2F1 (c.14_15insC), CYP2J2 (-76G>T), CYP2S1 (13106C>T, 13255A>G) и генетическая предрасположенность к хроническим заболеваниям органов дыхания, вызванных курением и профессиональной деятельностью // Генетика. 2011. Т. 47. № 10. С. 1402- 1410.
СПИСОК СОКРАЩЕНИЙ
ПДРФ ПЦР ХБЛЗ ХОБЛ ХПБ ADAM33 AhR CAT CYP EPHX1 GOLD GPX1 GST IL IL1RN LTA MMP NQO1 NAT2 OR PY SERPINA SOD TGFB1 TIMP TLR4 TNFA UGT2B7 VDBP | Полиморфизм длины рестрикционных фрагментов Полимеразная цепная реакция Хронические бронхолегочные заболевания Хроническая обструктивная болезнь легких Хронический профессиональный бронхит A disintegrin and metalloproteinase domain 33 (дизинтегриновая металлопротеаза 33) Aryl hydrocarbon receptor (арилгидрокарбоновый рецептор) Catalase (каталаза) Сytochrome P450 (цитохром P450) epoxide hydrolase 1, microsomal (эпоксидгидролаза 1 типа, микросомальная) Global Initiative for Chronic Obstructive Pulmonary Disease Glutathione peroxidase 1 (глутатионпероксидаза 1 типа, Se-зависимая) Glutathione S-transferase (глутатион-S-трансферазы) Interleukin (интерлейкин) Interleukin 1 receptor antagonist (рецепторный антагонист интерлейкина 1) Lymphotoxin alpha (TNF superfamily, member 1) имфотоксин- Multifactor-Dimensionality Reduction Matrix metallopeptidase (матриксные металлопротеиназы) NAD(P)H dehydrogenase, quinone 1 (НАД(Ф)Н-хинон оксидоредуктаза 1) N-acetyltransferase 2 (arylamine N-acetyltransferase) ариламин N ацетилтрансфераза 2 Odds Ratio (отношение шансов) Pack-years (индекс курения, выраженный в пачках/лет) Serpin peptidase inhibitor, clade A (сериновые ингибиторы протеаз кластера А) Superoxide dismutase (супероксиддисмутаза) Transforming growth factor, beta 1 (трансформирующий фактор роста, 1-субъединица) Tissue inhibitor of metalloproteinases (тканевые ингибиторы металлопротеиназ) Toll-like receptor 4 (толл - подобный рецептор 4) Tumor necrosis factor (фактор некроза опухоли, -субъединица) UDP glucuronosyltransferase 2 family, polypeptide B7 (УДФ-глюкуронозилтрансфераза) Group-specific component, vitamin D binding protein (витамин Д-связывающий белок) |