Polymorphism of LYPLAL1 and TGFA genes associated with progression of knee osteoarthritis in residents Central Chernozem Region of Russia



Cite item

Abstract

Background. Кnee Osteoarthritis (OA) is a multifactorial disease in which genetic factors play an important role. The share of the hereditary component in the development of OA, according to various literature sources, ranges from 40 to 65%. Кnee OA is a progressive disease that leads to a decrease in the quality of life and disability. The study aimed to evaluate the role of polymorphic markers of candidate genes rs2820436 and rs2820443 LYPLAL1, rs3771501 TGFA, rs11177 GNL3, rs6976 GLT8D1 in the progression of knee OA in the population of the Central Chernozem Region of Russia. Materials and methods. The study was performed in a case-control design on a sample of 500 patients with knee OA. Case – patients with III-IV stages of the disease according to J. Kellgren - J. Lawrence (n=325), control (individuals who do not have the analyzed sign - III-IV stages of the disease) - patients with stage II (n=175). Genotyping of five single nucleotide polymorphisms (SNPs) of candidate genes was performed using the polymerase chain reaction method for DNA synthesis. The study of the associations of the studied polymorphic loci, the calculation of haplotype frequencies and the analysis of their relationship with the progression of knee OA was carried out by the method of logistic regression in the program PLINK v 2.050. Results. Significant associations with the progression of OA of the knee were established for allelic variant A rs2820436 of LYPLAL1 gene according to allelic (OR=1.48, p=0.010, pperm=0.012), additive (OR=1.58, p=0.009, pperm=0.010), dominant (OR=1.61, p=0.024, pperm=0.030) genetic models and A/A genotype of the same polymorphism (OR=2.53, p=0.041). The genotypes C/C rs2820436 LYPLAL1 (OR=0.67, p=0.043), A/G rs3771501 TGFA (OR=0.67, p=0.042) have a protective role in the progression of the disease. It was found that the frequency of the AC haplotype of haploblock rs2820436-rs2820443 in the group of patients with III-IV stages of the disease was significantly higher than in patients with stage II (OR=1.83, p=0.002, pperm=0.002). The identified molecular genetic markers rs2820436 and rs2820443 of LYPLAL1 gene, rs3771501 of TGFA gene are associated both with the risk of developing OA according to previous genome-wide studies and, according to our data, are associated with the progression of knee OA. Conclusions. Genetic risk factors for the development of knee OA of III-IV radiological stages are allelic variant A and genotype A/A rs2820436 of LYPLAL1 gene, haplotype AC of haploblock rs2820436-rs2820443 in the population of the Central Chernozem Region of Russia. Genotypes C/C rs2820436 of LYPLAL1 gene and A/G rs3771501 of TGFA gene have a protective value in the progression of this disease.

Full Text

Полиморфизм генов LYPLAL1 и TGFA ассоциирован с прогрессированием остеоартроза коленного сустава у жителей Центрального Черноземья России

Введение

Под остеоартрозом (ОА) понимают хроническое заболевание опорно-двигательного аппарата, которое характеризуется разрушением суставного хряща, воспалением синовиальной оболочки, ремоделированием кости [1]. Остеоартрозом страдает более 240 миллионов человек во всем мире [2]. Распространённость ОА коленного сустава колеблется в достаточно широких диапазонах, зависит от популяции, пола, возраста [3, 4]. В работе Л.И. Алексеевой и др. (2019) приведены обобщённые данные нескольких эпидемиологических исследований гонартроза, которые свидетельствуют о том, что частота встречаемости данного недуга может варьировать в диапазоне от 2,0 до 42,4%, при условии, что для диагностики применялись исключительно клинические критерии, от 16,3 до 33,0% – только рентгенологические критерии, и от 1,5 до 15,9% – комбинация клинических и рентгенологических критериев [3]. В России ОА коленного и/или тазобедренного суставов страдает 13% населения старше 18 лет [5]. ОА коленного сустава приводит к снижению качества жизни, потере трудоспособности, ограничению переносимости физической нагрузки, а прогрессирование заболевания может стать причиной инвалидности [6].

ОА коленного сустава является хроническим, медленно прогрессирующим, с нарастанием симптомов заболеванием [7]. Данные по факторам риска прогрессирования ОА коленного сустава остаются весьма противоречивыми [8]. Согласно литературным данным, в развитии/прогрессировании ОА значимая роль отводится наследственным (генетическим) факторам [1, 3, 9-11]. Однако проведённые до настоящего времени молекулярно-генетические исследования ОА коленного сустава позволили установить сравнительно небольшое количество полиморфных вариантов генов-кандидатов (ADAM12, CILP, OPG, TNA, IL-1, GDF5), ассоциированных с прогрессированием данной патологии [12-16]. Стоит отметить малочисленность проведённых исследований по данной проблеме, а в отечественной литературе вообще отсутствуют исследования генетических факторов риска прогрессирования гонартроза. Выше изложенные факты и диктуют необходимость проведения дальнейших исследований, направленных на поиск полиморфных генетических маркеров прогрессирования ОА коленного сустава в различных популяциях мира, в том числе и России.

Цель исследования – оценить роль полиморфных маркеров генов-кандидатов rs2820436 и rs2820443 LYPLAL1, rs3771501 TGFA, rs11177 GNL3, rs6976 GLT8D1 в прогрессировании ОА коленного сустава у населения Центрального Черноземья РФ.

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

Дизайн исследования

Исследование выполнено в дизайне «случай-контроль» на выборке из 500 больных ОА коленного сустава. «Группа случаев» – больные с III-IV рентгенологическими стадиями заболевания по J. Kellgren – J. Lawrence, «группа контроля» (индивидуумы, не имеющие анализируемого признака – III-IV стадий заболевания) – пациенты со II стадией заболевания. Выборка для исследования была сформирована в период с февраля 2016 по декабрь 2018 гг.

Критерии включения: 1) русская национальность, территория рождения и проживания – Центральное Черноземье РФ, 2) возраст от 40 лет, 3) диагностированный первичный ОА коленного сустава II-IV стадий по J. Kellgren – J. Lawrence (1957), 4) наличие болевого синдрома по визуально-аналоговой шкале (ВАШ) более 40 баллов, 5) наличие добровольного информированного согласия на исследование.

