$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Źródła danych FAERS
Dane rzeczywiste do tego badania pochodzą z bazy danych FAERS (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). Jest to publicznie dostępna i zanonimizowana baza danych, więc nie było wymagane etyczne zatwierdzenie do tego badania. Informacje o pięciu zawartych lekach przedstawiono w Tabeli 1. Przeszukanie przeprowadzono poprzez ograniczenie początkowo do złamań biodra, a czas poszukiwań obejmował okres od I kwartału 2014 do IV kwartału 2023. Aby zapewnić wiarygodność i stabilność danych, badanie ustandaryzowało terminologię zgłaszanych zdarzeń niepożądanych za pomocą słownika MedDRA w wersji 26.142. Kwartalne pliki ASCII FAERS od I kwartału 2014 do IV kwartału 2023 zostały pobrane i zaimportowane do analizy. Wyodrębnione tabele FAERS obejmowały DEMO, DRUG, REAC, THER, RPSR oraz OUTC. Stoły te zostały połączone we wszystkich miejscach przed pokazem. Raporty były powiązane za pomocą CASEID i PRIMARYID, aby zapewnić spójność w zakresie demografii, leków, reakcji, terapii, raportujących i wyników wyników. Do leków docelowych należały etanercept, adalimumab, infliximab, certolizumab pegol oraz golimumab. Nazwy leków w tabeli DRUG były standaryzowane przez konwersję tekstu na wielkie litery, usunięcie zbędnych spacji i sprawdzanie wariantów pisowni w razie potrzeby. Leki docelowe zostały zidentyfikowane za pomocą standaryzowanych nazw generycznych w tabeli LEKÓW, a ograniczenie roli leku przeprowadzono przy użyciu ROLE_COD = "PS", wskazując główny podejrzany lek43. Konkretna mapa drogowa techniki przesiewowej badania przedstawiona jest na Rysunku 2. Po oczyszczeniu i przesiewaniu danych uzyskano unikalny i analizowalny zestaw kwalifikujących się raportów o złamaniach biodra do kolejnych analiz.
Zduplikowane raporty FAERS zostały usunięte przed wykryciem sygnału. Duplikaty zostały zidentyfikowane według CASEID i PRIMARYID. Gdy wiele raportów zawierało ten sam CASEID, według FDA_DT zachowano najnowszy raport. Jeśli wiele raportów miało ten sam CASEID i FDA_DT, zatrzymywano ten z najwyższym PRIMARYID. Po deduplikacji każdy CASEID dostarczał tylko jeden rekord do końcowego zbioru danych analitycznych. Raporty były uwzględniane, jeśli spełniały wszystkie następujące kryteria: data raportu między I kwartałem 2014 a IV kwartałem 2023; zdarzenie niepożądane zostało zakodowane jako "złamanie biodra"; co najmniej jeden z pięciu inhibitorów TNF został odnotowany w tabeli LEKÓW; a rola narkotykowa została zakodowana jako główny podejrzany. Raporty były wykluczane, jeśli były duplikatywnymi dokumentami, brakowało prawidłowych informacji CASEID lub PRIMARYID, nie zawierały odpowiadającego wpisu LEK lub REAC, nie zawierały docelowego zdarzenia niepożądanego lub wskazywały docelowy inhibitor TNF wyłącznie jako lek towarzyszący lub wtórny podejrzany.
Źródła danych GWAS dla randomizacji mendelowskiej
Dane ekspozycji dla TNF-α w randomizacji mendelowskiej dla tego badania pochodzą z bazy danych IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), z identyfikatorem GWAS prot-c-3722_49_2 z badania Suhre K i in. Populacja badana miała pochodzenie europejskie, a liczba SNP wynosiła 501 42844.
Dane ekspozycji dla sTNFR1 w randomizacji mendelowskiej dla tego badania pochodzą z bazy danych IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), z identyfikatorem GWAS prot-c-c-2654_19_1 z badania Suhre K i in. Populacja badana miała pochodzenie europejskie, a liczba SNP wynosiła 501 42844.
Dane ekspozycji dla sTNFR2 w randomizacji mendelowskiej dla tego badania pochodzą z bazy danych IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), z identyfikatorem GWAS prot-c-c-3152_57_1 z badania Suhre K i in. Populacja badana miała pochodzenie europejskie, a liczba SNP wynosiła 501 42844.
Dane dotyczące złamania biodra, GWAS ID GCST90161240, złożone w katalogu GWAS (https://www.ebi.ac.uk/gwas/studies/GCST90161240), pochodzą z metaanalizy dużej skali GWAS, która obejmowała 11 516 przypadków złamań biodra oraz 723 838 osób kontrolnych45. Klasyfikacja chorób jest zgodna z Międzynarodową Klasyfikacją Chorób (ICD; kody ICD-10 S72.0–S72.2 oraz kod ICD-9 820).
W oryginalnych badaniach GWAS uzyskano etyczną zgodę i świadomą zgodę. Ponieważ obecne badanie wykorzystywało publicznie dostępne, zanonimizowane dane FAERS oraz publicznie dostępne statystyki podsumowujące GWAS, nie było wymagane dodatkowe zatwierdzenie etyczne.
Środowisko programistyczne i implementacja przepływu pracy
Wszystkie analizy przeprowadzono w wersji R 4.3.2. Import, czyszczenie, łączenie i tabulacja danych FAERS były wykonywane przy użyciu workflow zarządzania danymi opartych na R. Tabele danych były importowane za pomocą funkcji takich jak data.table::fread() lub readr::read_delim(), scalane za pomocą CASEID i PRIMARYID oraz przetwarzane za pomocą funkcji dplyr. Statystyki opisowe oraz 2 × 2 tabele awaryjne zostały wygenerowane przy użyciu niestandardowych skryptów R.
Analizy randomizacji Mendelowskiej przeprowadzono przy użyciu wersji TwoSampleMR 0.5.6. Instrumenty ekspozycji były wyodrębniane z użyciem progu istotności P < 1 × 10⁻5 lub formatowane na podstawie statystyk podsumowujących GWAS, wykorzystując struktury wejściowe zgodne z TwoSampleMR. Zgrupowanie instrumentów wykonano za pomocą clump_data() przy clump_r2 = 0,001 i clump_kb = 10 000. Dane końcowe były wyodrębniane lub formatowane za pomocą extract_outcome_data() lub read_outcome_data(), w zależności od formatu źródłowego. Zestawy danych ekspozycji i wyników zostały zharmonizowane za pomocą harmonise_data(). Estymacje przyczynowe zostały wygenerowane za pomocą mr() przy użyciu następujących metod randomizacji mendelowskiej: MR-Egger, mediana ważona, odwrotna wariancja ważona, tryb prosty oraz tryb ważony. Heterogeniczność oceniano za pomocą mr_heterogeneity(), a poziomą plejotropię za pomocą mr_pleiotropy_test(). Wszystkie zbiory danych zostały zaimportowane, oczyszczone, ujednolicone i przeanalizowane w tym środowisku programowym, aby zapewnić spójny i powtarzalny przepływ pracy analitycznej.
Analiza farmakoczujności
Analizy opisowe zostały wykorzystane do podsumowania działań niepożądanych związanych ze złamaniem biodra związanego z pięcioma lekami. Następnie przeprowadzono analizy detekcji sygnału przy użyciu czterech algorytmów dysproporcji, w tym wskaźnika raportowania szans (ROR), proporcjonalnego wskaźnika raportowania (PRR), wieloelementowego zmniejszacza gamma Poissona (MGPS) oraz sieci neuronowej propagacji pewności bayesowskiej (BCPNN). Kryteria dla czterech głównych algorytmów przedstawiono w Tabeli 246.
Analiza randomizacji mendelowskiej
Statystyki podsumowujące TNF-α, sTNFR1 i sTNFR2 zostały wyodrębnione jako zbiory danych o ekspozycji, a statystyki podsumowujące złamania biodra zostały wyodrębnione jako zbiór wyników. Analizy ograniczały się do europejskich zbiorów danych o pochodzeniu, gdy były dostępne, aby zmniejszyć błędy stratyfikacyjne populacji.
Aby zminimalizować błędy spowodowane zaburzeniami sprzężenia i słabymi instrumentami, zastosowano następujące kryteria: próg istotności w całym genomie P < 1 × 10⁻5, próg nierównowagi sprzężeń r2 < 0,001, okno zgrupowania 10 000 kb oraz statystyka F > 20. Statystyka F została obliczona dla każdej zachowanej zmiennej instrumentalnej jako beta2/se 2 w celu oceny siły instrumentu. SNP o F-statystyce ≤ 20 zostały wykluczone z analiz końcowych.
Po selekcji SNP zestawy danych ekspozycji i wyników zostały zharmonizowane, aby wyrównać allele efektu. Podczas harmonizacji allele efektowe i inne allele były wyrównane pomiędzy zbiorami danych ekspozycji i wyników. SNP z niekompatybilnymi allelami zostały usunięte, a palindromiczne SNP o niejednoznacznych częstotliwościach alleli wykluczono, gdy nie udało się określić orientacji nici. Po harmonizacji sprawdzono zachowane SNP, aby potwierdzić, że współczynniki beta odpowiadają temu samemu allelowi efektu w obu zbiorach danych. Liczba SNP zachowanych po zgrupowaniu i harmonizacji była rejestrowana dla każdego ekspozycji jako punkt kontrolny pośredniej powtarzalności.
Zastosowano pięć metod randomizacji Mendelowskiej, w tym MR-Eggera, ważoną medianę, odwrotną wariancję ważoną, tryb prosty oraz ważony. Potencjalna heterogeniczność zmiennych instrumentalnych oceniano za pomocą testu Q Cochrana, a P < 0,05 uznano za oznakę istotnej heterogeniczności. Potencjalna pozioma plejotropia została oceniona za pomocą przecięcia MR-Eggera, a P < 0,05 wskazywało na plejotropię, co sugerowało obniżoną wiarygodność szacunku przyczynowego47. Analizy te generowały szacunki efektów przyczynowych wraz ze statystykami heterogeniczności i plejotropii dla każdej ekspozycji.
Pośrednie punkty kontrolne dla powtarzalności
Po każdym głównym etapie przetwarzania rejestrowano punkty kontrolne pośrednie, aby zapewnić powtarzalność przepływu pracy. W przepływie pracy FAERS punkty kontrolne obejmowały liczbę importowanych rekordów DEMO, liczbę unikalnych rekordów po deduplikacji, liczbę raportów zawierających złamanie biodra jako docelowe zdarzenie niepożądane, liczbę zgłoszeń dotyczących pięciu inhibitorów TNF oraz końcową liczbę kwalifikujących się zgłoszeń, w których inhibitory TNF zostały zarejestrowane jako główne podejrzane leki. W procesie randomizacji Mendelowskiej punkty kontrolne obejmowały liczbę wyodrębnionych SNP dla każdej ekspozycji, liczbę SNP zachowanych po zgrupowaniu nierównowagi sprzężenia, liczbę dostępnych SNP w zbiorze danych, liczbę SNP zachowanych po harmonizacji oraz końcową liczbę zmiennych instrumentalnych używanych w każdej analizie randomizacji mendelowskiej.
Raportowanie statystyczne
Wyniki ciągłe zostały przedstawione wraz z odpowiadającymi szacunkami efektów, 95% przedziałami ufności (95% CI) oraz wartościami P . O ile nie określono inaczej, istotność statystyczna była definiowana jako dwustronna P < 0,05. Do analizy farmakowigilancji podawano opisowe liczby i szacunki nieproporcjonalności dla każdego indywidualnego inhibitora TNF oraz dla grupy grupowej inhibitorów TNF. Dla analizy randomizacji mendelowskiej podano szacunki przyczynowe, błędy standardowe, przedziały ufności 95%, wartości P , statystyki heterogeniczności, wyniki testów pleiotropii oraz liczbę zachowanych SNP dla każdej ekspozycji.