$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
FAERS veri kaynakları
Bu çalışmanın gerçek veri verileri FAERS veritabanından (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html) alınmıştır. Bu, kamuya açık ve anonimleştirilmiş bir veritabanıdır, bu nedenle bu çalışma için etik onay gerekmemiştir. Dahil edilen beş ilaç hakkında bilgiler Tablo 1'de gösterilmiştir. Arama, önce yan etkileri kalça kırığıyla sınırlandırarak gerçekleştirildi ve arama süresi 2014'ün 1. çeyreğinden 2023'ün 4. çeyreğine kadar uzandı. Güvenilir ve istikrarlı veri sağlamak için çalışma, bildirilen yan olayların terimlerini MedDRA Sözlüğü 26.142 versiyonu üzerinden standartlaştırdı. FAERS'in 2014'ün birinci çeyreğinden 2023'ün 4. çeyreğine kadar olan ASCII dosyaları analiz için indirildi ve içe aktarıldı. Çıkarılan FAERS tabloları DEMO, DRUG, REAC, THER, RPSR ve OUTC'yi içeriyordu. Bu masalar tüm bölgelerde birleştirilerek eleme yapıldı. Raporlar, demografik, ilaç, reaksiyon, tedavi, raporlayıcı ve sonuç bilgileri arasında tutarlılık sağlamak için CASEID ve PRIMARYID kullanılarak bağlandı. Hedef ilaçlar arasında etanercept, adalimumab, infliksimab, certolizumab pegol ve golimumab bulunuyordu. DRUG tablosunda ilaç isimleri, metni büyük harflere dönüştürerek, fazladan boşlukları kaldırarak ve gerektiğinde yazım varyantlarını kontrol ederek standartlaştırıldı. Hedef ilaçlar, DRUG tablosunda standart jenerik isimlerle tanımlandı ve ilaç rolü kısıtlaması ROLE_COD = "PS" kullanılarak gerçekleştirildi; bu da birincil şüpheliilaç 43'ü gösteriyordu. Özel çalışma tarama tekniği yol haritası Şekil 2'de gösterilmiştir. Veri temizliği ve taramasından sonra, uygun kalça kırığı raporlarının benzersiz ve analiz edilebilir bir veri seti elde edildi ve sonraki analizler için yapıldı.
Tekrarlanan FAERS raporları sinyal tespitinden önce kaldırıldı. Kopyalar CASEID ve PRIMARYID verilerine göre belirlendi. Birden fazla rapor aynı CASEID'i paylaştığında, en son rapor FDA_DT'ye göre korundu. Birden fazla raporda aynı CASEID ve FDA_DT varsa, en yüksek PRIMARYID içeren rapor korunurdu. Deduplication'dan sonra, her CASEID nihai analitik veri setine yalnızca bir kayıt katkı sağladı. Raporlar, aşağıdaki kriterlerin tamamını karşılıyorsa dahil edildi: raporlama tarihi 2014 birinci çeyreği ile 2023 4. çeyreği arasında; yan etki "Kalça kırığı" olarak kodlandı; beş TNF inhibitöründen en az biri DRUG tablosunda kaydedildi; ve uyuşturucu rolü birincil şüpheli olarak kodlanmıştı. Raporlar, tekrarlanmış kayıtlar, geçerli CASEID veya PRIMARYID bilgisi yoksa, karşılık gelen DRUG veya REAC girişi yoksa, hedef yan etkinliği içermiyorsa veya hedef TNF inhibitörü yalnızca eşzamanlı veya ikincil şüpheli ilaç olarak listeleniyorsa hariç tutuluyordu.
Mendel rastgeleleştirmesi için GWAS veri kaynakları
Bu çalışma için Mendel randomizasyonunda TNF-α maruziyet verileri, IEU OpenGWAS veritabanından (https://gwas.mrcieu.ac.uk/) alındı ve Suhre K ve arkadaşlarının çalışmasından GWAS ID prot-c-3722_49_2 ile alındı. Çalışma popülasyonu Avrupa kökenliydi ve SNP sayısı 501.42844 idi.
Bu çalışma için mendel randomizasyonunda sTNFR1 için maruziyet verileri, IEU OpenGWAS veritabanından (https://gwas.mrcieu.ac.uk/) alındı ve Suhre K ve diğerlerinin çalışmasından GWAS ID prot-c-2654_19_1 ile alındı. Çalışma popülasyonu Avrupa kökenliydi ve SNP sayısı 501.42844 idi.
Bu çalışma için Mendel randomizasyonunda sTNFR2 için maruziyet verileri, IEU OpenGWAS veritabanından (https://gwas.mrcieu.ac.uk/) alındı ve Suhre K ve ark. çalışmasından GWAS ID prot-c-3152_57_1 olarak alındı. Çalışma popülasyonu Avrupa kökenliydi ve SNP sayısı 501.42844 idi.
Kalça kırığı için GWAS ID GCST90161240 sonuç verileri, GWAS Kataloğu'na (https://www.ebi.ac.uk/gwas/studies/GCST90161240) yerleştirilmiştir ve 11.516 kalça kırığı vakası ile 723.838 kontrol45'i içeren büyük ölçekli bir GWAS'ın meta-analizinden elde edilen verilerdir. Hastalık sınıflandırması, Uluslararası Hastalıklar Sınıflandırması (ICD; ICD-10 kodları S72.0–S72.2 ve ICD-9 kodu 820).
Orijinal GWAS çalışmalarında etik onay ve bilgilendirilmiş onay alınmıştı. Mevcut çalışma, kamuya açık ve anonimleştirilmiş FAERS verileri ile kamuya açık GWAS özet istatistiklerini kullandığı için ek etik onay gerekmemiştir.
Yazılım ortamı ve iş akışı uygulaması
Tüm analizler R sürüm 4.3.2 kullanılarak gerçekleştirildi. FAERS veri ithalatı, temizliği, birleştirme ve tablolama işlemleri R tabanlı veri yönetim iş akışları kullanılarak gerçekleştirildi. Veri tabloları data.table::fread() veya readr::read_delim(gibi fonksiyonlarla içe aktarıldı, CASEID ve PRIMARYID kullanılarak birleştirildi ve dplyr fonksiyonları kullanılarak işlendi. Betimleyici istatistikler ve 2 × 2 kontendens tablosu özel R betikleri kullanılarak oluşturuldu.
Mendel randomizasyon analizleri TwoSampleMR sürüm 0.5.6 kullanılarak gerçekleştirildi. Pozlama cihazları, P < 1 × 10⁻5 anlamlılık eşiği kullanılarak çıkarıldı veya TwoSampleMR uyumlu giriş yapılarıyla GWAS özet istatistiklerinden biçimlendirildi. Enstrüman topaklama, clump_r2 = 0.001 ve clump_kb = 10.000 olan clump_data() kullanılarak gerçekleştirildi. Sonuç verileri, kaynak formatına bağlı olarak extract_outcome_data() veya read_outcome_data() kullanılarak çıkarıldı veya biçimlendirildi. Maruz kalma ve sonuç veri setleri harmonise_data() kullanılarak uyumlaştırıldı. Nedensel tahminler, mr() kullanılarak aşağıdaki Mendel rastgeleleştirme yöntemleriyle oluşturuldu: MR-Egger, ağırlıklı medyan, ters varyans ağırlıklı, basit mod ve ağırlıklı mod. Heterojenlik mr_heterogeneity(), yatay pleiotropi ise mr_pleiotropy_test() kullanılarak değerlendirildi. Tüm veri setleri, tutarlı ve tekrarlanabilir bir analitik iş akışı sağlamak için bu yazılım ortamında ithal edildi, temizlendi, uyumlaştırıldı ve analiz edildi.
Farmakovigilyans analizi
Beş ilaçla ilişkili kalça kırığına bağlı yan etkileri özetlemek için betimleyici analizler kullanıldı. Sinyal algılama analizleri, raporlama oranı (ROR), orantılı raporlama oranı (PRR), çok öğeli gamma Poisson küçültücü (MGPS) ve Bayes güven yayılım sinir ağı (BCPNN) dahil olmak üzere dört orantısızlık algoritması kullanılarak gerçekleştirildi. Dört ana algoritmanın kriterleri Tablo 246'da gösterilmiştir.
Mendel rastgeleleştirme analizi
TNF-α, sTNFR1 ve sTNFR2 için özet istatistikler maruziyet veri setleri olarak çıkarıldı, kalça kırığı özet istatistikleri ise sonuç veri seti olarak çıkarıldı. Analizler, popülasyon tabakalaşma yanlılığını azaltmak için mevcut olduğunda Avrupa atalarının veri setleriyle sınırlandırıldı.
Bağlantı dengesizliği ve zayıf aletlerden kaynaklanan yanlılığı en aza indirmek için şu kriterler uygulandı: genom genelinde anlamlılık eşiği P < 1 × 10⁻5, bağlantı dengesizliği eşiği r2 < 0.001, 10.000 kb topaklanma penceresi ve F-istatistiği > 20. F-istatistik, korunan her enstrüman değişkeni için beta2/se2 olarak hesaplanarak enstrüman gücünü değerlendirdi. F-istatistikli ≤ 20'ye sahip SNP'ler aşağı akış analizlerinden dışlandı.
SNP seçiminden sonra, poz ve sonuç veri setleri etki allellerini hizalamak için uyumlandırıldı. Uyumlandırma sırasında, etki allelleri ve diğer alleller maruz kalma ve sonuç veri setleri arasında hizalanmıştır. Uyumsuz allellere sahip SNP'ler çıkarıldı ve belirsiz allel frekanslarına sahip palindromik SNP'ler iplik yönelimi belirlenemediğinde hariç tutuldu. Uyumlaştırmadan sonra, korunan SNP'ler, beta katsayılarının her iki veri setinde de aynı etki alleliyle eşleştiği doğrulanmak için kontrol edildi. Topaklanma ve uyumlaştırmadan sonra tutulan SNP sayısı, her bir poz için ara tekrarlanabilirlik kontrol noktası olarak kaydedildi.
MR-Egger, ağırlıklı medyan, ters varyans ağırlıklı, basit mod ve ağırlıklı mod dahil olmak üzere beş Mendel rastgeleleştirme yöntemi uygulandı. Enstrümantal değişkenlerin potansiyel heterojenliği Cochran'ın Q testi kullanılarak değerlendirildi ve P < 0.05 anlamlı heterojenlik göstergesi olarak değerlendirildi. Potansiyel yatay pleiotropi MR-Egger kesimiyle değerlendirildi ve P < 0.05 pleiotropiyi gösterdi; bu da nedenseltahminin 47 güvenilirliğinin azaldığını gösteriyor. Bu analizler, her maruziyet için nedensel etki tahminleri, heterojenlik ve pleiotropi istatistikleriyle birlikte üretti.
Tekrarlanabilirlik için ara kontrol noktaları
İş akışının tekrarlanabilirliğini sağlamak için her ana işlem adımından sonra ara kontrol noktaları kaydedildi. FAERS iş akışı için, kontrol noktaları arasında ithal edilen DEMO kayıt sayısı, deduplication sonrası benzersiz kayıt sayısı, kalça kırığını hedef olumsuz olay olarak içeren rapor sayısı, beş TNF inhibitörü ile ilgili rapor sayısı ve TNF inhibitörlerinin birincil şüpheli ilaç olarak kaydedildiği uygun raporların nihai sayısı yer alıyordu. Mendel rastgeleleştirme iş akışı için, kontrol noktaları arasında her maruziyet için çıkarılan SNP sayısı, bağlantı dengesizliği topaklandıktan sonra tutulan SNP sayısı, sonuç veri setinde mevcut SNP sayısı, uyumlaştırmadan sonra tutulan SNP sayısı ve her Mendel rastgeleleştirme analizinde kullanılan son araçsal değişken sayısı yer alıyordu.
İstatistiksel raporlama
Sürekli sonuçlar, karşılık gelen etki tahminleri, %95 güven aralıkları (%95 CI) ve P değerleri ile rapor edilmiştir. Aksi belirtilmedikçe, istatistiksel anlamlılık iki taraflı bir P < 0.05 olarak tanımlanmıştır. Farmakovijilans analizi için, her bireysel TNF inhibitörü ve birleştirilmiş TNF inhibitörü grubu için betimleyici sayılar ve orantısızlık tahminleri rapor edildi. Mendel randomizasyon analizi için, nedensel tahminler, standart hatalar, %95 güven aralıkları, P değerleri, heterojenlik istatistikleri, pleiotropi test sonuçları ve her maruziyet için korunan SNP sayısı rapor edilmiştir.