Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

ASSOCIATION OF TYPE 1 DIABETES WITH POLYMORPHIC ALLIES OF HLA CLASS II GENES IN POPULATION OF YAKUTS

Pavlova N.I. 1 Kurtanov K.A. 1 Dyakonova A.T. 1 Soloveva N.A. 1 Sydykova L.A. 2, 3 Aleksandrova T.N. 1 Filippova N.P. 1 Nikiforova M.E. 4 Dodokhov V.V. 1, 5
1 Yakut science center of complex medical problems
2 Federal State Autonomous Educational Institution of Higher Education "M. K. Ammosov North-Eastern Federal University"
3 State budgetary institution of the Republic of Sakha (Yakutia) "Yakut Republican Clinical Hospital"
4 State Autonomous Institution of the Republic of Sakha (Yakutia) Republican Hospital No. 1 - National Center of Medicine Department of Pediatric Endocrinology and Gastroenterology.
5 Federal State Budgetary Educational Institution of Higher Education Yakutsk State Agricultural Academy
For the first time, a molecular genetic analysis of the frequency distribution of alleles and genotypes of polymorphisms rs3104413, rs2854275 and rs9273363 of the HLA-DQA1 and HLA-DQB1 genes, respectively, was carried out in the Yakut population. 92 patients diagnosed with type 1 diabetes were tested, 44 of which belonged to men and 48 women. The average age of patients was 23.04 ± 0.27 years (from 4 to 56 years). A sample of 210 healthy volunteers was formed as a comparison group. Amplification of the HLA-DQA1 gene region containing the single nucleotide polymorphism rs3104413 was carried out in the process of real-time PCR, the HLA-DQB1 gene containing the rs2854275 and rs9273363 polymorphisms - allele specific PCR. Analysis of the distribution of genotypes showed that the most common genotype for the rs3104413 polymorphism in the group of patients with type 1 diabetes is C / C (69.6%), whereas in the control group the homozygous G / G variant was more common (41.1%). It was revealed that the most common genotype of the polymorphic variant rs2854275 in the control group is the homozygous genotype G / G (88.5%), and in the group of patients with type 1 diabetes the genotype G / T (56.3%). The polymorphic variant rs9273363 was characterized by a significantly low frequency of allele A (1.5%) in the group of healthy individuals compared with the group of patients with type 1 diabetes (53.6%). The highest relative risk indicators (RR) among the genotypes were found for the combination of the C / C genotype of the rs3104413 polymorphism with the heterozygous G / T genotype of the rs2854275 - RR polymorphism of 4.234 (p <0.001); for a combination of the C / C genotype of the rs3104413 polymorphism with the T / T genotype of the rs2854275 polymorphism - RR = 3.692 (p <0.001); for a combination of the G / G polymorphism rs3104413 genotype and the A / C genotype, the rs9273363 polymorphism corresponded to RR = 3.625 (p <0.001). The obtained data can be used as biological predictors of the development of type 1 diabetes in order to conduct timely personalized preventive measures.
hla class ii genes
type 1 diabetes mellitus
polymorphism
yakuts

Сахарный диабет 1-го типа (СД 1) – метаболическое (обменное) заболевание, характеризующееся гипергликемией, в основе которого лежит деструкция β-клеток, приводящая к абсолютному дефициту инсулина.

Сахарный диабет представляет собой важнейшую медико-социальную проблему во всем мире. Это объясняется его широким распространением, тяжестью поздних осложнений, дороговизной средств диагностики и лечения, которые необходимы больным в течение всей жизни. По данным Росстата на сегодняшний день в мире насчитывается более 400 млн человек, больных СД, и распространенность заболевания продолжает неуклонно расти. На начало 2017 г. общая численность пациентов, страдающих СД, в РФ составила 4,348 млн человек, из них: 4 млн – СД 2-го типа, 255 тыс. – СД 1-го типа и 75 тыс. – СД других типов. По последним данным в Республике Саха (Якутия) с диагнозом СД проживают 21 677 человек, из них: 20 508 – это пациенты с СД 2-го типа, 1099 – больные СД 1-го типа, 70 человек страдают СД других типов [1].

Согласно современным данным в развитии СД 1-го типа принимают участие большое число генов [2], более половины генетических рисков обусловлено участием полиморфных вариантов генов HLA, расположенных на коротком плече хромосомы 6 (6p21). Полагая, что вклад генетических факторов в риск развития СД 1-го типа составляет более 50%, полный вклад полиморфных аллелей локуса HLA может быть оценен более чем 25%.

