$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
FAERS-gegevensbronnen
De praktijkgegevens voor deze studie zijn verkregen uit de FAERS-database (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). Dit is een openbaar toegankelijke en geanonimiseerde database, dus ethische goedkeuring was niet vereist voor deze studie. Informatie over de vijf opgenomen geneesmiddelen is weergegeven in Tabel 1. De zoektocht werd uitgevoerd door eerst de bijwerkingen te beperken tot een heupfractuur, en de zoektijd liep van Q1 2014 tot Q4 2023. Om betrouwbare en stabiele gegevens te garanderen, standaardiseerde de studie de terminologie van gerapporteerde bijwerkingen via het MedDRA Dictionary versie 26.142. FAERS kwartaalbestanden van ASCII van Q1 2014 tot Q4 2023 werden gedownload en geïmporteerd voor analyse. De geëxtraheerde FAERS-tabellen omvatten DEMO, DRUG, REAC, THER, RPSR en OUTC. Deze tafels werden in alle vertrekken samengevoegd vóór de vertoning. Meldingen werden gekoppeld met behulp van CASEID en PRIMARYID om consistentie te waarborgen tussen demografische, geneesmiddel-, reactie-, therapie-, rapportage- en uitkomstinformatie. De doelmedicijnen waren onder andere etanercept, adalimumab, infliximab, certolizumab, pegol en golimumab. Geneesmiddelnamen in de DRUGS-tabel werden gestandaardiseerd door tekst om te zetten naar hoofdletters, extra spaties te verwijderen en spellingsvarianten te controleren wanneer nodig. Doelgeneesmiddelen werden geïdentificeerd met gestandaardiseerde generieke namen in de GENEESMIDDEL-tabel, en de restrictie van de rol van geneesmiddelen werd uitgevoerd met ROLE_COD = "PS", wat het primaire verdachte geneesmiddel43 aangeeft. De specifieke routekaart voor de studiescreeningstechniek is te zien in Figuur 2. Na het opschonen en screenen van de gegevens werd een unieke en analyseerbare dataset van geschikte rapporten over heupfracturen verkregen voor latere analyses.
Dubbele FAERS-rapporten werden verwijderd vóór signaaldetectie. Duplicaten werden geïdentificeerd volgens CASEID en PRIMARYID. Wanneer meerdere rapporten dezelfde CASEID deelden, werd het meest recente rapport volgens FDA_DT behouden. Als meerdere rapporten dezelfde CASEID en FDA_DT hadden, werd het rapport met de hoogste PRIMARYID behouden. Na deduplicatie droeg elke CASEID slechts één record bij aan de uiteindelijke analytische dataset. Meldingen werden opgenomen als ze aan alle volgende criteria voldeden: rapportagedatum tussen Q1 2014 en Q4 2023; het bijeffect werd gecodeerd als "Heupfractuur"; ten minste één van de vijf TNF-remmers werd vermeld in de GENEESMIDDEL-tabel; en de rol van drugs werd als hoofdverdachte gecodeerd. Meldingen werden uitgesloten als het dubbele records waren, het ontbraken aan geldige CASEID- of PRIMARYID-informatie, geen overeenkomstige GENEESMIDDEL- of REAC-invoer bevatten, het doelbijwerkingsincident niet bevatten, of de doel-TNF-remmer alleen als een gelijktijdig of secundair verdacht middel werden vermeld.
GWAS-databronnen voor Mendeliaanse randomisatie
De blootstellingsgegevens voor TNF-α in de Mendeliaanse randomisatie voor deze studie zijn verkregen uit de IEU OpenGWAS-database (https://gwas.mrcieu.ac.uk/), met GWAS ID prot-c-3722_49_2 uit de studie van Suhre K et al. De onderzoekspopulatie was van Europese afkomst en het aantal SNP's bedroeg 501.42844.
De blootstellingsgegevens voor sTNFR1 in de mendelische randomisatie voor deze studie zijn verkregen uit de IEU OpenGWAS-database (https://gwas.mrcieu.ac.uk/), met GWAS ID prot-c-2654_19_1 uit de studie van Suhre K et al. De onderzoekspopulatie was van Europese afkomst en het aantal SNP's bedroeg 501.42844.
De blootstellingsgegevens voor sTNFR2 in de Mendeliaanse randomisatie voor deze studie zijn verkregen uit de IEU OpenGWAS-database (https://gwas.mrcieu.ac.uk/), met GWAS ID prot-c-3152_57_1 uit de studie van Suhre K et al. De onderzoekspopulatie was van Europese afkomst en het aantal SNP's bedroeg 501.42844.
De uitkomstgegevens voor heupfractuur, GWAS ID GCST90161240, gedeponeerd in de GWAS Catalog (https://www.ebi.ac.uk/gwas/studies/GCST90161240), zijn gegevens van een meta-analyse van een grootschalige GWAS die 11.516 gevallen van heupfractuur en 723.838 controlesomvatte 45. De ziekteclassificatie komt overeen met de Internationale Ziekteclassificatie (ICD; ICD-10-codes S72.0–S72.2 en ICD-9-code 820).
Ethische goedkeuring en geïnformeerde toestemming waren verkregen in de oorspronkelijke GWAS-studies. Omdat de huidige studie gebruikmaakte van publiek toegankelijke, geanonimiseerde FAERS-gegevens en openbaar beschikbare GWAS-samenvattende statistieken, was geen aanvullende ethische goedkeuring vereist.
Softwareomgeving en workflowimplementatie
Alle analyses werden uitgevoerd met R versie 4.3.2. FAERS-gegevensimport, schoonmaak, samenvoegen en tabuleren werden uitgevoerd met behulp van R-gebaseerde databeheerworkflows. Datatabellen werden geïmporteerd met functies zoals data.table::fread() of readr::read_delim(), samengevoegd met CASEID en PRIMARYID, en verwerkt met dplyr-functies. Beschrijvende statistieken en 2 × 2 contingentietabellen werden gegenereerd met behulp van aangepaste R-scripts.
Mendeliaanse randomisatieanalyses werden uitgevoerd met TwoSampleMR versie 0.5.6. Blootstellingsinstrumenten werden geëxtraheerd met een significantiedrempel van P < 1 × 10⁻5 of geformatteerd uit GWAS-samenvattingsstatistieken met behulp van TwoSampleMR-compatibele invoerstructuren. Instrumentklonteren werd uitgevoerd met clump_data() met clump_r2 = 0,001 en clump_kb = 10.000. Uitkomstgegevens werden uitgepakt of opgemaakt met extract_outcome_data() of read_outcome_data(), afhankelijk van het bronformaat. Blootstellings- en uitkomstdatasets werden geharmoniseerd met behulp van harmonise_data(). Causale schattingen werden gegenereerd met mr() met de volgende mendeliaanse randomisatiemethoden: MR-Egger, gewogen mediaan, inverse variantiegewicht, eenvoudige modus en gewogen modus. Heterogeniteit werd beoordeeld met mr_heterogeneity(), en horizontale pleiotropie met behulp van mr_pleiotropy_test(). Alle datasets werden geïmporteerd, opgeschoond, geharmoniseerd en geanalyseerd binnen deze softwareomgeving om een consistente en reproduceerbare analytische workflow te garanderen.
Farmacovigilantie-analyse
Beschrijvende analyses werden gebruikt om bijwerkingen gerelateerd aan heupfracturen die met de vijf geneesmiddelen geassocieerd zijn, samen te vatten. Signaaldetectieanalyses werden vervolgens uitgevoerd met vier disproportionaliteitsalgoritmen, waaronder de rapportageodds ratio (ROR), proportionele rapportageratio (PRR), multi-item gamma Poisson-shrinker (MGPS) en Bayesian confidence propagation neural network (BCPNN). De criteria voor de vier belangrijkste algoritmen worden weergegeven in Tabel 246.
Mendeliaanse randomisatieanalyse
Samenvattende statistieken voor TNF-α, sTNFR1 en sTNFR2 werden geëxtraheerd als blootstellingsdatasets, en samenvattende statistieken van heupfracturen werden als uitkomstdataset geëxtraheerd. Analyses waren beperkt tot Europese afkomstdatasets wanneer beschikbaar om populatiestratificatie-bias te verminderen.
Om bias veroorzaakt door linkage-disequilibrium en zwakke instrumenten te minimaliseren, werden de volgende criteria toegepast: genoom-brede significantie-drempel P < 1 × 10⁻5, linkage-disequilibriumdrempel r2 < 0,001, klontervenster van 10.000 kb en F-statistiek > 20. De F-statistiek werd voor elke behouden instrumentele variabele berekend als beta2/se2 om de instrumentsterkte te evalueren. SNP's met F-statistiek ≤ 20 werden uitgesloten van downstream-analyses.
Na SNP-selectie werden blootstellings- en uitkomstdatasets geharmoniseerd om effectallelen uit te lijnen. Tijdens de harmonisatie werden effectallelen en andere allelen uitgelijnd tussen de blootstellings- en uitkomstdatasets. SNP's met incompatibele allelen werden verwijderd, en palindromische SNP's met ambigue allelfrequenties werden uitgesloten wanneer de oriëntatie van de streng niet kon worden bepaald. Na harmonisatie werden de behouden SNP's gecontroleerd om te bevestigen dat bètacecoëfficiënten overeenkwamen met hetzelfde effectallel in beide datasets. Het aantal SNP's dat na het samenvoegen en harmoniseren werd vastgelegd voor elke blootstelling als een tussentijdse reproduceerbaarheidscontrole.
Vijf Mendeliaanse randomisatiemethoden werden toegepast, waaronder MR-Egger, gewogen mediaan, inverse variantiegewicht, eenvoudige modus en gewogen modus. De potentiële heterogeniteit van instrumentele variabelen werd beoordeeld met behulp van de Q-test van Cochran, en P < 0,05 werd als indicatief voor significante heterogeniteit beschouwd. Potentiële horizontale pleiotropie werd geëvalueerd met behulp van het MR-Egger-intercept, en P < 0,05 gaf pleiotropie aan, wat wijst op een verminderde betrouwbaarheid van de causale schatting47. Deze analyses genereerden schattingen van causale effecten samen met heterogeniteits- en pleiotropiestatistieken voor elke blootstelling.
Tussenliggende controlepunten voor reproduceerbaarheid
Tussenliggende controlepunten werden na elke belangrijke verwerkingsstap geregistreerd om reproduceerbaarheid van de workflow te waarborgen. Voor de FAERS-workflow omvatten checkpoints het aantal geïmporteerde DEMO-records, het aantal unieke records na deduplicatie, het aantal meldingen met heupfractuur als doelbijwerking, het aantal meldingen met de vijf TNF-remmers, en het uiteindelijke aantal in aanmerking komende meldingen waarin TNF-remmers als primaire verdachte geneesmiddelen werden geregistreerd. Voor de Mendeliaanse randomisatieworkflow omvatten de checkpoints het aantal SNP's dat per blootstelling werd geëxtraheerd, het aantal SNP's dat behouden bleef na koppelingsdisequilibrium-clumping, het aantal SNP's beschikbaar in de uitkomstdataset, het aantal SNP's dat na harmonisatie werd behouden, en het uiteindelijke aantal instrumentele variabelen dat in elke Mendeliaanse randomisatieanalyse werd gebruikt.
Statistische rapportage
Continue resultaten werden gerapporteerd met overeenkomstige effectschattingen, 95% betrouwbaarheidsintervallen (95% BI) en P-waarden . Tenzij anders vermeld, werd statistische significantie gedefinieerd als een tweezijdige P < 0,05. Voor farmacovigilantie-analyse werden beschrijvende tellingen en disproportionaliteitsschattingen gerapporteerd voor elke individuele TNF-inhibitor en voor de gepoolde TNF-inhibitorgroep. Voor Mendeliaanse randomisatieanalyse werden causale schattingen, standaardfouten, 95% betrouwbaarheidsintervallen, P-waarden , heterogeniteitsstatistieken, pleiotropietestresultaten en het aantal behouden SNP's gerapporteerd voor elke blootstelling.