$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Abbiamo ottenuto l'approvazione etica e il consenso informato dal Comitato Etico della Ricerca Biomedica del Primo Ospedale Affiliato dell'Università di Nanchang. Numero etico: (2025)CDYFYYLK(08-007).
Analisi MR
Recupero dei dati
I dati plasmatici, pQTL, sono stati ottenuti dallo studio di Zheng et al.14, che ha integrato cinque dataset GWAS15, 16, 17, 18, 19, e dallo studio di Ferkingstad et al. I criteri di inclusione dei dati erano i seguenti: (i) associazioni significative a livello genomico (p < 5 × 10⁻⁸); e (ii) le proteine plasmatiche come potenziali bersagli terapeutici per l'OA. Il disegno dello studio è riassunto nella Figura 1. Per prima cosa, abbiamo identificato i target terapeutici candidati utilizzando dati GWAS dall'IEU OpenGWAS e dati plasmamici pQTL dagli studi di Zheng14 e Ferkingstad20 (Tabella Supplementare S1 e Tabella Supplementare S2). Sono stati quindi condotti filtri Steiger e scansioni fenotipiche per validare la robustezza dei risultati. L'OpenGWAS dell'IEU (https://gwas.mrcieu.ac.uk/) è stato utilizzato per ottenere statistiche riassuntive per l'OA dell'anca o del ginocchio (n = 417.596), l'OA del ginocchio (n = 403.124) e l'OA dell'anca (n = 393.873)21.
Comandi di filtraggio SNP
Gli SNP con significatività genomica (p < 5 × 10⁻⁸) sono stati sottoposti a un processo di aggregazione (r² < 0,001, statistica F > 10, dimensione finestra = 10.000 kb) prima dell'analisi RM.
Analisi MR
Per indagare potenziali bersagli farmacologici, è stata effettuata un'analisi RM utilizzando proteine plasmatiche come esposizioni e OA come risultato, implementata tramite il package "TwoSampleMR" in R (v4.3.1). Quando era disponibile un solo pQTL per una proteina, veniva utilizzato il rapporto di Wald; altrimenti, è stata applicata la MR (MR-IVW) ponderata per varianza inversa, seguita da valutazioni di eterogeneità e pleiotropia. La correzione di Bonferroni è stata utilizzata per tenere conto di test multipli, con una soglia di p < 5,63 × 10⁻⁵ per la priorità delle proteine.
Filtraggio di Steiger e scansione fenotipica
Per valutare la causalità inversa, abbiamo effettuato un filtraggio Steiger. Un risultato di "TRUE" con p < 0,05 non indicava causalità inversa. La scansione fenotipica è stata condotta utilizzando LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 per esaminare le associazioni dei pQTL con altri caratteristi. Le soglie erano R² = 0,1 e una finestra di ±500.000 coppie di base. Gli effetti pleiotropi sono stati assegnati ai pQTL che soddisfano sia i seguenti requisiti: (i) associazione significativa a livello genomico (p < 5 × 10⁻⁸), sia (ii) associazione con fattori di rischio noti per l'OA.
Studio di associazione a livello di fenomeno
Per spiegare la pleiotropia genica e gli effetti fuori bersaglio, abbiamo condotto uno studio di associazione a livello di fenomeno (PheWAS) utilizzando il Portale PheWAS di AstraZeneca (https://azphewas.com/), che contiene 15.500 fenotipi binari e 1.500 fenotipi continui da ~450.000 partecipanti al BiobankUK 23. Le soglie sono state impostate ai valori predefiniti per minimizzare i falsi positivi.
Rete di interazione proteina-proteina (PPI)
Per visualizzare le interazioni tra potenziali bersagli proteici identificati tramite RM, abbiamo utilizzato GeneMANIA (https://genemania.org/) per l'analisi delle interazioni proteina-proteina e la visualizzazionedei risultati 24.
Analisi dell'arricchimento
Per indagare la rilevanza biologica, abbiamo condotto un'analisi di arricchimento utilizzando strumenti bioinformatici di https://www.bioinformatics.com.cn per l'analisi e visualizzazione dei dati.
Flusso di lavoro trascrittomico
L'RNA totale è stato estratto utilizzando il kit di reagenti per estrazione di RNA seguendo le linee guida del produttore. La qualità dell'RNA veniva valutata utilizzando un sistema automatizzato di valutazione della qualità dell'RNA; sono stati utilizzati solo campioni con RIN ≥7.0. La qualità è stata confermata dall'elettroforesi del gel di agarosio senza RNasi (gel all'1,5%). L'mRNA eucariotico è stato arricchito utilizzando perle di Oligo(dT); l'mRNA procariotico è stato arricchito utilizzando il Magnetic Kit di Eliminazione dell'RNA. L'mRNA è stato frammentato (200-700 nt) e convertito in cDNA utilizzando il RNA Library Prep Kit. La libreria di cDNA veniva riparata all'estremità, con coda A, legata agli adattatori, purificata usando sfere magnetiche purificatrici del DNA (1,0×) e amplificata tramite PCR. Il sequenziamento veniva eseguito su una piattaforma di sequenziamento di nuova generazione ad alta produttività. I geni espressi differenzialmente sono stati definiti da log₂FC > 1 e aggiustati p < 0,05.
Farmacologia delle reti
Per identificare potenziali farmaci come proteine bersaglio, abbiamo utilizzato BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Un cutoff di punteggio di 0,74 (LR = 32,5) è stato utilizzato per selezionare i composti noti e previsti. I componenti erboristici sono stati prelevati dal TCMSP (https://old.tcmsp-e.com/index.php) e filtrati con OB > 30% e DL > 0,1826.
Aggancio molecolare
Il docking molecolare è stato utilizzato per valutare le interazioni di legame. Le strutture proteiche sono state prelevate dal PDB (https://www.rcsb.org/). UCSF Chimera veniva utilizzata per preprocessare strutture rimuovendo liganti e solventi. AutoDock Tools veniva utilizzato per calcolare le cariche di Gasteiger e definire i centri e le dimensioni delle scatole. Le strutture dei farmaci sono state ottenute da PubChem (https://pubchem.ncbi.nlm.nih.gov/) e pre-elaborate in modo simile. L'aggancio è stato effettuato utilizzando AutoDock Vina. Le dimensioni delle scatole variavano a seconda del bersaglio. Le affinità di legame sono state calcolate e i risultati sono stati visualizzati in UCSF Chimera.