Определение аллелей HLA помогает в определении риска развития СД 1-го типа. Это особенно полезно в профилактике у пациентов с высокой степенью риска [3]. HLA-типирование также требуется в генетических исследованиях, направленных на определение молекулярной основы СД 1-го типа [4]. Тем не менее высокая стоимость HLA генотипирования является одной из основных проблем, препятствующих проведению масштабных исследований. В связи с этим ведутся исследования по определению HLA аллелей с использованием однонуклеотидных полиморфизмов (SNP) (Walsh E.C. et al., 2003; Leslie S. et al., 2008; Dilthey A.T. et al., 2011; Ferreira R.C. et al., 2012; Barker J.M. etal. 2008; Cao Ngueyen et al., 2013).

В России исследования по определению HLA-аллелей с использованием однонуклеотидных полиморфизмов не проводятся. В связи с этим на сегодняшний день существует необходимость в проведении исследований, направленных на разработку регионально-адаптированного метода HLA-типирования СД 1-го типа с помощью однонуклеотидных полиморфизмов (SNP).

Целью данного исследования являлось выявление ассоциации полиморфных вариантов генов HLA II класса с СД 1-го типа среди якутской популяции.

Материалы и методы исследования

Экспериментальная часть работ по генотипированию полиморфизмов rs3104413, rs2854275, rs9273363 была проведена в лаборатории наследственной патологии отдела молекулярной генетики Якутского научного центра комплексных медицинских проблем (ЯНЦ КМП). Для исследования использованы образцы ДНК из коллекции биоматериала ЯНЦ КМП с применением УНУ «Геном Якутии» (рег. № USU_507512). Выборка состояла из 92 пациентов больницы ФГБНУ «Якутский научный центр комплексных медицинских проблем», ГАУ РС(Я) РБ № 1 НЦМ Педиатрический центр и эндокринологического отделения ГБУ РС (Я) «Якутская городская клиническая больница» г. Якутска. В состав выборки вошли 92 пациента с диагнозом СД 1-го типа в возрасте от 4 лет до 56 лет, проживающих в РС (Я), якутов по этнической принадлежности. По половому признаку пациенты разделились на 44 мужчин (47,8%) и 48 женщин (52,2%). Средний возраст пациентов составил 23,04±0,27 года (от 4 до 56 лет), средний возраст пациентов мужского пола – 20,5±2,3 года (от 5 до 40 лет), женского пола – 25,11±2,59 года (от 4 до 56 лет). Популяционную выборку составили 210 якутов, не страдающих СД 1-го типа. Этническая принадлежность учитывалась до третьего поколения.

Амплификация области гена HLA-DQA1, содержащего однонуклеотидный полиморфизм rs3104413, осуществлялась в процессе ПЦР в реальном времени с использованием пар праймеров и аллельспецифичных зондов для амплификации ДНК, описанных в журнале «NATURE COMMUNICATIONS» (Isabelle Serr et al., 2015). Праймеры и зонды синтезированы компанией ООО «Биотех-Индустрия (Lumiprobe)» (Москва, Россия). Последовательность праймеров: Форвард праймер 5'-CAGCTGAGCACTGAGTAG-3', реверс праймер 5'-GCAGTTGAGAAGTGAGAG-3'. Структура зондов: FAM – Probe rs3104413 LPC [6FAM]CAGCCT[+G]CT[+C]TC[+C]TA[+T]TGG[BHQ1], HEX – Probe rs3104413 LPG [HEX]CAGCCT[+G]CT[+G]TC[+C]TA[+T]TGG[BHQ1].

Амплификация проводилась согласно температурной программе, приведенной ниже (табл. 1).

Таблица 1

Температурная программа амплификации полиморфизма rs3104413

Стадии

Температура

Время

Циклы

1

Первая денатурация

95°С

0:10

1 цикл

2

Денатурация

95°С

0:30

50 циклов

3

Отжиг

55°С

1:00

 

Измерение сигнала флуоресценции проводилось на втором этапе реакции (55°С – 1:00). Детекция флуоресценции осуществлялась «по конечной точке» согласно протоколу прибора «Real-time CFX 96 Touch» («Biorad», США). Пример распределения облаков генотипов проведения ПЦР и детекции флуоресценции «по конечной точке» представлен на рисунке 1. Соответствие флуоресцентных красителей: аллель C – канал FAM, аллель G – канал HEX.

