Características basais
Um total de 8.218 participantes do NHANES 2011–2018 foi incluído, dos quais 623 foram classificados como tendo DRC e 7.595 como Não-DRC. Diferenças significativas entre os grupos DRC e Não-DRC foram observadas em idade, sexo, raça/etnia, estado civil, escolaridade, PIR, IMC, pressão arterial, hipertensão, consumo de álcool, tabagismo, ingestão de vitamina B12 e ingestão alimentar (todos p < 0,05). Diferenças significativas também foram observadas em ACR, TFG, Scr, BUN, CT e TG (Tabela 1).
Foram observadas diferenças significativas entre os quartis de ingestão de vitamina B12 em relação ao nível educacional, hipertensão, estado de DNT e ingestão energética. Em contraste, nenhuma diferença significativa foi observada entre os quartis em relação à idade, estado civil, IMC, pressão arterial, consumo de álcool, TFG, RAC, ácido úrico sérico, Scr ou TG (Table 2).
Associação entre vitamina B12 e DNT
A ingestão de vitamina B12 foi categorizada em quartis (Q1–Q4), e os resultados das análises de regressão logística multivariável são apresentados na Tabela 3. No Modelo 3 completamente ajustado, os participantes no quartil mais alto de ingestão de vitamina B12 (Q4) apresentaram odds significativamente menores de DNT do que aqueles no quartil mais baixo (Q1) (OR = 0,70, IC 95%: 0,50–0,97, p = 0,034), enquanto nenhuma associação significativa foi observada para Q2 (OR = 1,14, IC 95%: 0,85–1,53, p = 0,386) ou Q3 (OR = 0,94, IC 95%: 0,69–1,28, p = 0,684). Além disso, a tendência linear no modelo completamente ajustado não atingiu significância estatística (p para tendência = 0,086), embora a direção geral da associação fosse inversa. Esses achados sugerem que uma ingestão mais alta de vitamina B12, particularmente no quartil mais elevado, esteve associada a menores odds de DNT, mas os resultados não apoiam uma relação dose-resposta linear estável entre a ingestão de vitamina B12 e a DNT.
Análises de subgrupo da associação entre vitamina B12 e DKD
Análises de subgrupos foram realizadas para verificar se a associação entre a ingestão de vitamina B12 e DNT era consistente entre diferentes características dos participantes. No Modelo 1, participantes no quartil mais alto de ingestão de vitamina B12 (Q4) apresentaram odds significativamente menores de DNT do que aqueles no Q1, tanto no subgrupo de homens quanto no de mulheres (p < 0,05; Tabela 4). As associações nos subgrupos de TRA foram mais heterogêneas. Entre participantes com TRA < 30 mg/g, associações significativas foram observadas em alguns quartis intermediários de ingestão de vitamina B12, enquanto o padrão foi atenuado após estratificação por níveis mais altos de TRA. Nenhuma associação significativa foi observada nos subgrupos de TFGe.
Relação não linear entre a vitamina B12 e a DNT
A análise de RCS completamente ajustada mostrou uma associação não linear significativa entre a ingestão de vitamina B12 e DTD (Figura 2, p para não linearidade = 0,037). O modelo RCS utilizou quatro nós localizados nos percentis 5, 35, 65 e 95 da distribuição da ingestão de vitamina B12, com 3,24 μg/dia usado como valor de referência. A análise subsequente de pontos de quebra identificou dois pontos de quebra baseados nos dados em 3,24 μg/dia (IC 95%: 2,86–3,62) e 8,36 μg/dia (IC 95%: 7,52–9,20), dividindo a ingestão de vitamina B12 em três faixas. Para ingestão de vitamina B12 < 3,24 μg/dia, nenhuma associação significativa com DTD foi observada (OR = 1,07, p = 0,575). Na faixa de 3,24–8,36 μg/dia, a ingestão de vitamina B12 foi inversamente associada à DTD (OR = 0,87, p = 0,025). Para ingestão > 8,36 μg/dia, nenhuma associação significativa foi observada (OR = 1,00, p = 0,995) (Tabela 5). Esses achados foram consistentes com um padrão não linear, com uma associação inversa significativa observada apenas na faixa intermediária de ingestão. Para avaliar a robustez desses achados, foram realizadas análises de sensibilidade utilizando especificações alternativas de RCS com três e cinco nós. Padrões de associação gerais semelhantes foram observados, com associações gerais estatisticamente significativas em todos os Modelos 1–3 sob ambas as especificações de nós, enquanto os modelos com cinco nós forneceram apoio adicional a uma relação não linear (Figura Suplementar 1).
Identificação de genes candidatos relacionados à vitamina B12 na DNT
No conjunto de dados GSE30528, identificamos 493 genes diferencialmente expressos (GDEs) entre amostras de DKD e normais, compreendendo 126 genes superexpressos e 367 genes subexpressos (|log₂FC| > 1, FDR < 0,05; Figura 3A,B, Tabela Suplementar 1). A análise de vias KEGG revelou enriquecimento significativo desses GDEs nas vias de interação entre matriz extracelular (MEC) e receptor, adesão focal, sinalização por integrina e regulação do citoesqueleto de actina (Figura 3C, Tabela Suplementar 2). A análise de enriquecimento de GO indicou ainda que os GDEs estavam predominantemente associados à matriz extracelular contendo colágeno, desenvolvimento do glomérulo, adesão focal e desenvolvimento do néfron (Figura 3D, Tabela Suplementar 2).
A análise de interseção entre os genes diferencialmente expressos (DEGs) e o conjunto de genes relacionados à vitamina B12 identificou 20 genes candidatos com sobreposição (Figura 3E). Esses genes foram posteriormente utilizados para construir uma rede de interação proteína-proteína (PPI) utilizando o banco de dados STRING, com as interações filtradas por um escore mínimo de confiança de 0,4. A rede resultante consistiu em 11 nós e 20 arestas, representando proteínas e suas interações, respectivamente (Figura 3F). As relações entre os 20 genes candidatos foram ainda visualizadas utilizando o GeneMANIA (Figura 3G).
Caracterização exploratória de P2RY14 como um gene candidato relacionado à vitamina B12 na DN
Com base na associação observada entre a ingestão de vitamina B12 e o estado de DNT, em seguida buscamos identificar genes-chave relacionados à vitamina B12 implicados na patogênese da DNT por meio de análise transcriptômica. Dentre os 20 genes candidatos sobrepostos, comparamos ainda seus padrões de expressão em diferentes conjuntos de dados de DNT. Nos conjuntos de dados GSE30528 e GSE30529, a expressão de P2RY14 foi significativamente aumentada em amostras de DNT em comparação com controles normais (Figura 4A,B). P2RY14 codifica um receptor acoplado à proteína G e já foi relatado como envolvido na quimiotaxia de células imunes, reorganização do citoesqueleto, respostas inflamatórias e estresse oxidativo28,29. No entanto, seu papel na DNT ainda não foi diretamente investigado. Considerando sua regulação positiva consistente em diferentes conjuntos de dados independentes de DNT e sua relevância biológica potencial, P2RY14 foi priorizado para validação subsequente e análises funcionais. A análise da curva ROC foi posteriormente realizada para avaliar preliminarmente a capacidade discriminatória de P2RY14. No conjunto de dados de descoberta GSE30528, a AUC foi de 0,923 (IC 95%: 0,807–1,000), enquanto no conjunto de dados de validação externa preliminar GSE30529, a AUC foi de 0,883 (IC 95%: 0,744–1,000). No conjunto de dados adicional de validação GSE142025, a AUC foi de 0,663 (IC 95%: 0,478–0,847), indicando capacidade discriminatória limitada (Figura 4C–E). Em conjunto, esses achados fornecem evidências preliminares sobre o valor discriminatório potencial de P2RY14, embora seu desempenho diagnóstico exija validação adicional em coortes independentes maiores.
Para explorar funções biológicas e vias associadas a P2RY14 expressão em DKD, as amostras de GSE30528 foram divididas em P2RY14-alto e P2RY14-baixos de acordo com o nível mediano de expressão de P2RY14. A análise de expressão diferencial foi realizada entre os dois grupos. Em seguida, foram realizadas análises de enriquecimento GO e KEGG para caracterizar as funções biológicas e vias associadas à expressão de P2RY14. A análise de enriquecimento GO revelou que P2RY14 foi associado ao ECM contendo colágeno, desenvolvimento do glomérulo, componente estrutural do ECM e desenvolvimento renal (Figura 4F, Tabela Suplementar 3). A análise KEGG implicou ainda P2RY14 na interação entre matriz extracelular e receptor, adesão focal e cascatas do complemento e da coagulação (Figura 4G, Tabela Suplementar 3).
Além disso, a GSEA mostrou que a via de sinalização de NF-κB e a via de sinalização de TNF estavam significativamente enriquecidas no grupo com alta expressão de P2RY14 , sugerindo que a alta expressão de P2RY14 pode estar associada à sinalização inflamatória na DNT (Figura 4H, Tabela Suplementar 3).
Ensaio exploratório de docking molecular de MC-1 com P2RY14
Cinco compostos candidatos associados ao P2RY14 foram identificados por meio da triagem no GeneCards e no DGIdb, e MC-1 foi selecionado para análise exploratória de docking molecular (Figura 5A). Utilizando MC-1 como ligante e P2RY14 como receptor, a análise de docking previu uma possível conformação de ligação de MC-1 dentro de um sítio candidato do P2RY14, com um escore Vina de −6,0 kcal/mol. A conformação prevista envolveu interações com resíduos de Tyr, Lys, Asn, Asp e Glu e a formação de sete ligações de hidrogênio (Figura 5B–D).
Declaração sobre Disponibilidade de Dados
Os dados utilizados neste estudo estão disponíveis publicamente. Os dados do NHANES dos ciclos de pesquisa de 2011–2018 estão disponíveis no National Center for Health Statistics. Os conjuntos de dados transcriptômicos GSE30528, GSE142025 e GSE30529 estão disponíveis publicamente no banco de dados Gene Expression Omnibus (GEO). Os resultados processados de expressão diferencial, enriquecimento funcional e GSEA gerados neste estudo são fornecidos nos Arquivos Suplementares.