Критерии исключения: 1) не русская национальность, проживание и/или рождение вне Центрально-Чернозёмного региона России, 2) наличие тяжёлых форм артериальной гипертензии, ишемической болезни сердца, сахарного диабета, почечно-печёночной недостаточности, онкологических заболеваний, системных заболеваний соединительной ткани, травм суставов в анамнезе, воспалительные заболевания суставов, врождённых пороков развития опорно-двигательной системы, 3) отказ от участия в исследовании.

У всех больных, включённых в исследование, был диагностирован первичный ОА коленного сустава на основе клинического осмотра, лабораторного и рентгенологического исследований согласно Федеральным клиническим рекомендациям по диагностике и лечению остеоартроза. У 175 больных (39,0%) была установлена II рентгенологическая стадия ОА коленного сустава по шкале J. Kellgren – J. Lawrence, 230 пациентов (46,0%) имели III рентгенологическую стадию, и у 95 больных (19,0%) установлена IV стадия заболевания.

Для исследования ассоциаций полиморфных маркеров генов-кандидатов с прогрессированием ОА коленного сустава изучаемая выборка больных была разделена на две группы. Первая группа пациентов (n=325) представляла собой объединённую выборку больных с III и IV рентгенологическими стадиями, во вторую группу (n=175) были включены больные со II рентгенологической стадией. Объединение групп пациентов с III-IV стадиями заболевания было выполнено с целью повышения мощности проводимого анализа ассоциаций.

В качестве объекта генетического исследования использовалась геномная ДНК, выделенная из периферической венозной крови методом фенол-хлороформной экстракции. Полиморфные маркеры генов-кандидатов были отобраны для исследования на основании следующих критериев: 1) наличие ассоциаций с ОА по данным ранее проведенных полногеномных исследований (GWAS) [17]; 2) значимый регуляторный потенциал согласно данным онлайн-ресурса HaploReg (v.4.1) [18]; 3) частота полиморфизма не менее 5%.

В соответствии с данными критериями в исследование были включены пять SNPs четырёх генов-кандидатов: rs2820436 и rs2820443 LYPLAL1, rs3771501 TGFA, rs11177 GNL3, rs6976 GLT8D1. Все полиморфные локусы ассоциированы с ОА по данным GWAS [19-23], имеют выраженный регуляторный потенциал (rs2820443, rs3771501, rs11177 расположены в эволюционно-консервативном регионе ДНК; rs2820443, rs11177, rs6976 находятся в регионе связывания с регуляторными белками и в участках повышенной чувствительности к ДНКазе; rs2820443, rs11177 локализуются в области гистоновых белков, маркирующих промоторы; все 5 SNPs расположены в регионе гистонов, маркирующих энхансеры), частоты минорных аллелей составили более 5%.

Образцы ДНК были генотипированы на амплификаторе Real-Time CFX96 (США) методом TaqMan зондов с использованием соответствующих праймеров и зондов.

Статистический анализ

Для количественных показателей (возраст, рост, ИМТ) рассчитывали среднее арифметическое и стандартное отклонение (M±m). Для качественных показателей в изученных группах определяли их абсолютные значения и проценты. Для сравнения количественных показателей применяли критерий Манна – Уитни, для качественных признаков – χ2. Статистически значимым считали результат при p<0,05.

Проводили оценку наблюдаемого распределения генотипов для пяти SNPs, включенным в исследование, и его соответствие ожидаемому распределению, согласно равновесия Харди–Вайнберга (отклонения принимались за статистически значимые при p≤0,05). При сравнительном анализе частот аллелей и генотипов между двумя группами больных ОА коленного сустава применяли критерий χ2 с поправкой Йетса на непрерывность, использовали таблицы сопряженности 2х2. Все вычисления выполнялись в программе STATISTICA for Windows 10.0. За статистически значимый уровень принимали р<0,05. Изучение ассоциаций исследуемых полиморфных локусов с прогрессированием ОА коленного сустава проводили с использованием метода логистической регрессии (рассматривали четыре генетические модели: аллельная, аддитивная, рецессивная и доминантная) с поправками на ковариаты – возраст, ИМТ в программе PLINK v. 2.050 [24]. Для оценки ассоциаций полиморфных маркеров генов-кандидатов с прогрессированием заболевания использовали ОR (отношение шансов) и 95%CI (95% доверительный интервал для OR). Проводилась коррекция на множественные сравнения с использованием пермутационных процедур. За статистически значимый уровень принимали рperm<0,05.

Изучение неравновесия по сцеплению между рядом расположенными SNPs, расчёт частот гаплотипов и анализ их ассоциаций с прогрессированием ОА коленного сустава проводили в программе PLINK v. 2.050 [24]. Блочная структура определялась посредством алгоритма «Confidence intervals» (r2≥0,2) [25]. Частоты гаплотипов определялись с помощью EM-алгоритма. Оценку характера ассоциации гаплотипов с прогрессированием заболевания осуществляли с использованием показателя отношения шансов (OR). В исследование включались ковариаты (возраст, ИМТ, распределение по полу), и проводилась коррекция на множественные сравнения при помощи пермутационного теста [26]. За статистически значимый уровень принимали рperm<0,05.

Для полиморфных вариантов генов-кандидатов, показавших значимые ассоциации с прогрессирование ОА коленного сустава, был изучен регуляторный потенциал с использованием онлайн ресурса HaploReg (v4.1) [18, 27] и их связь с экспрессией генов, альтернативным сплайсингом с применением онлайн программы GTExportal [28].

Результаты

Медико-биологические и клинико-анамнестические характеристики изучаемых групп больных представлены в таблице 1. Сформированные группы больных не отличались по росту и наличию соматической сопутствующей патологии (p>0,05) (Таблица 1). Однако в объединённой выборке больных с III-IV стадиями наблюдался достоверно более высокий показатель индекса массы тела (ИМТ) по сравнению с группой больных, имеющих II рентгенологическую стадию (p<0,05), также данные группы пациентов отличались по возрасту (p<0,001) (Таблица 1).

Таблица 1

Характеристика исследуемых групп

Показатели

Больные ОА

III-IV стадий

Больные ОА

II стадии

p

(n = 325)

(n = 175)

Мужчины/женщины, n (%)

133/192

(40,92/59,08)

75/100

(42,86/57,14)

0,75

Средний возраст (М±m), лет

59,71±6,27

51,53±4,31

<0,001

Рост (М±m), см

168,94±7,88

169,68±7,85

0,34

ИМТ (М±m), кг/м2

32,39±5,49

31,06±4,63

0,01

