CC BY
25
© Р.Р. Фурина, Я.А. Фурман, С.А. Казанкова, С.А. Казанков, Е.А. Голодюк, М.О. Гремякина, Н.Н. Митракова, 2018
УДК 54.061:543.06:577.121:612.015.3:616.15:616.633:616.24-008.8-074
СТАТИСТИЧЕСКИЙ АНАЛИЗ КОЛИЧЕСТВЕННОГО ПРЕДСТАВЛЕНИЯ ЛЕТУЧИХ МЕТАБОЛИТОВ МОЧИ КАК ВОЗМОЖНЫХ ОНКОМАРКЕРОВ
Р.Р. Фурина1,2, Я.А. Фурман1, С.А. Казанкова1, С.А. Казанков1, Е.А. Голодюк1,2, М.О. Гремякина1, Н.Н. Митракова1,2
1ГБУ РМЭ «Республиканская клиническая больница», г. Йошкар-Ола
2ФГБОУ ВПО «Государственный Поволжский технический университет», г. Йошкар-Ола
STATISTICAL ANALYSIS OF QUANTITATIVE PRESENTATION OF VOLATILE METABOLITIS OF URINE AS POSSIBLE ONCOMMARKERS
R.R. Furina1,2, Ya.A. Furman1, S.A. Kazankova1, S.A. Kazankov1, E.A. Golodyuk1,2, M.O. Gremyakina1, N.N. Mitrakova1,2
Republican Clinical Hospital, Yoshkar-Ola 2Volga State Technological University, Yoshkar-Ola
Митракова Нина Николаевна - доктор медицинских наук, профессор кафедры радио-технических и медико-биологических систем ФГБОУ ВПО «Государственный Поволжский технический университет», заведующая эндоскопическим отделением ГБУ РМЭ «Республиканская клиническая больница»
424030, г. Йошкар-Ола, ул. Осипенко, д. 33, тел.: (362) 42-64-54, +7-927-871-90-22, e-mail: endomitrakova@mail.ru Mitrakova N.N. - D. Med. Sc., Professor of the Department of Radio-Technical and Biomedical Systems of Volga State Technical University, Head of the Endoscopic Department of the Republican Clinical Hospital
33 Osipenko Str., Yoshkar-Ola, Russian Federation, 424030, tel.: (362) 42-64-54, +7-927-871-90-22, e-mail: endomitrakova@mail.ru
Реферат. Решается задача выбора информативных признаков летучих метаболитов мочи пациентов как онкомаркеров для диагностики рака легких. На хромато-масс-спектрограммах, полученных на хроматографе «Хроматэк-Кристалл-5000», определены характеристики порядка 150 метаболитов, выделенных из мочи пациентов. Выполнен предварительный отбор 16-ти метаболитов для дальнейших исследований. Проведен двухэтапный математический анализ хроматографиче-ских данных. В качестве информативных признаков были приняты значения высоты h и площади s пиков каждого метаболита. На первом этапе анализа проверялась справедливость нулевой гипотезы для выборочных данных параметров h и s. Показано, что только для метаболита Disulfide dimethyl нулевая гипотеза уверенно отвергается. На втором этапе анализа параметров данного метаболита методом максимального правдоподобия получены оценки вероятностей правильной диагностики рака легких для контрольной группы и группы, больных раком легкого, соответственно равные 0,733 и 0,677. Ключевые слова: онкомаркеры, метаболиты, диагностика рака легких.
Abstract. The problem of the choice of informative signs of the urine volatile metabolites of patients as oncomarkers for the diagnosis of lung cancer is being solved. The chromatography-mass spectrograms obtained on chromatograph «Chromatech-Crystal-5000» determined the characteristics of about 150 metabolites isolated from urine of patients. Preliminary selection of 16 metabolites for further studies was performed. A two-stage mathematical analysis of chromatographic data was carried out. As information signs, the heights and peak areas of each metabolite were taken. At the first stage of the analysis, the validity of the null hypothesis for the sample data of the parameters and was verified. It is shown that only for the metabolite Disulfide dimethyl the null hypothesis is surely rejected. The second stage of the analysis of the parameters of this metabolite using the maximum likelihood method was used to estimate the probabilities of correct diagnosis of lung cancer for the control group and the group with lung cancer, respectively, equal to 0.733 and 0.677. Кеу words: tumor markers, metabolites, lung cancer diagnosis.
Введение
На сегодня во всем мире проводятся активные исследования по обнаружению онкомаркеров. Большие надежды в данном вопросе возлагаются на молекулярную биологию, активно развивающуюся в последнее время. Результатом расшифровки генома человека является возникновене ряда научных направлений принципиально нового характера, получивших собирательное название «постгеномные технологии». К ним относятся метаболомика, изучающая совокупность метаболитов (метаболический профиль) у пациента при данной патологии [1, 5]. Поэтому на метаболомику возлагаются большие надежды в поиске метаболитов - биомаркеров заболеваний, и в первую очередь онкологических.
Одним из направлений подобных исследований служит разработка простого, неинвазивного, обладающего низкой стоимостью метода диагностики онкопатологии при диспансеризации населения и, в перспективе, для самостоятельного обследования с использованием тест-систем [1, 2, 5, 6]. В последнем
десятилетии активно обсуждаются методы, основанные на анализе метаболизма человека при наличии онкологических заболеваний [3, 4, 12]. Среди них существенную роль играют методы, основанные на хроматографии с масс-спектрометрическим детектированием [9, 10, 11]. Применение хромато-масс-спектрометрии в обнаружении рака получила широкое применение благодаря работам Филлипса и его коллег, разработавших «электронный нос» для автоматизированного анализа дыхания и обнаружения летучих органических соединений для различных видов онкопатологии [7, 8].
Цель данной работы состоит в разработке двух последовательно выполняемых процедур. Первая из них является тест-пробой для разделения сформированных выборок с количественными данными метаболита в контрольной и рабочей группах на два класса: группа с информативным метаболитом и группа с нонинформативным метаболитом. Вторая процедура является классическим статистическим анализом, в результате которого для группы с информативным метаболитом принимается диагностическое решение.
Рис. 1. Шприц ТФМЭ
Методика подготовки материалов для анализа состава летучих метаболитов мочи
Исследования проводились в лаборатории компании «Хроматэк» на базе газового хроматографа «Хро-матэк-Кристалл 5000» отечественного производства. Для подготовки пробы использовался метод твердофазной микроэкстракции (ТФМЭ). Метод основан на сорбции компонентов газовой фазы образца на нить с полимерным покрытием и ее последующей термической десорбцией в предварительно нагретом инжекторе газового хроматографа (рис. 1).
Суть газовой хроматографии заключается в разделении смеси образца на компоненты, перемещающиеся по хроматографической колонке с потоком газа-носителя. По мере движения разделяемая смесь многократно распределяется между газом-носителем и нелетучей неподвижной жидкой фазой. Масс-спектрометрический детектор обнаруживает в потоке газа-носителя анализируемые вещества по какому-либо физико-химическому свойству. Отклик осуществляется за счет преобразования свойств в электрический сигнал. В масс-спектрометрическом детекторе ионизированные компоненты разделяются в соответствии с их m/z (отношение массы иона к заряду), и регистрируемый ионный ток для каждого значения m/z преобразуется в масс-спектр.
£ Mjtfbiait EIH.UI Ak'kiCL'v kg 'ЕШЙ
иА) <frtfr г-ама ВС авм Фшнт СЕМГ Aswit.e jf.s-. Qfuou ' IjEfftwdJiTn^nnut " - ж
^'nTjJRLr.^ Iff * -ЛЯ1 ЛЛШТ1., .А , A - 1
13СЭ v ft Mi
A -L. в. . с О U I F С - if: : -
i 4чПл (иниршшии: F:\temp\l\ArialiCve I.D\Prol»cii\M9020)\6Го\Ноя6рь\1й.ID. 12\Aletv ks .uiidn
I Ная*ня« илнала: M 44,5 - 300 j I
3 Nt 'Ншанщ Время Высота Площадь Плоивдь, "i в(рСШТНОСГЬ -1
a 1 0,0 J! D,0'J Cl.iKi o.ao -
5 L_ * Carbon cii rle 4,2? 41G144.2S 833140,12 0,96 *937 ■ " ........- '
6 3 Carbonyl sulfide 4,34 143700.99 204492.02 0,31
J 4 IVimil h^ttilyl ftuorfd 4.46 210396,« 411742,11 0,03 47,79
f 3 hexane 4,67 M36S42.S5 14455008,20 1,»
A 4 Cyclopropane, etfiylii J,il 37M76.12 193557,50 t,Ot 12,
u 7 Butane, 2,2, i,3-litrjn \U 676366.64 1734513,61 C.L2 36,
11 ijMeihinetMQf 3fsa «437173,51. 4,52 94,65 34,97 ...........-■
Li — s Heptane № 211752,74 403354,74 O.DJ
id сзгЬлг disulfide Mi BISU.IH !ч'| III,' 4.- 613794,47 0.04 90,6
— U 13 p]methyl sulfide ь.ьг 600660.57 0,04 72,72 ........
Sljn-uniplriitliyHph 7.02 72121.33 173Н£5,6Д 0,01 10.47
Ш 13 Heptane, 2.4-Arwthi 7,61 439544.43 1753057,03 0.12 16,1
m 14 huran 7,42 li3Ui«i.Yj 947113S, 1С QjM м.сг -MS
m IS Atetont M2 2241232,44 9959975,47 0,69 »2,0 a
19 16 Cydotrlilloiane, hei< в,17 «775267.1? 35340455,53 2,44 95.26
20 17 12,14,14-Trimithy!-3r 9,S4 2E1S 72,43 0.04 4,75 i
11 IE F uran, 2iri#lhyl* 10,37 «03101,27 3594533,11 0J5 57,7 i
U 15 J,44>tJTlitlYv 1 -1 -h epte 10,55 284543,72 JC9J&11,44 0,01 27,94!
13 JO CMoram{1li4nesulfor 10,74 169J79,43 703711,32 0.A5 65,03: 4'
W 31 Furai, 3-metliyl- 11,44; НИН71 U321«S.60[ 0,14 50,53. ---------------.
15 11 Mutmani ц,ц] -. ц -.м ,'1 J 505728,72 0.W:
lb 23 euunal, 2-m«thyl- _ Ш»Г. = 478391,27 229167L22 0,1*1 4», 44. *
17 24 Morphinui-AJ-dloi-q 12,42] 140325, W 343042,93 0.02! 4.35
IS 25|№епгепе »jMelhylcne Chi pride Ц.1ШM ioaatjij. /1 D,06 65,93
JS 34530492' 1531414,36: 0,11 $3>0<j
3D H 1 Z7 2A6-fodoheDtatrJen > K SIJibhH f _13,35 euiS.85 566255 0.04 > 31Ы ГШМttllf ■—ГГВ ■ HI ИТ L 1
Рис. 2. Пример файла хромато-масс-спектра (матрица)
Формирование метаболических профилей больных онкопатологией
В процессе проведения исследования была сформирована группа из 94 пациентов с клинически подтвержденным раком легкого (РЛ) и контрольная группа (КГ) из 134 клинически здоровых людей.
Результатом каждого хромато-масс-спектро-метрического анализа является набор летучих метаболитов мочи пациента, представляющий собой таблицу с перечнем детектированных системой ГХ-МС химических соединений, с соответствующей информацией о времени выхода из колонки и достоверностью интерпретации состава соединений. Пик характеризовался следующими количественными параметрами: высота Ь пика и его площадь б. Из по-
лученного профиля исключались вещества, выделяющиеся с поверхности колонки и других частей хроматографа. Пример файла с результатами, полученными методом ГХ-МС приведен на рисунке 2.
Отбор данных для проведения диагностических исследований
При проведении исследований биологического материала на газовом хроматографе «Хроматэк-Кри-сталл 5000» зарегистрированы 114 летучих метаболитов. Из них для пациентов контрольной группы КГ - 46, и для пациентов группы РЛ - 68 метаболитов. Как уже было отмечено выше, каждый метаболит в обеих группах был количественно охарактеризован следующими параметрами: время выхода t , высота пика Ь,
площадь пика s и достоверность сравнения Рр с соединениями из банка эталонов. В качестве информативных параметров были приняты параметры h и s - соответственно высота и площадь пика на хромато-грамме по выбранному иону при условии достаточно высоких значений параметра Рпр (Рпр>0,8). Параметры h и s рассматриваются в качестве количественных концентраторов информации о наличии и отсутствии у пациента рака легких. Из большого количества полученных с помощью газового хроматографа летучих метаболитов большая часть была исключена из дальнейшего анализа по следующим причинам:
• отсутствия значительного количества значений параметров h и s; для пациентов обеих групп;
• малая частота встречаемости метаболита в группе (например, метаболит 3-methylfuran был получен только для семи пациентов КГ);
• низкое значение параметра Рр.
С учетом этих ограничений диагностические исследования далее проводились для следующих 16-ти метаболитов: methanethiol; hexanal; 2-ethyl-methylfuran 2,5dimethylfuran; 2,3,5-triethylfuran; dimethyldisulfide dimethyltrisulfide; acetone; 4-heptanone; 2-pentanone 2-methoxythiophene; 2h-1-benzopyran; 2-butanone; 2,6,10,10-tetramethyl-1-oxaspiro4.5dec-6-ene; 2-ethyl-1-hexanol; 1-bromo-2-methyl-2propanol.
Методика тест-пробы полученных выборок
В соответствии со сформулированной целью данной статьи, для быстрого определения из полученных на хроматографе профилей 16-ти метаболитов была выполнена проверка для каждого из них нулевой гипотезы. Нулевая гипотеза в рамках данного исследования для конкретного метаболита заключалась в следующем утверждении: выборки значений параметров h и s для групп КГ и РЛ относятся к одной и той же генеральной совокупности. Очевидно, если нулевая гипотеза имеет место, то эти выборки в плане вынесения диагноза являются нонинформативными. Следует отметить, что отрицание нулевой гипотезы для конкретного метаболита еще не является достаточным условием для принятия решения о наличии или отсутствия рака легкого у пациента. Нулевая гипотеза используется для оценки генеральных параметров случайной величины по их выборочным данным. Если в результате проверки эта гипотеза принимается, то это значит, что изменчивость генеральных параметров случайных величин вызвана второстепенными случайными причинами. Поэтому обе анализируемые
выборки относятся к одной и той же генеральной совокупности и, следовательно, пациенты групп КГ и РЛ либо оба здоровы, либо страдают патологией. Критерием оценки служит стандартная величина нормированного отклонения , с которой сравнивается фактическое значение 1ф этого критерия. Критерием достоверности различия между средними т1 и т2 выборок групп КГ и РЛ служит отношения разности средних т1 - т2 к своей статистической ошибке. Вычисление этого соотношения выполняется по формуле 1.
V
(П\ - 1)с12 + (n2 - п\ + п
(1)
n2 - 2
В этом выражении через п 2 и а12 обозначены объем выборки и ее дисперсия для группы КГ, а через п2 и а22 аналогичные параметры выборки для группы РЛ. Предполагая нормальный закон распределения вероятностей в этих выборках можно задать в рамках нулевой гипотезы вероятность того, что выборочная средняя т не отклонится от средней М генеральной совокупности больше, чем на 21, где 1=(т-М)/а. Заранее установленный процент случаев, когда такое отклонение будет иметь место, называемый уровнем значимости гипотезы, позволяет принять или отвергнуть нулевую гипотезу. Обычно в исследовательской практике приняты три уровня значимости гипотезы, равные соответственно 5%, 1% и 0,1%. Им соответствуют следующие величины к 1=1,96; 1=2,58 и 1=3,29. Если из соображений значимости проводимого эксперимента выбран конкретный уровень значимости то при
ф st
(2)
нулевая гипотеза принимается. Здесь 1ф - величина, вычисляемая по формуле (1). Следует отметить, что в целом ряде случаев диагностическое решение принимается не по одному, а по нескольким информативным признакам. При этом для отклонения нулевой гипотезы надо учитывать вклад каждого из таких признаков и для каждого из них устанавливать отдельный уровень значимости.
В таблице 1 приведены результаты тест-пробы 16-ти отобранных для диагностических исследований метаболитов. Для каждого из них вычислялись два значения параметра 1ф - одно значение для 1ф(Ь)
Таблица 1. Результаты тест-пробы 16-ти метаболитов
n Метаболит t* (s) h (s)
1 Methanethil 0,212 0,930
2 Hexanal 0,419 0,354
3 Furan,2-ethyl-5methyl 1,24 1,256
4 Furan,2,5dimethyl 0,153 0,179
5 Furan,2,3,5-triethyl 0,209 0,983
6 Disulfide dimethyl 2,69 3,175
7 Dimethyl-trisulfide 1,584 1,556
8 Acetone 1,096 0,557
9 4-Heptanone 0,532 0,588
10 2-Pentanone 0,617 0,22
11 2-methoxythiophene 0,313 0,77
12 2H-1-Benzopyran 0,483 0,22
13 2-Butanone 0,627 0,78
14 1-Oxaspiro4 1,313 1,2
15 1-Hexanol,2-ethyl 0,97 0,94
16 1-Bromo-2-methyl-2propanol 0,89 0,87
выборки высоты пиков, другое значение tф(s) - для выборки площади этих пиков. Как следует из данных этой таблицы перспективным метаболитом для диагностики рака легкого является Disulfide dimethyl со значениями t(h)=2,69 и tф(s)=3,175. Среди других метаболитов, позволяющих получить полезную информацию для принятия диагностического решения, отметим следующие: Furan, 2-ethyl-5-methyl (tj(h)=1,211; tys)=1,256), Dimethyltrisulfide (t$(h)=1,584; t^s)=1,556) и 1 -Oxaspiro4 (t/h)=1,313; t/s.)=1,2).
Необходимо отметить, что эти результаты имеют предварительное, ориентировочное значение, так как получены при условии нормального закона распределения вероятностей в исследуемых выборках. Более точное решение будет найдено ниже по результатам статистического анализа полученных выборок для метаболита Disulfide dimethyl.
Методика уточненного анализа выборочных значений метаболита Disulfide dimethyl на основе значений функции правдоподобия
Как видно из данных таблицы 1 для метаболита Disulfide dimethyl были получены значения параметров tф(h) и tф(s), с высокой вероятностью отвергающие нулевую гипотезу для выборок высоты и площади пиков данного метаболита на хроматограммах.
Поэтому он является достаточно перспективным для эффективной диагностики рака легких. Целью уточненного анализа будет получение оценок вероятности правильной диагностики Рпр по выборочным значениям параметров И и 5 данного метаболита для групп КГ и РЛ. При проведении анализа примем следующие условия: 1) выборочные значения взаимно независимы и 2) случайные величины в выборке И и 5 групп КГ и РЛ подчинены нормальному закону распределения вероятностей. Хотя справедливость этих предположений достаточно сомнительна, но это не приведет к получению окончательных ошибочных результатов, так как вероятности Рпр будут получены опытным путем использования объективных данных, полученных по метаболиту Р^ШАСе-сИте^у! методом ГХ-МС для обеих групп КГ и РЛ. Использование при расчетах Рпр приведенных выше предположений дает возможность использовать хорошо известный математический аппарат для нормального закона распределения, а полезность полученных значений Рпр далее будет подтверждена экспериментальным путем. В качестве критерия принятия решения об отнесении выборочного значения параметра И или 5 метаболита к здоровому или больному пациентам будет использован критерий максимального правдоподобия. Широкое применение данного критерия на
практике объясняется не только получением результатов сравнительно простыми расчетами, но и возможностью использования для этого минимального количества данных. Решение принимается на основании сравнения отношения правдоподобия Л с пороговым значением Л0 , т.е.
Л = L >Л0
(3)
А0=-|.
(4)
Параметры m и о в этих формулах для каждой из групп КГ или РЛ берутся из таблицы 2. Например, если нормируется отсчет 10,94 высоты пика метаболита Disulfide dimethyl, в предположении, что он относится к группе КГ, то
10,94 -8,9
5,44
■■ 0,375.
где и 12 - отсчеты функций правдоподобия, полученные при подстановке в выражение для плотности нормального распределения вероятностей выборочных значений. В целом ряде случаев пороговое значение Л0 принимается равным единице, т.е.
Воспользовавшись таблицей для функции (5) получаем отсчет функции правдоподобия 1(И/КГ)=0,372 .
Если же предполагается, что этот отсчет относится к группе РЛ, то
^ = ^^ = - 0,574.
h/РЛ
9,36
Отсчет функции правдоподобия для этого случая имеет вид 1(Ь/РЛ)=0,338. Отношение правдоподобия для рассмотренного примера равно
При использовании выражений (3) и (4) предварительно должны быть найдены параметры распределения законов вероятностей для выборочных значений. В нашем случае был принят нормальный закон распределения плотности вероятностей, задаваемый параметрами т и о, найденных по экспериментально полученным выборкам. В таблице 2 приведены значения этих параметров.
Для нахождения отсчета функции правдоподобия I используется выражение для плотности распределения вероятностей нормально распределенной случайной величины
л = ¡-(h/КГ) КГ/РЛ -(h/РЛ)'
0,372 0,338
= 1,101.
С учетом выражений (3) и (4) принимается решение о том, что данный отсчет высоты пика, относящийся к анализируемому метаболиту, принадлежит здоровому пациенту. Если взять отсчет площади s метаболита Disulfide dimethyl, равный 130, то аналогичный расчет приводит к следующим результатам:
130-394 = 13039,4 = 3,44;
x,
130-83,3
w 26,3 ' ' */РЛ 46,5
-{s/КГ) = 1,001; -(s/РЛ) = 2,42;
1,003;
f (x)--
V2nexp I 2
(5)
где х - нормированное значение выборочной случайной величины, получаемое следующим образом:
И - т s - т
а
а
Таблица 2. Параметры, выбранные для рабочей (РЛ) и контрольной (КГ) групп
Группа Параметры пика
Высота h Площадь s
КГ m = 8,9 m = 39,44
о = 5,44 о = 26,3
РЛ m = 16,31 m = 83,3
о = 9,36 о = 46,5
Л ,
-(s/КГ)
0,001
■■ 0,0005<1.
-(s/РЛ) 2,42
Поскольку отношение отсчетов функции правдоподобия -(s/КГ) и -(s/РЛ) меньше единицы, то принимается решение о том, что пациент с метаболитом Disulfide dimethyl, у которого площадь пика равна 130, болен. Аналогичным образом были проанализированы все данные по метаболиту Disulfide dimethyl, полученные для пациентов групп КГ и РЛ. Итоговые результаты приведены в таблицах 3-6.
Данные, приведенные в таблицах 3-6, дают возможность оценить условные вероятности правильных решений Pnp при использовании только одного информативного признака, высоту h или площадь s на выходе хроматографа для метаболита Disulfide dimethyl (табл. 7). В качестве оценки вероятности P берется значение относительной частоты пра-
h/КГ
1
Таблица 3. Результаты автоматической диагностики пациентов контрольной группы по высоте И пика
Номер пациента 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
h 10,94 3,75 4,01 10,03 17,03 10,55 6,83 13,42 18,25 6,8 3,17 13,28 1,95 9,4 1,86
Л КГ/ РЛ 1,1 1,52 1,58 1,18 0,44 1,15 1,55 0,49 0,24 1,55 1,54 0,74 1,43 1,31 0,78
Решение прав./ош. прав прав прав прав ош прав прав ош ош прав прав ош прав прав ош
Таблица 4. Результаты автоматической диагностики пациентов контрольной группы по площади б пика
Номер пациента 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
s 47,1 12,6 17,1 43,8 48,5 44,1 31,55 65,93 104,4 30,38 14,23 64,19 9,3 50,4 8,04
Л КГ/ РЛ 1,29 1,92 1,9 1,4 1,25 1,4 1,76 0,64 0,05 1,79 1,89 0,7 1,82 1,18 1,86
Решение прав./ош. прав прав прав прав прав прав прав ош ош прав прав ош прав прав прав
Таблица 5. Результаты автоматической диагностики пациентов рабочей группы по высоте И пика
Номер пациента 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
h 24 24,6 33 5,28 6,88 22 10,24 15 1,38 33,18 17,27 11,4 8,8 3,82 11,75
Л РЛ/ КГ 33,2 44,3 82,0 0,63 0,64 15,1 0,83 1,85 1,76 79,0 1,98 0,97 0,74 0,64 2,4
Решение прав./ош. прав прав прав ош ош прав ош прав прав прав прав ош ош ош прав
Таблица 6. Результаты автоматической диагностики пациентов рабочей группы по площади б пика
Номер пациента 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
s 130 117 157 22,8 28,8 91,16 50,91 74,63 91,31 173,7 99,57 70,3 53,31 21,01 65,94
Л РЛ/ КГ 220 60,2 112 0,52 0,55 76 0,81 2,42 7,3 610 13,4 1,91 0,95 0,67 1,54
Решение прав./ош. прав прав прав ош ош прав ош прав прав прав прав прав ош ош прав
Таблица 7. Оценки вероятностей правильных решений по величине параметров h и s метаболита Disulfide dimethyl
вильных и ошибочных решений, приведенных в последних строках таблиц 3-6.
Получим величины оценок вероятностей Рпр для групп КГ и РЛ при совместном влиянии параметров И и 5 на принятие диагностического решения. Для этого найдем значения функций правдоподобия при одновременном учете этих параметров и подсчитаем значения компонент относительной частоты правильных и ошибочных решений для каждой из групп. Результа-
ты расчетов дают следующие значения оценок вероятностей Р :
пр „
для группы КГ Рп =0,773; а для группы
рл Р =0,666. "Р
пр
Заключение
В процессе проведения исследования летучих метаболитов мочи клинически здоровых пациентов (контрольная группа КГ) и пациентов с клинически установленным диагнозом рака легких (рабочая группа РЛ) установлено статистическое различие параметров метаболита □¡БиЖСе-сПте^у! для этих групп. Для данного метаболита не подтвердилась нулевая гипотеза в отношении выборочных значений высот и площадей пиков на хроматограмме соответствующих
Группа По параметру h По параметру s
КГ 0,667 0,8
РЛ 0,6 0,667
параметров метаболита Disulfide-dimethyl. Дальнейший анализ выборочных данных на основе аппарата функций правдоподобия позволил оценить вероятности правильной постановки диагноза дл я группы КГ (пациент здоров), равную 0,666 и для РЛ группы (пациент болен), равную 0,733. Данные оценки постановки правильного диагноза охватывают достаточно высокий 70-ти процентный уровень. Они представляют практический интерес и стимулируют проведение дальнейших исследований в поисках других метаболитов, позволяющих увеличить этот уровень до 95%.
Литература
1. Фурина Р.Р., Рыжков В.Л., Митракова Н.Н., и др. Метод ранней диагностики рака желудка на основе матаболомических исследований // Экспериментальная и клиническая гастроэнтерология. -2014. - вып. 110, №10. - С. 14-17.
2. Фурина Р.Р., Митракова Н.Н., Рыжков В.Л., Сафиул-лин И.К. Метаболомические исследования в медицине // Казанский медицинский журнал. - 2014. -XCV (1). - С. 1-6.
3. Царев Н.И., Царев В.И., Катраков И.Б. Практическая газовая хроматография. - Барнаул: Изд-во Алт. Унта, 2000. - 156 с.
4. Kouremenos К.А., Pitt J., Marriott P.J. Metabolic profiling of infant urine using comprehensive two-dimensional gas chromatography: Application to the diagnosis of organic acidurias and biomarker discovery // Journal of Chromatography A. - 2010. -1217. - P. 104-111.
5. Oliveira P.A., Colaco A., Chaves H.R. et al. Chemical carcinogenesis // An Acad. Bras. Cienc. - 2007. -79. - P. 593-616.
6. Patti G.J., Yanes O., Siuzdak G. Innovation: Metabolomics: the apogee of the omics trilogy // Nat. Rev. Mol. Cell. Biol. - 2012. - 13. - P. 263-269.
7. Phillips M., Cataneo R.N., Cummin A.R.C. et al. Detection of lung cancer with volatile markers in the breath // Chest. - 2003. - 123 (6). - P. 2115-2123.
8. Phillips M., Gleeson K., Hughes J.M., et al. Volatile organic compounds in breath as markers of lung cancer: a crosssectional study // Lancet. - 1999. -353. - P. 1930-1933.
9. Pauling L., Robinson A.B., Teranishi R., Cary P. Quantitative analysis of urine vapor and breath by gas-liquid partition chromatography // Proc. Natl. Acad. Sci. USA. - 1971. - 68. - P. 2374-2376.
10. Silva C.L., Passos M., Camara J.S. Investigation of urinary volatile organic metabolites as potential cancer biomarkers by solid-phase microextraction in combination with gas chromatography-mass spectrometry // British Journal of Cancer. - 2011. -105. - P. 1894-1904.
11. Zimmermann D., Hartmann M., Moyer M.P. et al. Determination of volatile products of human colon cell line metabolism by GC/MS analysis // Metabolomics. - 2007. - 31. - P. 13-17.
12. Woo H.M., Kim K.M., Choi M.H. et al. Mass spectrometry based metabolomic approaches in urinary biomarker study of women's cancers // Clin. Chim. Acta. -2009. - 400. - P. 63-69.
