Витамин D и смертность: метаанализ стандартизированного 25-гидроксивитамина D по отдельным участникам

Витамин D участвует в регуляции гомеостаза кальция и оказывает благотворное влияние на здоровье скелета. Уровни 25-гидроксивитамина D (25[OH]D) в сыворотке измеряют для оценки статуса витамина D, который в основном определяется выработкой витамина D в коже под действием солнечного света (ультрафиолета-B) и, в меньшей степени, пищей или дополнительный прием витамина D. Плохой статус витамина D стал фактором риска различных неблагоприятных последствий для здоровья, включая смертность, но существуют разногласия по поводу классификации статуса витамина D и формы связи между концентрациями 25(OH)D и смертностью.

Существующие знания о статусе витамина D и смертности ограничены отсутствием стандартизации лабораторных методов. Предыдущие исследования показали, что различия в анализе и лабораторных условиях оказывают значительное влияние на сообщаемые концентрации 25(OH)D и, следовательно, на связь между 25(OH)D и последствиями для здоровья. Поэтому Программа стандартизации витамина D (VDSP), совместная инициатива, возглавляемая Национальным институтом здравоохранения и Управлением пищевых добавок (NIH-ODS), разработала протоколы для стандартизации данных по 25(OH)D из текущих и предыдущих исследований.

В этой работе, которая является частью проекта ЕС «Пищевые решения для ликвидации дефицита витамина D и укрепления здоровья на протяжении всего жизненного цикла» (ODIN), мы стремимся устранить пробел в знаниях о связи между стандартизированным 25(OH )D концентрации и смертность в совместном метаанализе с использованием данных отдельных участников (IPD) из восьми когорт исследований по всей Европе. Подробно мы изучим связь 25(OH)D со смертностью от всех причин, сердечно-сосудистых заболеваний и рака. Мы используем одноэтапный подход для этого метаанализа, который имеет то преимущество, что ИПЗ всех исследований моделируется одновременно, в то время как традиционные двухэтапные подходы основаны на совокупных данных каждого отдельного исследования. Учитывая, что в большинстве отдельных исследований этого метаанализа уже сообщалось об исходных данных по 25(OH)D, мы стремимся сравнить исходные и стандартизированные концентрации 25(OH)D в отношении различий в зарегистрированных концентрациях и их связи со смертностью.

Методы ИССЛЕДОВАНИЕ ИДЕНТИФИКАЦИЯ И ОТБОР Мы установили сотрудничество для проведения этого мета-анализа IPD. Потенциальные участники рабочего пакета Седьмой рамочной программы ЕС ODIN были приглашены на однодневный семинар в Амстердаме в ноябре 2012 года для обсуждения целей, реализации и развития задачи. Приглашенные европейские участники были отобраны на основе недавно опубликованных данных о больших предполагаемых когортах витамина D со смертностью и сердечно-сосудистыми исходами. Был ряд предпосылок для включения индивидуальных когортных исследований в ODIN: необходимо было иметь качественные биобанки образцов для единообразного отбора и анализа отобранных образцов, проверенные проспективные данные о клинических исходах, готовность к сотрудничеству и опыт в этой области.

ИССЛЕДОВАНИЯ И УЧАСТНИКИ Все отдельные когорты являются частью консорциума ODIN. Мы будем использовать IPD из восьми независимых проспективных когортных исследований из Норвегии, Германии, Дании, Нидерландов и Исландии. Включены исследования: 4-й обзор исследования Тромсё, Людвигсхафенское исследование риска и здоровья сердечно-сосудистой системы (LURIC), Рейкьявикское исследование возраста, генетической и окружающей среды (AGES), Нью-Хорнское исследование (NHS), Орхусское маммографическое когортное исследование, немецкое обследование состояния здоровья взрослых (DEGS) и когорта пожилых и молодых людей продольного исследования старения в Амстердаме (LASA). Все когортные исследования проводились в соответствии с Хельсинкской декларацией, и от всех участников исследования было получено письменное информированное согласие.

ИЗМЕРЕНИЕ 25(OH)D В соответствии с протоколом VDSP мы проведем процедуру отбора проб, чтобы выбрать подмножество из 100-150 биобанков образцов сыворотки из каждого отдельного когортного исследования для повторного анализа 25(OH)D с помощью стандартизированного и сертифицированный метод жидкостной хроматографии-тандемной масс-спектрометрии (LC-tandem MS), который восходит к эталонной процедуре измерения высшего порядка Национального института стандартов и технологий США (NIST). Повторно проанализированные значения 25(OH)D будут использоваться для разработки основных уравнений регрессии для каждого существующего когортного исследования и для повторной калибровки существующих измерений 25(OH)D. NHS не проводила предыдущих измерений 25(OH)D и будет полностью проанализирована.

