$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
FAERS-Datenquellen
Die realen Daten für diese Studie stammen aus der FAERS-Datenbank (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). Dies ist eine öffentlich zugängliche und anonymisierte Datenbank, daher war für diese Studie keine ethische Genehmigung erforderlich. Informationen zu den fünf enthaltenen Medikamenten sind in Tabelle 1 dargestellt. Die Suche wurde durchgeführt, indem zunächst die unerwünschten Ereignisse auf einen Hüftbruch beschränkt wurden, und die Suchzeit reichte vom ersten bis zum vierten Quartal 2023. Um verlässliche und stabile Daten sicherzustellen, standardisierte die Studie die Terminologie der gemeldeten unerwünschten Ereignisse im MedDRA-Wörterbuch Version 26.142. Die vierteljährlichen ASCII-Dateien von FAERS vom ersten Quartal 2014 bis zum vierten Quartal 2023 wurden heruntergeladen und zur Analyse importiert. Die extrahierten FAERS-Tabellen umfassten DEMO, DRUG, REAC, THER, RPSR und OUTC. Diese Tische wurden vor der Sichtung in allen Bereichen zusammengelegt. Die Berichte wurden mittels CASEID und PRIMARYID verknüpft, um Konsistenz zwischen demografischen, Medikamenten-, Reaktions-, Therapie-, Berichts- und Ergebnisinformationen sicherzustellen. Zu den Zielmedikamenten gehörten Etanercept, Adalimumab, Infliximab, Certolizumab, Pegol und Golimumab. Die Arzneimittelnamen in der DRUG Tabelle wurden standardisiert, indem Text in Großbuchstaben umgewandelt, zusätzliche Leerzeichen entfernt und bei Bedarf Rechtschreibvarianten überprüft wurden. Zielmedikamente wurden mit standardisierten Generikanamen in der ARZNEIMITTEL-Tabelle identifiziert, und die Arzneimittelrolle-Einschränkung wurde mit ROLE_COD = "PS" durchgeführt, was das primäre verdächtigeMedikament 43 anzeigt. Die konkrete Roadmap für die Screening-Technik der Studien ist in Abbildung 2 dargestellt. Nach der Datenbereinigung und Screening wurde ein einzigartiger und analysierbarer Datensatz mit geeigneten Hüftfrakturberichten für spätere Analysen gewonnen.
Duplikate FAERS-Berichte wurden vor der Signalerkennung entfernt. Duplikate wurden gemäß CASEID und PRIMARYID identifiziert. Als mehrere Berichte dasselbe CASEID teilten, wurde laut FDA_DT der jüngste Bericht beibehalten. Wenn mehrere Berichte dieselbe CASEID und FDA_DT hatten, wurde der Bericht mit der höchsten PRIMARYID behalten. Nach der Deduplizierung trug jede CASEID nur einen Datensatz zum finalen analytischen Datensatz bei. Die Berichte wurden aufgenommen, wenn sie alle folgenden Kriterien erfüllten: Berichtsdatum zwischen dem ersten Quartal 2014 und dem vierten Quartal 2023; das unerwünschte Ereignis wurde als "Hüftfraktur" kodiert; mindestens einer der fünf TNF-Inhibitoren wurde in der Arzneimitteltabelle verzeichnet; und die Drogenrolle wurde als Hauptverdächtiger kodiert. Berichte wurden ausgeschlossen, wenn es sich um doppelte Einträge handelte, keine gültigen CASEID- oder PRIMARYID-Informationen enthielten, keine entsprechenden ARZNEIMITTEL- oder REAC-Einträge enthielten, das Ziel-Nebenereignis nicht enthielten oder der Ziel-TNF-Inhibitor nur als begleitendes oder sekundäres verdächtiges Medikament aufgeführt wurden.
GWAS-Datenquellen für Mendelsche Randomisierung
Die Expositionsdaten für TNF-α in der Mendelschen Randomisierung dieser Studie stammen aus der IEU OpenGWAS-Datenbank (https://gwas.mrcieu.ac.uk/) mit GWAS ID prot-c-3722_49_2 aus der Studie von Suhre K et al. Die Studienpopulation war europäischer Abstammung, und die Zahl der SNPs betrug 501.42844.
Die Expositionsdaten für sTNFR1 in der mendelschen Randomisierung dieser Studie wurden aus der IEU OpenGWAS-Datenbank (https://gwas.mrcieu.ac.uk/) mit GWAS ID prot-c-2654_19_1 aus der Studie von Suhre K et al. gewonnen. Die Studienpopulation war europäischer Abstammung, und die Zahl der SNPs betrug 501.42844.
Die Expositionsdaten für sTNFR2 in der Mendelschen Randomisierung dieser Studie stammen aus der IEU OpenGWAS-Datenbank (https://gwas.mrcieu.ac.uk/) mit GWAS ID prot-c-3152_57_1 aus der Studie von Suhre K et al. Die Studienpopulation war europäischer Abstammung, und die Zahl der SNPs betrug 501.42844.
Die Ergebnisdaten für Hüftfraktur, GWAS ID GCST90161240, die im GWAS-Katalog (https://www.ebi.ac.uk/gwas/studies/GCST90161240) abgelegt sind, stammen aus einer Metaanalyse eines groß angelegten GWAS, das 11.516 Hüftfrakturfälle und 723.838 Kontrollgruppenumfasste. Die Krankheitsklassifikation entspricht der Internationalen Krankheitsklassifikation (ICD; ICD-10-Codes S72.0–S72.2 und ICD-9-Code 820).
In den ursprünglichen GWAS-Studien wurden ethische Genehmigung und informierte Zustimmung eingeholt. Da die vorliegende Studie öffentlich zugängliche, anonymisierte FAERS-Daten und öffentlich verfügbare GWAS-Zusammenfassungsstatistiken verwendete, war keine zusätzliche ethische Genehmigung erforderlich.
Softwareumgebung und Workflow-Implementierung
Alle Analysen wurden mit R-Version 4.3.2 durchgeführt. Die Import, Reinigung, Zusammenführung und Tabellierung von FAERS-Daten wurden mithilfe von R-basierten Datenmanagement-Workflows durchgeführt. Datentabellen wurden mit Funktionen wie data.table::fread() oder readr::read_delim() importiert, mit CASEID und PRIMARYID zusammengeführt und mit dplyr-Funktionen verarbeitet. Beschreibende Statistiken und 2 × 2 Kontingenztabellen wurden mit benutzerdefinierten R-Skripten erstellt.
Mendelsche Randomisierungsanalysen wurden mit TwoSampleMR Version 0.5.6 durchgeführt. Expositionsinstrumente wurden mit einer Signifikanzschwelle von P < 1 × 10⁻5 extrahiert oder aus GWAS-Zusammenfassungsstatistiken mit TwoSampleMR-kompatiblen Eingabestrukturen formatiert. Das Klumpen der Instrumente erfolgte mit clump_data() mit clump_r2 = 0,001 und clump_kb = 10.000. Die Ergebnisdaten wurden je nach Quellformat mit extract_outcome_data() oder read_outcome_data() extrahiert oder formatiert. Expositions- und Ergebnisdatensätze wurden mittels harmonise_data() harmonisiert. Kausale Schätzungen wurden mit mr() mit den folgenden mendelschen Randomisierungsmethoden erstellt: MR-Egger, gewichteter Median, inverse Varian-gewichtet, einfacher Modus und gewichteter Modus. Die Heterogenität wurde mit mr_heterogeneity() bewertet, die horizontale Pleiotropie mit mr_pleiotropy_test(). Alle Datensätze wurden innerhalb dieser Softwareumgebung importiert, bereinigt, harmonisiert und analysiert, um einen konsistenten und reproduzierbaren analytischen Workflow sicherzustellen.
Pharmakovigilanzanalyse
Deskriptive Analysen wurden verwendet, um hüftfrakturbedingte Nebenwirkungen der fünf Medikamente zusammenzufassen. Signalerkennungsanalysen wurden anschließend mit vier Disproportionalitätsalgorithmen durchgeführt, darunter das Reporting Odds Ratio (ROR), das proportionale Reporting Ratio (PRR), den Multi-Item-Gamma-Poisson-Schrumpfer (MGPS) und das Bayessche Confidence Propagation Neural Network (BCPNN). Die Kriterien für die vier Hauptalgorithmen sind in Tabelle 246 dargestellt.
Mendelsche Randomisierungsanalyse
Zusammenfassungsstatistiken für TNF-α, sTNFR1 und sTNFR2 wurden als Expositionsdatensätze extrahiert, und Hüftfraktur-Zusammenfassungsstatistiken wurden als Ergebnisdatensatz extrahiert. Die Analysen beschränkten sich auf europäische Abstammungsdatensätze, sofern verfügbar, um den Populationsstratifizierungs-Bias zu verringern.
Um den durch Verknüpfungsungleichgewicht und schwache Instrumente verursachte Verzerrung zu minimieren, wurden folgende Kriterien angewandt: genomweite Signifikanzschwelle P < 1 × 10⁻5, Kopplungsungleichgewichtsschwelle r2 < 0,001, Klumpenfenster von 10.000 kb und F-Statistik > 20. Die F-Statistik wurde für jede erhaltene instrumentelle Variable als Beta2/se2 berechnet, um die Instrumentenstärke zu bewerten. SNPs mit F-Statistik ≤ 20 wurden von nachgelagerten Analysen ausgeschlossen.
Nach der Auswahl des SNP wurden die Expositions- und Ergebnisdatensätze harmonisiert, um Effektallele auszurichten. Während der Harmonisierung wurden Effektallele und andere Allele zwischen Expositions- und Ergebnisdatensätzen ausgerichtet. SNPs mit inkompatiblen Allelen wurden entfernt, und palindromische SNPs mit mehrdeutigen Allelfrequenzen wurden ausgeschlossen, wenn die Strangorientierung nicht bestimmt werden konnte. Nach der Harmonisierung wurden die erhaltenen SNPs überprüft, um zu bestätigen, dass die Beta-Koeffizienten demselben Effektallel in beiden Datensätzen entsprechen. Die Anzahl der nach der Klumpung und Harmonisierung erhaltenen SNPs wurde für jede Exposition als zwischenliegender Reproduzierbarkeitskontrollpunkt erfasst.
Fünf Mendelsche Randomisierungsmethoden wurden angewandt, darunter MR-Egger, gewichteter Median, inverse Varianzgewichtet, einfacher Modus und gewichteter Modus. Die potenzielle Heterogenität instrumenteller Variablen wurde mit dem Cochranschen Q-Test bewertet, und P < 0,05 wurde als Hinweis auf signifikante Heterogenität angesehen. Die potenzielle horizontale Pleiotropie wurde mit dem MR-Egger-Interzept bewertet, und P < 0,05 deutete auf Pleiotropie hin, was auf eine verminderte Zuverlässigkeit der kausalen Schätzung47 hindeutet. Diese Analysen erzeugten kausale Effektschätzungen sowie Heterogenitäts- und Pleiotropiestatistiken für jede Exposition.
Zwischenkontrollpunkte für die Reproduzierbarkeit
Zwischenkontrollpunkte wurden nach jedem wichtigen Verarbeitungsschritt aufgezeichnet, um die Reproduzierbarkeit des Arbeitsablaufs sicherzustellen. Für den FAERS-Workflow umfassten die Checkpoints die Anzahl der importierten DEMO-Datensätze, die Anzahl der eindeutigen Datensätze nach Deduplizierung, die Anzahl der Berichte mit Hüftfraktur als Ziel-Nebenwirkung, die Anzahl der Meldungen mit den fünf TNF-Inhibitoren und die endgültige Anzahl der zulässigen Berichte, in denen TNF-Inhibitoren als primäre Verdächtige verzeichnet wurden. Für den Mendelschen Randomisierungs-Workflow umfassten die Checkpoints die Anzahl der für jede Exposition extrahierten SNPs, die Anzahl der nach der Klumpung der Linkage Disequilibrium erhaltenen SNPs, die Anzahl der im Ergebnisdatensatz verfügbaren SNPs, die Anzahl der nach der Harmonisierung beibehaltenen SNPs und die endgültige Anzahl der in jeder Mendelschen Randomisierungsanalyse verwendeten Instrumentvariablen.
Statistische Berichterstattung
Kontinuierliche Ergebnisse wurden mit entsprechenden Effektschätzungen, 95%-Konfidenzintervallen (95%-KI) und P-Werten berichtet. Sofern nicht anders angegeben, wurde die statistische Signifikanz als zweiseitiges P < 0,05 definiert. Für die Pharmakovigilanzanalyse wurden deskriptive Zählungen und Disproportionalitätsschätzungen für jeden einzelnen TNF-Inhibitor und für die gepoolte TNF-Inhibitor-Gruppe berichtet. Für die Mendelsche Randomisierungsanalyse wurden kausale Schätzungen, Standardfehler, 95%-Konfidenzintervalle, P-Werte , Heterogenitätsstatistiken, Pleiotropietestergebnisse und die Anzahl der erhaltenen SNPs für jede Exposition berichtet.