Bazal özellikler
NHANES 2011–2018'den toplam 8.218 katılımcı dahil edilmiş olup, bunlardan 623'ü DKD olan ve 7.595'i Non-DKD olarak sınıflandırılmıştır. DKD ve Non-DKD grupları arasında yaş, cinsiyet, ırk/etnisite, medeni durum, eğitim düzeyi, PIR, BMI, kan basıncı, hipertansiyon, alkol tüketimi, sigara kullanımı, vitamin B12 alımı ve diyetle alım açısından anlamlı farklar gözlemlenmiştir (tüm p < 0.05). Ayrıca ACR, eGFR, Scr, BUN, TC ve TG değerlerinde de anlamlı farklar saptanmıştır (Tablo 1).
Vitamin B12 alımı kuartilleri arasında eğitim düzeyi, hipertansiyon, DKD durumu ve enerji alımı açısından anlamlı farklılıklar gözlemlenmiştir. Buna karşın; yaş, medeni durum, BMI, kan basıncı, alkol tüketimi, eGFR, ACR, serum ürik asit, Scr veya TG değerleri açısından kuartiller arasında anlamlı bir fark görülmemiştir (Tablo 2).
Vitamin B12 ile DKD arasındaki ilişki
Vitamin B12 alımı kuartillere (Q1–Q4) ayrılmış olup, çok değişkenli lojistik regresyon analizlerinin sonuçları Tablo 3'te sunulmuştur. Tam olarak düzeltilmiş Model 3'te, en yüksek vitamin B12 alımı kuartilindeki (Q4) katılımcıların DKD olasılığı, en düşük kuartildeki (Q1) katılımcılara göre anlamlı derecede daha düşüktü (OR = 0.70, 95% CI: 0.50–0.97, p = 0.034); buna karşın Q2 (OR = 1.14, 95% CI: 0.85–1.53, p = 0.386) veya Q3 (OR = 0.94, 95% CI: 0.69–1.28, p = 0.684) için anlamlı bir ilişki gözlemlenmemiştir. Ayrıca, ilişkinin genel yönü ters olsa da, tam olarak düzeltilmiş modeldeki doğrusal eğilim istatistiksel anlamlılığa ulaşmamıştır (p for trend = 0.086). Bu bulgular, özellikle en yüksek kuartildeki daha yüksek vitamin B12 alımının daha düşük DKD olasılığı ile ilişkili olduğunu ancak sonuçların vitamin B12 alımı ile DKD arasında sabit bir doğrusal doz-yanıt ilişkisini desteklemediğini göstermektedir.
Vitamin B12 ile DKD arasındaki ilişkinin alt grup analizleri
Vitamin B12 alımı ile DKD arasındaki ilişkinin farklı katılımcı özellikleri genelinde tutarlı olup olmadığını incelemek için alt grup analizleri yapılmıştır. Model 1'de, en yüksek vitamin B12 alımı kuartilindeki (Q4) katılımcıların, hem erkek hem de kadın alt gruplarında Q1'dekilere göre DKD olasılıkları anlamlı derecede daha düşük bulunmuştur (p < 0.05; Tablo 4). ACR alt grupları arasındaki ilişkiler daha heterojen seyretmiştir. ACR < 30 mg/g olan katılımcılar arasında, bazı orta düzey vitamin B12 alımı kuartilleri için anlamlı ilişkiler gözlemlenirken, daha yüksek ACR seviyelerine göre tabakalaştırma yapıldığında bu örüntü zayıflamıştır. eGFR alt grupları genelinde anlamlı bir ilişme gözlemlenmemiştir.
Vitamin B12 ile DKD arasındaki doğrusal olmayan ilişki
Tam düzeltilmiş RCS analizi, vitamin B ile arasında anlamlı, doğrusal olmayan bir ilişki olduğunu göstermiştir12 alımı ve DKD (Şekil 2, p doğrusallık dışı değer = 0,037). RCS modeli; vitamin B'nin 5., 35., 65. ve 95. persentillerinde yer alan dört düğüm kullanmıştır12 giriş dağılımı, ile 3.24 μg/day referans değer olarak kullanılmıştır. Sonraki kırılma noktası analizi, şu noktalarda veri güdümlü iki kırılma noktası belirlemiştir: 3.24 μg/day (%95 GA: 2,86–3,62) ve 8.36 μg/day (%95 GA: 7,52–9,20), vitamin B'yi bölen12 üç aralığa ayrılmıştır. B vitamini için12 alım < 3.24 μg/day, DKD ile anlamlı bir ilişki gözlemlenmemiştir (OR = 1,07, p = 0,575). 3,24– aralığında8.36 μg/day, vitamin B12 alım, DKD ile ters ilişkiliydi (OR = 0,87, p = 0,025). Alım için > 8.36 μg/day, anlamlı bir ilişki gözlemlenmemiştir (OR = 1.00, p = 0.995) (Tablo 5). Bu bulgular, yalnızca orta düzeydeki alım aralığında anlamlı bir ters ilişkiyle gözlemlenen doğrusal olmayan bir örüntü ile tutarlıdır. Bu bulguların sağlamlığını değerlendirmek için üç ve beş düğümlü alternatif RCS spesifikasyonları kullanılarak duyarlılık analizleri gerçekleştirilmiştir. Her iki düğüm spesifikasyonu altında Model 1–3 genelinde istatistiksel olarak anlamlı genel ilişkilerle birlikte benzer genel ilişki örüntüleri gözlemlenmiş; beş düğümlü modeller ise doğrusal olmayan bir ilişkiyi daha fazla desteklemiştir (Ek Şekil 1).
DKD'de aday vitamin B12-ilişkili genlerin belirlenmesi
GSE30528 veri setinde, DKD ve normal örnekler arasında 126'sı yukarı yönlü ve 367'si aşağı yönlü düzenlenen genlerden oluşan 493 DEG saptadık (|log₂FC| > 1, FDR < 0.05; Şekil 3A,B, Ek Tablo 1). KEGG yolak analizi, bu DEG'lerin hücre dışı matris (ECM)-reseptör etkileşimi, fokal adezyon, integrin sinyalizasyonu ve aktin hücre iskeleti düzenlemesi yolaklarında anlamlı şekilde zenginleştiğini ortaya koydu (Şekil 3C, Ek Tablo 2). GO zenginleştirme analizi, DEG'lerin ağırlıklı olarak kolajen içeren ECM, glomerulus gelişimi, fokal adezyon ve nefron gelişimi ile ilişkili olduğunu gösterdi (Şekil 3D, Ek Tablo 2).
DEGs ile vitamin B12 ilişkili gen setinin kesişim analizi, örtüşen 20 aday gen belirledi (Şekil 3E). Bu genler, daha sonra etkileşimleri minimum 0,4 güven puanına göre filtrelenerek STRING veritabanı kullanılarak bir PPI ağı oluşturmak için kullanıldı. Ortaya çıkan ağ, sırasıyla proteinleri ve etkileşimlerini temsil eden 11 düğüm ve 20 kenardan oluşmuştur (Şekil 3F). 20 aday gen arasındaki ilişkiler GeneMANIA kullanılarak ayrıca görselleştirilmiştir (Şekil 3G).
DKD'de aday bir vitamin B12 ilişkili gen olarak P2RY14'nin keşifsel karakterizasyonu
Vitamin B ile gözlemlenen ilişkiye dayanarak12 alım ve DKD durumuyla birlikte, bir sonraki adımda temel B vitamini12transkriptomik analiz kullanarak DKD patogenezinde rol oynayan -ilişkili genler belirlenmiştir. Örtüşen 20 aday gen arasından, bu genlerin farklı DKD veri setlerindeki ekspresyon modelleri ayrıca karşılaştırılmıştır. GSE30528 ve GSE30529 veri setlerinde, P2RY14 ekspresyonu, normal kontrollerle karşılaştırıldığında DKD örneklerinde anlamlı düzeyde artmıştır (Şekil 4A,B). P2RY14 bir G proteini kenetli reseptörü kodlar; immün hücre kemotaksisi, hücre iskeleti yeniden düzenlenmesi, inflamatuar yanıtlar ve oksidatif stres ile ilişkili olduğu bildirilmiştir28,29Ancak, DKD üzerindeki rolü henüz doğrudan araştırılmamıştır. Bağımsız DKD veri setleri genelindeki tutarlı artışı ve potansiyel biyolojik önemi göz önüne alındığında, P2RY14 sonraki doğrulama ve fonksiyonel analizler için önceliklendirildi. Ardından, ...'nın ayırt edici yeteneğini ön değerlendirmek amacıyla ROC eğrisi analizi gerçekleştirildi. P2RY14GSE30528 keşif veri setinde AUC 0,923 (%95 GA: 0,807–1,000) iken, GSE30529 ön dış doğrulama veri setinde AUC 0,883 (%95 GA: 0,744–1,000) idi. GSE142025 ek doğrulama veri setinde ise AUC 0,663 (%95 GA: 0,478–0,847) olarak belirlendi ve bu durum sınırlı bir ayırt edici yeteneğe işaret etti (Şekil 4C–E). Bir bütün olarak bu bulgular, ...'nın potansiyel ayırıcı değerine ilişkin ön kanıtlar sunmaktadır P2RY14ancak tanısal performansının daha büyük bağımsız kohortlarda ayrıca doğrulanması gerekmektedir.
DKD'deki P2RY14 ekspresyonu ile ilişkili biyolojik fonksiyonları ve yolakları araştırmak için, GSE30528 örnekleri P2RY14'nin medyan ekspresyon düzeyine göre P2RY14-yüksek ve P2RY14-düşük gruplarına ayrılmıştır. İki grup arasında diferansiyel ekspresyon analizi gerçekleştirilmiştir. Ardından, P2RY14 ekspresyonu ile ilişkili biyolojik fonksiyonları ve yolakları karakterize etmek için GO ve KEGG zenginleştirme analizleri yürütülmüştür. GO zenginleştirme analizi, P2RY14'nün kolajen içeren ECM, glomerulus gelişimi, ECM yapısal bileşeni ve böbrek gelişimi ile ilişkili olduğunu ortaya koymuştur (Şekil 4F, Ek Tablo 3). KEGG analizi, P2RY14'nün ECM-reseptör etkileşimi, fokal adezyon ile kompleman ve koagülasyon kaskadlarında rol oynadığını daha ileri düzeyde göstermiştir (Şekil 4G, Ek Tablo 3).
Ek olarak, GSEA, NF-κB sinyal yolunun ve TNF sinyal yolunun P2RY14-yüksek grubunda anlamlı şekilde zenginleştiğini göstermiştir; bu durum, yüksek P2RY14 ekspresyonunun DKD'deki enflamatuar sinyalizasyon ile ilişkili olabileceğini düşündürmektedir (Şekil 4H, Ek Tablo 3).
MC-1'in P2RY14 ile keşifsel moleküler kı husbandırması
GeneCards ve DGIdb taramaları aracılığıyla P2RY14 ile ilişkili beş aday bileşik tanımlanmış ve keşif amaçlı moleküler docking analizi için MC-1 seçilmiştir (Şekil 5A). Ligand olarak MC-1 ve reseptör olarak P2RY14 kullanılarak yapılan docking analizi, MC-1'in P2RY14'nin bir aday cebi içinde −6.0 kcal/mol Vina skoru ile potansiyel bir bağlanma pozuna sahip olduğunu öngörmüştür. Öngörülen poz; Tyr, Lys, Asn, Asp ve Glu kalıntıları ile etkileşimleri ve yedi hidrojen bağının oluşumunu içermektedir (Şekil 5B–D).
Veri Kullanılabilirliği Beyanı
Bu çalışmada kullanılan veriler halka açıktır. 2011–2018 anket döngülerine ait NHANES verilerine Ulusal Sağlık İstatistikleri Merkezi'nden (National Center for Health Statistics) ulaşılabilir. GSE30528, GSE142025 ve GSE30529 transkriptomik veri setleri Gene Expression Omnibus (GEO) veri tabanında halka açık olarak sunulmaktadır. Bu çalışmada oluşturulan işlenmiş diferansiyel ekspresyon, fonksiyonel zenginleştirme ve GSEA sonuçları Ek Dosyalarda sunulmuştur.

