Исходные характеристики
В исследование было включено в общей сложности 8 218 участников NHANES за 2011–2018 гг., из которых 623 были классифицированы как пациенты с ДБП и 7 595 как пациенты без ДБП. Значимые различия между группами с ДБП и без ДБП наблюдались по возрасту, полу, расе/этнической принадлежности, семейному положению, уровню образования, PIR, ИМТ, артериальному давлению, наличию гипертонии, потреблению алкоголя, статусу курения, потреблению витамина B12 и рациону питания (все p < 0,05). Также были отмечены значимые различия в показателях ACR, рСКФ, Scr, BUN, TC и TG (Таблица 1).
Значимые различия между квартилями потребления витамина B12 наблюдались в уровне образования, наличии гипертензии, статусе DKD и потреблении энергии. Напротив, значимых различий между квартилями не было выявлено по возрасту, семейному положению, ИМТ, артериальному давлению, потреблению алкоголя, рСКФ, соотношению ACR, уровню мочевой кислоты в сыворотке, Scr или TG (Таблица 2).
Связь между витамином B12 и ДКП
Потребление витамина B12 было разделено на квартили (Q1–Q4), и результаты многофакторного логистического регрессионного анализа представлены в Таблице 3. В полностью скорректированной Модели 3 у участников с самым высоким уровнем потребления витамина B12 (Q4) шансы развития DKD были значительно ниже, чем у участников с самым низким уровнем потребления (Q1) (OR = 0.70, 95% CI: 0.50–0.97, p = 0.034), в то время как для Q2 (OR = 1.14, 95% CI: 0.85–1.53, p = 0.386) или Q3 (OR = 0.94, 95% CI: 0.69–1.28, p = 0.684) значимых ассоциаций обнаружено не было. Кроме того, линейный тренд в полностью скорректированной модели не достиг статистической значимости (p для тренда = 0.086), хотя общее направление ассоциации было обратным. Эти результаты позволяют предположить, что более высокое потребление витамина B12, особенно в верхнем квартиле, было связано с более низкими шансами развития DKD, однако полученные данные не подтверждают наличие стабильной линейной зависимости «доза — ответ» между потреблением витамина B12 и DKD.
Анализ в подгруппах взаимосвязи между витамином B12 и ДКП
Был проведен анализ в подгруппах, чтобы проверить, является ли связь между потреблением витамина B12 и ДБП однородной при различных характеристиках участников. В модели 1 участники в самом высоком квартиле потребления витамина B12 (Q4) имели значимо более низкую вероятность развития ДБП, чем участники в Q1, как в подгруппе мужчин, так и в подгруппе женщин (p < 0.05; Таблица 4). Связи в подгруппах по уровню ACR были более неоднородными. Среди участников с ACR < 30 mg/g значимые связи наблюдались для некоторых промежуточных квартилей потребления витамина B12, тогда как при стратификации по более высоким уровням ACR эта закономерность сглаживалась. В подгруппах по уровню eGFR значимых связей выявлено не было.
Нелинейная связь между витамином B12 и ДКП
Полностью скорректированный анализ RCS выявил значимую нелинейную связь между потреблением витамина B12 и ДКП (Рисунок 2, p для нелинейности = 0,037). В модели RCS использовались четыре узла, расположенные на 5-м, 35-м, 65-м и 95-м процентилях распределения потребления витамина B12, при этом в качестве референтного значения было принято 3,24 μg/day. Последующий анализ точек перегиба выявил две определенные данными точки перегиба при 3,24 μg/day (95% CI: 2,86–3,62) и 8,36 μg/day (95% CI: 7,52–9,20), что разделило потребление витамина B12 на три диапазона. При потреблении витамина B12 < 3,24 μg/day значимой связи с ДКП не наблюдалось (OR = 1,07, p = 0,575). В диапазоне 3,24–8,36 μg/day потребление витамина B12 было обратно связано с ДКП (OR = 0,87, p = 0,025). При потреблении > 8,36 μg/day значимой связи не наблюдалось (OR = 1,00, p = 0,995) (Таблица 5). Эти результаты соответствуют нелинейному паттерну, при котором значимая обратная связь наблюдалась только в диапазоне промежуточного потребления. Для оценки надежности этих результатов был проведен анализ чувствительности с использованием альтернативных спецификаций RCS с тремя и пятью узлами. Были отмечены схожие общие паттерны связи со статистически значимыми общими ассоциациями в Моделях 1–3 при обеих спецификациях узлов, в то время как модели с пятью узлами дополнительно подтвердили нелинейный характер связи (Дополнительный рисунок 1).
Идентификация генов-кандидатов, связанных с витамином B12, при ДКП
В наборе данных GSE30528 мы выявили 493 дифференциально экспрессируемых гена (ДЭГ) между образцами с диабетической болезнью почек (ДБП) и нормальными образцами, включая 126 гнов с повышенной экспрессией и 367 генов с пониженной экспрессией (|log₂FC| > 1, FDR < 0.05; Рисунок 3A,B, Дополнительная таблица 1). Анализ путей KEGG выявил значительное обогащение этих ДЭГ в путях взаимодействия внеклеточного матрикса (ВКМ) с рецепторами, фокальной адгезии, интегринового сигналинга и регуляции актинового цитоскелета (Рисунок 3C, Дополнительная таблица 2). Анализ обогащения GO дополнительно показал, что ДЭГ преимущественно связаны с коллагенсодержащим ВКМ, развитием клубочков, фокальной адгезией и развитием нефронов (Рисунок 3D, Дополнительная таблица 2).
Анализ пересечения дифференциально экспрессируемых генов (DEG) и набора генов, связанных с витамином B12, позволил выявить 20 перекрывающихся генов-кандидатов (Рисунок 3E). Эти гены были впоследствии использованы для построения сети белок-белковых взаимодействий (PPI) с помощью базы данных STRING, при этом взаимодействия фильтровались по минимальному показателю достоверности 0,4. Полученная сеть состояла из 11 узлов и 20 ребер, представляющих белки и их взаимодействия соответственно (Рисунок 3F). Связи между 20 генами-кандидатами были дополнительно визуализированы с помощью GeneMANIA (Рисунок 3G).
Предварительная характеристика P2RY14 как потенциального гена, связанного с витамином B12, при диабетической нефропатии
Основываясь на наблюдаемой связи между потреблением витамина B12 и статусом диабетической болезни почек (ДБП), мы приступили к поиску ключевых генов, связанных с витамином B12 и участвующих в патогенезе ДБП, с помощью транскриптомного анализа. Среди 20 перекрывающихся генов-кандидатов мы дополнительно сравнили паттерны их экспрессии в различных наборах данных по ДБП. В наборах данных GSE30528 и GSE30529 экспрессия P2RY14 была значительно повышена в образцах с ДБП по сравнению с нормальным контролем (Рисунок 4A,B). P2RY14 кодирует G-белоксопряженный рецептор; сообщалось о его участии в хемотаксисе иммунных клеток, перестройке цитоскелета, воспалительных реакциях и окислительном стрессе28,29. Однако его роль при ДБП еще не была напрямую исследована. Учитывая стабильное повышение экспрессии в независимых наборах данных по ДБП и его потенциальную биологическую значимость, ген P2RY14 был выбран в качестве приоритетного для последующей валидации и функционального анализа. Затем был проведен анализ ROC-кривых для предварительной оценки дискриминационной способности P2RY14. В открывающем наборе данных GSE30528 значение AUC составило 0,923 (95% ДИ: 0,807–1,000), тогда как в предварительном внешнем валидационном наборе GSE30529 AUC составило 0,883 (95% ДИ: 0,744–1,000). В дополнительном валидационном наборе GSE142025 AUC составило 0,663 (95% ДИ: 0,478–0,847), что указывает на ограниченную дискриминационную способность (Рисунок 4C–E). В совокупности эти результаты предоставляют предварительные доказательства потенциальной диагностической ценности P2RY14, хотя его диагностическая эффективность требует дальнейшего подтверждения на более крупных независимых когортах.
Для изучения биологических функций и путей, связанных с P2RY14 экспрессия при ДКП, образцы GSE30528 были разделены на P2RY14высокий и P2RY14—группы с низким уровнем экспрессии в соответствии с медианным уровнем экспрессии P2RY14Был проведен анализ дифференциальной экспрессии между двумя группами. Затем были выполнены анализы обогащения GO и KEGG для характеристики биологических функций и путей, связанных с экспрессией P2RY14. Анализ обогащения GO выявил, что P2RY14 был связан с коллагенсодержащим внеклеточным матриксом (ВКМ), развитием клубочков, структурным компонентом ВКМ и развитием почек (Рисунок 4F, Дополнительная таблица 3). Анализ KEGG дополнительно указал на участие P2RY14 во взаимодействии внеклеточного матрикса с рецепторами, фокальной адгезии, а также в каскадах активации комплемента и коагуляции (Рисунок 4G, Дополнительная таблица 3).
Кроме того, анализ GSEA показал, что сигнальный путь NF-κB и сигнальный путь TNF были значительно обогащены в группе с высоким уровнем экспрессии P2RY14 , что позволяет предположить, что высокая экспрессия P2RY14 может быть связана с воспалительной сигнализацией при ДКП (Рисунок 4H, Дополнительная таблица 3).
Исследовательский молекулярный докинг MC-1 с P2RY14
С помощью скрининга GeneCards и DGIdb было выявлено пять соединений-кандидатов, связанных с P2RY14, из которых MC-1 было выбрано для предварительного анализа молекулярного докинга (Рисунок 5A). При использовании MC-1 в качестве лиганда и P2RY14 в качестве рецептора анализ докинга предсказал потенциальную позу связывания MC-1 в одном из предполагаемых карманов P2RY14 со значением Vina score −6.0 kcal/mol. Предсказанная поза включала взаимодействия с остатками Tyr, Lys, Asn, Asp и Glu, а также образование семи водородных связей (Рисунок 5B–D).
Заявление о доступности данных
Данные, использованные в этом исследовании, находятся в открытом доступе. Данные NHANES за циклы обследования 2011–2018 гг. доступны в Национальном центре статистики здравоохранения. Транскриптомные наборы данных GSE30528, GSE142025 и GSE30529 общедоступны в базе данных Gene Expression Omnibus (GEO). Обработанные результаты анализа дифференциальной экспрессии, функционального обогащения и GSEA, полученные в ходе данного исследования, представлены в дополнительных файлах.

