В данное ретроспективное одномоментное ассоциативное исследование были включены последовательно поступавшие пациенты, госпитализированные в отделение неврологии во 2-ю народную больницу города Чанчжоу с жалобами на головокружение, головную боль, для проведения скрининга на инсульт или когнитивной оценки в период с января 2022 года по декабрь 2024 года. Информированное согласие не запрашивалось в связи с ретроспективным, неинтервенционным характером исследования и использованием полностью анонимизированных данных. Исследование было одобрено Комитетом по этике клинических медицинских технологий 2-й народной больницы города Чанчжоу (идентификатор одобрения IRB: [2023] YLJSA069) и проводилось в соответствии с внутренними правилами учреждения.
Дизайн исследования и популяция пациентов
Первоначально было обследовано 850 пациентов. Критерии исключения были строго определены следующим образом: (1) отсутствие полных серий МРТ головного мозга (в частности, отсутствие взвешенного по магнитной восприимчивости изображения [SWI], что делает невозможной точную оценку микрокровоизлияний); (2) отсутствие подробных записей зондирования периодонта или данных о количестве зубов во всей полости рта; (3) неполные данные общего анализа крови, препятствующие расчету SII; (4) наличие в анамнезе обширного инсульта, опухоли головного мозга, черепно-мозговой травмы или инфекции центральной нервной системы; и (5) наличие активной острой инфекции, тяжелого аутоиммунного заболевания или злокачественного новообразования в течение последнего месяца для минимизации влияния смешивающих факторов на системные маркеры воспаления. Для обеспечения целостности данных использовался метод анализа полных случаев, при котором исключались лица с отсутствием любого из основных клинико-радиологических параметров. После скрининга в анализ было включено 234 подходящих участника. Эта выборка, включающая 129 событий с высокой нагрузкой, соответствовала эмпирическому правилу количества событий на переменную (EPV) >10 для стабильного многофакторного моделирования.
Оценка состояния периодонта и воспаления
Двое калиброванных пародонтологов провели все клинические обследования полости рта с использованием стандартизированного пародонтологического зонда.
Перед началом исследования оба эксперта прошли стандартизирующее обучение; межисследовательская надежность была отличной (κ Коэна = 0.82). Глубина зондирования кармана (PPD) и уровень клинической прикрепленности (CAL) измерялись в шести точках каждого зуба. На основании классификации 2018 года, описанной Tonetti et al.15, заболевания пародонта были разделены на легкую/отсутствующую, умеренную или тяжелую формы. В частности, участники были распределены по следующим категориям: (1) легкий/отсутствующий пародонтит (включая здоровье пародонта и пародонтит I стадии), определяемый как интердентальный уровень клинической прикрепленности (CAL) в месте наибольшей потери прикрепленности ≤2 mm и глубина зондирования кармана (PPD) ≤4 mm при отсутствии потери зубов, связанной с пародонтитом; (2) умеренный пародонтит (пародонтит II стадии), определяемый как интердентальный CAL 3–4 mm, максимальный PPD ≤5 mm и потеря не более четырех зубов вследствие пародонтита; и (3) тяжелый пародонтит (пародонтит III/IV стадии), определяемый как интердентальный CAL ≥5 mm, PPD ≥6 mm и/или потеря четырех и более зубов вследствие деструкции пародонта15.
Для каждого участника также было задокументировано количество сохранившихся зубов. Поскольку сохранение минимум 20 функциональных зубов является общепризнанным клиническим критерием сохранения базовых функций полости рта и успешного старения ротовой полости16,17,18, общее количество зубов было преобразовано в бинарную переменную: выраженная потеря зубов (<20 зубов) в противовес отсутствию выраженной потери зубов (≥20 зубов).
Что касается системного воспаления, показатель SII был рассчитан на основе образцов крови, взятых натощак при поступлении. Образцы венозной крови собирали в пробирки с ЭДТА и обрабатывали в соответствии со стандартным протоколом клинической лаборатории больницы перед анализом с использованием автоматического гематологического анализатора для определения количества нейтрофилов, лимфоцитов и тромбоцитов. Индекс рассчитывали следующим образом:

