$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Nanchang Üniversitesi'nin İlk Bağlı Hastanesi'nin Biyomedikal Araştırma Etik Komitesi'nden etik onay ve bilgilendirilmiş onay aldık. Etik Numarası: (2025) CDYFYYLK(08-007).
MR analizi
Veri alımı
Plazma pQTL verileri, Zheng ve ark.14 tarafından yapılan çalışmadan elde edilmiştir; bu çalışmada beş GWASveri seti 15,16,17,18,19 entegre edilmiştir ve Ferkingstad ve diğerlerinin çalışmasından elde edilmiştir. Verilerin dahil etme kriterleri şu şekildeydi: (i) genom çapında anlamlı ilişkiler (p < 5 × 10⁻⁸); ve (ii) OA için potansiyel terapötik hedefler olarak plazma proteinleri. Çalışma tasarımı Şekil 1'de özetlenmiştir. İlk olarak, IEU OpenGWAS'tan alınan GWAS verileri ve Zheng14 ileFerkingstad 20 (Ek Tablo S1 ve Ek Tablo S2) çalışmalarından plazma pQTL verilerini kullanarak aday tedavi hedefleri belirledik. Sonrasında Steiger filtreleme ve fenotip taraması yapılarak sonuçların sağlamlığını doğrulamak için yapıldı. IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), kalça veya diz OA (n = 417.596), diz OA (n = 403.124) ve kalça OA (n = 393.873)21 için özet istatistikler elde etmek için kullanıldı.
SNP filtreleme komutları
Genom çapında anlamlılığa sahip SNP'ler (p < 5 × 10⁻⁸) MR analizinden önce bir kümelenme sürecine maruz kalmıştır (r² < 0.001, F-istatistikleri > 10, pencere boyutu = 10.000 kb).
MR analizi
Potansiyel ilaç hedeflerini araştırmak için, MR analizi plazma proteinleri maruziyet olarak ve OA olarak kullanılarak, R (v4.3.1) içindeki "TwoSampleMR" paketi üzerinden uygulandı. Bir protein için yalnızca tek bir pQTL bulunduğunda, Wald oranı kullanıldı; aksi takdirde, ters varyans ağırlıklı MR (MR-IVW) uygulandı, ardından heterojenlik ve pleiotropi değerlendirmeleri yapıldı. Bonferroni düzeltmesi, çoklu testleri hesaba katmak için kullanıldı; proteinlerin önceliklendirilmesi için p < 5.63 × 10⁻⁵ eşik olarak belirlendi.
Steiger filtreleme ve fenotip taraması
Ters nedenselliği değerlendirmek için Steiger filtreleme yaptık. "TRUE" sonucu p < 0.05 ile ters nedensellik olmadığını göstermiştir. Fenotip taraması, pQTL'lerin diğer özelliklerle ilişkisini incelemek için LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 kullanılarak gerçekleştirilmiştir. Eşikler R² = 0.1 ve ±500.000 baz çifti penceresi idi. Pleiotropik etkiler, aşağıdaki iki şeyi karşılayan pQTL'lere atandı: (i) genom çapında anlamlı ilişki (p < 5 × 10⁻⁸) ve (ii) bilinen OA risk faktörleriyle ilişki.
Fenome çapında ilişki çalışması
Gen pleiotropisi ve hedef dışı etkileri hesaba katmak için, ~450.000 İngiliz Biobank katılımcısından 15.500 ikili fenotip ve 1.500 sürekli fenotip içeren AstraZeneca PheWAS Portalı (https://azphewas.com/) kullanılarak fenom çapında bir ilişki çalışması (PheWAS)yürüttük. Eşikler, yanlış pozitifleri en aza indirmek için varsayılan değerlere ayarlandı.
Protein - protein etkileşimi (PPI) ağı
MR ile tanımlanan potansiyel protein hedefleri arasındaki etkileşimleri görselleştirmek için protein-protein etkileşim analizi ve sonuç görselleştirme24 için GeneMANIA (https://genemania.org/) kullandık.
Zenginleştirme analizi
Biyolojik önemi araştırmak için, veri analizi ve görselleştirme için https://www.bioinformatics.com.cn biyoinformatik araçlarıyla zenginleştirme analizi yaptık.
Transkriptomik iş akışı
Toplam RNA, üreticinin yönergelerine uygun olarak RNA ekstraksiyon reaktif kiti kullanılarak çıkarıldı. RNA kalitesi otomatik bir RNA kalite değerlendirme sistemi kullanılarak değerlendirildi; yalnızca RIN ≥7.0 içeren örnekler kullanıldı. Kalite, RNaz-free agaroz jel elektroforezi (%1,5 jel) ile doğrulandı. Ökaryotik mRNA, Oligo(dT) boncukları kullanılarak zenginleştirildi; prokaryotik mRNA, RNA eliminasyon Manyetik Kiti kullanılarak zenginleştirildi. mRNA, parçalanmış (200-700 nt) ve RNA Library Prep Kit kullanılarak cDNA'ya dönüştürülmüştür. cDNA kütüphanesi uçtan tamir edilmiş, A-kuyruk yapılmış, adaptörlere bağlanmış, DNA arındırıcı manyetik boncuklarla (1.0×) arıtılmış ve PCR ile amplifikasyon yapılmıştır. Dizileme yüksek verimli bir yeni nesil dizileme platformunda gerçekleştirildi. Farklı olarak eksprese edilen genler log₂FC > 1 ile tanımlandı ve p < 0.05 ile ayarlandı.
Ağ farmakolojisi
Hedef proteinler için potansiyel ilaçları belirlemek için BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25 kullandık. Bilinen ve tahmin edilen bileşikler seçilmek için 0,74 (LR = 32,5) bir skor kesme noktası kullanıldı. Bitkisel bileşenler TCMSP'den (https://old.tcmsp-e.com/index.php) alındı ve OB >% 30 ile ve DL > 0.1826 ile filtrelendi.
Moleküler kenetleme
Moleküler kenetlenme bağlanma etkileşimlerini değerlendirmek için kullanıldı. Protein yapıları PDB'den alındı (https://www.rcsb.org/). UCSF Chimera, ligandlar ve çözücüler çıkararak yapıları ön işlemede kullanıldı. AutoDock Tools, Gasteiger yüklerini hesaplamak ve kutu merkezleri ile boyutlarını tanımlamak için kullanıldı. İlaç yapıları PubChem (https://pubchem.ncbi.nlm.nih.gov/)'dan alındı ve benzer şekilde ön işlenmiş oldu. Kenetlenme AutoDock Vina ile gerçekleştirildi. Kutu boyutları hedefe göre değişiyordu. Bağlanma eğilimleri hesaplandı ve sonuçlar UCSF Kimera'da görselleştirildi.