Figura 1 riassume il flusso di lavoro complessivo del protocollo di randomizzazione mendeliana a due campioni, dalla selezione dei dataset e dallo screening degli strumenti fino alle analisi di sensibilità e alla visualizzazione dei risultati. Di seguito sono presentati gli output rappresentativi di ciascuna fase.
Variabili strumentali genetiche e controllo di qualità dell'armonizzazione
Un totale di 52 SNP associati all'IS sono stati inizialmente estratti dal dataset GWAS dell'esposizione utilizzando una soglia di significatività di p < 5 × 10-6. Le associazioni con l'outcome sono state quindi estratte dal dataset GWAS della demenza vascolare (VaD). Dopo l'abbinamento dell'outcome, la ricerca di SNP proxy e l'armonizzazione degli alleli, sono stati mantenuti 51 SNP per l'analisi finale di MR. Due varianti sono state abbinate utilizzando SNP proxy. Gli strumenti genetici mantenuti sono elencati nella Tabella 1. I valori di F-statistica dei 51 strumenti mantenuti variavano da 20,91 a 46,06, con una mediana di 25,08 e un IQR di 22,95–30,71. Nessuno dei SNP mantenuti presentava un valore di F-statistica inferiore a 10, indicando una bassa probabilità di bias dovuto a strumenti deboli.
Stime MR per IS e VaD
Il diagramma di dispersione ha mostrato la direzione e l'entità delle stime specifiche per SNP ottenute con diversi metodi di MR (Figura 2). I metodi IVW e della mediana pesata hanno evidenziato associazioni positive, mentre la stima ottenuta con MR-Egger era coerente nella direzione ma non ha raggiunto la significatività statistica.
Il grafico a foreste delle stime individuali dei SNP ha mostrato le associazioni specifiche per SNP con il rischio di demenza vascolare (VaD) (Figura 3). Le stime complessive di MR sono riassunte nella Figura 4. Il metodo IVW ha evidenziato un'associazione positiva statisticamente significativa tra ictus ischemico (IS) predetto geneticamente e rischio di VaD (OR = 1.63, 95% CI: 1.21–2.21, p = 0.0013). Il metodo della mediana pesata ha prodotto una stima significativa concorde (OR = 1.60, 95% CI: 1.06–2.42, p = 0.0240). La stima MR-Egger era direzionalmente coerente ma non ha raggiunto la significatività statistica (OR = 2.21, 95% CI: 0.97–5.01, p = 0.0645).
Lo studio GWAS sull'esito di VaD di FinnGen includeva 881 casi e 211.508 controlli. In un calcolo post hoc del test di Wald bilaterale condizionato alla stima osservata IVW, la potenza stimata per l'analisi primaria IVW era dell'89,5% con α = 0,05. L'effetto minimo rilevabile corrispondente per una potenza dell'80% era un OR di 1,54. Al contrario, l'effetto minimo rilevabile per l'analisi MR-Egger era un OR di 3,23, valore superiore alla stima MR-Egger osservata con direzione concorde (OR = 2,21). Pertanto, il risultato non significativo dell'analisi MR-Egger deve essere interpretato come riflesso della limitata precisione di questo estimatore di sensibilità meno efficiente, piuttosto che, di per sé, come evidenza contraria alla direzione della stima primaria IVW.
Poiché è stata valutata una sola associazione prespecificata tra esposizione e outcome, la soglia corretta con il metodo di Bonferroni per l'inferenza primaria di MR era pari a 0,05/1 = 0,05. Pertanto, il risultato IVW ha soddisfatto la soglia di significatività corretta. Nel complesso, questi risultati forniscono evidenze suggestive di un possibile effetto positivo dell'IS predetta geneticamente sul rischio di demenza vascolare (VaD), principalmente basate sulla stima IVW e sostenute dall'analisi di sensibilità per la mediana pesata. Tuttavia, la stima MR-Egger era coerente nella direzione, ma non ha raggiunto la significatività statistica; pertanto, i risultati non devono essere interpretati come evidenza definitiva di causalità.
Analisi dell'eterogeneità e della pleiotropia orizzontale
L'eterogeneità tra SNP è stata valutata utilizzando la statistica Q di Cochran. Il test di eterogeneità IVW ha prodotto un valore Q = 57,46 con 50 gradi di libertà (p = 0,218), mentre il test di eterogeneità MR-Egger ha prodotto un valore Q = 56,77 con 49 gradi di libertà (p = 0,208). Questi risultati non hanno indicato un'eterogeneità sostanziale tra le stime specifiche per SNP. La pleiotropia orizzontale direzionale è stata valutata mediante il test dell'intercetta MR-Egger. L'intercetta MR-Egger è risultata pari a -0,0195 (SE = 0,0253, p = 0,444), indicando l'assenza di evidenza statistica di pleiotropia orizzontale direzionale. L'analisi MR-PRESSO è stata condotta utilizzando i 51 SNP armonizzati conservati, con 10.000 simulazioni. Il test globale MR-PRESSO non ha mostrato evidenza di pleiotropia orizzontale globale (RSSobs = 59,61; p empirico = 0,2388). Poiché il test globale non è risultato statisticamente significativo, i test individuali di outlier e di distorsione non erano applicabili e non è stata generata alcuna stima corretta per gli outlier. Il grafico a imbuto ha fornito una valutazione visiva della simmetria delle stime specifiche per SNP (Figura 5). Le impostazioni e i risultati dettagliati di MR-PRESSO sono riportati in Tabella Supplementare 1.
Analisi di sensibilità leave-one-out e convalida dell'output del protocollo
È stata eseguita un'analisi di tipo leave-one-out per valutare se la stima complessiva di MR fosse influenzata da un singolo SNP. Il grafico leave-one-out ha mostrato che la rimozione sequenziale di singoli SNP non modificava sostanzialmente la stima complessiva (Figura 6), suggerendo che nessuno strumento genetico singolo dominasse l'associazione. Nel complesso, questi risultati rappresentativi dimostrano l'output pratico del flusso di lavoro di MR a due campioni descritto in questo protocollo. Tabella 1 mostra gli strumenti genetici conservati dopo la selezione e l'armonizzazione degli SNP. Figura 2 illustra la direzione delle stime di MR attraverso i diversi metodi, Figura 3 presenta le stime specifiche per SNP, Figura 4 riassume le stime complessive di MR, Figura 5 supporta la valutazione visiva della pleiotropia o asimmetria, e Figura 6 valuta l'influenza dei singoli SNP. Tabella Supplementare 2 collega ogni passaggio principale del protocollo al relativo output di validazione.
Tutte le statistiche riassuntive grezze GWAS analizzate in questo studio sono disponibili pubblicamente. Il dataset relativo all'esposizione all'ictus ischemico è stato ottenuto dal database IEU OpenGWAS con l'ID del dataset ebi-a-GCST90018864. Il dataset relativo all'esito demenza vascolare è stato ottenuto dal dataset GWAS di FinnGen con l'ID del dataset finn-b-F5_VASCDEM. La tabella delle variabili strumentali estratte, il dataset dell'analisi armonizzata, le tabelle dei risultati MR e i file di output pertinenti sono forniti come file supplementari.