Figura 1: Fluxograma da seleção de participantes para a análise do NHANES. O fluxograma mostra os critérios de exclusão sequenciais aplicados aos participantes dos ciclos do NHANES de 2011 a 2018 e a população final do estudo incluída na análise. Clique aqui para visualizar uma versão maior desta figura.

Figura 2: Associação não linear entre a ingestão de vitamina B12 e DNT utilizando um modelo de spline cúbico restrito. A linha sólida representa a razão de chances estimada para DNT ao longo da ingestão dietética de vitamina B12, e a área sombreada representa o intervalo de confiança de 95%. O modelo foi ajustado utilizando as covariáveis incluídas no Modelo 3. Clique aqui para visualizar uma versão maior desta figura.

Figura 3: Identificação de genes candidatos relacionados à vitamina B12 na doença renal diabética (DKD). (A,B) Gráfico de dispersão e mapa de calor mostrando os perfis de expressão dos genes diferencialmente expressos (DEGs) entre amostras de doença renal diabética (DKD) e amostras normais no conjunto de dados GSE30528. (C,D) As 15 vias KEGG e termos GO mais enriquecidos entre os DEGs. (E) Análise cruzada entre DEGs e genes relacionados à vitamina B12. (F) Rede de interação proteína-proteína (PPI). (G) Rede GeneMANIA dos 20 genes relacionados à vitamina B12. Clique aqui para visualizar uma versão maior desta figura.

