Progettazione dello studio
Questa ricerca ha utilizzato un approccio MR a due campioni per esplorare la relazione causale tra tratti infiammatori associati alle NETs e il diabete di tipo 2 (T2DM) e le sue complicanze. Il modello MR si basa su tre assunzioni fondamentali20: 1) la variante genetica deve essere significativamente associata all'esposizione, 2) la variante genetica deve influenzare l'esito unicamente attraverso l'esposizione (esclusione della pleiotropia) e 3) la variante genetica non deve essere associata a fattori di confondimento. Per una panoramica dettagliata del processo dello studio, si prega di consultare Figura 1. Il manoscritto rispetta le linee guida MR-STROBE per la presentazione degli studi MR, garantendo standard rigorosi di trasparenza e riproducibilità21.
Fonti dei dati
I dati delle analisi di associazione sull'intero genoma (GWAS) utilizzati in questa analisi MR provengono da database pubblici GWAS. Per quanto riguarda gli studi GWAS di livello riassuntivo relativi al diabete di tipo 2 (T2DM), sono stati inclusi la più ampia meta-analisi su una popolazione europea, composta da 74.124 casi e 824.006 controlli, insieme a uno studio prospettico caso-coorte annidato condotto in Europa, comprendente 9.978 casi e 12.348 controlli22. Per le complicanze diabetiche, tra cui T2DM con complicanze renali e T2DM con complicanze circolatorie periferiche, questo studio ha recuperato le statistiche riassuntive GWAS da studi in cui i casi erano pazienti con T2DM affetti dalla specifica complicanza e i controlli erano individui senza T2DM. Informazioni dettagliate sulle fonti dei dati sono disponibili nella Tabella Supplementare 1.
I dati GWAS per i tratti infiammatori associati alle NET sono stati ottenuti dal catalogo GWAS (vedere Tabella Supplementare 2). Per affrontare l'eterogeneità nella specificità biologica, questo studio ha suddiviso le esposizioni incluse in due gruppi distinti in base ai loro ruoli funzionali nella biologia dei neutrofili: (1) Fattori principali associati alla NETosi: questa categoria comprende marcatori direttamente coinvolti nella formazione strutturale delle NET o nel processo enzimatico della decondensazione della cromatina. In particolare, questo studio ha incluso le NET stesse, la mieloperossidasi (MPO), l'elastasi neutrofila (NE) e il complesso MPO-DNA. MPO e NE sono enzimi essenziali per la degradazione delle istoni e la decondensazione della cromatina, mentre il complesso MPO-DNA è un marcatore surrogato specifico per le NET23. (2) Mediatori infiammatori associati alle NET: questa categoria comprende citochine e mediatori che agiscono come regolatori iniziali o effettori finali strettamente legati alla NETosi, ma coinvolti anche in vie infiammatorie più ampie. Questo gruppo include l'interleuchina-6 (IL-6)24, il fattore di necrosi tumorale-alfa (TNF-α)25, i livelli di lipocalina associata alla gelatinasi dei neutrofili (NGAL)26 e il fattore della rete di comunicazione cellulare 1 (CCN1)27.
I dati utilizzati in questo studio sono stati ottenuti da database ad accesso aperto o da studi GWAS pubblicati in ricerche precedenti; pertanto, non è stata richiesta un'approvazione etica per questa ricerca.
Selezione di strumenti genetici associati ai NET
Poiché la formazione delle NET rappresenta un processo biologico dinamico non direttamente misurato negli studi GWAS convenzionali, questo studio ha adottato un approccio basato sui geni per identificare strumenti genetici associati a tratti infiammatori legati alle NET. È stato selezionato un insieme completo di 257 geni noti essere crucialmente coinvolti nella formazione e regolazione delle NET, come MPO, insieme ai geni che codificano per le istoni e le proteine dei granuli neutrofili, sulla base di studi meccanicistici consolidati e della letteratura pubblicata28.
Selezione della variabile strumentale
Per identificare strumenti genetici robusti per i tratti infiammatori associati alle NETs, questo studio ha implementato un flusso di lavoro sequenziale di filtraggio utilizzando statistiche riassuntive di GWAS. I passaggi operativi specifici sono stati i seguenti:
Screening iniziale: sono state estratte le varianti genetiche significativamente associate ai tratti di esposizione29 in base a una soglia di significatività di P < 5 × 10⁻6.
Filtraggio per frequenza dell'allele minore (MAF): per garantire un'adeguata potenza statistica, sono state escluse le polimorfismi a singolo nucleotide (SNPs) con una MAF ≤ 0,0130.
Raggruppamento per squilibrio di associazione (LD): per ridurre il confondimento dovuto allo squilibrio di associazione (LD), sono stati selezionati SNP indipendenti utilizzando la funzione di raggruppamento con parametri impostati su r2 < 0,001 all'interno di una finestra di 10.000 kb31.
Valutazione della forza dello strumento: La forza di ciascuna IV rimanente è stata quantificata utilizzando la statistica F, calcolata come F = R2 × (N-2) / (1-R2). Sono stati mantenuti solo gli SNP con una statistica F > 10 per ridurre al minimo il bias dovuto a strumenti deboli32.
Identificazione di SNP proxy e armonizzazione dei dati
Per affrontare la mancanza di SNP nel set di dati dell'esito, questo studio ha effettuato una ricerca di proxy e un processo di armonizzazione dei dati:
Sostituzione con proxy: quando un SNP bersaglio proveniente dal GWAS dell'esposizione non era disponibile nel GWAS dell'esito, questo studio ha utilizzato la funzione LDproxy() basata sul panel di riferimento europeo del Progetto 1,000 Genomi. Un SNP proxy candidato è stato selezionato solo se mostrava un'elevata associazione di linkage disequilibrium con il SNP originale (r2> 0,8). Se non veniva trovato alcun proxy idoneo, il SNP veniva escluso.
Armonizzazione: i set di dati relativi all'esposizione e all'esito sono stati allineati utilizzando la funzione harmonise_data() del pacchetto TwoSampleMR (R versione 4.0.5). In questo studio è stato impostato il parametro action = 2 per allineare automaticamente tutti gli SNP al filamento diretto e rimuovere gli SNP palindromi con orientamento del filamento ambiguo.
Verifica: dopo l'armonizzazione, questo studio ha esaminato manualmente il dataset armonizzato per confermare che le frequenze degli alleli effetto fossero coerenti tra i dati di esposizione e quelli di esito.
Analisi MR e test di sensibilità
L'inferenza causale è stata condotta utilizzando il pacchetto TwoSampleMR (R versione 4.0.5).
Analisi primarie e secondarie: questo studio ha applicato il metodo del peso inverso della varianza (IVW) come approccio principale33. Sono state effettuate analisi complementari utilizzando i metodi MR-Egger34, mediana pesata e moda pesata per garantire la robustezza35.
Valutazione della sensibilità: l'eterogeneità tra le VI è stata valutata mediante il test Q di Cochran36 tramite la funzione mr_heterogeneity(). La pleiotropia orizzontale è stata valutata utilizzando il test dell'intercetta MR-Egger34 (mr_pleiotropy_test()).
Rilevamento degli outlier: questo studio ha utilizzato il pacchetto MR-PRESSO per rilevare possibili outlier37. La funzione mr_presso() è stata eseguita con 1.000 simulazioni. Gli SNP outlier identificati (P < 0,05) sono stati rimossi e le stime causali sono state ricalcolate per verificare la stabilità dei risultati. Inoltre, è stata effettuata un'analisi leave-one-out per assicurare che l'associazione causale non fosse determinata da un singolo SNP38.
Correzione Statistica
Per tener conto dei test multipli, i valori P ottenuti dalle analisi MR sono stati corretti utilizzando il metodo del False Discovery Rate (FDR). Questo è stato applicato mediante la funzione P.adjust() in R con il parametro method = "fdr". Le associazioni con un valore P corretto (PFDR) < 0,05 sono state considerate statisticamente significative.