Basischarakteristika
Insgesamt wurden 8.218 Teilnehmer aus NHANES 2011–2018 eingeschlossen, von denen 623 als an DKD erkrankt und 7.595 als Nicht-DKD klassifiziert wurden. Zwischen den Gruppen DKD und Nicht-DKD zeigten sich signifikante Unterschiede hinsichtlich Alter, Geschlecht, Rasse/Ethnizität, Familienstand, Bildungsabschluss, PIR, BMI, Blutdruck, Hypertonie, Alkoholkonsum, Raucherstatus, Vitamin-B12-Aufnahme und Ernährungsaufnahme (alle p < 0,05). Signifikante Unterschiede zeigten sich ebenfalls bei ACR, eGFR, Scr, BUN, TC und TG (Tabelle 1).
Signifikante Unterschiede zwischen den Quartilen der Vitamin-B12-Aufnahme wurden hinsichtlich des Bildungsabschlusses, der Hypertonie, des DKD-Status und der Energiezufuhr beobachtet. Im Gegensatz dazu zeigten sich zwischen den Quartilen keine signifikanten Unterschiede hinsichtlich Alter, Familienstand, BMI, Blutdruck, Alkoholkonsum, eGFR, ACR, serumuremischer Säure, Scr oder TG (Tabelle 2).
Zusammenhang zwischen Vitamin B12 und DKD
Die Zufuhr von Vitamin B12 wurde in Quartile (Q1–Q4) eingeteilt, und die Ergebnisse der multivariablen logistischen Regressionsanalysen sind in Tabelle 3 dargestellt. Im vollständig adjustierten Modell 3 wiesen Teilnehmer im höchsten Quartil der Vitamin-B12-Zufuhr (Q4) signifikant niedrigere Chancen für DKD auf als jene im niedrigsten Quartil (Q1) (OR = 0,70, 95 % KI: 0,50–0,97, p = 0,034), während für Q2 (OR = 1,14, 95 % KI: 0,85–1,53, p = 0,386) oder Q3 (OR = 0,94, 95 % KI: 0,69–1,28, p = 0,684) keine signifikanten Assoziationen beobachtet wurden. Zudem erreichte der lineare Trend im vollständig adjustierten Modell keine statistische Signifikanz (p für Trend = 0,086), obwohl die allgemeine Richtung der Assoziation invers war. Diese Ergebnisse deuten darauf hin, dass eine höhere Zufuhr von Vitamin B12, insbesondere im höchsten Quartil, mit geringeren Chancen für DKD assoziiert ist, die Ergebnisse stützen jedoch keine stabile lineare Dosis-Wirkungs-Beziehung zwischen der Zufuhr von Vitamin B12 und DKD.
Untergruppenanalysen der Assoziation zwischen Vitamin B12 und DKD
Subgruppenanalysen wurden durchgeführt, um zu untersuchen, ob die Assoziation zwischen der Aufnahme von Vitamin B12 und DKD über verschiedene Teilnehmermerkmale hinweg konsistent war. In Modell 1 wiesen Teilnehmer im höchsten Quartil der Vitamin-B12-Aufnahme (Q4) signifikant niedrigere Chancen für DKD auf als jene in Q1, und zwar sowohl in den Untergruppen der Männer als auch der Frauen (p < 0,05; Tabelle 4). Die Assoziationen über die ACR-Subgruppen hinweg waren heterogener. Bei Teilnehmern mit ACR < 30 mg/g wurden signifikante Assoziationen für einige mittlere Quartile der Vitamin-B12-Aufnahme beobachtet, wohingegen der Muster nach der Stratifikation nach höheren ACR-Werten abgeschwächt war. Über die eGFR-Subgruppen hinweg wurden keine signifikanten Assoziationen beobachtet.
Nichtlineare Beziehung zwischen Vitamin B12 und DKD
Die vollständig adjustierte RCS-Analyse zeigte eine signifikante nichtlineare Assoziation zwischen der Zufuhr von Vitamin B12 und DKD (Abbildung 2, p für Nichtlinearität = 0,037). Im RCS-Modell wurden vier Knoten an den 5., 35., 65. und 95. Perzentilen der Verteilung der Vitamin-B12-Zufuhr verwendet, wobei 3,24 μg/Tag als Referenzwert diente. Eine anschließende Bruchpunktanalyse identifizierte zwei datenbasierte Bruchpunkte bei 3,24 μg/Tag (95 % KI: 2,86–3,62) und 8,36 μg/Tag (95 % KI: 7,52–9,20), die die Vitamin-B12-Zufuhr in drei Bereiche unterteilten. Bei einer Vitamin-B12-Zufuhr < 3,24 μg/Tag wurde keine signifikante Assoziation mit DKD beobachtet (OR = 1,07, p = 0,575). Im Bereich von 3,24–8,36 μg/Tag war die Zufuhr von Vitamin B12 umgekehrt mit DKD assoziiert (OR = 0,87, p = 0,025). Bei einer Zufuhr > 8,36 μg/Tag wurde keine signifikante Assoziation festgestellt (OR = 1,00, p = 0,995) (Tabelle 5). Diese Ergebnisse entsprachen einem nichtlinearen Muster, wobei eine signifikante inverse Assoziation nur im mittleren Zufuhrbereich beobachtet wurde. Um die Robustheit dieser Ergebnisse zu überprüfen, wurden Sensitivitätsanalysen mit alternativen RCS-Spezifikationen mit drei und fünf Knoten durchgeführt. Ähnliche allgemeine Assoziationsmuster wurden beobachtet, wobei unter beiden Knotenspezifikationen statistisch signifikante Gesamtassoziationen in den Modellen 1–3 bestanden; die Modelle mit fünf Knoten stützten darüber hinaus einen nichtlinearen Zusammenhang (Ergänzende Abbildung 1).
Identifizierung von potenziellen vitamin B12-bezogenen Genen bei DKD
Im GSE30528-Datensatz identifizierten wir 493 differentiell exprimierte Gene (DEGs) zwischen DKD- und Normalproben, bestehend aus 126 hochregulierten und 367 herunterregulierten Genen (|log₂FC| > 1, FDR < 0,05; Abbildung 3A,B, Zusätzliche Tabelle 1). Die KEGG-Pfadanalyse zeigte eine signifikante Anreicherung dieser DEGs in den Wegen der Wechselwirkung zwischen extrazellulärer Matrix (ECM) und Rezeptor, fokalen Adhäsion, Integrin-Signaltransduktion und Regulation des Aktin-Zytoskeletts (Abbildung 3C, Zusätzliche Tabelle 2). Die GO-Anreicherungsanalyse ergab zudem, dass die DEGs vorwiegend mit der kollagenhaltigen extrazellulären Matrix, der Glomerulusentwicklung, der fokalen Adhäsion und der Nephronentwicklung assoziiert waren (Abbildung 3D, Zusätzliche Tabelle 2).
Die Schnittmengenanalyse der DEGs und des vitamin-B12-bezogenen Gen-Sets ergab 20 überlappende Kandidatengene (Abbildung 3E). Anschließend wurden diese Gene verwendet, um mithilfe der STRING-Datenbank ein PPI-Netzwerk zu erstellen, wobei die Interaktionen auf einen Mindestvertrauenswert von 0,4 gefiltert wurden. Das resultierende Netzwerk bestand aus 11 Knoten und 20 Kanten, die jeweils Proteine und ihre Interaktionen darstellen (Abbildung 3F). Die Beziehungen zwischen den 20 Kandidatengenen wurden anschließend mit GeneMANIA visualisiert (Abbildung 3G).
Explorative Charakterisierung von P2RY14 als einem potenziellen vitamin B12-assoziierten Gen bei DKD
Basierend auf der beobachteten Assoziation zwischen der Zufuhr von Vitamin B12 und dem DKD-Status suchten wir anschließend mithilfe einer transkriptomischen Analyse nach wesentlichen, mit Vitamin B12 assoziierten Genen, die an der DKD-Pathogenese beteiligt sind. Unter den 20 überschneidenden Kandidatengenen verglichen wir ihre Expressionsmuster zusätzlich in verschiedenen DKD-Datensätzen. In den Datensätzen GSE30528 und GSE30529 war die Expression von P2RY14 in DKD-Proben im Vergleich zu normalen Kontrollen signifikant hochreguliert (Abbildung 4A,B). P2RY14 kodiert einen G-Protein-gekoppelten Rezeptor und wurde bereits mit der Chemotaxis von Immunzellen, der Umorganisation des Cytoskeletts, entzündlichen Reaktionen und oxidativem Stress in Verbindung gebracht28,29. Seine Rolle bei DKD wurde jedoch bisher nicht direkt untersucht. Angesichts der konsistenten Hochregulation in unabhängigen DKD-Datensätzen und seiner potenziellen biologischen Relevanz wurde P2RY14 für nachfolgende Validierungs- und Funktionsanalysen priorisiert. Anschließend wurde eine ROC-Kurvenanalyse durchgeführt, um die Differenzierungsfähigkeit von P2RY14 vorläufig einzuschätzen. Im Entdeckungsdatensatz GSE30528 betrug die AUC 0,923 (95 % KI: 0,807–1,000), während sie im vorläufigen externen Validierungsdatensatz GSE30529 bei 0,883 (95 % KI: 0,744–1,000) lag. Im zusätzlichen Validierungsdatensatz GSE142025 betrug die AUC 0,663 (95 % KI: 0,478–0,847), was auf eine begrenzte Differenzierungsfähigkeit hindeutet (Abbildung 4C–E). Insgesamt liefern diese Ergebnisse erste Hinweise auf den potenziellen Differenzierungswert von P2RY14, wobei die diagnostische Leistungsfähigkeit in größeren, unabhängigen Kohorten weiter validiert werden muss.
Um biologische Funktionen und Signalwege zu untersuchen, die mit P2RY14 Expression in DKD, die GSE30528-Proben wurden aufgeteilt in P2RY14-hoch und P2RY14-niedrigen Gruppen entsprechend dem medianen Expressionsniveau von P2RY14. Eine Differenzialexpressionsanalyse wurde zwischen den beiden Gruppen durchgeführt. Anschließend wurden GO- und KEGG-Anreicherungsanalysen durchgeführt, um die biologischen Funktionen und Signalwege zu charakterisieren, die mit der Expression von P2RY14 assoziiert sind. Die GO-Anreicherungsanalyse zeigte, dass P2RY14 war mit kollagenhaltiger ECM, Glomerularentwicklung, struktureller Bestandteile der ECM und Nierenentwicklung assoziiert (Abbildung 4F, Ergänzungstabelle 3). Die KEGG-Analyse deutete darüber hinaus darauf hin P2RY14 in der Interaktion zwischen ECM und Rezeptor, der fokalen Adhäsion sowie den Komplement- und Gerinnungskaskaden (Abbildung 4G, Ergänzungstabelle 3).
Darüber hinaus zeigte die GSEA, dass der NF-κB-Signalweg und der TNF-Signalweg in der P2RY14-hohen Gruppe signifikant angereichert waren, was darauf hindeutet, dass ein hoher P2RY14-Expressionsgrad mit entzündlichen Signalwegen bei DKD assoziiert sein könnte (Abbildung 4H, Ergänzende Tabelle 3).
Exploratives molekulares Docking von MC-1 mit P2RY14
Fünf Kandidatenverbindungen, die mit P2RY14 assoziiert sind, wurden durch eine Screening-Analyse mit GeneCards und DGIdb identifiziert, und MC-1 wurde für eine explorative molekulare Docking-Analyse ausgewählt (Abbildung 5A). Unter Verwendung von MC-1 als Ligand und P2RY14 als Rezeptor prognostizierte die Docking-Analyse eine mögliche Bindungsposition von MC-1 in einer Kandidaten-Pocket von P2RY14 mit einem Vina-Score von −6,0 kcal/mol. Die vorhergesagte Pose beinhaltete Wechselwirkungen mit Tyr-, Lys-, Asn-, Asp- und Glu-Resten sowie die Bildung von sieben Wasserstoffbrückenbindungen (Abbildung 5B–D).
Datennutzungsrichtlinie
Die in dieser Studie verwendeten Daten sind öffentlich zugänglich. NHANES-Daten aus den Erhebungszyklen 2011–2018 sind beim National Center for Health Statistics erhältlich. Transkriptomische Datensätze (GSE30528, GSE142025 und GSE30529) sind öffentlich im Gene Expression Omnibus (GEO)-Datenbank abrufbar. Die in dieser Studie erstellten verarbeiteten Ergebnisse zur differentiellen Expression, funktionellen Anreicherung und GSEA sind in den ergänzenden Dateien enthalten.