Figura 4: Caracterização exploratória de P2RY14 como gene candidato relacionado à vitamina B12 na DN. (A,B) Gráficos de caixa mostrando os níveis de expressão dos genes candidatos sobrepostos em amostras de DN e controles normais nos conjuntos de dados GSE30528 e GSE30529, respectivamente. (C–E) Análises da curva ROC (característica operacional do receptor) de P2RY14 nos conjuntos de dados GSE30528, GSE30529 e GSE142025, respectivamente. (F) Análise de enriquecimento da Ontologia Genética (GO). (G) Análise de enriquecimento de vias do Kyoto Encyclopedia of Genes and Genomes (KEGG). (H) GSEA dos grupos com alta e baixa expressão de P2RY14, mostrando enriquecimento das vias de sinalização NF-κB e TNF em GSE30528. Clique aqui para visualizar uma versão maior desta figura.

Figura 5: Triagem de fármaco-alvo e análise exploratória de encaixe molecular de P2RY14. (A) Rede de interação entre fármacos e alvos dos medicamentos associados ao P2RY14. (B–D) Encaixe molecular da MC-1 com a proteína P2RY14. Clique aqui para visualizar uma versão maior desta figura.
Tabela 1: Características básicas da população do estudo. Apresentam-se as características demográficas, clínicas, bioquímicas e dietéticas dos participantes de acordo com o status de DNT. Variáveis contínuas e categóricas são resumidas utilizando análise ponderada por inquérito. Clique aqui para baixar esta tabela.
Tabela 2: Características básicas dos participantes de acordo com os quartis de ingestão de vitamina B12. As características demográficas, clínicas, bioquímicas e dietéticas da população do estudo são apresentadas conforme os quartis de ingestão dietética de vitamina B12. Clique aqui para baixar esta tabela.
Tabela 3: Associação entre a ingestão de vitamina B12 e o estado de DNT em adultos dos Estados Unidos. São apresentadas as razões de chances (ORs) e os intervalos de confiança de 95% (ICs) para a associação entre a ingestão dietética de vitamina B12 e o estado de DNT em modelos de regressão logística progressivamente ajustados, sendo o quartil mais baixo de ingestão de vitamina B12 (Q1) o grupo de referência. Clique aqui para baixar esta tabela.
Tabela 4: Associação entre a ingestão de vitamina B12 e o estado de DNT em análises de subgrupos. São apresentadas as razões de chances (ORs) e os intervalos de confiança de 95% (ICs) para DNT nos quartis de ingestão de vitamina B12, estratificados por subgrupos de idade, sexo, IMC, TFGe, RCA, hipertensão e hiperlipidemia, e ajustados para os modelos especificados. Clique aqui para baixar esta tabela.
Tabela 5: Análise de regressão segmentada em três partes da associação entre a ingestão de vitamina B12 e DNT. As associações entre a ingestão dietética de vitamina B12 e DNT são apresentadas para as três faixas de ingestão definidas pelos dois pontos de corte determinados pelos dados, em 3,24 e 8,36 μg/dia. Clique aqui para baixar esta tabela.
Figura Suplementar 1: Análises de sensibilidade da associação entre a ingestão dietética de vitamina B12 e DNT utilizando modelos de spline cúbico restrito. (A–C) Modelos com três nós. (D–F) Modelos com cinco nós. A linha sólida representa a OR estimada, a área sombreada representa o IC de 95% e a linha tracejada indica OR = 1. Valores de p para a associação geral e não linearidade são apresentados em cada painel.Clique aqui para baixar este arquivo.
Tabela Suplementar 1: Genes diferencialmente expressos entre doenças renais diabéticas (DKD) e amostras normais no conjunto de dados GSE30528. Os resultados de expressão diferencial para genes identificados com |log₂FC| > 1 e FDR < 0,05 são apresentados juntamente com o conjunto de genes relacionados à cobalamina utilizado para triagem de genes candidatos.Clique aqui para baixar este arquivo.
Tabela Suplementar 2: Análise de enriquecimento funcional dos genes diferencialmente expressos na DN. Apresentam-se os resultados do enriquecimento da Ontologia Genética (GO) e da Enciclopédia de Quioto de Genes e Genomas (KEGG) para os genes diferencialmente expressos identificados no conjunto de dados GSE30528.Clique aqui para baixar este arquivo.
Tabela Suplementar 3: Resultados da análise de enriquecimento funcional e da análise de enriquecimento de conjuntos de genes associados ao P2RY14 no conjunto de dados GSE30528. Apresentam-se os resultados do GO, KEGG e GSEA utilizados para caracterizar as funções biológicas e vias associadas ao P2RY14.Clique aqui para baixar este arquivo.