Сопутствующая патология, n (%)

Сердечно-сосудистой системы, n (%)

123 (37,84%)

61 (34,85%)

0,57

Пищеварительной системы, n (%)

39 (12,00%)

21 (12,00%)

1

Мочеполовой системы, n (%)

23 (7,08%)

6 (3,43%)

0,14

Центральной нервной системы, n (%)

37 (11,38%)

15 (8,57%)

0,41

Эндокринных органов, n (%)

35 (10,77%)

16 (9,14%)

0,66

Дыхательной системы, n (%)

42 (11,20%)

17 (9,71%)

0,36

Прочие, n (%)

22 (6,77%)

7 (4,00%)

0,29

Примечание: М — среднее арифметическое; m — стандартная ошибка; p – уровень статистической значимости

 

Популяционно-генетический анализ показал, что для всех рассматриваемых полиморфных локусов генов-кандидатов (rs2820436 и rs2820443 LYPLAL1, rs3771501 TGFA, rs11177 GNL3, rs6976 GLT8D1) в группе больных со II стадией и объединённой группе пациентов с III-IV стадиями наблюдаемое распределение генотипов соответствовало ожидаемому согласно равновесию Харди-Вайнберга (pHWE>0,05) (Таблица 2).

 

Таблица 2

Частоты аллелей и генотипов полиморфных локусов генов-кандидатов

в исследуемых группах больных

Полиморфизм

Аллели, генотипы

Больные ОА

III-IV стадий

n (%)

Больные ОА

II стадии

n (%)

OR (95%CI)

p

rs2820436

С

447 (68,98%)

267 (76,72%)

0,68 (0,50-0,92)

0,013

А

201 (31,02%)

81 (23,28%)

1,48 (1,09-2,02)

А/А

31 (9,57%)

7 (4,02%)

2,53 (1,03-6,44)

0,041

А/С

139 (42,90%)

67 (38,51%)

1,20 (0,81-1,78)

0,394

С/С

154 (47,53%)

100 (57,47%)

0,67 (0,45-0,99)

0,043

Ho/He

(PHWE)

0,429/0,428

(1,000)

0,385/0,357

(0,397)

 

 

Минорный аллель А

(модель аллельная)

 

 

1,48 (1,10-2,00)

0,010

A/A vs. А/С vs. С/С (модель аддитивная)

 

 

1,58 (1,12-2,23)

0,009

C/C vs. А/С + А/А (модель доминантная)

 

 

1,61 (1,06-2,44)

0,024

А/C + С/С vs. А/А (модель рецессивная)

 

 

2,50 (0,99-6,26)

0,052

rs2820443

T

466 (71,69%)

266 (76,00%)

0,80 (0,59-1,09)

0,164

C

184 (28,31%)

84 (24,00%)

1,25 (0,92-1,71)

C/C

27 (8,30%)

13 (7,43%)

1,13 (0,54-2,38)

0,864

C/T

130 (40,00%)

58 (33,14%)

1,35 (0,90-2,01)

0,158

T/T

168 (51,70%)

104 (59,43%)

0,73 (0,50-1,08)

0,118

Ho/He

(PHWE)

0,400/0,406

(0,785)

0,331/0,365

 (0,218)

 

 

Минорный аллель C

(модель аллельная)

 

 

1,25 (0,93-1,68)

0,142

C/C vs. C/T vs. T/T (модель аддитивная)

 

 

1,37 (0,98-1,09)

0,064

T/T vs. C/C + C/T (модель доминантная)

 

 

1,54 (1,01-2,33)

0,043

C/T + T/T vs. C/C (модель рецессивная)

 

 

1,28 (0,59-2,81)

0,533

rs3771501

G

367 (56,46%)

199 (57,18%)

0,97 (0,74-1,27)

0,880

A

283 (43,54%)

149 (42,82%)

1,03 (0,79-1,35)

A/A

69 (21,23%)

27(15,52%)

1,47 (0,88-2,47)

0,155

A/G

145 (44,62%)

95 (54,59%)

0,67 (0,46-0,99)

0,042

G/G

111 (34,15%)

52 (29,89%)

1,22 (0,80-1,85)

0,386

Ho/He

(PHWE)

0,446/0,492

 (0,092)

0,546/0,490

 (0,164)

 

 

Минорный аллель A

(модель аллельная)

 

 

1,03 (0,79-1,34)

0,826

A/A vs. A/G vs. G/G (модель аддитивная)

 

 

0,95 (0,71-1,28)

0,749

G/G vs. A/A + A/G (модель доминантная)

 

 

0,73 (0,47-1,13)

0,160

A/G + G/G vs. A/A (модель рецессивная)

 

 

1,37 (0,79-2,36)

0,258

rs11177

G

375 (57,69%)

189 (54,00%)

1,16 (0,89-1,52)

0,292

A

275 (42,31%)

161 (46,00%)

0,86 (0,66-1,13)

A/A

61 (18,77%)

36 (20,57%)

0,89 (0,55-1,45)

0,713

A/G

153 (47,08%)

89 (50,86%)

0,86 (0,59-1,26)

0,476

G/G

111 (34,15%)

50 (28,57%)

1,30 (0,85-1,97)

0,241

Ho/He

(PHWE)

0,471/0,488

(0,570)

0,509/0,497

(0,879)

 

 

Минорный аллель A

(модель аллельная)

 

 

0,86 (0,66-1,12)

0,261

A/A vs. A/G vs. G/G (модель аддитивная)

 

 

0,83 (0,62-1,12)

0,226

G/G vs. A/A + A/G (модель доминантная)

 

 

0,72 (0,46-1,12)

0,147

A/G + G/G vs. A/A (модель рецессивная)

 

 

0,89 (0,53-1,43)

0,655

rs6976

C

262 (57,28%)

187 (54,68%)

0,80 (0,61-1,07)

0,134

T

270 (42,72%)

155 (45,32%)

1,24 (0,94-1,65)

T/T

62 (19,62%)

35 (20,47%)

0,95 (0,58-1,55)

0,918

C/T

146 (46,20%)

85 (49,71%)

0,87 (0,59-1,28)

0,520

C/C

108 (34,18%)

51 (29,82%)

1,22 (0,80-1,86)

0,382

Ho/He

(PHWE)

0,462/0,489

(0,357)

0,497/0,496

(1,000)

 

 

Минорный аллель T

(модель аллельная)

 

 

0,90 (0,69-1,17)

0,435