Abbildung 1: Flussdiagramm zur Teilnehmerauswahl für die NHANES-Analyse. Das Flussdiagramm zeigt die schrittweise angewendeten Ausschlusskriterien für Teilnehmer der NHANES-Zyklen 2011–2018 sowie die endgültige Studienpopulation, die in die Analyse eingeschlossen wurde. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 2: Nichtlinearer Zusammenhang zwischen der Aufnahme von Vitamin B12 und diabetischer Nierenerkrankung (DKD) unter Verwendung eines eingeschränkten kubischen Spline-Modells. Die durchgezogene Linie stellt die geschätzte Odds-Ratio für DKD über verschiedene Werte der Ernährungsaufnahme von Vitamin B12 dar, und der schattierte Bereich zeigt das 95-%-Konfidenzintervall an. Das Modell wurde unter Berücksichtigung der in Modell 3 enthaltenen Kovariaten angepasst. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 3: Identifizierung von potenziellen vitamin-B12-bezogenen Genen bei diabetischer Nierenerkrankung (DKD). (A,B) Volcano-Plot und Heatmap zur Darstellung der Expressionsprofile differentiell exprimierter Gene (DEGs) zwischen Proben mit diabetischer Nierenerkrankung (DKD) und normalen Proben im GSE30528-Datensatz. (C,D) Die 15 am stärksten angereicherten KEGG-Wege und GO-Terme der DEGs. (E) Kreuzanalyse zwischen DEGs und vitamin-B12-bezogenen Genen. (F) Protein-Protein-Interaktions-(PPI-)Netzwerk. (G) GeneMANIA-Netzwerk der 20 vitamin-B12‑bezogenen Gene. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 4: Erkundende Charakterisierung von P2RY14 als kandidierendem vitamin B12-assoziierten Gen bei DKD. (A,B) Boxplots, die die Expressionsniveaus der überlappenden Kandidatengene in DKD-Proben und normalen Kontrollen in den Datensätzen GSE30528 und GSE30529 zeigen. (C–E) ROC-Kurvenanalysen (Receiver Operating Characteristic) von P2RY14 in den Datensätzen GSE30528, GSE30529 und GSE142025. (F) Gene Ontology (GO)-Anreicherungsanalyse. (G) Kyoto Encyclopedia of Genes and Genomes (KEGG)-Pfadanreicherungsanalyse. (H) GSEA der P2RY14-hohen und P2RY14-niedrigen Gruppen, die eine Anreicherung der NF-κB- und TNF-Signalwege in GSE30528 zeigt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 5: Arzneimittel-Zielmolekül-Screening und explorative molekulare Docking-Analyse von P2RY14. (A) Netzwerk der Wechselwirkungen zwischen Arzneimitteln und Zielmolekülen für P2RY14-assoziierte Arzneimittel. (B–D) Molekulares Docking von MC-1 mit dem P2RY14-Protein. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Tabelle 1: Grundlegende Merkmale der Studienpopulation. Demografische, klinische, biochemische und Ernährungsmerkmale der Teilnehmer entsprechend dem DKD-Status werden dargestellt. Stetige und kategoriale Variablen werden mithilfe einer analysegewichteten Umfragezusammenfassung zusammengefasst. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 2: Grundlegende Merkmale der Teilnehmer entsprechend der Quartile der Vitamin-B12-Aufnahme. Demografische, klinische, biochemische und Ernährungsmerkmale der Studienpopulation werden über die Quartile der diätetischen Vitamin-B12-Aufnahme dargestellt. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 3: Zusammenhang zwischen der Aufnahme von Vitamin B12 und dem DKD-Status bei Erwachsenen in den USA. Es werden die Odds-Ratios (ORs) und 95 %-Konfidenzintervalle (CIs) für den Zusammenhang zwischen der Ernährungsaufnahme von Vitamin B12 und dem DKD-Status über schrittweise angepasste logistische Regressionsmodelle dargestellt, wobei das niedrigste Quartil der Vitamin-B12-Aufnahme (Q1) als Referenzgruppe dient. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 4: Zusammenhang zwischen der Aufnahme von Vitamin B12 und dem DKD-Status in Subgruppenanalysen. Dargestellt sind die Odds-Ratios (ORs) und 95 %-Konfidenzintervalle (CIs) für DKD über die Quartile der Vitamin-B12-Aufnahme hinweg, stratifiziert nach Alter, Geschlecht, BMI, eGFR, ACR, Hypertonie und Hyperlipidämie-Subgruppen, und adjustiert für die angegebenen Modelle. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 5: Dreigliedrige stückweise Regressionsanalyse der Assoziation zwischen Vitamin-B12-Zufuhr und DKD. Die Assoziationen zwischen der diätetischen Zufuhr von Vitamin B12 und DKD werden für die drei durch die zwei datenbasierten Bruchpunkte bei 3,24 und 8,36 μg/Tag definierten Zufuhrbereiche dargestellt. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Zusatzabbildung 1: Sensitivitätsanalysen der Assoziation zwischen diätetischem Vitamin B12 Aufnahme und diabetische Nierenerkrankung (DKD) mithilfe von Modellen mit eingeschränkten kubischen Splines (A–C) Modelle mit drei Knoten. (D–F) Modelle mit fünf Knoten. Die durchgezogene Linie stellt die geschätzte OR dar, der schattierte Bereich das 95%-KI und die gestrichelte Linie zeigt OR = 1 an. p Werte für die Gesamtassoziation und Nichtlinearität sind in jedem Diagramm dargestellt.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Zusatz-Tabelle 1: Differenziell exprimierte Gene zwischen diabetischer Nierenerkrankung (DKD) und normalen Proben im GSE30528-Datensatz. Ergebnisse zur differentiellen Expression von Genen mit einem |log₂FC| > 1 und einem FDR < 0,05 werden zusammen mit dem Cobalamin-bezogenen Gen-Satz, der für die Kandidatengensuche verwendet wurde, dargestellt.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 2: Funktionsanreicherungsanalyse der differentiell exprimierten Gene bei DKD. Dargestellt sind die Ergebnisse der Anreicherungsanalyse für die genetische Ontologie (GO) und die Kyoto-Enzyklopädie für Gene und Genome (KEGG) für die in dem Datensatz GSE30528 identifizierten differentiell exprimierten Gene.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Zusatz tabelle 3: Ergebnisse der funktionalen Anreicherung und der Gen-Set-Anreicherungsanalyse in Verbindung mit P2RY14 im Datensatz GSE30528. Dargestellt sind GO-, KEGG- und GSEA-Ergebnisse, die zur Charakterisierung der mit P2RY14 assoziierten biologischen Funktionen und Signalwege verwendet wurden.Bitte klicken Sie hier, um diese Datei herunterzuladen.