50 2.2.3. Генетико – статистические методы анализа

Распределение генотипов по исследованным полиморфным локусам проверяли на соответствие равновесию Харди-Вайнберга (РХВ) с помощью точного теста Фишера [Вейр Б., 1995]. Рассчитывали ожидаемую гетерозиготность полиморфизма генов NRAMP1, IL12B, VDR, IL1B, IL1RN [Nei M., 1975].

D=(hobs–hexp)/hexp,

где hobs и hexp – ожидаемая и наблюдаемая гетерозиготность соответственно.

Для анализа ассоциации маркеров исследуемых генов с туберкулезом, а также с качественными патогенетически важными признаками заболевания, сравнивали частоты аллелей и генотипов в группах больных и здоровых индивидов, используя критерий χ2 с поправкой Йетса на непрерывность. При численностях генотипов менее пяти использовали точный тест Фишера. В дополнение к этому об ассоциации разных генотипов (или их комбинаций) с заболеванием судили по величине отношения шансов (odds ratio (OR)), которая показывает, во сколько раз выше вероятность заболеть для индивида с определенным генотипом (или комбинацией генотипов) [Pearce N., 1993].

OR= (A/B)/(C/D), где

А – число (процент) людей с данным генотипом (комбинацией генотипов) в группе больных;

С - число (процент) людей с данным генотипом (комбинацией генотипов) в группе здоровых;

В – число (процент) индивидов, не имеющих данного генотипа (комбинации генотипов) в группе больных;

D - число (процент) индивидов, не имеющих данного генотипа (комбинации генотипов) в группе здоровых.

Значения OR>1 указывают на возможную положительную ассоциацию с заболеванием. Обсуждение величин OR проводили при уровне значимости не более 5%.

На материале семейной выборки больных изучение ассоциаций полиморфизма исследованных генов с туберкулезом проводили с использованием теста на неравновесие при переносе (Transmission/Disequilibrium Test, TDT), который в случае диаллельного маркерного локуса М сводится к анализу таблицы сопряженности 2´2, где в ячейках матрицы суммированы случаи наследования и не наследования от родителей больными детьми маркерных аллелей [Spielman R. S.

a – число случаев наследования аллеля М1 от родителей М1М1;

b – число случаев наследования аллеля М1 от родителей М1М2;

c – число случаев наследования аллеля М2 от родителей М1М2;

d – число случаев наследования аллеля М2 от родителей М2М2;

Используются данные только от гетерозиготных родителей. Статистика теста рассчитывается по формуле:

TDT=(b-c)2/(b+c)

и в случае верной нулевой гипотезы (Н0: нет ассоциации) асимптотически распределена как χ2 с 1 степенью свободы.

С целью выявления ассоциации маркеров исследуемых генов с количественными, патогенетически важными признаками туберкулеза, проводили сравнение средних значений уровней метрических показателей у носителей разных генотипов с помощью однофакторного дисперсионного анализа по Фишеру и теста LSD. При наличии зависимости признака от пола показатели анализировались отдельно в группе мужчин и женщин. В случае влияния возраста на количественный параметр проводилась его корректировка, которая осуществлялась с помощью уровня линейной регрессии и рассчитывалась по формуле [Лильин Е.Т. и др., 1984]:

y=x+b(t0 -t),

где y – коррегированное значение исходной величины (х) признака;

t – возраст индивида

t0 – определенный возраст, к которому приводятся все значения;

b – коэффициент линейной регрессии признака по возрасту, который рассчитывается по формуле:

b=rxt/st2

где rxt – коэффициент корреляции признака с возрастом;

st - стандартное отклонение возраста в выборке.

Проверку на нормальность распределений осуществляли с помощью критерия Колмогорова-Смирнова и Лилифорса. В случае неравных дисперсий использовали непараметрические тесты Манна-Уитни, Краскела-Уоллиса и медианный тест [Лакин Г.Ф., 1990]. Сравнение дисперсий проводили по критерию Левене.

Расчеты гаметического неравновесия между парами молекулярно-генетических маркеров проводили по Hill W. G. (1974). Все расчеты осуществляли с помощью программ «STATISTICA for Windows 6.0» и «Microsoft Excel 7.0».