T/T vs. C/T vs. C/C (модель аддитивная)

 

 

0,88 (0,65-1,19)

0,401

C/C vs. T/T + C/T (модель доминантная)

 

 

0,76 (0,48-1,21)

0,247

C/T + C/C vs. T/T (модель рецессивная)

 

 

0,96 (0,57-1,63)

0,888

Примечание: OR – отношение шансов, 95% CI – 95% доверительный интервал отношения щансов, p – уровень значимости, Ho/He – наблюдаемая/ожидаемая гетерозиготность, PHWE – уровень значимости отклонения от закона Харди-Вайнберга

 

При сравнительном анализе частот аллелей и генотипов изучаемых полиморфных маркеров генов-кандидатов выявлено, что среди больных III-IV стадий частота аллельного варианта А rs2820436 достоверно больше, чем в группе больных со II стадией заболевания (p=0,013) (Таблица 2). Также установлено, что встречаемость генотипа А/А rs2820436 в группе больных III-IV стадиями в 2,38 раза больше и частота генотипа С/C этого же локуса в 1,21 раза меньше, чем среди пациентов со II рентгенологической стадией ОА коленного сустава (p=0,041 и p=0,043 соответственно). Достоверные различия по частоте генотипа A/G выявлены для локуса rs3771501 гена TGFA: встречаемость данного генотипа в объединённой группе больных III-IV рентгенологических стадий была меньше, чем среди пациентов со II стадией (в 1,22 раза, p=0,042). Таким образом аллельный вариант А и генотип А/А rs2820436 являются факторами риска прогрессирования ОА коленного сустава (OR=1,48 и OR=2,53 соответственно), генотипы С/С rs2820436 и A/G rs3771501 определяют протективную направленность в прогрессировании заболевания (OR=0,67 и OR=0,67 соответственно).

Полученные в рамках четырёх генетических моделей показатели ассоциации изучаемых SNPs с прогрессированием ОА представлены в таблице 2. В рамках аллельной генетической модели выявлено, что частота минорного аллеля A rs2820436 в группе больных с III-IV стадиями достоверно выше, чем в группе пациентов со II стадией (OR=1,48, 95%СI 1,10-2,00, р=0,010, рperm=0,012). Значимые ассоциации с прогрессированием ОА коленного сустава установлены для аллеля А rs2820436 гена LYPLAL1 также согласно аддитивной (OR=1,58, 95%СI 1,12-2,23, р=0,009, рperm=0,010) и доминантной (OR=1,61, 95%СI 1,06-2,44, р=0,024, рperm=0,030) генетическим моделям. Стоит отметить, что в рамках доминантной генетической модели установлена связь аллельного варианта С rs2820443 гена LYPLAL1 c прогрессированием заболевания (OR=1,54, 95%СI 1,01-2,33, р=0,043), однако при проведении пермутационных процедур выявленный уровень ассоциаций не достиг статистически значимого (pperm=0,076).

Установлено неравновесие по сцеплению между парами SNPs, расположенными на одной хромосоме: rs2820436-rs2820443 (r2=0,47) и rs11177 и rs6976 (r2=0,93). Эти пары SNPs образуют соответственно два гаплоблока – H1 и H2 (Таблица 3). Поэтому на следующем этапе работы проведено изучение ассоциаций гаплотипов двух установленных нами гаплоблоков полиморфных маркеров генов-кандидатов с прогрессированием ОА коленного сустава (Таблица 3). Выявлено, что частота гаплотипа AC гаплоблока H1 в группе больных с III-IV стадиями заболевания достоверно выше (в 1,49 раза), чем у пациентов со II стадией (OR=1,83, р=0,002, рperm=0,002). По гаплоблоку H2 значимых ассоциаций с прогрессированием ОА коленного сустава не обнаружено.

Таблица 3

Частоты гаплотипов полиморфных локусов генов-кандидатов в исследуемых группах больных

Гаплоблок

(хромосома)

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

Гаплотип

Частота гаплотипа

OR

P

Больные ОА

III-IV стадий

 (n=325)

Больные ОА

II стадии

 (n=175)

H1

(1 хромосома)

rs2820436-rs2820443

AC

0,240

0,161

1,83

0,002

CC

0,041

0,078

0,59

0,094

AT

0,071

0,072

0,89

0,684

CT

0,649

0,690

0,77

0,109

H2

(3 хромосома)

rs11177-rs6976

AT

0,413

0,450

0,85

0,284

GT

0,014

0,003

2,06

0,342

AC

0,008

0,015

0,43

0,118

GC

0,565

0,532

1,21

0,214

Примечание: OR – отношение шансов; p – уровень статистической значимости

 

Обсуждение

В настоящем исследовании установлены ассоциации полиморфных маркеров генов-кандидатов (rs2820436 и rs2820443 LYPLAL1, rs3771501 TGFA) с прогрессированием ОА коленного сустава у населения Центрального Чернозёмного региона России. Нами выявлено, что аллельный вариант А rs2820436 гена LYPLAL1 является рисковым в прогрессировании заболевания (OR=1,48-1,61), как и генотип A/A данного полиморфизма (OR=2,53). Однако генотип С/С rs2820436 имеет протективную направленность в прогрессировании ОА коленного сустава (OR=0,67). Полученные нами результаты согласуются с литературными данными по этой проблеме. Согласно результатам GWAS исследования [21], выполненном на европейцах, противоположный аллель С rs2820436 гена LYPLAL1 определяет протективную направленность в развитии ОА (OR=0,93, р=2,01x10-09). В GWAS исследовании U. Styrkarsdottir et al. (2018) также показано, что аллель С полиморфного маркера rs2820436 связан с низким риском развития ОА тазобедренного сустава у европейцев (OR=0,93, р=9,4x10-09) [20]. Однако ассоциация данного полиморфизма с ОА коленного сустава в GWAS исследовании [20] не достигла полногеномного уровня значимости у европейцев (OR=0,96, р=1,1x10-05).

