$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Fonti dati FAERS
I dati reali per questo studio sono stati ottenuti dal database FAERS (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). Si tratta di un database pubblico e anonimizzato, quindi non era richiesta un'approvazione etica per questo studio. Le informazioni sui cinque farmaci inclusi sono mostrate nella Tabella 1. La ricerca è stata effettuata limitando innanzitutto gli eventi avversi alla frattura dell'anca, e il periodo di ricerca è andato dal primo trimestre 2014 al quarto trimestre 2023. Per garantire dati affidabili e stabili, lo studio ha standardizzato la terminologia degli eventi avversi segnalati tramite la versione 26.142 del Dizionario MedDRA. I file trimestrali ASCII di FAERS dal primo trimestre 2014 al quarto trimestre 2023 sono stati scaricati e importati per l'analisi. Le tabelle estratte di FAERS includevano DEMO, DRUG, REAC, THER, RPSR e OUTC. Questi tavoli venivano uniti in tutti gli ambienti prima della proiezione. I rapporti sono stati collegati tramite CASEID e PRIMARYID per garantire coerenza tra dati demografici, farmacologici, reazionali, terapeutici, segnalatori e risultati. I farmaci target includevano etanercet, adalimumab, infliximab, certolizumab pegol e golimumab. I nomi dei farmaci nella tabella DRUG sono stati standardizzati convertendo il testo in maiuscolo, eliminando spazi extra e controllando le varianti ortografiche quando necessario. I farmaci target sono stati identificati usando nomi generici standardizzati nella tabella DRUG e la restrizione del ruolo del farmaco è stata eseguita usando ROLE_COD = "PS", indicando il principale farmacosospetto 43. La roadmap specifica della tecnica di screening dello studio è mostrata nella Figura 2. Dopo la pulizia e lo screening dei dati, è stato ottenuto un dataset unico e analizzabile di segnalazioni di fratture dell'anca idonee per analisi successive.
I rapporti duplicati del FAERS venivano rimossi prima del rilevamento del segnale. I duplicati sono stati identificati secondo CASEID e PRIMARYID. Quando più segnalazioni condividevano lo stesso CASEID, secondo FDA_DT il rapporto più recente veniva mantenuto. Se più segnalazioni avevano lo stesso CASEID e FDA_DT, il rapporto con il PRIMARYID più alto veniva trattenuto. Dopo la deduplicazione, ogni CASEID ha contribuito con un solo record al dataset analitico finale. I rapporti sono stati inclusi se soddisfacevano tutti i seguenti criteri: data di segnalazione tra il primo trimestre 2014 e il quarto trimestre 2023; l'evento avverso fu codificato come "frattura dell'anca"; almeno uno dei cinque inibitori del TNF era registrato nella tabella DRUGS; e il ruolo della droga era indicato come sospetto principale. I referti sono stati esclusi se erano record duplicati, privi di informazioni valide su CASEID o PRIMARYID, senza una corrispondente voce DRUG o REAC, non includevano l'evento avverso bersaglio, o indicavano l'inibitore del TNF come farmaco sospetto concomitante o secondario.
Fonti di dati GWAS per la randomizzazione mendeliana
I dati di esposizione per TNF-α nella randomizzazione mendeliana di questo studio sono stati ottenuti dal database IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), con GWAS ID prot-c-3722_49_2 dallo studio di Suhre K et al. La popolazione studiata era di origine europea e il numero di SNP era di 501.42844.
I dati di esposizione per sTNFR1 nella randomizzazione mendeliana di questo studio sono stati ottenuti dal database IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), con GWAS ID prot-c-2654_19_1 dallo studio di Suhre K et al. La popolazione studiata era di origine europea e il numero di SNP era di 501.42844.
I dati di esposizione per sTNFR2 nella randomizzazione mendeliana di questo studio sono stati ottenuti dal database IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), con GWAS ID prot-c-3152_57_1 dallo studio di Suhre K et al. La popolazione studiata era di origine europea e il numero di SNP era di 501.42844.
I dati sugli esiti per la frattura dell'anca, GWAS ID GCST90161240, depositati nel GWAS Catalog (https://www.ebi.ac.uk/gwas/studies/GCST90161240), sono dati provenienti da una meta-analisi di un GWAS su larga scala che ha incluso 11.516 casi di frattura dell'anca e 723.838controlli 45. La classificazione della malattia è in linea con la Classificazione Internazionale delle Malattie (ICD; codici ICD-10 S72.0–S72.2 e codice ICD-9 820).
L'approvazione etica e il consenso informato erano stati ottenuti negli studi GWAS originali. Poiché il presente studio utilizzava dati FAERS anonimizzati e pubblici e statistiche GWAS riassuntive disponibili pubblicamente, non era necessaria un'ulteriore approvazione etica.
Ambiente software e implementazione dei flussi di lavoro
Tutte le analisi sono state effettuate utilizzando la versione 4.3.2 di R. L'importazione, pulizia, fusione e tabulazione dei dati FAERS venivano effettuate utilizzando flussi di lavoro di gestione dati basati su R. Le tabelle dati venivano importate usando funzioni come data.table::fread() o readr::read_delim(), unite usando CASEID e PRIMARYID, e processate con funzioni dplyr. Statistiche descrittive e 2 × 2 tabelle di contingenza sono state generate utilizzando script R personalizzati.
Le analisi di randomizzazione mendeliane sono state condotte utilizzando TwoSampleMR versione 0.5.6. Gli strumenti di esposizione sono stati estratti utilizzando una soglia di significatività di P < 1 × 10⁻5 oppure formattati da statistiche riassuntive GWAS utilizzando strutture di input compatibili con TwoSampleMR. Il raggruppamento degli strumenti è stato eseguito usando clump_data() con clump_r2 = 0,001 e clump_kb = 10.000. I dati di risultato sono stati estratti o formattati usando extract_outcome_data() o read_outcome_data(), a seconda del formato di origine. I dataset di esposizione ed esito sono stati armonizzati usando harmonise_data(). Le stime causali sono state generate utilizzando mr() con i seguenti metodi di randomizzazione mendeliana: MR-Egger, mediana ponderata, ponderata per varianza inversa, modalità semplice e modalità ponderata. L'eterogeneità è stata valutata usando mr_heterogeneity(), e la pleiotropia orizzontale con mr_pleiotropy_test(). Tutti i dataset sono stati importati, puliti, armonizzati e analizzati all'interno di questo ambiente software per garantire un flusso di lavoro analitico coerente e riproducibile.
Analisi della farmacovigilanza
Sono state utilizzate analisi descrittive per riassumere gli eventi avversi correlati alle fratture dell'anca associati ai cinque farmaci. Le analisi di rilevamento del segnale sono state quindi eseguite utilizzando quattro algoritmi di sproporzionalità, tra cui il rapporto di segnalazione (ROR), il rapporto proporzionale (PRR), il restringente di Poisson gamma multi-item (MGPS) e la rete neurale bayesiana di propagazione della fiducia (BCPNN). I criteri per i quattro principali algoritmi sono mostrati nella Tabella 246.
Analisi di randomizzazione mendeliana
Le statistiche riassuntive per TNF-α, sTNFR1 e sTNFR2 sono state estratte come dataset di esposizione, mentre le statistiche riassuntive sulle fratture dell'anca sono state estratte come dataset di risultati. Le analisi sono state limitate a dataset di ascendenza europei quando disponibili per ridurre il bias di stratificazione della popolazione.
Per minimizzare il bias causato da disequilibrio di collegamento e strumenti deboli, sono stati applicati i seguenti criteri: soglia di significatività genomica P < 1 × 10⁻5, soglia di disequilibrio di collegamento r2 < 0,001, finestra di aggregazione di 10.000 kb e F-statistica > 20. La F-statistica è stata calcolata per ciascuna variabile strumentale mantenuta come beta2/se2 per valutare la forza dello strumento. Gli SNP con F-statistico ≤ 20 sono stati esclusi dalle analisi a valle.
Dopo la selezione degli SNP, i dataset di esposizione ed esito sono stati armonizzati per allineare gli alleli di effetto. Durante l'armonizzazione, gli alleli di effetto e gli altri alleli sono stati allineati tra i dataset di esposizione e di esame. Gli SNP con alleli incompatibili sono stati rimossi e gli SNP palindromici con frequenze alleli ambigue sono stati esclusi quando non è stato possibile determinare l'orientamento del filo. Dopo l'armonizzazione, gli SNP trattenuti sono stati controllati per confermare che i coefficienti beta corrispondevano allo stesso allele di effetto in entrambi i dataset. Il numero di SNP mantenuti dopo aggregazione e armonizzazione è stato registrato per ogni esposizione come checkpoint di riproducibilità intermedia.
Sono stati applicati cinque metodi di randomizzazione mendeliana, tra cui MR-Egger, mediana ponderata, ponderata per varianza inversa, modalità semplice e modalità pesata. La potenziale eterogeneità delle variabili strumentali è stata valutata utilizzando il test Q di Cochran, e P < 0,05 è stata considerata indicativa di eterogeneità significativa. La potenziale pleiotropia orizzontale è stata valutata utilizzando l'intercettazione MR-Egger, e P < 0,05 ha indicato pleiotropia, suggerendo una riduzione dell'affidabilità della stimacausale 47. Queste analisi hanno generato stime degli effetti causali insieme a statistiche di eterogeneità e pleiotropia per ogni esposizione.
Checkpoint intermedi per la riproducibilità
Punti di controllo intermedi venivano registrati dopo ogni fase principale di elaborazione per garantire la riproducibilità del flusso di lavoro. Per il flusso di lavoro FAERS, i checkpoint includevano il numero di record DEMO importati, il numero di record unici dopo la deduplicazione, il numero di segnalazioni che contenevano fratture dell'anca come evento avverso bersaglio, il numero di segnalazioni che coinvolgevano i cinque inibitori del TNF e il numero finale di segnalazioni idonee in cui gli inibitori del TNF sono stati registrati come farmaci sospetti principali. Per il flusso di lavoro di randomizzazione mendeliana, i checkpoint includevano il numero di SNP estratti per ogni esposizione, il numero di SNP mantenuti dopo l'aggregazione di disequilibrio di collegamento, il numero di SNP disponibili nel dataset di risultati, il numero di SNP mantenuti dopo l'armonizzazione e il numero finale di variabili strumentali utilizzate in ogni analisi di randomizzazione mendeliana.
Reportistica statistica
Sono stati riportati risultati continui con stime di effetto corrispondenti, intervalli di confidenza al 95% (IC 95%) e valori P . Salvo diversa specificazione, la significanza statistica è stata definita come una P a due facce < 0,05. Per l'analisi della farmacovigilanza, sono stati riportati conteggi descrittivi e stime di sproporzionalità per ciascun singolo inibitore del TNF e per il gruppo aggregato degli inibitori del TNF. Per l'analisi di randomizzazione mendeliana, sono state riportate stime causali, errori standard, intervalli di confidenza al 95%, valori P , statistiche di eterogeneità, risultati del test di pleiotropia e numero di SNP trattenuti per ogni esposizione.