Рис. 1. Распределение облаков генотипов полиморфизма rs3104413 гена HLA-DQA1. Примечание: «К–» – отрицательный контроль, C/C – гомозигота по предковому аллелю C. C/G – гетерозигота, G/G – гомозигота по мутантному аллелю G

Анализ однонуклеотидных полиморфизмов (SNP) rs2854275, rs9273363 проводился с помощью аллельспецифичных праймеров (АС-ПЦР). Праймеры были синтезированы компанией ООО «Сибэнзим» (г. Новосибирск, Россия). В рамках настоящей работы были созданы оригинальные аллельспецифичные праймеры на основе референсного генома hg38/Human по базе UCSC Genome Browser (http://genome.ucsc.edu/index.html). Последовательность праймеров и условия ПЦР представлены в таблице 2.

Таблица 2

Последовательность праймеров и условия АС-ПЦР*

Ген/

Антиген

Метод генотипирования

Структура праймеров

tоС

отжига праймеров

,

Мутант- ный аллель

rs2854275

АС-ПЦР*

Внутренний форвард праймер:

5’-CTTAACTTTGGTGGCATCTTCTTGTG-3’

Внутренний реверс праймер:

5’- AAGCTGTGGTTCTGGCTCCACATTTA-3’

Наружный форвард праймер:

5’-TAAGAGGGAAGAGCATGAGCTGAGT-3’

Наружный реверс праймер:

5’-TCTTTCAGTCACTGGAAAATGCTTACA-3’

64

134, 192 п.о.**

rs9273363

АС-ПЦР*

Внутренний форвард праймер:

5’-TGCTTTCAGGGTCATGGCCTTC-3’

Внутренний реверс праймер:

5’-TGGGAGCTTCTGCGAGGTTATGT -3’

Наружный форвард праймер:

5’-TAGTTCTCCACTCTGACCTCATC-3’

Наружный реверс праймер:

5’-GTTATGAGGGAGGTCACTACACA -3’

67

146, 141 п.о.**

β-globin***

АС-ПЦР* (Вутрен-ний контроль)

F: 5’-CAACTTCATCCACGTTCACC-3’

R: 5’-GAAGAGCCAAGGACAGGTAC-3’

57

268 п.о.**

 

Примечание. * – АС-ПЦР – аллель-специфическая ПЦР, ** – п.о. – пара оснований. *** – Dominquez O et al., 1992; Michelle A.M. et al., 2000; Steffens Nakken et al., 1995

Результаты амплификации и рестрикции фракционировали в 3–4%-ном агарозном геле с бромистым этидием при напряжении 120–300 Вт в течение 30–45 мин. Анализ результатов электрофореза проводили под UV-лучами в гель-документирующем приборе Vilbеr Lourmat (рис. 2, 3).

Рис. 2. Электрофореграмма продукта амплификации rs2854275 в 2%-ном агарозном геле Примечание: п.н. – пар нуклеотидов. М – ДНК маркер pUC19/MspI. 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 – нумерация образцов. 270 п.н. – внутренний контроль. 192 п.н. – аллель T, 134 п.н. – аллель G

Рис. 3. Электрофореграмма продукта амплификации rs9273363 в 2%-ном агарозном геле Примечание: п.н. – пар нуклеотидов. М – ДНК маркер pUC19/MspI. 1, 2, 3, 4, 5 – нумерация образцов. 270 п.н. – внутренний контроль. 141 п.н. – аллель A, 146 п.н. – аллель C

Статистический анализ полученных результатов медико-генетического исследования был проведен с помощью программы: «Office Microsoft Excel 2010», «Statistica 8.0». Частоты аллелей и генотипов rs3104413, rs2854275 и rs9273363 определяли путем прямого подсчета. Результаты считались значимыми, когда значение «р» было меньше 0,05 (p<0,05).

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

В результате проведенного анализа установлено, что только для полиморфного варианта rs2854275 в популяции якутов наблюдаемое распределение генотипов соответствовало ожидаемому при равновесии Харди–Вайнберга (РХВ). Отклонение установлено для полиморфного варианта rs3104413 (Х²=65,966, р<0,001), которое характеризовалось уменьшением уровня наблюдаемой гетерозиготности (Но=0,219) по сравнению с ее ожидаемым уровнем (Не=0,498). Для полиморфного варианта rs9273363 (Х² = 134, р<0,001), которое характеризовалось уменьшением уровня наблюдаемой гетерозиготности (Но=0,0) по сравнению с ее ожидаемым уровнем (Не=0,029) (табл. 3).

Таблица 3

Распределение генотипов и аллелей полиморфизмов rs3104413, rs2854275, rs9273363 в популяционной выборке якутов

Полимор-физм

Кол-во человек, n

Генотип, %

Распределение

аллелей

Но

Не

Х2

р

rs3104413

210

С/С

C/G

G/G

С

G

0,219

0,498

65,966

<0,001

36,2

21,9

41,9

0,471

0,529

rs2854275

122

Т/Т

Т/G

G/G

Т

G

0,157

0,145

1,297

0,255

0,0

15,7

84,3

0,079

0,921

rs9273363

134

А/А

А/С

С/С

А

С

0,0

0,029

134

<0,001

1,5

0,0

98,5

0,015

0,985

Примечание: n – численность выборки, Ho, He – наблюдаемая и ожидаемая гетерозиготность, Х2 по Харди–Вайнбергу, р – уровень достоверности

Согласно базе данных проекта «1000 геномов» для полиморфизма rs3104413 во всех исследованных выборках популяций мира установлены высокая частота мутантного аллеля С (78–92%) и низкая частота предкового аллеля G (8–22%) [5]. В исследованной выборке якутов частота предкового аллеля G была выше и составила 52,9%, тогда как частота мутантного аллеля С характеризовалась более низкими значениями и составила 47,1%, что может быть следствием как недостаточного объема выборки, так и этнической особенности.

Для полиморфизма rs2854275 в исследованной выборке якутов частота предкового аллеля G (92,1%) преобладала над частотой мутантного аллеля Т (7,9%), что сопоставимо с результатами ранее проведенных исследований, в которых отмечаются высокая частота предкового аллеля G (90–97%) и низкая частота мутантного аллеля Т (3–10%) [5].

Полиморфизм rs9273363 у якутов характеризовался преобладанием частоты предкового аллеля С (98,5%) по сравнению с частотой мутантного аллеля А (1,5%), что согласуется с результатами ранее проведенных исследований, в которых отмечаются высокая частота предкового аллеля С (67–91%) и низкая частота мутантного аллеля А (9–33%) [5].

Результаты анализа распределения частот аллелей и генотипов полиморфизмов rs3104413, rs2854275 и rs9273363 среди больных СД 1-го типа и контрольной выборки представлены в таблице 4.

Таблица 4

Частота встречаемости генотипов и аллелей полиморфизмов rs3104413, rs2854275 и rs9273363 в группе больных СД 1-го типа и контрольной выборки

Генотипы, аллели

Больные СД 1-го типа

Контрольная выборка

Х2

OR

(95% CI)

для аллелей

Значимость, p

rs3104413

 

n = 92, абс. (%)

n = 210, абс.(%)

     

С/С

64 (69,6)

76 (36,2)

32,099

4,036

(2,708-6,017)

<0,001*

С/G

16 (17,4)

46 (21,9)

G/G

12 (13,0)

88 (41,9)

С

144 (0,783)

198 (0,471)

49,184

<0,001**

G

40 (0,217)

222 (0,529)

rs2854275

 

n = 92, абс. (%)

n = 122, абс. (%)

     

Т/Т

16 (16,7)

0 (0,0)

88,7

13,8

(7,48-25,46)

<0,001*

G/T

52 (56,3)

14 (11,5)

G/G

24 (27,1)

108 (88,5)

T

84 (44,8)

14 (5,7)

92,41

<0,001**

G

100 (55,2)

230 (94,3)

rs9273363

 

n = 20, абс.(%)

n = 134, абс.(%)

     

А/А

4 (14,3)

2 (1,5)

138,89

76,15

(24,89-233,02)

<0,001*

А/С

22 (78,6)

0 (0,0)

С/С

2 (7,1)

132 (98,5)

А

30 (53,6)

4 (1,5)

128,28

<0,001**

С

26 (46,4)

264 (98,5)

Примечание. Достигнутый уровень значимости при сравнении распределения генотипов (*) и частоты аллелей (**) в группах сравнения 1 и 2. n – численность выборок, Х2 с поправкой Йейтса.

Анализ распределения генотипов показал, что наиболее распространенным генотипом по полиморфизму rs3104413 в группе больных СД 1-го типа является С/С (69,6%), тогда как в контрольной группе чаще встречался гомозиготный вариант G/G (41,1%). В ходе исследования установлено, что преобладающим аллелем полиморфизма rs3104413 в популяции якутов является аллель G (52,9%), при том что в ранее исследованных выборках популяций мира установлено преобладание частоты мутантного аллеля С (78–92%) [5]. Таким образом, мутантный аллель С для якутов может использоваться в качестве маркера повышенного риска развития СД 1 (OR – 4,036; 95% CI: 2,71–6,02; р<0,001).

Выявлено, что самым распространенным генотипом полиморфного варианта rs2854275 в группе контроля является гомозиготный генотип G/G (88,5%), а в группе больных СД 1-го типа – генотип G/Т (56,3%). Анализ отношения шансов показал, что частота встречаемости аллеля Т в группе людей с СД 1-го типа статистически достоверно выше (44,8%), чем среди людей, не страдающих сахарным диабетом 1-го типа (5,7%), что позволяет считать носительство аллеля Т предрасполагающим фактором к развитию СД 1-го типа (OR= 13,8; 95% CI: 7,48–25,46; р<0,001).

Полиморфный вариант rs9273363 характеризовался достоверно низкой частотой аллеля А (1,5%) в группе здоровых индивидов по сравнению с группой больных СД 1-го типа (53,6%). Отношение шансов (OR) риска развития заболевания для носителей аллеля А составило 76,15 (95% CI:24,89–233,02; р<0,001).

Результаты анализа распределения частот генотипов rs3104413 в сочетании с генотипами полиморфизмов rs2854275, rs9273363 среди больных СД 1-го типа и контрольной выборки представлены в таблице 5.

Таблица 5

Частота встречаемости генотипов полиморфизмов rs3104413, rs2854275, rs9273363 среди больных СД 1-го типа и контрольной выборки

Вариант сочетания генотипов, №

rs3104413

rs2854275

rs9273363

Частота в группе СД 1-го типа, абс. (%)

Частота в контрольной группе, абс. (%)

RR

(95% CI)

р

1

С/С

T/T

14 (15,2)

0

3,692

(3,055 – 4,463)

<0,001

2

С/С

G/G

8 (8,7)

68 (32,4)

0,283

(0,144 – 0,557)

<0,001

3

С/С

G/T

42 (45,7)

8 (3,8)

4,234

(3,212 – 5,58)

<0,001

4

G/G

А/С

12 (13)

0

3,625

(3,008 – 4,368)

<0,001

5

G/G

C/C

0

88 (41,9)

0,000

<0,001

6

С/G

G/G

А/А

2 (2,2)

0

3,333

(2,804 – 3,962)

0.170

7

С/G

G/G

А/С

6 (6,5)

0

3,442

(2,881 – 4,113)

0.002

8

С/G

G/G

C/C

0

40 (19)

0,000

<0,001

9

С/G

G/T

C/C

2 (2,2)

4 (1,9)

1,096

(0,349 – 3,444)

0.769

10

С/G

G/T

А/С

4 (4,3)

0

3,386

(2,842 – 4,036)

0.013

11

С/G

G/T

А/А

2 (2,2)

2 (1)

1,656

(0,612 – 4,478)

0.759

Всего выявлено 11 вариантов сочетаний генотипов полиморфизмов rs3104413, rs2854275, rs9273363. В группе больных с СД 1-го типа наиболее часто встречались варианты № 3, 1 и 4, тогда как в группе контроля чаще отмечались № 5, 2 и 8.

Заключение

Наиболее высокими показателями относительного риска среди сочетаний генотипов являлись С/С полиморфизма rs3104413 в сочетании с гетерозиготным генотипом G/T полиморфизма rs2854275 – RR = 4,234 (р<0,001); генотип С/С полиморфизма rs3104413 в сочетании с генотипом Т/Т полиморфизма rs2854275 – RR = 3,692 (р<0,001); генотип G/G полиморфизма rs3104413 в сочетании с генотипом А/С полиморфизма rs9273363 – RR = 3,625 (р<0,001).

Таким образом, полученные данные могут быть использованы в качестве биологических предикторов развития СД 1-го типа с целью проведения своевременных персонализированных профилактических мероприятий.

Исследование было проведено в рамках НИР Изучение генетической структуры и груза наследственной патологии популяций Республики Саха (Якутия).