Стоит отметить, что имеется ряд исследований, которые демонстрируют связь rs2820436 с признаками, имеющими значимую роль в патогенезе ОА (например, ИМТ). В трансэтническом мета-анализе J. Gong et al. (2018) (изучена выборка из афроамериканцев, латиноамериканцев, азиатов, коренных американцев и американцев европеоидного происхождения) выявлена связь аллельного варианта А rs2820436 с ИМТ, причём данный аллель определяет повышенный ИМТ (β=0,0049, р=3,79×10−08) [29]. А как известно, повышенный ИМТ и ожирение, являются ведущими фактором риска развития и прогрессирования ОА [8, 30-32]. В работе L. Fernández-Rhodes et al. (2017) на выборке из более 102 000 взрослых африканцев, латиноамериканцев, азиатов и европейцев, напротив, установлена ассоциация аллеля А rs2820436 с более низким индексом массы тела (OR=0,50, р=3,2x10-08) [33]. В исследовании [34] показана связь экспрессии rs2820436 гена LYPLAL1 в жировой ткани с ИМТ (р=1,2x10-04) и уровнем базального липолиза (р=0,026). Стоит отметить, что в доступной литературе не обнаружены исследования, посвящённые изучению роли полиморфного генетическогого маркера rs2820436 гена LYPLAL1 в прогрессировании ОА.

Медико-биологической основой выявленных ассоциаций могут быть выраженные функциональные эффекты rs2820436 гена LYPLAL1. C помощью онлайн-ресурса HaploReg (v.4.1) [18] установлено, что данный полиморфный маркер расположен в области ДНК, связанной с гистоновыми белками H3K4me1, маркирующими неактивные энхансеры в хондроцитах, фибробластах, остеобластах, жировой ткани, находится в регионе гистона H3K27ac, маркирующего активные энхансеры в лимфобластах, хондроцитах, а также в области гистона H3K9ac, определяющего активные промоторы в жировой ткани, регионе трёх регуляторных мотивов ДНК (MZF1, PLAG1, TFIIA). Стоит отметить, что согласно нашим данным минорный аллель А rs2820436 LYPLAL1 является фактором риска прогрессирования ОА коленного сустава, снижает аффинность к факторам транскрипции PLAG1, TFIIA (различие между LOD scores аллелей А (alt) и С (ref) составляет -11,9 и -11,7 соответственно). Также с применением онлайн программы GTExportal [28] выявлено, что аллельный вариант А rs2820436 связан с повышенной экспрессией гена RP11-392O17.1 в жировой ткани (p=1,9х10-06) и низкой экспрессией этого же гена в надпочечниках (p=4,0х10-06), определяет пониженный уровень альтернативного сплайсинга транскрипта гена RP11-392O17.1 в жировой ткани (p=4,8 Е-21).

Аллельный вариант С rs2820443 гена LYPLAL1, который, согласно нашим данным, в составе гаплотипа АС гаплоблока rs2820436-rs2820443 гена LYPLAL1 является фактором риска прогрессирования ОА коленного сустава у населения европейской части России, в GWAS исследовании [22] имеет рисковое значение для развития ОА коленного и/или тазобедренного сустава у европейцев (OR=1,06, р=6x10-11). Стоит отметить, что в каталоге GWAS представлены исследования, демонстрирующие на полногеномном уровне связь rs2820443 гена LYPLAL1 с патогенетически значимыми для ОА коленного сустава признаками (распределение абдоминального жира, окружность бедра). Согласно данным GWAS [35-36], аллель С rs2820436 связан с более высоким коэффициентом распределения абдоминальной жировой ткани в организме (β=0,025; р=8х10-20 для европейцев) [35] и большей окружностью бедра с поправкой на ИМТ (β=0,060; р=3х10-83 для европейской популяции) [36]. В полногеномном исследовании [37] показана ассоциация противоположного аллеля Т rs2820436 гена LYPLAL1 с меньшей окружностью бедра с поправкой на ИМТ (β=0,048; р=1х10-35 в европейской выборке).

Выявленные в рамках нашего исследования ассоциации rs2820443 гена LYPLAL1 с прогрессированием ОА коленного сустава могут быть обусловлены выраженными регуляторными эффектами данного SNP [18, 28]. Полиморфный локус rs2820443, находящийся на расстоянии 22kb от 3' конца гена RP11-95P13.2, расположен в эволюционно консервативном регионе, в регионе гистонов, маркирующих промоторы в жировой ткани и энхансеры в остеобластах, мышечной ткани, фибробластах, в области гиперчувствительности к ДНКазе-1, регионе связывания с двумя регуляторными белками (STAT3, FOXA1) и регионе 8 регуляторных мотивов ДНК [18]. Минорный аллель С rs2820443, который, согласно нашим данным, в составе гаплотипа АС гаплоблока rs2820436-rs2820443 является фактором риска прогрессирования заболевания, повышает аффинность к трансформационному фактору DBP (различия между LOD scores аллелей C(alt) и T(ref) составляют 11,8) и понижает аффинность к транскрипционным факторам: Fox, Foxc1, Foxd1, Foxj1, Foxk1, Foxo, Sox (различия между LOD scores аллелей C(alt) и T(ref) = -4 – -12) [18]. Данный аллельный вариант связан с низкой экспрессией гена RIMKLBP2 в жировой ткани (β=-0,17, р=3,6x10-05) и низким уровнем альтернативного сплайсинга транскрипта гена RP11-392O17.1 в жировой ткани (β=-0,47, р=2,1x10-11) [28].

Нами установлена связь генотипа A/G rs3771501 гена TGFA с прогрессированием ОА коленного сустава (OR=0,67). Ассоциация аллеля G rs3771501 гена TGFA с ОА на GWAS уровне впервые была установлена E. Zengini et al. в 2018 году в европейской популяции [21]. Авторами показано, что данный аллельный вариант rs3771501 является протективным фактором развития ОА (OR=0,94, p=1,66x10-08). В GWAS исследовании C.G. Boer et al. (2021) установлена связь rs3771501 гена TGFA с ОА суставов руки у европейцев (p<5,8x10–4) [38]. Ещё в двух GWAS исследованиях [22, 23] выявлены ассоциации аллельного варианта А rs3771501 гена TGFA с ОА в европейской популяции (OR=1,05, p=4,27x10–16) [22] и на смешанных выборках европейцев, азиатов, американцев европейского происхождения (OR=1,04, p=4,05x10–15) [23].

По данным базы генетической информации GeneCards: The Human Gene Database [39] ген TGFA (трансформирующий фактор роста альфа) кодирует фактор роста, являющегося лигандом для рецептора эпидермального фактора роста, который активирует сигнальный путь для пролиферации, дифференцировки и развития клеток. Ещё в работах C.T. Appleton et al. (2007) сообщалось о повышенном уровне TGFA в хондроцитах больных ОА коленного сустава [40]. Известно также, что TGFA подавляет анаболические и способствует катаболическим процессам в суставном хряще. TGFA является мощным стимулятором деградации хряща за счёт активации сигнальных путей Rho/ROCK и MEK/ERK [41]. C.T. Appleton et al. (2015) установлено, что ингибирование сигнального пути TGFA-CCL2 снижает риск прогрессирования экспериментального посттравматического ОА in vivo [42].