Для учета сильного перекоса распределения значений SII перед включением в статистический анализ было применено натуральное логарифмическое преобразование (Log_SII).
Сбор данных МРТ и общая нагрузка CSVD
Всем испытуемым была проведена стандартизированная МРТ головного мозга на томографе 3.0 T, оснащенном стандартной катушкой для головы. Протокол визуализации включал Т1-взвешенные, Т2-взвешенные изображения, режим инверсии-восстановления с подавлением сигнала от свободной жидкости (FLAIR) и магнитно-резонансную томографию, взвешенную по магнитной восприимчивости (SWI); толщина среза составляла 5 mm. Стандартные параметры визуализации были следующими: Т1-взвешенные изображения (время повторения [TR] = 2,000 ms, время эха [TE] = 9 ms, поле обзора [FOV] = 230 × 230 mm2, размер матрицы = 256 × 256); Т2-взвешенные изображения (TR = 4,500 ms, TE = 85 ms, FOV = 230 × 230 mm2, размер матрицы = 256 × 256); изображения в режиме инверсии-восстановления с подавлением сигнала от свободной жидкости (FLAIR) (TR = 8,500 ms, TE = 120 ms, время инверсии [TI] = 2,400 ms, FOV = 230 × 230 mm2, размер матрицы = 256 × 256); и изображения, взвешенные по магнитной восприимчивости (SWI) (TR = 28 ms, TE = 20 ms, угол отклонения = 15°, FOV = 230 × 230 mm2, размер матрицы = 256 × 256). Все последовательности были получены с толщиной среза 5 mm и межсрезовым интервалом 1.0 mm. Модель сканера указана в таблице материалов. Маркеры визуализации болезни мелких сосудов мозга (CSVD) независимо оценивались двумя врачами-рентгенологами, которые не имели доступа к клиническим данным участников. Межэкспертная надежность общей оценки степени выраженности CSVD была высокой (κ Коэна = 0.85). Любые расхождения разрешались путем консенсуса с третьим, старшим врачом.
В соответствии с установленным международным консенсусом1,3, совокупная степень выраженности церебральных заболеваний малых сосудов (CSVD) оценивалась по шкале от 0 до 4 на основе анализа четырех нейровизуализационных признаков. Один балл присваивался за каждый из следующих результатов МРТ: (1) наличие как минимум одного лакунарного инфаркта; (2) один или более церебральных микрокровоизлияний (CMBs); (3) расширение периваскулярных пространств (EPVS) в базальных ганглиях умеренной или тяжелой степени (степень ≥2); и (4) выраженные гиперинтенсивности белого вещества (WMH), определяемые как балл по шкале Fazekas ≥2 в глубоком белом веществе или 3 в перивентрикулярной области. Диагностические определения строго соответствовали международному консенсусу «Стандарты отчетности об изменениях сосудов при нейровизуализации» (STRIVE-1)1. Гиперинтенсивности белого вещества (WMH) оценивались по шкале Fazekas19, а расширение периваскулярных пространств (EPVS) в базальных ганглиях — с использованием валидированной 4-балльной визуальной оценочной шкалы, описанной Potter et al.20. Для разработки модели оценки рисков пациенты с общим баллом CSVD ≥2 классифицировались как лица с высокой степенью выраженности CSVD, поскольку данный порог последовательно ассоциируется с ускоренным когнитивным снижением и смертностью21,22,23.
Сбор данных о ковариатах
Демографические характеристики пациентов и клинический анамнез были извлечены из институциональных электронных медицинских карт. Собранные переменные включали возраст, пол, индекс массы тела (ИМТ), статус курения и потребление алкоголя. Кроме того, для каждого участника были задокументированы кардиометаболические и сосудистые сопутствующие заболевания, включая дислипидемию, сахарный диабет, гипертензию, ишемическую болезнь сердца (ИБС), перенесенный инфаркт миокарда (ИМ) и перенесенный ишемический инсульт.
Все сопутствующие заболевания были подтверждены путем анализа диагнозов врачей, зафиксированных в историях болезни стационарных пациентов, истории приема препаратов, а также результатов рутинных лабораторных и инструментальных исследований при поступлении в соответствии с установленными клиническими рекомендациями. Гипертензия определялась как систолическое артериальное давление ≥140 mmHg, диастолическое артериальное давление ≥90 mmHg или текущая антигипертензивная терапия. Сахарный диабет определялся как концентрация глюкозы в плазме натощак ≥7,0 mmol/L, уровень гликированного гемоглобина (HbA1c) ≥6,5% или применение препаратов, снижающих уровень глюкозы. Дислипидемия определялась по отклонениям показателей липидов сыворотки крови натощак или текущей гиполипидемической терапией. ИБС, инфаркт миокарда в анамнезе и перенесенный ишемический инсульт подтверждались на основании документально зафиксированного клинического анамнеза, а также архивных отчетов о нейроваскулярной или кардиологической визуализации, доступных в системе электронных медицинских карт.
Статистический анализ и разработка модели
Перед проведением статистического анализа нормальность распределения непрерывных переменных оценивали с помощью критерия Шапиро–Уилка. Описательная статистика для категориальных переменных представлена в виде абсолютного числа (в процентах), их сравнение проводили с использованием критерия хи-квадрат Пирсона или точного теста Фишера, в зависимости от ситуации. Непрерывные переменные с нормальным распределением выражали как средние значения ± стандартные отклонения и сравнивали с помощью t-критерия Стьюдента, в то время как переменные с ненормальным распределением выражали как медианы с интерквартильным размахом и сравнивали с помощью U-критерия Манна–Уитни. Для изучения возможных нелинейных зависимостей «доза-эффект» между количеством сохранившихся зубов и риском высокой степени выраженности ЦСБ (CSVD) была построена модель ограниченных кубических сплайнов (RCS), скорректированная по возрасту, гипертензии и Log_SII, с использованием четырех узлов, расположенных на 5-м, 35-м, 65-м и 95-м процентилях распределения количества зубов. Формальная проверка на нелинейность проводилась с помощью дисперсионного анализа (ANOVA), чтобы определить, что лучше объясняет наблюдаемую связь: пороговый эффект или непрерывный линейный тренд.
Для отбора переменных при наличии коррелирующих предикторов все исходные переменные были включены в модель регрессии LASSO. Оптимальный параметр настройки (λ) определяли с помощью 10-кратной перекрестной проверки; его выбирали по критерию минимальной биномиальной девиации (λmin). Для последующего анализа были оставлены переменные с ненулевыми коэффициентами. Отобранные переменные затем вводились в модель многофакторной логистической регрессии для оценки отношений шансов (OR) и 95% доверительных интервалов (CI) их связи с высокой степенью выраженности CSVD. На основе итоговой многофакторной модели была построена исследовательская клиническая номограмма для индивидуальной оценки риска.
Для оценки дополнительной значимости периодонтальных параметров была проведена сравнительная analiza базовой модели (включающей возраст, гипертензию и Log_SII) и расширенной модели, учитывающей степень периодонтита и количество зубов. Эффективность моделей оценивалась с помощью нескольких взаимодополняющих показателей. Дискриминационная способность оценивалась по площади под ROC-кривой (AUC), а различия между моделями сравнивались с помощью теста Делонга24. При оптимальном значении отсечки, определенном по максимальному индексу Юдена, также рассчитывались чувствительность, специфичность, прогностическая значимость положительного результата (PPV) и прогностическая значимость отрицательного результата (NPV). Для количественной оценки улучшения реклассификации риска, обеспечиваемого расширенной моделью, рассчитывались непрерывный индекс улучшения реклассификации (NRI) и интегрированный индекс улучшения дискриминации (IDI)25. Калибровка оценивалась с помощью калибровочных кривых, построенных на основе 1000 бутстреп-повторных выборок с использованием пакета rms; также рассчитывалась средняя абсолютная ошибка между прогнозируемым и наблюдаемым рисками. Чтобы учесть возможный переобучение и обеспечить несмещенную оценку эффективности модели, бутстреп-валидация включала коррекцию оптимизма, при этом сообщались как скорректированная по смещению, так и видимая эффективность модели. Анализ кривых принятия решений (DCA) был проведен для оценки потенциальной клинической полезности путем расчета чистого выигрыша в диапазоне пороговых вероятностей. Все статистические анализы и визуализации были выполнены с использованием статистического программного обеспечения и пакетов rms, glmnet, pROC, PredictABEL и dcurves. Статистическая значимость определялась как двустороннее значение P <0,05.