ОРИГИНАЛЬНАЯ СТАТЬЯ ..........ФУНДАМЕНТАЛЬНАЯ
4»llHni<ILIIin<IUIinn ТПШ1 1 МО "5
И КЛИНИЧЕСКАЯ МЕДИЦИНА lum '» |ч- л
КОМПЬЮТЕРИЗИРОВАННАЯ СИСТЕМА ОЦЕНКИ РЕЗУЛЬТАТОВ КТ-ОСТЕОДЕНСИТОМЕТРИИ ПОЗВОНКОВ
КОЛПИНСКИЙ Г.И.1, ЗАХАРОВ И.С.1, УШАКОВА Г.А.1, УШАКОВ А.В.2, ИВАНОВ В.И.3, КАГАН Е.С.3
ФГБОУ ВО «Кемеровский государственный медицинский университет» Минздрава России
2ГАУЗ КО «Кемеровская областная клиническая больница им. С.В. Беляева» 3ФГБОУ ВО «Кемеровский государственный университет»
ORIGINAL ARTICLE
A SOFTWARE FOR EVALUATION OF LUMBAR SPINE COMPUTED TOMOGRAPHY DENSITOMETRY
GLEB I. KOLPINSKIY1, IGOR S. ZAKHAROV1, GALINA A. USHAKOVA1, ALEXANDER V. USHAKOV2, VADIM I. IVANOV3, ELENA S. KAGAN3
1Kemerovo State Medical University (22a, Voroshilova Street, Kemerovo, 650056), Kemerovo, Russian Federation
2Kemerovo Regional Clinical Hospital (22a, Oktyabr'skiy Prospekt, Kemerovo, 650066), Kemerovo, Russian Federation
3Kemerovo State University (6, Krasnaya Street, Kemerovo, 650043), Kemerovo, Russian Federation
Резюме
Цель. Разработка компьютеризированной системы оценки результатов КТ-остеоденси-тометрии позвонков, позволяющей проводить комплексное прогнозирование риска остеопо-ротических переломов у женщин в постмено-паузальном периоде.
Материалы и методы. Созданию прогностической системы предшествовало формирование модели, способной на основании результатов количественной компьютерной томографии осуществлять оценку вероятности возникновения риска остеопоротических переломов. При разработке модели 282 женщинам постменопаузально-го периода было выполнено изучение минеральной плотности кости поясничных позвонков.
Результаты. Чувствительность модели составила 77,8 %, специфичность - 86,7 %. Площадь под ROC-кривой (АиС) - 0,894 (0,8550,932). Для удобства использования модели в практической работе врача была написана прикладная программа для ЭВМ.
Заключение. Созданная компьютерная программа будет способствовать оптимизации диагностики остеопороза и прогнозирования остеопоротических переломов позвонков у женщин в по-стменопаузальном периоде, что в итоге позволит осуществлять своевременную профилактику и коррекцию указанных инволюционных нарушений.
Ключевые слова: постменопаузальный остеопороз, количественная компьютерная томография, прогностические модели.
English ► Abstract
Aim: To develop a software for evaluation of lumbar spine computed tomography densitometry to predict the risk of osteoporotic fractures in postmenopausal women.
Materials and Methods: We evaluated the lumbar spine computed tomography densitometry data in 282 postmenopausal women following a computerized modeling.
VOL. 1, № 3
FUNDAMENTAL
AND CLINICAL MEDICINE
ОРИГИНАЛЬНАЯ СТАТЬЯ
Results: The sensitivity and specificity of the model were 77.8% and 86.7%, respectively, with the area under the ROC curve of 0.894 (95% confidence interval 0.855-0.932). We have also developed an original software for prediction of osteoporotic fractures in postmenopausal women.
Conclusion: Our software can assist the clinicians in predicting the risk of osteoporotic fractures in postmenopausal women.
Keywords: postmenopausal osteoporosis, computed tomography, predictive models.
Введение
В постменопаузальном периоде у женщины происходят существенные метаболические изменения, затрагивающие различные системы организма [1]. Одним из нарушений, обусловленных инволюционными процессами, является постменопаузальный остеопороз, распространённость которого неуклонно увеличивается с возрастом. Ведущее место в диагностике данной патологии занимают методы, основанные на оценке уровня минеральной плотности кости (МПК) - двухэнергетическая рентгеновская абсорбциометрия и количественная компьютерная томография [2]. Считается, что показатели МПК на 75-80 % определяют прочность кости [3].
Остеоденситометрия, проводимая методом количественной компьютерной томографии, имеет ряд преимуществ, в частности возможность раздельного исследования трабекуляр-ной и кортикальной костной ткани [4].
Характерными локализациями остеопоро-тических переломов являются позвонки, кости предплечья и шейка бедренной кости [5]. По данным литературы, распространённость осте-опоротических переломов позвонков составляет от 2,9-25,3 % [6], повторение которых способствует снижению качества жизни и росту инвалидизации пациенток.
Существуют различные способы оценки риска переломов. Широко распространён метод изучения вероятности возникновения переломов в ближайшие 10 лет ^ЯАХ). Указанная система производит оценку факторов риска обследуемого: возраста, пола, индекса массы тела, наличия переломов в анамнезе и у ближайших родственников, курения, приёма глюкокортико-идов, наличия ревматоидного артрита, прочих причин вторичного остеопороза, злоупотребления алкоголем и уровня минеральной плотности шейки бедренной кости [7]. В то же время, некоторые авторы продемонстрировали невысокую чувствительность данного способа. При ретроспективном десятилетнем исследовании чувствительность российской версии модели
FRAX составила 42 %, а специфичность - 74 % [8]. По данным F. A. Tremollieres (2010) [9], площадь под ROC-кривой модели FRAX равна 0,63 (0,56-0,69), что свидетельствует о среднем качестве модели. Учитывая представленные данные литературы, имеется необходимость дальнейшего поиска оптимальных моделей прогнозирования остеопоротических переломов, обладающих достаточно высокой чувствительностью и специфичностью. Это, в свою очередь, позволит проводить своевременную профилактику указанных осложнений.
В настоящее время, с целью оптимизации работы практического врача, всё больше находит применение использование прогностических систем, основанных на современных компьютерных технологиях.
Цель исследования
Разработка компьютеризированной системы оценки результатов КТ-остеоденситометрии позвонков, позволяющей проводить комплексное прогнозирование риска остеопоротических переломов у женщин в постменопаузальном периоде.
Материалы и методы
Созданию компьютеризированной прогностической системы предшествовало формирование модели, способной, на основании результатов количественной компьютерной томографии (ККТ) выполнять оценку вероятности возникновения риска остеопоротических переломов [10]. Для разработки модели использовался метод бинарной логистической регрессии с применением программы IBM SPSS Statistics 21.
Основой прогностической модели явились результаты остеоденситометрии, провёденной 282 женщинам постменопаузального периода методом количественной компьютерной томографии. Для создания модели были отобраны пациентки с остеопоротическими компрессионными переломами XII грудного - IV поясничного позвонков (n = 72). Факт наличия компрес-
ОРИГИНАЛЬНАЯ СТАТЬЯ
ФУНДАМЕНТАЛЬНАЯ И КЛИНИЧЕСКАЯ МЕДИЦИНА
ТОМ 1, № 3
Таблица 1. Коэффициенты регрессии, рассчитанные при построении прогностической модели
Table 1. Regression coefficients used for calculation of the predictive model
сионного перелома позвонка устанавливали на основании полуколичественного метода (H. K. Genant, 1993) [11]. Носителями противопоставленного признака явились 210 женщин постме-нопаузального периода, не имеющих компрессионных переломов позвонков. Были выбраны показатели минеральной плотности трабеку-лярной, кортикальной костной ткани и индексы билатеральной асимметрии МПК II-IV поясничных позвонков. Вычислялись коэффициенты регрессии. С применением ROC-анали-за выносилось суждение об уровне качества созданной модели, при этом использовалось значение величины площади под ROC-кривой (AUC, area under the curve).
При моделировании применялось регрессивное уравнение:
y=a+b xX+bxX+-b.xX.
J 112 2 ii
где y - зависимая переменная, принимающая два значения: 0 - нет перелома, 1 - есть пере-
лом; а - константа; Ы - коэффициенты регрессии; Хi - переменные.
Рассчитывалась вероятность возникновения переломов:
(1+е-у)
где Р - прогностическая вероятность, е -экспонента, приближенное значение которой равно 2,718.
Валидность модели оценивалась на основании процента верно переклассифицированных случаев и критерия Somers'D. Проверка общей согласованности прогностической модели с реальными данными осуществлена по критерию согласия Хосмера-Лемешова. Граница критического уровня значимости р соответствовала 0,05.
Результаты и их обсуждение
При разработке прогностической модели получены следующие показатели (таблица 1).
Показатели бинарной Показатели КТ-остеоденситометрии
логистической регрессии МПК траб. (Xi) МПК корт. (X2) ИА МПК траб. (Хз) ИА МПК корт. (Х4) Константа
B, коэффициент регрессии -0,044 -0,014 3,443 2,395 -2,551
Стандартная ошибка 0,020 0,007 1,361 1,187 3,337
Примечание:
МПК траб. и МПК корт. - минеральная плотность трабекулярной и кортикальной костной ткани;
ИА МПК траб. и ИА МПК корт. - индексы билатеральной асимметрии минеральной плотности трабекулярной и кортикальной костной ткани поясничных позвонков.
Indices are for trabecular and cortical bone tissue mineral density;
Indices are for bilateral asymmetry indices of trabecular and cortical bone tissue mineral density
Формула прогностической вероятности возникновения переломов позвонков имеет следующий вид:
1 + е(2,555-0,044хХ1-0,014хХ2+3,443хХ3+2,395хХ,1)
При проведении ЯОС—-анализа площадь под ЯОС-кривой (АиС) созданной модели составила 0,894 (0,855-0,932), что свидетельствует о высокой прогностической способности. Чувствительность модели равна 77,8 %, специфичность - 86,7 %. Процент верной переклассификации (конкордантности) составил 84,4 %. Коэффициент Somers'D равен 0,778 (р
= 0,001). Изучение общей согласованности прогностической модели с реальными данными выполнялось на основании критерия согласия Хосмера-Лемешова (х2 = 0,61, р = 0,644). Приведённые показатели демонстрируют высокую валидность модели. Получено решение Роспатента о выдаче патента на изобретение «Способ прогнозирования риска остеопоротических переломов позвонков у женщин постменопа-узального периода» от 15.09.2016; заявка № 2015116979/14 (026400) от 05.05.2015.
После формирования модели были рассчитаны диапазоны качественной оценки прогностической вероятности возникновения переломов:
VOL. 1, № 3
FUNDAMENTAL
AND CLINICAL MEDICINE
ОРИГИНАЛЬНАЯ СТАТЬЯ
Р выше 0,5 свидетельствует о высоком риске переломов; прогностическая вероятность, находящаяся в пределах от 0,5 до 0,371 говорит о среднем риске, ниже 0,371 - о низком.
Обсуждение
Созданная модель прогнозирования позволила сформировать алгоритм определения вероятности возникновения остеопоротических переломов позвонков у женщин в постменопаузаль-ном периоде на основе КТ-остеоденситометрии.
Алгоритм включает несколько этапов. Первый этап: «на входе» системы происходит определение показателей минеральной плотности трабе-кулярной и кортикальной костной ткани и индексов билатеральной асимметрии МПК поясничных позвонков методом количественной компьютерной томографии. Единицы измерения мг/см3 или mgCa-Ha/ml. Второй этап - обработка информации, с использованием формулы прогностической вероятности. Третий этап: «на выходе» получаются значения прогностической вероятности возникновения переломов позвонков. Формируется заключение о степени риска переломов.
Для удобства использования прогностической модели в практической деятельности врача была написана прикладная компьютерная программа «Прогнозирование риска переломов позвонков» (номер свидетельства о государственной регистрации программы для ЭВМ: 2015618028; опубликовано 20.08.2015), которая может быть использована совместно с основными операционными системами: Microsoft Windows 9x/ NT/2000/Vista, 7, 8. Объём программы составляет 697 Кб. Алгоритм программы также состоит из трёх этапов: «на входе» вводится фамилия, имя, отчество, возраст обследуемой, затем - абсолютные показатели минеральной плотности кости трабекулярной и кортикальной костной ткани и индексов билатеральной асимметрии МПК поясничных позвонков (рисунок 1).
После компьютерной обработки, «на выходе» выводятся показатели прогностической вероятности. Даётся заключение о риске остеопоротических переломов позвонков (рисунок 2).
Заключение
1. Разработанная прогностическая модель позволяет проводить оценку результатов КТ-о-стеоденситометрии поясничных позвонков с целью прогнозирования риска остеопоротиче-ских переломов у женщин в постменопаузаль-ном периоде.
Рисунок 1. Внешний вид программы «Прогнозирование риска переломов позвонков»
Figure 1. Our original software for prognosis of vertebral fractures
Рисунок 2. Результаты обработки показателей КТ-осте-оденситометрии, программой «Прогнозирование риска переломов позвонков»
Figure 2. An output obtained by our software
2. Созданная прогностическая модель обладает достаточно высокой чувствительностью, специфичностью и адекватной валидностью.
3. Компьютерная программа «Прогнозирование риска переломов позвонков» будет способствовать повышению качества комплексного прогнозирования остеопоротических переломов позвонков.
Таким образом, созданная компьютерная программа будет способствовать оптимизации диагностики остеопороза и прогнозирования остеопоротических переломов позвонков у женщин в постменопаузальном периоде, что в итоге позволит осуществлять своевременную профилактику и коррекцию указанных инволюционных нарушений.
ОРИГИНАЛЬНАЯ СТАТЬЯ
ФУНДАМЕНТАЛЬНАЯ И КЛИНИЧЕСКАЯ МЕДИЦИНА
ТОМ 1, № 3
Литература / References:
1. Medical menopause / edited by Smetnik VP. Yaroslavl: OOO «Izdatel'stvo Litern», 2006. 848 p. Russian (Медицина климактерия / под. ред. В.П. Сметник. Ярославль: ООО «Издательство Литера», 2006. 848 с.)
2. Zakharov IS. Features of bone densitometry in postmenopausal women. Nuclear Medicine and Radiation Safety. 2015; 60 (2): 56-59. Russian (Захаров И.С. Особенности костной денситометрии у женщин в постменопаузальном периоде // Медицинская радиология и радиационная безопасность. 2015. Т.60, №2. С. 56-59.)
3. Sveshnikov AA, Astapenkov DS. Threshold values of skeletal mineral density and fracture incidence. Genius of оrthopaedics. 2010; (2): 85-90. Russian (Свешников А.А., Астапенков Д.С. Пороговые величины минеральной плотности костей скелета и частота переломов // Гений ортопедии. 2010. № 2. С. 85-90.)
4. Kolpinsky GI, Zaharov IS. Diagnosis and prognosis of postmenopausal osteoporosis. Kemerovo: "Int", 2015. 202 p. Russian (Колпинский Г.И., Захаров И.С. Диагностика и прогнозирование постменопаузального остеопороза. Кемерово: «Инт», 2015. 202 с.)
5. Averkieva YV, Raskina TA. The prevalence of osteoporotic fractures among residents of the city of Kemerovo in older age groups. Bulletin Kuzbass scientific center. 2010; (11): 14. Russian (Аверкиева Ю.В., Раскина Т.А. Распространённость остеопороти-ческих переломов среди жителей г. Кемерово в старших возрастных группах // Вестник Кузбасского научного центра. 2010. № 11. С. 14.)
6. Benevolenskaja LI. Guidelines for osteoporosis. М.: BINOM, 2003. 524 p. Russian (Беневоленская Л.И. Руководство по осте-опорозу. М.: БИНОМ, 2003. 524 c.).
7. Kanis JA, Oden A, Johansson H., Borgstrom F, Ström O, McCloskey EV. FRAX®, a new tool for assessing fracture risk: clinical applications and intervention thresholds. Medicographia. 2010; 32 (1): 33-40.
8. Nikitinskaya OA, Toroptsova NV. Assessment of fractures risk using the FRAX tool (A ten-year retrospective study). Almanac of Clinical Medicine. 2014; (32): 50-55. Russian (Никитинская О.А., Торопцова Н.В. Оценка риска переломов с использованием модели FRAX (ретроспективное десятилетнее исследование) // Альманах клинической медицины. 2014. № 32. С. 50-55.)
9. Tremollieres FA, Pouilles JM, Drewniak N, Laparra J, Ribot CA, Dargent-Molina P. Fracture risk prediction using BMD and clinical risk factors in early postmenopausal women: sensitivity of the WHO FRAX tool. J. Bone Miner. Res. 2010; 25 (5): 1002-1009.
10. Zakharov IS, Kolpinsky GI, Ushakova GA, Kagan ES. Model prediction of risk of osteoporotic vertebral fractures in women using quantitative computed tomography. Radiology - practice. 2015; (4): 19-27. (Захаров И.С., Колпинский Г.И., Ушакова Г.А., Каган Е.С. Модель прогнозирования риска остеопоротических переломов позвонков у женщин с использованием количественной компьютерной томографии // Радиология - практика. 2015. №4. С. 19-27.)
11. Genant HK, Wu CY, van Kuijk C, Nevitt MC. Vertebral fracture assessment using a semiquantitative technique. J. Bone Min. Res. 1993; 8 (9): 1137-1148.
Сведения об авторах
Колпинский Глеб Иванович, доктор медицинских наук, профессор кафедры лучевой диагностики, лучевой терапии и онкологии, ФГБОУ ВО «Кемеровский государственный медицинский университет» Минздрава России
Захаров Игорь Сергеевич, кандидат медицинских наук, доцент кафедры акушерства и гинекологии № 1, ФГБОУ ВО «Кемеровский государственный медицинский университет» Минздрава России Ушакова Галина Александровна, профессор, доктор медицинских наук, заведующая кафедрой акушерства и гинекологии № 1, ФГБОУ ВО «Кемеровский государственный медицинский университет» Минздрава России
Ушаков Александр Владимирович, кандидат медицинских наук, заведующий рентгенологическим отделением № 2 ГАУЗ КО «Кемеровская областная клиническая больница им. С. В. Беляева» Иванов Вадим Иванович, кандидат биологических наук, доцент кафедры физиологии человека и безопасности жизнедеятельности ФГБОУ ВО «Кемеровский государственный университет» Каган Елена Сергеевна, кандидат технических наук, доцент кафедры прикладной математики ФГБОУ ВО «Кемеровский государственный университет»
Корреспонденцию адресовать:
Захаров Игорь Сергеевич Адрес: 650029, г. Кемерово, ул. Ворошилова, 22а. Тел.: 8 (3842) 46-51-62. E-mail: [email protected]
VOL 1 № 3 FUNDAMENTAL
' AND CLINICAL MEDICINE
ОРИГИНАЛЬНАЯ СТАТЬЯ
Authors
Prof. Gleb I. Kolpinskiy, MD, PhD, Professor, Department of Radiation Diagnosis, Radiotherapy and Oncology, Kemerovo State Medical University, Kemerovo, Russian Federation
Dr. Igor S. Zakharov, MD, PhD, Associate Professor, Department of Obstetrics and Gynecology № 1, Kemerovo State Medical University, Kemerovo, Russian Federation
Prof. Galina A. Ushakova, MD, PhD, Head of the Department of Obstetrics and Gynecology № 1, Kemerovo State Medical University, Kemerovo, Russian Federation
Dr. Aleksandr V. Ushakov, MD, PhD, Head of the Radiology Unit №2, Kemerovo Regional Clinical Hospital, Kemerovo, Russian Federation Dr. Vadim I. Ivanov, MD, PhD, Associate Professor, Department of Human Physiology and Life Safety, Kemerovo State University, Kemerovo, Russian Federation
Dr. Elena S. Kagan, PhD, Associate Professor, Department of Applied Mathematics, Kemerovo State University, Kemerovo, Russian Federation
Corresponding author:
Dr. Igor S. Zakharov,
Voroshilova Street 22a, Kemerovo, 650029, Russian Federation E-mail: [email protected]
Acknowledgements: There was no funding for this article.