В литературе имеются данные о связи rs2862851 гена TGFA, находящегося в сильном сцеплении с изучаемым полиморфным маркером rs3771501 (r2=1), с ОА. В исследовании [43] установлена ассоциация аллеля T rs2862851 гена TGFA с риском развития ОА коленного сустава в китайской популяции Хань (OR =1,41, p =1,68x10-4). Кроме того, показано, что концентрации TGFA в крови у больных ОА коленного сустава были выше, чем у здоровых людей (p=8,15x10-13), а также установлено увеличение концентрации TGFA в крови с возрастанием тяжести заболевания (p =1,36x10-72).

Полиморфный маркер rs3771501 гена TGFA обладает выраженными регуляторными эффектами [18]: локализован в функционально активном регионе генома, в области ДНК, связанной с гистоновыми белками (H3K4me1), маркирующими неактивные энхансеры в клетках хондроцитов, гиппокампе, черной субстанции, височной доле мозга, коре. С помощью онлайн ресурса GTExportal [28] выявлено, что аллель G rs3771501 ассоциирован с повышенной экспрессией гена TGFA в большеберцовых нервах (β=0,09, р=1,6x10-04), в различных отделах головного мозга (β=0,21-0,31, р<3,1x10-06).

Таким образом, молекулярно-генетические маркеры rs2820436 и rs2820443 гена LYPLAL1, rs3771501 гена TGFA связаны как с риском развития ОА по данным ранее проведённых полногеномных исследований, так и, согласно нашим данным, ассоциированы с прогрессированием ОА коленного сустава у жителей Центрального Черноземья России.

Заключение

В рамках настоящего исследования установлена связь полиморфных маркеров rs2820436 и rs2820443 гена LYPLAL1, rs3771501 TGFA с прогрессированием ОА коленного сустава у населения Центрального Черноземья России. Аллельный вариант А и генотип A/A rs2820436 гена LYPLAL1 (OR=1,48-2,53), гаплотип АС гаплоблока rs2820436-rs2820443 (OR=1,83) являются факторами риска прогрессирования заболевания. Гепотипы C/C rs2820436 гена LYPLAL1 (OR=0,67) и A/G rs3771501 гена TGFA (OR=0,67) имеют протективное значение в прогрессировании ОА коленного сустава.

Заявленный вклад авторов.

Новаков В.Б. – литературный обзор публикаций по теме статьи, редактирование.

Новакова О.Н. –сбор и обработка материала, литературный обзор публикаций по теме статьи и написание текста рукописи, этапное и заключительное редактирование рукописи.

Чурносов М.И. – разработка дизайна исследования, статистическая обработка полученного материала, анализ полученных результатов, утверждение окончательного текста статьи.

Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования.

Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.

Этическая экспертиза. Проведение исследования одобрено локальным этическим комитетом при Медицинском институте ФГАОУ ВО «Белгородский государственный национальный исследовательский университет», протокол № 2 (7) от 02.02.2016 г.

Информированное согласие на публикацию. Авторы получили письменное согласие пациентов на участие в исследовании и публикацию медицинских данных.

Authors’ contributions.

Vitaly B. Novakov – literary review of publications on the topic of the article, editing.

Olga N. Novakova – collection and processing of material, literary review of publications on the topic of the article and writing the text of the manuscript, stage and final editing of the manuscript.

Mikhail I. Churnosov – development of the study design, statistical processing of the obtained material, analysis of the results, approval of the final text of the article.

Funding source. The authors declare no external funding for the study.

Competing interests The authors declare the absence of obvious and potential conflicts of interest related to the publication of this article.

Ethics approval. The study was approved by the local ethics committee at the Medical Institute of the Belgorod State National Research University, protocol No. 2 (7) dated February 2, 2016.

Consent for publication. The authors obtained written consent from patients to participate in the study and publish medical data.

СВЕДЕНИЯ ОБ АВТОРАХ:

  1. Новаков Виталий Борисович – врач травматолог-ортопед, ОГБУЗ «Городская больница №2 г. Белгорода», Белгород, Россия; аспирант медицинского института, ФГАОУ ВО «Белгородский государственный национальный исследовательский университет» Министерства образования и науки Российской Федерации, Белгород, Россия
    Email: v.novakov@bk.ru

ORCID: https://orcid.org/0000-0001-5337-2156

  1. Новакова Ольга Николаевна - к. б. н., доцент кафедры медико-биологических дисциплин медицинского института, ФГАОУ ВО «Белгородский государственный национальный исследовательский университет» Министерства образования и науки Российской Федерации, Белгород, Россия

Email: litovkina@bsu.edu.ru

ORCID: https://orcid.org/0000-0003-2700-1284

  1. Чурносов Михаил Иванович - д. м. н., профессор, заведующий кафедрой медико-биологических дисциплин Медицинского института, ФГАОУ ВО «Белгородский государственный национальный исследовательский университет» Министерства образования и науки Российской Федерации, Белгород, Россия

Email: churnosov@bsu.edu.ru

ORCID: https://orcid.org/0000-0003-1254-6134

AUTHORSʼ INFORMATION:

  1. Vitaly B. Novakov - Orthopaedic Surgeon, City Hospital №2 of Belgorod, Belgorod, Russia; PhD Student, Belgorod National Research University, Belgorod, Russia

Email: v.novakov@bk.ru

ORCID: https://orcid.org/0000-0001-5337-2156

  1. Olga N. Novakova - Cand. Sci. (Biol.), Associate Professor of the Department of Biomedical Disciplines of the Medical Institute, Belgorod National Research University, Belgorod, Russia
    Email: litovkina@bsu.edu.ru.

ORCID: https://orcid.org/0000-0003-2700-1284

  1. Mikhail I. Churnosov - Dr. Sci. (Med.), Professor, Head of the Department of Medical and Biological Disciplines of the Medical Institute, Belgorod National Research University, Belgorod, Russiа
    Email: churnosov@bsu.edu.ru.

ORCID: https://orcid.org/0000-0003-1254-6134

×

About the authors

Vitaly Novakov

