$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Uzyskaliśmy etyczną zgodę i świadomą zgodę od Komitetu Etyki Badań Biomedycznych Pierwszego Szpitala Afiliowanego Uniwersytetu Nanchang. Numer etyki: (2025)CDYFYYLK(08-007).
Analiza MR
Pobieranie danych
Dane pQTL plazmy zostały uzyskane z badania Zheng i in.14, które zintegrowało pięć zestawów danych GWAS 15,16,17,18,19, oraz z badania Ferkingstada i in. Kryteria włączenia do danych były następujące: (i) istotne powiązania w całym genomie (str. < 5 × 10⁻⁸); oraz (ii) białka osocza jako potencjalne cele terapeutyczne dla OA. Projekt badania jest podsumowany na Rysunku 1. Najpierw zidentyfikowaliśmy kandydatów na cele terapeutyczne, korzystając z danych GWAS z IEU OpenGWAS oraz danych pQTL osocza z badań przeprowadzonych przez Zheng14 i Ferkingstad20 (Tabela Uzupełniająca S1 i Tabela Uzupełniająca S2). Następnie przeprowadzono filtrowanie Steigera i skanowanie fenotypowe, aby potwierdzić odporność wyników. IEU OpenGWAS (https://gwas.mrcieu.ac.uk/) został użyty do uzyskania podsumowania statystyk dotyczących odrzewczenia stawu stawowego lub biodrowego (n = 417 596), odrzewczenia kolan (n = 403 124) oraz oddechu stawu biodrowego (n = 393 873)21.
Polecenia filtrowania SNP
SNP o istotności ogólnogenomowej (p < 5 × 10⁻⁸) zostały poddane procesowi zgrudkowania (r² < 0,001, F-statystyka > 10, rozmiar okna = 10 000 kb) przed analizą MR.
Analiza MR
Aby zbadać potencjalne cele leków, przeprowadzono analizę MR z wykorzystaniem białek osocza jako ekspozycji oraz OA jako wyniku, zaimplementowanej za pomocą pakietu "TwoSampleMR" w języku R (v4.3.1). Gdy dla białka dostępny był tylko jeden pQTL, stosowano stosunek Walda; w przeciwnym razie stosowano odwrotną wariancję ważoną MR (MR-IVW), a następnie oceniano heterogeniczność i plejotropię. Korekta Bonferroniego została użyta do uwzględnienia wielu testów, z progiem p < 5,63 × 10⁻⁵ dla priorytetyzacji białek.
Filtrowanie Steigera i skanowanie fenotypowe
Aby ocenić odwrotną przyczynowość, przeprowadziliśmy filtrowanie Steigera. Wynik "PRAWDZIWY" z p < 0,05 wskazywał na brak odwrotnej przyczynowości. Skanowanie fenotypowe przeprowadzono za pomocą LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 w celu zbadania powiązania pQTL z innymi cechami. Progi wynosiły R² = 0,1 oraz okno par bazowych ±500 000. Efekty plejotropowe przypisano pQTL spełniającym oba następujące kryteria: (i) istotne genomowe powiązanie (p < 5 × 10⁻⁸) oraz (ii) powiązanie z znanymi czynnikami ryzyka OA.
Badanie skojarzeń na poziomie fenomowym
Aby uwzględnić pleiotropię genów i efekty poza celem, przeprowadziliśmy badanie asocjacji fenomowej (PheWAS) z wykorzystaniem portalu AstraZeneca PheWAS (https://azphewas.com/), który zawiera 15 500 fenotypów binarnych i 1 500 ciągłych fenotypów od ~450 000 uczestników UK Biobank(23). Progi zostały ustawione na wartości domyślne, aby zminimalizować fałszywe alarmy.
Sieć interakcji białko-białko (PPI)
Aby zwizualizować interakcje między potencjalnymi celami białkowymi zidentyfikowanymi za pomocą MR, użyliśmy GeneMANIA (https://genemania.org/) do analizy interakcji białko-białko oraz wizualizacji wyników24.
Analiza wzbogacenia
Aby zbadać istotność biologiczną, przeprowadziliśmy analizę wzbogacenia z wykorzystaniem narzędzi bioinformatycznych z https://www.bioinformatics.com.cn do analizy i wizualizacji danych.
Przepływ pracy transkryptomicznej
Całkowite RNA zostało wyekstrahowane za pomocą zestawu odczynnika do ekstrakcji RNA, zgodnie z wytycznymi producenta. Jakość RNA oceniano za pomocą zautomatyzowanego systemu oceny jakości RNA; użyto tylko próbek o RIN ≥7.0. Jakość potwierdzono elektroforezą agarozy wolną od RNazy (żel 1,5%). Eukariotyczne mRNA zostało wzbogacone za pomocą kulek oligo(dT); prokariotyczne mRNA zostało wzbogacone za pomocą zestawu magnetycznego eliminacyjnego RNA. mRNA zostało fragmentowane (200-700 nt) i przekształcone w cDNA za pomocą zestawu RNA Library Prep Kit. Biblioteka cDNA była naprawiana na końcach, była ogoniasta w kształcie A, podwiązywała się z adapterami, oczyszczała za pomocą magnetycznych kulek oczyszczających DNA (1,0×) i wzmacniała PCR. Sekwencjonowanie wykonywano na wysokoprzepustowej platformie sekwencjonowania nowej generacji. Geny różnicowo ekspresyjne zostały zdefiniowane przez log₂FC > 1 i skorygowane p < 0,05.
Farmakologia sieciowa
Aby zidentyfikować potencjalne leki dla białek docelowych, użyliśmy BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Do wyboru znanych i przewidywanych związków używano progu wyniku 0,74 (LR = 32,5). Składniki ziołowe pobrano z TCMSP (https://old.tcmsp-e.com/index.php) i przefiltrowano OB > 30% oraz DL > 0,1826.
Dokowanie molekularne
Do oceny interakcji wiązań zastosowano dokowanie molekularne. Struktury białkowe zostały pobrane z PDB (https://www.rcsb.org/). Chimera UCSF była używana do wstępnego przetwarzania struktur poprzez usuwanie ligandów i rozpuszczalników. AutoDock Tools był używany do obliczania ładunków Gasteigera oraz definiowania środkow i rozmiarów pudełek. Struktury leków zostały uzyskane z PubChem (https://pubchem.ncbi.nlm.nih.gov/) i przetworzone w podobny sposób. Dokowanie odbywało się za pomocą AutoDock Vina. Wymiary pudełka różniły się w zależności od celu. Obliczono powinowactwa wiązania, a wyniki zobrazowano w chimerze UCSF.