Рисунок 1: Блок-схема отбора участников для анализа NHANES. На блок-схеме показаны последовательные критерии исключения, примененные к участникам циклов NHANES 2011–2018 годов, и итоговая исследуемая популяция, включенная в анализ. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 2: Нелинейная связь между потреблением витамина B12 и DKD с использованием модели ограниченного кубического сплайна. Сплошная линия представляет расчетное отношение шансов развития DKD в зависимости от потребления витамина B12 с пищей, а заштрихованная область — 95% доверительный интервал. Модель была скорректирована с использованием ковариат, включенных в Модель 3. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 3: Идентификация потенциальных генов, связанных с витамином B12, при ДКП. (A,B) График «вулкан» (volcano plot) и тепловая карта, демонстрирующие профили экспрессии дифференциально экспрессируемых генов (ДЭГ) в образцах с диабетической болезнью почек (ДКП) и нормальных образцах из набора данных GSE30528. (C,D) 15 наиболее обогащенных путей KEGG и терминов GO для ДЭГ. (E) Перекрестный анализ ДЭГ и генов, связанных с витамином B12. (F) Сеть белок-белковых взаимодействий (PPI). (G) Сеть GeneMANIA для 20 генов, связанных с витамином B12. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 4: Исследовательская характеристика P2RY14 как ген-кандидата, связанного с витамином B12, при диабетической болезни почек (ДБП). (A,B) Диаграммы размаха, демонстрирующие уровни экспрессии перекрывающихся генов-кандидатов в образцах с ДБП и в нормальном контроле в наборах данных GSE30528 и GSE30529 соответственно. (C–E) Анализ рабочих характеристик приемника (ROC-кривые) для P2RY14 в наборах данных GSE30528, GSE30529 и GSE142025 соответственно. (F) Анализ обогащения генных онтологий (GO). (G) Анализ обогащения путей Киотской энциклопедии генов и геномов (KEGG). (H) GSEA групп с высокой и низкой экспрессией P2RY14, демонстрирующий обогащение сигнальных путей NF-κB и TNF в GSE30528. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 5: Скрининг взаимодействия препаратов с мишенями и эксплораторный анализ молекулярного докинга P2RY14. (A) Сеть взаимодействий между лекарственными средствами и мишенями для препаратов, связанных с P2RY14. (B–D) Молекулярный докинг MC-1 с белком P2RY14. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Таблица 1: Основные характеристики исследуемой популяции. Представлены демографические, клинические, биохимические и диетические характеристики участников в зависимости от статуса DKD. Непрерывные и категориальные переменные обобщены с использованием анализа с учетом весов выборки. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 2: Основные характеристики участников в зависимости от квартилей потребления витамина B12. Демографические, клинические, биохимические и диетические характеристики исследуемой популяции представлены по квартилям потребления витамина B12 с пищей. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 3: Связь между потреблением витамина B12 и статусом DKD среди взрослых жителей США. Отношения шансов (OR) и 95% доверительные интервалы (CI) для связи между потреблением витамина B12 с пищей и статусом DKD представлены в последовательно скорректированных моделях логистической регрессии, где в качестве референтной группы использовался нижний квартиль потребления витамина B12 (Q1). Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 4: Связь между потреблением витамина B12 и статусом DKD в анализах подгрупп. Отношения шансов (OR) и 95% доверительные интервалы (CI) для DKD представлены по квартилям потребления витамина B12, стратифицированные по подгруппам возраста, пола, ИМТ, СКФ, соотношения альбумин/креатинин, артериальной гипертензии и гиперлипидемии, и скорректированные согласно указанным моделям. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 5: Трехсегментный кусочно-линейный регрессионный анализ связи между потреблением витамина B12 и ДКП. Связи между потреблением витамина B12 с пищей и ДКП представлены для трех диапазонов потребления, определенных двумя точками излома на основе данных при 3.24 и 8.36 μg/day. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Дополнительный рисунок 1: Анализ чувствительности связи между потреблением витамина B12 с пищей и ДБП с использованием моделей ограниченных кубических сплайнов. (A–C) Модели с тремя узлами. (D–F) Модели с пятью узлами. Сплошная линия представляет расчетный ОШ, заштрихованная область — 95% ДИ, а пунктирная линия указывает на ОШ = 1. Значения p для общей связи и нелинейности приведены на каждой панели.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительная таблица 1: Дифференциально экспрессируемые гены при сравнении образцов с диабетической болезнью почек (ДБП) и нормальных образцов в наборе данных GSE30528.Результаты дифференциальной экспрессии для генов, определенных по критериям |log₂FC| > 1 и FDR < 0.05, представлены вместе с набором кобаламин-зависимых генов, использованным для скрининга генов-кандидатов.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительная таблица 2: Анализ функционального обогащения дифференциально экспрессируемых генов при ДКП.** Представлены результаты анализа обогащения по Gene Ontology (GO) и Киотской энциклопедии генов и геномов (KEGG) для дифференциально экспрессируемых генов, выявленных в наборе данных GSE30528.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительные материалы Table 3: Результаты функционального обогащения и анализа обогащения наборов генов, связанных с P2RY14, в наборе данных GSE30528. Представлены результаты анализа GO, KEGG и GSEA, использованные для характеристики биологических функций и путей, связанных с P2RY14.Пожалуйста, нажмите здесь, чтобы скачать этот файл.