Belgorod National Research University;
City Hospital No. 2 of Belgorod

Email: v.novakov@bk.ru
ORCID iD: 0000-0001-5337-2156

Orthopedic Surgeon of the Belgorod City Hospital No. 2, postgraduate student of the Medical Institute of the Belgorod National Research University

Russian Federation, 308015, Belgorod, Pobeda St., 85; 308036, Belgorod, Gubkina St., 46

Olga Novakova

Belgorod National Research University

Author for correspondence.
Email: litovkina@bsu.edu.ru
ORCID iD: 0000-0003-2700-1284
SPIN-code: 1001-8369
Scopus Author ID: 57193996396

Cand. Sci. (Biol.), Associate Professor of the Department of Biomedical Disciplines

Russian Federation, 308015, Belgorod, Pobeda St., 85

Mikhail Churnosov

Belgorod National Research University

Email: churnosov@bsu.edu.ru
ORCID iD: 0000-0003-1254-6134
SPIN-code: 7407-9649
Scopus Author ID: 6601948788

Dr. Sci. (Med.), Professor, Head of the Department of Biomedical Disciplines

308015, Belgorod, Pobeda St., 85

References

  1. Reynard L.N., Barter M.J. Osteoarthritis year in review 2019: genetics, genomics and epigenetics. Osteoarthritis Cartilage. 2020;28(3):275-284. doi: 10.1016/j.joca.2019.11.010.
  2. Katz J.N., Arant KR, Loeser RF. Diagnosis and Treatment of Hip and Knee Osteoarthritis: A Review. JAMA. 2021;325(6):568-578. doi: 10.1001/jama.2020.22171.
  3. Alekseeva L.I., Taskina E.A., Kashevarova N.G. [Osteoarthritis: epidemiology, classification, risk factors and progression, clinic, diagnosis, treatment]. Sovremennaya revmatologiya. 2019;13(2): 9-21. (In Russian). doi: 10.14412/1996-7012-2019-2-9-21.
  4. Panikar V.I., Shcherban E.A., Pavlova I.A. [Complex geriatric assessment of osteoarthrosis of knee joints in the senior age]. Nauchnyye rezulʹtaty biomeditsinskikh issledovaniy [Research Results in Biomedicine]. 2019;5(1):131-139. (In Russian). doi: 10.18413/2313-8955-2019-5-1-0-10.
  5. Galushko E.A., Nasonov E.L. [Prevalence of rheumatic diseases in Russia]. Alʹmanakh klinicheskoy meditsiny. 2018;46(1):32-39. (In Russian). doi: 10.18786/2072-0505-2018-46-1-32-39.
  6. Kavalersky G., Smetanin S., Lychagin A., Moisov A. [Faktory riska razvitiya osteoartroza kolennogo sustava]. Vrach. 2017;3:22-24. Available from: https://vrachjournal.ru/ru/25877305-2017-03-06 (15.06.2022). (In Russian).
  7. Matveev R.P., Bragina S.V. [Knee osteoarthritis: problems and social significance]. Ekologiya cheloveka. 2012;9:53-62. (In Russian). doi: 10.33396/1728-0869-2012-9-53-62.
  8. Kashevarova N.G., Alekseeva L.I. [Faktory riska progressirovaniya osteoartroza kolennykh sustavov]. Nauchno-prakticheskaya revmatologiya [Rheumatology Science and Practice]. 2014;52(5):553-561. (In Russian). doi: 10.14412/1995-4484-2014-553-561.
  9. Van Meurs J.B. Osteoarthritis year in review 2016: genetics, genomics and epigenetics. Osteoarthritis Cartilage. 2017;25(2):181-189. doi: 10.1016/j.joca.2016.11.011.
  10. Primorac D., Molnar V., Rod E., Jeleč Ž., Čukelj F., Matišić V. et al. Knee Osteoarthritis: A Review of Pathogenesis and State-Of-The-Art Non-Operative Therapeutic Considerations. Genes (Basel). 2020;11(8):854. doi: 10.3390/genes11080854.
  11. Novakov V.B., Novakova O.N., Churnosov M.I. [Risk factors and molecular entities of the etiopathogenesis of the knee osteoarthritis (literature review)]. Genij Ortopedii. 2021;27(1):112-120. (In Russian). doi: 10.18019/1028-4427-2021-27-1-112-120.
  12. Valdes A.M., Hart D.J., Jones K.A., Surdulescu G., Swarbrick P., Doyle D.V., Schafer A.J., Spector T.D. Association study of candidate genes for the prevalence and progression of knee osteoarthritis. Arthritis Rheum. 2004;50(8):2497-507. doi: 10.1002/art.20443. PMID: 15334463.
  13. Abd Elazeem M.I., Abdelaleem E.A., Mohamed R.A. Genetic influence of growth and differentiation factor 5 gene polymorphism (+104T/C) on the development of knee osteoarthritis and its association with disease severity. Eur J Rheumatol. 2017;4(2):98-103. doi: 10.5152/eurjrheum.2017.160093.
  14. Gok K., Cemeroglu O., Cakirbay H., Gunduz E., Acar M., Cetin E.N. et al. Relationship between cytosine-adenine repeat polymorphism of ADAMTS9 gene and clinical and radiologic severity of knee osteoarthritis. Int J Rheum Dis. 2018;21(4):821-827. doi: 10.1111/1756-185X.12849.
  15. Wu X., Kondragunta V., Kornman K.S., Wang H.Y., Duff G.W., Renner J.B. et al. IL-1 receptor antagonist gene as a predictive biomarker of progression of knee osteoarthritis in a population cohort. Osteoarthritis Cartilage. 2013;21(7):930-8. doi: 10.1016/j.joca.2013.04.003.
  16. Chen R., Zhang Y., Xu H., Hu H., Chen M., Shuai Z. Val109Asp Polymorphism of the Omentin-1 Gene and Incidence of Knee Osteoarthritis in a Chinese Han Population: A Correlation Analysis. Drug Des Devel Ther. 2021;15:5075-5086. doi: 10.2147/DDDT.S340410.
  17. Catalog GWAS (Genome-Wide Association Studies). Available from: http://www.genome.gov/gwastudies/.
  18. HaploReg v4.1. Available from: https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php.
  19. ArcOGEN Consortium; arcOGEN Collaborators, Zeggini E., Panoutsopoulou K., Southam L., Rayner N.W., Day-Williams A.G., Lopes M.C. et al. Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study. Lancet. 2012;380(9844):815-23. doi: 10.1016/S0140-6736(12)60681-3.
  20. Styrkarsdottir U., Lund S.H., Thorleifsson G., Zink F., Stefansson O.A., Sigurdsson J.K. et al. Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis. Nat Genet. 2018;50(12):1681-1687. doi: 10.1038/s41588-018-0247-0.
  21. Zengini E., Hatzikotoulas K., Tachmazidou I., Steinberg J., Hartwig F.P., Southam L. et al. Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis. Nat Genet. 2018;50(4):549-558. doi: 10.1038/s41588-018-0079-y.
  22. Tachmazidou I., Hatzikotoulas K., Southam L., Esparza-Gordillo J., Haberland V., Zheng J. et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet. 2019;51(2):230-236. doi: 10.1038/s41588-018-0327-1.
  23. Boer C.G., Hatzikotoulas K., Southam L., Stefánsdóttir L., Zhang Y., Coutinho de Almeida R. et al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. 2021;184(18):4784-4818.e17. doi: 10.1016/j.cell.2021.07.038.
  24. Purcell S., Neale B., Todd-Brown K., Thomas L., Ferreira M.A., Bender D. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559-75. doi: 10.1086/519795.
  25. Gabriel S.B., Schaffner S.F., Nguyen H., Moore J.M., Roy J., Blumenstiel B. et al. The structure of haplotype blocks in the human genome. Science. 2002;296(5576):2225-9. doi: 10.1126/science.1069424.
  26. Che R., Jack J.R., Motsinger-Reif A.A., Brown C.C. An adaptive permutation approach for genome-wide association study: evaluation and recommendations for use. BioData Min. 2014;7:9. doi: 10.1186/1756-0381-7-9.
  27. Polonikov A.V., Klyosova E.Yu., Azarova I.E. [Bioinformatic tools and internet resources for functional annotation of polymorphic loci detected by genome wide association studies of multifactorial diseases (review)]. Research Results in Biomedicine. 2021;7(1):15-31. (In Russian). doi: 10.18413/2658-6533-2020-7-1-0-2.
  28. GTEx (Genotype-Tissue Expression) portal. Available from: http://www.gtexportal.org/.
  29. Gong J., Nishimura K.K., Fernandez-Rhodes L., Haessler J., Bien S., Graff M. et al. Trans-ethnic analysis of metabochip data identifies two new loci associated with BMI. Int J Obes (Lond). 2018;42(3):384-390. doi: 10.1038/ijo.2017.304.
  30. Kulkarni K., Karssiens T., Kumar V., Pandit H. Obesity and osteoarthritis. Maturitas. 2016;89:22-8. doi: 10.1016/j.maturitas.2016.04.006.
  31. Belluzzi E., El Hadi H., Granzotto M., Rossato M., Ramonda R., Macchi V. et al. Systemic and Local Adipose Tissue in Knee Osteoarthritis. J Cell Physiol. 2017;232(8):1971-1978. doi: 10.1002/jcp.25716.
  32. Lespasio M.J., Piuzzi N.S., Husni M.E., Muschler G.F., Guarino A., Mont M.A. Knee Osteoarthritis: A Primer. Perm J. 2017;21:16-183. doi: 10.7812/TPP/16-183.
  33. Fernández-Rhodes L., Gong J., Haessler J., Franceschini N., Graff M., Nishimura K.K. et al. Trans-ethnic fine-mapping of genetic loci for body mass index in the diverse ancestral populations of the Population Architecture using Genomics and Epidemiology (PAGE) Study reveals evidence for multiple signals at established loci. Hum Genet. 2017;136(6):771-800. doi: 10.1007/s00439-017-1787-6.
  34. Lundbäck V., Kulyte A., Strawbridge R.J., Ryden M., Arner P., Marcus C. et al. FAM13A and POM121C are candidate genes for fasting insulin: functional follow-up analysis of a genome-wide association study. Diabetologia. 2018;61(5):1112-1123. doi: 10.1007/s00125-018-4572-8.
  35. Rask-Andersen M., Karlsson T., Ek W.E., Johansson Å. Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects. Nat Commun. 2019;10(1):339. doi: 10.1038/s41467-018-08000-4.
  36. Christakoudi S., Evangelou E., Riboli E., Tsilidis K.K. GWAS of allometric body-shape indices in UK Biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci Rep. 2021;11(1):10688. doi: 10.1038/s41598-021-89176-6.
  37. Shungin D., Winkler T.W., Croteau-Chonka D.C., Ferreira T., Locke A.E., Mägi R. et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518(7538):187-196. doi: 10.1038/nature14132.
  38. Boer C.G., Yau M.S., Rice S.J., Coutinho de Almeida R., Cheung K., Styrkarsdottir U. et al. Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene. Ann Rheum Dis. 2021;80(3):367-375. doi: 10.1136/annrheumdis-2020-217834.
  39. GeneCards: The Human Gene Database Available from: http://www.genecards.org/.
  40. Appleton C.T., Usmani S.E., Bernier S.M., Aigner T., Beier F. Transforming growth factor alpha suppression of articular chondrocyte phenotype and Sox9 expression in a rat model of osteoarthritis. Arthritis Rheum. 2007;56(11):3693-705. doi: 10.1002/art.22968.
  41. Appleton C.T., Usmani S.E., Mort J.S., Beier F. Rho/ROCK and MEK/ERK activation by transforming growth factor-alpha induces articular cartilage degradation. Lab Invest. 2010;90(1):20-30. doi: 10.1038/labinvest.2009.111.
  42. Appleton C.T., Usmani S.E., Pest M.A., Pitelka V., Mort J.S., Beier F. Reduction in disease progression by inhibition of transforming growth factor α-CCL2 signaling in experimental posttraumatic osteoarthritis. Arthritis Rheumatol. 2015;67(10):2691-701. doi: 10.1002/art.39255.
  43. Cui G., Liu D., Wei R., Wu J., Liu R., Wang K. Association of rs2862851 in TGFA Gene with Peripheral TGFA Levels and the Severity of Knee Osteoarthritis in the Han Chinese Population. Genet Test Mol Biomarkers. 2020;24(12):771-776. doi: 10.1089/gtmb.2020.0119.

Supplementary files

There are no supplementary files to display.


Copyright (c) Eco-Vector

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: серия ПИ № ФС 77 - 82474 от 10.12.2021.


This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies