$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
We hebben ethische goedkeuring en geïnformeerde toestemming verkregen van de Biomedical Research Ethics Committee van het First Affiliated Hospital van de Nanchang University. Ethisch Nummer: (2025)CDYFYYLK(08-007).
MR-analyse
Data-opvraging
De plasma pQTL-gegevens zijn verkregen uit de studie van Zheng et al.14, die vijf GWAS-datasets 15,16,17,18,19 integreerde, en uit de studie van Ferkingstad et al. De inclusiecriteria voor de gegevens waren als volgt: (i) genoombrede significante associaties (p < 5 × 10⁻⁸); en (ii) plasma-eiwitten als potentiële therapeutische doelwitten voor OA. Het studieontwerp wordt samengevat in Figuur 1. Ten eerste identificeerden we kandidaat-therapeutische doelen met behulp van GWAS-gegevens van IEU OpenGWAS en plasma pQTL-gegevens uit de studies van Zheng14 en Ferkingstad20 (Supplemental Table S1 en Supplemental Table S2). Steiger-filtering en fenotype-scanning werden vervolgens uitgevoerd om de robuustheid van de resultaten te valideren. De IEU OpenGWAS (https://gwas.mrcieu.ac.uk/) werd gebruikt om samenvattende statistieken te verkrijgen voor heup- of knie-OA (n = 417.596), knie-OA (n = 403.124) en heup-OA (n = 393.873)21.
SNP-filtercommando's
SNP's met genoombrede significantie (p < 5 × 10⁻⁸) werden onderworpen aan een klonteringsproces (r² < 0,001, F-statistieken > 10, venstergrootte = 10.000 kb) vóór MR-analyse.
MR-analyse
Om potentiële geneesmiddeldoelen te onderzoeken, werd MR-analyse uitgevoerd met plasma-eiwitten als blootstelling en OA als uitkomst, geïmplementeerd via het "TwoSampleMR"-pakket in R (v4.3.1). Wanneer er slechts één enkele pQTL beschikbaar was voor een eiwit, werd de Wald-verhouding gebruikt; verder werd inverse variantiegewogen MR (MR-IVW) toegepast, gevolgd door heterogeniteits- en pleiotropiebeoordelingen. Bonferronicorrectie werd gebruikt om meerdere tests te kunnen toelichten, met een drempel van p < 5,63 × 10⁻⁵ voor het prioriteren van eiwitten.
Steigerfiltering en fenotype-scanning
Om omgekeerde causaliteit te beoordelen, voerden we Steiger-filtering uit. Een resultaat van "TRUE" met p < 0,05 gaf geen omgekeerde causaliteit aan. Fenotype-scanning werd uitgevoerd met behulp van LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 om associaties van pQTL's met andere kenmerken te onderzoeken. De drempels waren R² = 0,1 en een venster van ±500.000 basisparen. Pleiotrope effecten werden toegekend aan pQTL's die beide van de volgende punten vervulden: (i) genoombrede significante associatie (p < 5 × 10⁻⁸), en (ii) associatie met bekende OA-risicofactoren.
Fenoombrede associatiestudie
Om rekening te houden met genpleiotropie en off-target effecten, hebben we een fenoombrede associatiestudie (PheWAS) uitgevoerd met behulp van het AstraZeneca PheWAS Portal (https://azphewas.com/), dat 15.500 binaire fenotypen en 1.500 continue fenotypen bevat van ~450.000 UK Biobank-deelnemers(23). Drempels werden ingesteld op standaardwaarden om valse positieven te minimaliseren.
Eiwit-eiwitinteractie (PPI) netwerk
Om interacties tussen potentiële eiwitdoelen die door MR zijn geïdentificeerd te visualiseren, gebruikten we GeneMANIA (https://genemania.org/) voor eiwit-eiwitinteractie-analyse en resultaatvisualisatie24.
Verrijkingsanalyse
Om biologische relevantie te onderzoeken, voerden we verrijkingsanalyses uit met behulp van bio-informatica-instrumenten van https://www.bioinformatics.com.cn voor data-analyse en visualisatie.
Transcriptomische workflow
Totaal RNA werd gewonnen met behulp van de RNA-extractiereagenskit volgens de richtlijnen van de fabrikant. De kwaliteit van RNA werd beoordeeld met behulp van een geautomatiseerd RNA-kwaliteitsbeoordelingssysteem; alleen monsters met RIN ≥7,0 werden gebruikt. De kwaliteit werd bevestigd door RNase-vrije agarosegelelektroforese (1,5% gel). Eukaryote mRNA werd verrijkt met behulp van Oligo(dT)-kralen; prokaryot mRNA werd verrijkt met behulp van de RNA elimination Magnetic Kit. mRNA werd gefragmenteerd (200-700 nt) en omgezet naar cDNA met behulp van de RNA Library Prep Kit. De cDNA-bibliotheek werd eindgerepareerd, A-tailed, geligeerd aan adapters, gezuiverd met DNA-zuiverende magnetische kralen (1.0×) en PCR-amplificeerd. De sequencing werd uitgevoerd op een high-throughput next-generation sequencingplatform. Differentieel expressieve genen werden gedefinieerd door log₂FC > 1 en gecorrigeerd p < 0,05.
Netwerkfarmacologie
Om potentiële geneesmiddelen voor doeleiwitten te identificeren, gebruikten we BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Een scoregrens van 0,74 (LR = 32,5) werd gebruikt om bekende en voorspelde verbindingen te selecteren. Kruidencomponenten werden uit TCMSP (https://old.tcmsp-e.com/index.php) gehaald en gefilterd met OB > 30% en DL > 0,1826.
Moleculaire koppeling
Moleculaire koppeling werd gebruikt om bindingsinteracties te beoordelen. Eiwitstructuren werden uit de PDB gehaald (https://www.rcsb.org/). UCSF Chimera werd gebruikt om structuren voor te bewerken door liganden en oplosmiddelen te verwijderen. AutoDock Tools werd gebruikt om Gasteiger-ladingen te berekenen en de centra en afmetingen van de doos te definiëren. Geneesmiddelstructuren werden verkregen uit PubChem (https://pubchem.ncbi.nlm.nih.gov/) en op vergelijkbare wijze voorbewerkt. Het docken werd uitgevoerd met AutoDock Vina. Afmetingen van de doos varieerden per doelwit. Bindingsaffiniteiten werden berekend en de resultaten werden zichtbaar in UCSF Chimera.