ГРУППИРОВКА Население ИПЛ будет разделено на семь групп в соответствии с исходным статусом 25(OH)D. Распределение индивидуумов в одну из семи подгрупп будет выполнено для исходных и стандартизированных измерений 25(OH)D.

Согласно отчету Института медицины (IOM) за 2011 г., пороговые значения для групп 25(OH)D будут определены как с тяжелым дефицитом витамина D (≤29•99 нмоль/л), так как две группы пациентов с риском неадекватности (от 30 до ≤ 39•99 нмоль/л и от 40 до ≤ 49•99 нмоль/л), как достаточный уровень витамина D (от 50 до ≤ 75 нмоль/л), как две группы уровней витамина D, которые не всегда связаны с повышенной пользой (от 75 до ≤ 99•99 нмоль/л; от 100 до ≤ 124•99 нмоль/л) и высокий уровень витамина D, вызывающий беспокойство (≥ 125 нмоль/л; для преобразования нмоль/л в нг/мл разделите на 2•496).

СТАТИСТИЧЕСКИЙ АНАЛИЗ Различия между исходными и стандартизированными концентрациями 25(OH)D будут оцениваться с помощью t-критерия для парной выборки. Для сравнения других исходных характеристик в исходных стандартизированных группах витамина D мы будем использовать ANOVA для непрерывных и критерий χ2 для категорийных данных, в зависимости от ситуации.

Смертность от всех причин является первичным исходом и доступна во всех участвующих когортных исследованиях. Вторичными исходами являются смертность от сердечно-сосудистых заболеваний и смертность от рака, и они доступны во всех когортных исследованиях, кроме исследований NHS. Все конечные точки были найдены в мельчайших деталях из свидетельств о смерти, муниципальных реестров, медицинских карт и местных органов власти.

Все анализы исходов будут выполняться для исходных и стандартизированных значений 25(OH)D с помощью оценок метаанализа IPD и оценок, специфичных для исследования. Анализы будут основываться на лицах с полными данными о возрасте, поле, индексе массы тела (ИМТ), сезоне забора крови, уровне 25(OH)D, жизненном статусе при последующем наблюдении и во время наблюдения. Время наблюдения должно быть > 0 дней. Участники с отсутствующими данными будут исключены из анализа, и мы не будем выполнять вменение данных.

Связь между уровнями 25(OH)D и смертностью от всех причин будет оцениваться с использованием IPD в рамках одноэтапного подхода. В целом было показано, что метаанализ IPD после одноэтапной процедуры является более кратким подходом к бинарным результатам по сравнению с часто используемым двухэтапным подходом, в котором анализируются агрегированные данные. Мы будем использовать иерархическую параметрическую модель выживания, которая более приемлема по сравнению с моделью Кокса при анализе бинарных результатов. Единое уравнение будет обрабатываться в параметрической модели Вейбулла с ускоренным временем отказа (AFT), которая, как представляется, лучше всего соответствует базовым данным по сравнению с экспоненциальным, логарифмически-логистическим и логарифмически-нормальным распределениями. Модель будет построена с использованием SAS PROC NLMIXED (SAS Institute Inc., 100 SAS Campus Drive, Кэри, США) и случайного перехвата для учета случайных эффектов в когортных исследованиях.

Для анализа смертности 25(OH)D будет моделироваться с использованием 1) традиционного подхода категориальных переменных с группами в соответствии с ранее упомянутой классификацией IOM и 2) подхода ограниченных кубических сплайнов. Подход кубических сплайнов был выбран, чтобы сохранить непрерывный характер значений 25(OH)D и рассчитать отношения рисков (HR) с 95% доверительными интервалами (ДИ) для среднего значения каждой группы. В качестве эталона мы выбрали группу 25(OH)D с самым низким риском смертности. Наш анализ результатов будет кумулятивно скорректирован с учетом факторов риска смертности и детерминант статуса витамина D. В модели 1 мы делаем поправку на возраст (в годах), пол (мужской/женский) и сезон взятия крови (весна, лето, осень, зима). В модели 2, нашей основной статистической модели, мы дополнительно корректируем ИМТ (в кг/м²). В модели 3 мы дополнительно делаем поправку на сахарный диабет (да/нет) и артериальную гипертензию (да/нет), а в модели 4 добавляем онкологические заболевания в анамнезе (да/нет), сердечно-сосудистые заболевания в анамнезе (да/нет) и текущие статус курения как ковариант (да/нет).