Şekil 1: NHANES analizi için katılımcı seçim akış şeması. Akış şeması, 2011–2018 NHANES döngülerindeki katılımcılara uygulanan ardışık dışlama kriterlerini ve analize dahil edilen nihai çalışma popülasyonunu göstermektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayınız.

Şekil 2: Kısıtlanmış kübik spline modeli kullanılarak vitamin B12 alımı ile DKD arasındaki doğrusal olmayan ilişki. Düz çizgi, diyetle alınan vitamin B12 miktarına göre DKD için tahmin edilen olasılık oranını, gölgeli alan ise %95 güven aralığını temsil etmektedir. Model, Model 3'te yer alan kovaryantlar kullanılarak düzeltilmiştir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 3: DKD'de aday vitamin B12-ilişkili genlerin tanımlanması.(A,B) GSE30528 veri setindeki diyabetik böbrek hastalığı (DKD) ve normal örnekler arasındaki DEG'lerin ekspresyon profillerini gösteren yanardağ grafiği ve ısı haritası. (C,D) DEG'lerin en zenginleşmiş ilk 15 KEGG yolu ve GO terimleri. (E) DEG'ler ve vitamin B12-ilişkili genler arasındaki çapraz analiz. (F) Protein-protein etkileşim (PPI) ağı. (G) 20 vitamin B12‑ilişkili genin GeneMANIA ağı. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 4: DKD'de aday vitamin B12-ilişkili bir gen olarak P2RY14'ün keşifsel karakterizasyonu. (A,B) Sırasıyla GSE30528 ve GSE30529 veri setlerindeki DKD örneklerinde ve normal kontrollerde örtüşen aday genlerin ekspresyon seviyelerini gösteren kutu grafikleri. (C–E) Sırasıyla GSE30528, GSE30529 ve GSE142025 veri setlerinde P2RY14'ün alıcı işletim karakteristik (ROC) eğrisi analizleri. (F) Gen Ontolojisi (GO) zenginleştirme analizi. (G) Kyoto Gen ve Genom Ansiklopedisi (KEGG) yolak zenginleştirme analizi. (H) GSE30528'de NF-κB ve TNF sinyalizasyon yolaklarının zenginleşmesini gösteren P2RY14-yüksek ve P2RY14-düşük gruplarının GSEA'sı. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 5: P2RY14'nin ilaç-hedef taraması ve keşifsel moleküler docking analizi. (A) P2RY14 ile ilişkili ilaçların ilaç-hedef etkileşim ağı. (B–D) MC-1'in P2RY14 proteini ile moleküler docking çalışması. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Tablo 1: Çalışma popülasyonunun temel özellikleri. Katılımcıların DKD durumuna göre demografik, klinik, biyokimyasal ve diyetsel özellikleri sunulmuştur. Sürekli ve kategorik değişkenler, ankete göre ağırlıklandırılmış analiz kullanılarak özetlenmiştir. Bu tabloyu indirmek için lütfen buraya tıklayın.
Tablo 2: Vitamin B12 alımı kuartillerine göre katılımcıların temel özellikleri. Çalışma popülasyonunun demografik, klinik, biyokimyasal ve diyetle ilgili özellikleri, diyetle alınan vitamin B12 kuartillerine göre sunulmuştur. Bu Tabloyu indirmek için lütfen buraya tıklayın.
Tablo 3: Vitamin B ile ilişki12 ABD'li yetişkinlerde alım ve DKD durumu. Diyetle alınan vitamin B ile ilişki için olasılık oranları (OR'lar) ve %95 güven aralıkları (GA'lar)12 alım miktarları ve DKD durumu, en düşük vitamin B ile başlayan ve kademeli olarak düzeltilmiş lojistik regresyon modelleri üzerinden sunulmaktadır12 referans grup olarak alım çeyreği (Q1). Lütfen bu Tabloyu indirmek için buraya tıklayın.
Tablo 4: Alt grup analizlerinde vitamin B12 alımı ile DKD durumu arasındaki ilişki. DKD için odds oranları (OR'lar) ve %95 güven aralıkları (GA'lar); yaş, cinsiyet, BMI, eGFR, ACR, hipertansiyon ve hiperlipidemi alt gruplarına göre stratifiye edilmiş vitamin B12 alım kuartilleri üzerinden ve belirtilen modeller için ayarlanmış olarak sunulmuştur. Bu tabloyu indirmek için lütfen buraya tıklayın.
Tablo 5: Vitamin B12 alımı ile DKD arasındaki ilişkiye dair üç segmentli parçalı regresyon analizi. Diyetsel vitamin B12 alımı ile DKD arasındaki ilişkiler, 3,24 ve 8,36 μg/gün değerlerindeki iki veri tabanlı kırılma noktasıyla tanımlanan üç alım aralığı için sunulmuştur. Bu tabloyu indirmek için lütfen buraya tıklayın.
Ek Şekil 1: Diyetteki vitamin B ile olan ilişkiye dair duyarlılık analizleri12 kısıtlanmış kübik spline modelleri kullanılarak alım ve DKD. (A–C) Üç düğümlü modeller. (D–F) Beş düğümlü modeller. Düz çizgi tahmini OR'yi, gölgeli alan %95 GA'yı ve kesikli çizgi OR = 1'i temsil etmektedir. p genel ilişki ve doğrusal olmama değerleri her bir panelde gösterilmiştir.Lütfen bu dosyayı indirmek için buraya tıklayın.
Ek Tablo 1: GSE30528 veri setindeki diyabetik böbrek hastalığı (DKD) ve normal örnekler arasındaki diferansiyel ekspresyon gösteren genler. |log₂FC| > 1 ve FDR < 0.05 değerlerine göre belirlenen genlerin diferansiyel ekspresyon sonuçları, aday gen taraması için kullanılan kobalamin ilişkili gen setiyle birlikte sunulmuştur.Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Tablo 2: DKD'de diferansiyel ekspresyon gösteren genlerin fonksiyonel zenginleştirme analizi. GSE30528 veri setinde tanımlanan diferansiyel ekspresyon gösteren genlere ait Gen Ontolojisi (GO) ve Kyoto Gen ve Genom Ansiklopedisi (KEGG) zenginleştirme sonuçları sunulmuştur.Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Tablo 3: GSE30528 veri setinde P2RY14 ile ilişkili fonksiyonel zenginleştirme ve gen seti zenginleştirme analizi sonuçları. P2RY14 ile ilişkili biyolojik fonksiyonları ve yolakları tanımlamak için kullanılan GO, KEGG ve GSEA sonuçları sunulmuştur.Bu dosyayı indirmek için lütfen buraya tıklayın.