Figura 1. Flusso di lavoro del protocollo MR a due campioni. Questa figura riassume i passaggi principali del protocollo, inclusa la selezione del set di dati dell'esposizione, la selezione delle variabili strumentali, l'estrazione dei dati dell'esito, l'armonizzazione degli alleli, la stima MR, il test di eterogeneità, la valutazione della pleiotropia orizzontale, l'analisi leave-one-out e la visualizzazione. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 2. Grafico a dispersione delle stime di MR per l'associazione tra IS e VaD. Ogni punto rappresenta una stima specifica per SNP. Le linee tracciate rappresentano l'associazione stimata ottenuta utilizzando diversi metodi di MR. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3. Grafico a foreste delle stime MR per singolo SNP. Questa figura mostra le stime specifiche per singolo SNP relative all'associazione tra varianti genetiche associate all'ictus e il rischio di demenza vascolare. Le linee orizzontali rappresentano intervalli di confidenza al 95%. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 4. Forest plot delle stime complessive di MR ottenute con diversi metodi. Questa figura riassume le stime complessive di MR ottenute utilizzando i metodi IVW, mediana pesata e MR-Egger. Le linee orizzontali rappresentano intervalli di confidenza al 95%. MR, randomizzazione mendeliana; IVW, pesato per la varianza inversa. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 5. Grafico a imbuto delle stime di MR specifiche per SNP. Questa figura mostra la distribuzione delle stime di MR specifiche per SNP e consente una valutazione visiva della simmetria tra gli strumenti genetici. MR, randomizzazione mendeliana; IVW, ponderato per la varianza inversa. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 6. Analisi di sensibilità leave-one-out. Questa figura mostra la stima MR dopo la rimozione sequenziale di ciascun SNP. Il grafico è stato utilizzato per valutare se la stima complessiva fosse influenzata da un singolo strumento genetico. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| N. | SNP | Gene | Chr. | EA | OA | EAF.IS | EAF.VD | IS β (SE) | VD β (EE) | statistica F | R² |
| 1 | rs10886430 | GRK5 | 10 | G | A | 0.118244 | 0.09536 | 0.1147 (0.0220) | -0.0172 (0.0852) | 27.18 | 5,61E-05 |
| 2 | rs10936572 | LOC107986051 | 3 | T | C | 0.179643 | 0.1235 | -0.0512 (0.0104) | 0.1444 (0.0743) | 24.24 | 5,01E-05 |
| 3 | rs11045239 | PDE3A | 12 | A | G | 0.492943 | 0.4629 | 0.0604 (0.0089) | 0.0303 (0.0494) | 46.06 | 9,51E-05 |
| 4 | rs11047532 | LOC105369698 | 12 | G | C | 0.305214 | 0.1513 | 0.0499 (0.0105) | 0.0159 (0.0689) | 22.59 | 4,66E-05 |
| 5 | rs11065836 | CUX2 | 12 | A | G | 0.159123 | 0.07442 | -0.0608 (0.0103) | 0.0245 (0.0929) | 34.84 | 7,20E-05 |
| 6 | rs11105378 | ATP2B1 | 12 | T | C | 0.203466 | 0.07522 | -0.0540 (0.0098) | -0.1174 (0.0916) | 30.36 | 6,27E-05 |
| 7 | rs117140252 | - | 14 | A | G | 0.041216 | 0.06988 | 0.1568 (0.0311) | 0.1498 (0.0967) | 25.42 | 5,25E-05 |
| 8 | rs117343276 | - | 10 | G | A | 0.030533 | 0.0431 | -0.0905 (0.0196) | -0.1449 (0.1233) | 21.32 | 4,40E-05 |
| 9 | rs11831940 | HDAC7 | 12 | A | G | 0.251914 | 0.3526 | 0.0631 (0.0130) | 0.1226 (0.0513) | 23.56 | 4,87E-05 |
| 10 | rs11880613 | DNM2 | 19 | A | G | 0.184857 | 0.1804 | -0.0613 (0.0110) | 0.1499 (0.0640) | 31.06 | 6,41E-05 |
| 11 | rs12445022 | LOC124903748 | 16 | A | G | 0.271667 | 0.3193 | 0.0605 (0.0112) | 0.0218 (0.0527) | 29.18 | 6,03E-05 |
| 12 | rs12509595 | - | 4 | C | T | 0.297929 | 0.3124 | 0.0577 (0.0091) | 0.1462 (0.0533) | 40.2 | 8,30E-05 |
| 13 | rs12633109 | - | 3 | T | G | 0.360897 | 0.3263 | 0.0436 (0.0089) | 0.0395 (0.0526) | 24 | 4,96E-05 |
| 14 | rs1275980 | KCNK3 | 2 | T | C | 0.498545 | 0.5481 | -0.0582 (0.0094) | -0.0128 (0.0495) | 38.33 | 7,92E-05 |
| 15 | rs13123551 | - | 4 | A | T | 0.616442 | 0.5065 | 0.0552 (0.0097) | 0.0818 (0.0493) | 32.38 | 6,69E-05 |
| 16 | rs147871383 | MIR99AHG | 21 | A | G | 0.02368 | 0.01925 | 0.2259 (0.0494) | 0.2534 (0.1903) | 20.91 | 4,32E-05 |
| 17 | rs16918175 | - | 10 | C | T | 0.104037 | 0.0811 | -0.0649 (0.0131) | 0.0332 (0.0895) | 24.54 | 5,07E-05 |
| 18 | rs17182166 | ACVR1 | 2 | T | G | 0.12494 | 0.137 | 0.0751 (0.0156) | -0.0046 (0.0703) | 23.18 | 4,79E-05 |
| 19 | rs1906779 | - | 15 | A | G | 0.211386 | 0.1239 | 0.0438 (0.0095) | -0.0342 (0.0758) | 21.26 | 4,39E-05 |
| 20 | rs1948696 | ITGB5 | 3 | C | T | 0.650138 | 0.6041 | -0.0459 (0.0091) | 0.0219 (0.0501) | 25.44 | 5,25E-05 |
| 21 | rs1973765 | LSP1 | 11 | C | T | 0.486717 | 0.4124 | 0.0441 (0.0092) | -0.0229 (0.0505) | 22.98 | 4,75E-05 |
| 22 | rs2429123 | CACNA1C, DCP1B | 12 | C | T | 0.729926 | 0.725 | 0.0473 (0.0098) | 0.0190 (0.0549) | 23.3 | 4,81E-05 |
| 23 | rs2447561 | - | 8 | T | A | 0.873442 | 0.8382 | -0.0589 (0.0116) | -0.0570 (0.0674) | 25.78 | 5,33E-05 |
| 24 | rs245015 | MSH3 | 5 | A | G | 0.685061 | 0.6562 | -0.0456 (0.0089) | 0.0065 (0.0519) | 26.25 | 5,42E-05 |
| 25 | rs2501968 | CENPQ | 6 | G | A | 0.465302 | 0.4229 | -0.0490 (0.0086) | -0.0879 (0.0498) | 32.46 | 6,71E-05 |
| 26 | rs2526620 | - | 7 | G | A | 0.252144 | 0.2459 | 0.0500 (0.0092) | 0.0505 (0.0576) | 29.54 | 6,10E-05 |
| 27 | rs284160 | TGFBR3 | 1 | A | G | 0.192644 | 0.1186 | 0.0572 (0.0096) | 0.1573 (0.0768) | 35.5 | 7,33E-05 |
| 28 | rs2842870 | PMF1, PMF1-BGLAP | 1 | C | T | 0.360398 | 0.3434 | -0.0446 (0.0087) | -0.0727 (0.0518) | 26.28 | 5,43E-05 |
| 29 | rs2880492 | NCOR2 | 12 | C | T | 0.035841 | 0.01984 | -0.1293 (0.0266) | -0.1917 (0.1781) | 23.63 | 4,88E-05 |
| 30 | rs35790371 | RBFOX1 | 16 | A | G | 0.005125 | 0.001903 | 0.5587 (0.1176) | 0.1447 (0.5420) | 22.57 | 4,66E-05 |
| 31 | rs5752720 | TTC28 | 22 | T | C | 0.28441 | 0.1844 | 0.0422 (0.0091) | -0.0311 (0.0635) | 21.51 | 4,44E-05 |
| 32 | rs57694670 | SH3PXD2A | 10 | G | A | 0.433638 | 0.3834 | -0.0507 (0.0086) | 0.0080 (0.0507) | 34.76 | 7,18E-05 |
| 33 | rs6462001 | - | 7 | T | G | 0.881494 | 0.8329 | -0.0765 (0.0151) | -0.0567 (0.0662) | 25.67 | 5,30E-05 |
| 34 | rs6843082 | - | 4 | A | G | 0.65294 | 0.6907 | -0.0438 (0.0092) | -0.0568 (0.0534) | 22.67 | 4,68E-05 |
| 35 | rs7091346 | SH3PXD2A | 10 | T | C | 0.45273 | 0.2664 | -0.0582 (0.0098) | -0.0143 (0.0552) | 35.27 | 7,28E-05 |
| 36 | rs7194129 | CFDP1 | 16 | T | C | 0.573691 | 0.5538 | 0.0399 (0.0084) | 0.0293 (0.0495) | 22.56 | 4,66E-05 |
| 37 | rs7341574 | ZFPM2 | 8 | T | C | 0.313243 | 0.3895 | 0.0463 (0.0094) | -0.0774 (0.0505) | 24.26 | 5,01E-05 |
| 38 | rs7451833 | - | 6 | G | A | 0.103042 | 0.1166 | 0.1280 (0.0220) | 0.1253 (0.0777) | 33.85 | 6,99E-05 |
| 39 | rs74617384 | LPA | 6 | T | A | 0.070999 | 0.04573 | 0.1415 (0.0281) | -0.0532 (0.1162) | 25.36 | 5,24E-05 |
| 40 | rs74849463 | PIK3C2B | 1 | T | C | 0.240949 | 0.2043 | 0.0445 (0.0092) | -0.0319 (0.0616) | 23.4 | 4,83E-05 |
| 41 | rs757241 | AFAP1-AS1 | 4 | C | G | 0.673377 | 0.6704 | -0.0681 (0.0145) | -0.0274 (0.0523) | 22.06 | 4,56E-05 |
| 42 | rs76099321 | CNNM2 | 10 | A | G | 0.052085 | 0.02979 | -0.0691 (0.0150) | -0.2215 (0.1464) | 21.22 | 4,38E-05 |
| 43 | rs7670136 | - | 4 | C | T | 0.571088 | 0.6269 | -0.0473 (0.0101) | 0.0391 (0.0509) | 21.93 | 4,53E-05 |
| 44 | rs77455924 | NTM | 11 | T | C | 0.056457 | 0.02996 | 0.1632 (0.0331) | 0.0444 (0.1477) | 24.31 | 5,02E-05 |
| 45 | rs7820334 | - | 8 | T | C | 0.246549 | 0.2949 | -0.0643 (0.0127) | -0.1327 (0.0648) | 25.63 | 5,29E-05 |
| 46 | rs7859727 | CDKN2B-AS1 | 9 | T | C | 0.518845 | 0.4151 | 0.0569 (0.0087) | 0.0281 (0.0498) | 42.77 | 8,83E-05 |
| 47 | rs7989823 | COL4A1, COL4A2 | 13 | C | A | 0.591615 | 0.6188 | 0.0527 (0.0089) | 0.0693 (0.0518) | 35.06 | 7,24E-05 |
| 48 | rs79960344 | - | 17 | G | T | 0.102808 | 0.1031 | 0.0651 (0.0130) | -0.1596 (0.0805) | 25.08 | 5,18E-05 |
| 49 | rs880315 | CASZ1 | 1 | C | T | 0.445525 | 0.4139 | 0.0416 (0.0089) | 0.0444 (0.0501) | 21.85 | 4,51E-05 |
| 50 | rs9112 | LOC100505841 | 5 | A | G | 0.436595 | 0.3573 | 0.0407 (0.0085) | -0.0274 (0.0523) | 22.93 | 4,74E-05 |
| 51 | rs979380 | - | 17 | A | G | 0.531844 | 0.6197 | -0.0417 (0.0084) | -0.0196 (0.0506) | 24.64 | 5,09E-05 |
| EAF.IS e EAF.VD indicano le frequenze degli alleli d'effetto rispettivamente nei dati relativi all'esposizione all'ictus ischemico e all'esito relativo alla demenza vascolare. Statistica F = (β/SE)². R² per SNP = F/(F + N − 2), dove N = 484.121. Gli SNP sono stati selezionati utilizzando p < 5 × 10-6 e raggruppamento per disequilibrio di associazione con r² < 0,001 all'interno di una finestra di 10.000 kb. |
Tabella 1: Strumenti genetici per l'ictus ischemico e associazioni corrispondenti SNP–esito per la demenza vascolare (VaD). La tabella elenca i 51 SNP finali, inclusi l'ID del SNP, il gene mappato, il cromosoma, gli alleli, la FAE, le stime SNP–carattere, le statistiche F e l'R2 per ogni SNP. I SNP sono stati selezionati con p < 5 × 10-6 e raggruppati con r2 < 0,001 entro 10.000 kb. EAF.IS ed EAF.VD indicano le frequenze dell'allele effetto nei dataset IS e VD, rispettivamente. Chr., cromosoma; EA, allele effetto; OA, altro allele; EAF, frequenza dell'allele effetto; IS, ictus ischemico; VD, demenza vascolare. Cliccare qui per scaricare questa tabella.
Tabella supplementare 1. Valutazione della pleiotropia orizzontale mediante MR-PRESSO per l'analisi ictus ischemico–demenza vascolare. Questa tabella riassume l'analisi MR-PRESSO eseguita per valutare la pleiotropia orizzontale globale e le possibili varianti outlier utilizzando il set di dati finale armonizzato composto da 51 SNP.Cliccare qui per scaricare il file.
Tabella supplementare 2. Passaggi del protocollo e risultati corrispondenti di convalida. Questa tabella collega ogni passaggio principale del protocollo al relativo risultato rappresentativo e alla posizione riportata nel manoscritto, dimostrando l'implementazione e la convalida del flusso di lavoro analitico.Cliccare qui per scaricare il file.