Дополнительные настройки имеют модель два в качестве эталонной модели. Дополнительные корректирующие ковариаты доступны не в каждом когортном исследовании, поэтому дополнительные корректировки будут выполняться только в тех исследованиях, которые могут предоставить эти ковариаты.

Во-первых, во второй модели будет добавлен дополнительный прием кальция (да/нет). Первый дополнительный корректировочный анализ будет проводиться во всех исследованиях, кроме DEGS и LASA, молодой когорты, поскольку информация о дополнительном использовании отсутствует.

Во-вторых, дополнительное потребление витамина D (да/нет) будет добавлено ко второй модели во всех исследованиях, кроме LASA, когорты пожилых и молодых людей и DEGS, поскольку в этих исследованиях нет информации о дополнительном использовании витамина D.

В-третьих, дополнительная корректировка второй модели для физической активности (три фиктивные переменные для низкой, средней и высокой частоты физической активности) будет обработана во всех исследованиях, кроме когорты Орхусской маммографии. Для анализа чувствительности мы также опустим DEGS, поскольку среди участников DEGS высока доля молодых физически активных людей.

В-четвертых, во вторую модель будет добавлена ​​поправка на расчетную скорость клубочковой фильтрации (рСКФ; в мл/мин/1•73 м²). рСКФ будет рассчитываться на основе креатинина при исходном посещении в соответствии с уравнением исследования модификации диеты при заболеваниях почек (MDRD) с четырьмя переменными и будет добавлена ​​ко второй модели во всех исследованиях, кроме NHS, поскольку измерения креатинина недоступны в NHS.

В пятом анализе поправка на паратиреоидный гормон (в пмоль/л) будет добавлена ​​ко второй модели во всех исследованиях, кроме NHS, когорты маммографов в Орхусе, и LASA, когорты молодых.

В шестом анализе поправка на С-реактивный белок (в мг/л) будет добавлена ​​ко второй модели во всех исследованиях, кроме NHS, когорты маммографов в Орхусе, и LASA, когорты молодых.

В седьмом анализе ко второй модели будет добавлена ​​поправка на систолическое артериальное давление (в мм рт. ст.) во всех исследованиях, кроме когорты Орхусской маммографии и DEGS.

В восьмом анализе поправка на холестерин липопротеинов низкой плотности (в ммоль/л) будет добавлена ​​ко второй модели во всех исследованиях, кроме когорты Орхусской маммографии и DEGS.

В девятом анализе поправка на глюкозу (в ммоль/л) будет добавлена ​​ко второй модели во всех исследованиях, кроме исследования Тромсё, когорты маммографов в Орхусе и DEGS.

Все модели на исходном 25(OH)D будут выполняться без NHS, так как у NHS не было исходных измерений 25(OH)D. Для стандартизированного 25(OH)D мы вычислим все модели 1) с данными NHS и 2) без данных NHS, чтобы получить сопоставимые результаты между моделями исходных и стандартизированных измерений 25(OH)D.

Анализы чувствительности Анализы подгрупп будут проводиться для стратификации факторов риска дефицита витамина D и смертности. Подробно мы провели стратификацию по полу (женщины/мужчины), возрастным группам (<60 лет; от 60 до 69•9 лет; 70+ лет), группам ИМТ (<25 кг/м²; 25<30 кг/м²; ≥30 кг). / м²), добавки кальция (да/нет), добавки витамина D (да/нет), сердечно-сосудистые заболевания в анамнезе (да/нет) и рак в анамнезе (да/нет). Дальнейший анализ чувствительности будет ограничен лицами, умершими > 1 года и > 3 лет после исходного обследования, и ограничен когортами населения в целом (т. все когорты, кроме LURIC).

Вторичные исходы Для оценки вторичных исходов смертности от сердечно-сосудистых заболеваний и рака мы используем традиционные пропорциональные риски Кокса и модифицированную регрессию рисков в соответствии с методом Файна и Грея для учета конкурирующих рисков. Короче говоря, пропорциональные риски могут быть неудовлетворительными в многоцентровых условиях, поэтому допускается варьирование исходных рисков в исследованиях с одной когортой.