$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Wir haben ethische Zustimmung und informierte Zustimmung vom Biomedical Research Ethics Committee des First Affiliated Hospital der Nanchang University eingeholt. Ethiknummer: (2025)CDYFYYLK(08-007).
MR-Analyse
Datenabruf
Die Plasma-pQTL-Daten stammten aus der Studie von Zheng et al.14, die fünf GWAS-Datensätze 15,16,17,18,19 integrierte, sowie aus der Studie von Ferkingstad et al. Die Einschlusskriterien für die Daten waren wie folgt: (i) genomweite signifikante Assoziationen (p < 5 × 10⁻⁸); und (ii) Plasmaproteine als potenzielle therapeutische Ziele für OA. Das Studiendesign ist in Abbildung 1 zusammengefasst. Zunächst identifizierten wir anhand von GWAS-Daten aus der IEU OpenGWAS und Plasma-pQTL-Daten aus den Studien von Zheng14 und Ferkingstad20 (Supplemental Table S1 und Supplemental Table S2) potenzielle therapeutische Ziele. Anschließend wurden Steiger-Filterung und Phänotyp-Scans durchgeführt, um die Robustheit der Ergebnisse zu validieren. Das IEU OpenGWAS (https://gwas.mrcieu.ac.uk/) wurde verwendet, um zusammenfassende Statistiken für Hüft- oder Knie-OA (n = 417.596), Knie-OA (n = 403.124) und Hüft-OA (n = 393.873)21 zu erhalten.
SNP-Filterbefehle
SNPs mit genomweiter Bedeutung (p < 5 × 10⁻⁸) wurden vor der MR-Analyse einem Klumpungsprozess (r² < 0,001, F-Statistik > 10, Fenstergröße = 10.000 kb) unterzogen.
MR-Analyse
Zur Untersuchung potenzieller Arzneimittelziele wurde eine MR-Analyse durchgeführt, wobei Plasmaproteine als Exposition und OA als Ergebnis verwendet wurden, implementiert über das "TwoSampleMR"-Paket in R (v4.3.1). Wenn für ein Protein nur ein einziges pQTL verfügbar war, wurde das Wald-Verhältnis verwendet; ansonsten wurde das inverse Varianzgewichtete MR (MR-IVW) angewendet, gefolgt von Heterogenitäts- und Pleiotropiebewertungen. Die Bonferroni-Korrektur wurde verwendet, um mehrfache Tests zu berücksichtigen, mit einem Schwellenwert von p < 5,63 × 10⁻⁵ für die Priorisierung von Proteinen.
Steiger-Filterung und Phänotyp-Scanning
Um die umgekehrte Kausalität zu bewerten, führten wir eine Steiger-Filterung durch. Ein Ergebnis von "TRUE" mit p < 0,05 zeigte keine umgekehrte Kausalität. Die Phänotyp-Scanning wurde mit LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 durchgeführt, um Assoziationen von pQTLs mit anderen Merkmalen zu untersuchen. Die Schwellenwerte waren R² = 0,1 und ein Fenster von ±500.000 Basispaaren. Pleiotrope Effekte wurden pQTLs zugeordnet, die beide folgende Kriterien erfüllen: (i) genomweite signifikante Assoziation (p < 5 × 10⁻⁸) und (ii) Assoziation mit bekannten OA-Risikofaktoren.
Phenomweite Assoziationsstudie
Um Genpleiotropie und Off-Target-Effekte zu berücksichtigen, führten wir eine phänomenweite Assoziationsstudie (PheWAS) mit dem AstraZeneca PheWAS Portal (https://azphewas.com/) durch, das 15.500 binäre Phänotypen und 1.500 kontinuierliche Phänotypen von ~450.000 UK Biobank-Teilnehmern enthält,23. Die Schwellenwerte wurden auf Standardwerte gesetzt, um Fehlalarme zu minimieren.
Protein-Protein-Interaktionsnetzwerk (PPI)
Um Interaktionen zwischen potenziellen durch MR identifizierten Proteinzielen zu visualisieren, verwendeten wir GeneMANIA (https://genemania.org/) für die Analyse der Protein-Protein-Interaktionen und die Ergebnisvisualisierung24.
Anreicherungsanalyse
Um die biologische Relevanz zu untersuchen, führten wir eine Anreicherungsanalyse mit Bioinformatik-Tools von https://www.bioinformatics.com.cn für Datenanalyse und Visualisierung durch.
Transkriptomischer Arbeitsablauf
Die Gesamt-RNA wurde mit dem RNA-Extraktionsreagenz-Set nach den Richtlinien des Herstellers extrahiert. Die RNA-Qualität wurde mittels eines automatisierten RNA-Qualitätsbewertungssystems bewertet; es wurden nur Proben mit RIN ≥7.0 verwendet. Die Qualität wurde durch RNase-freie Agarosegel-Elektrophorese (1,5 % Gel) bestätigt. Eukaryotische mRNA wurde mit Oligo(dT)-Perlen angereichert; Die prokaryotische mRNA wurde mit dem RNA-Elimination Magnetic Kit angereichert. mRNA wurde fragmentiert (200–700 nt) und mit dem RNA Library Prep Kit in cDNA umgewandelt. Die cDNA-Bibliothek wurde endlich repariert, mit A-Tail versehen, an Adapter gebunden, mit DNA-reinigenden magnetischen Perlen (1,0×) gereinigt und PCR-amplifiziert. Die Sequenzierung erfolgte auf einer hochdurchsatzfähigen Next-Generation-Sequenzierungsplattform. Differenziell exprimierte Gene wurden durch log₂FC > 1 definiert und p < 0,05 angepasst.
Netzwerkpharmakologie
Um potenzielle Medikamente für Zielproteine zu identifizieren, verwendeten wir BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Ein Score-Cutoff von 0,74 (LR = 32,5) wurde verwendet, um bekannte und vorhergesagte Verbindungen auszuwählen. Kräuterkomponenten wurden aus TCMSP (https://old.tcmsp-e.com/index.php) gewonnen und mit OB > 30 % und DL > 0,1826 gefiltert.
Molekulare Andockung
Molekulares Andocken wurde verwendet, um Bindungsinteraktionen zu bewerten. Proteinstrukturen wurden aus dem PDB gewonnen (https://www.rcsb.org/). UCSF Chimera wurde verwendet, um Strukturen vorzuverarbeiten, indem Liganden und Lösungsmittel entfernt wurden. AutoDock Tools wurde verwendet, um Gasteiger-Ladungen zu berechnen und Kastenzentren sowie -größen zu definieren. Die Arzneimittelstrukturen wurden aus PubChem (https://pubchem.ncbi.nlm.nih.gov/) gewonnen und ähnlich vorverarbeitet. Das Andocken erfolgte mit AutoDock Vina. Die Boxmaße variierten je nach Ziel. Bindungsaffinitäten wurden berechnet und Ergebnisse wurden in UCSF Chimera visualisiert.