$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Obtuvimos la aprobación ética y el consentimiento informado del Comité de Ética en Investigación Biomédica del Primer Hospital Afiliado de la Universidad de Nanchang. Número de ética: (2025)CDYFYYLK(08-007).
Análisis de RM
Recuperación de datos
Los datos plasmáticos de pQTL se obtuvieron del estudio de Zheng et al.14, que integró cinco conjuntos de datosGWAS: 15, 16, 17, 18, 19, y del estudio de Ferkingstad et al. Los criterios de inclusión para los datos fueron los siguientes: (i) asociaciones significativas a nivel genómico (p < 5 × 10⁻⁸); y (ii) proteínas plasmáticas como posibles objetivos terapéuticos para la OA. El diseño del estudio se resume en la Figura 1. En primer lugar, identificamos objetivos terapéuticos candidatos utilizando datos GWAS del OpenGWAS del IEU y datos plasmáticos de pQTL de los estudios de Zheng14 yFerkingstad 20 (Tabla Suplementaria S1 y Tabla Suplementaria S2). A continuación, se realizaron filtrados Steiger y escaneos fenotípicos para validar la robustez de los resultados. El OpenGWAS del IEU (https://gwas.mrcieu.ac.uk/) se utilizó para obtener estadísticas resumen de la OA de cadera o rodilla (n = 417.596), OA de rodilla (n = 403.124) y OA de cadera (n = 393.873)21.
Comandos de filtrado SNP
Los SNPs con significación genómica (p < 5 × 10⁻⁸) fueron sometidos a un proceso de agrupamiento (r² < 0,001, estadística F > 10, tamaño de ventana = 10.000 kb) antes del análisis por RM.
Análisis de RM
Para investigar posibles objetivos farmacológicos, se realizó un análisis de RM utilizando proteínas plasmáticas como exposiciones y OA como resultado, implementado mediante el encapsulado "TwoSampleMR" en R (v4.3.1). Cuando solo había un único pQTL disponible para una proteína, se utilizaba la proporción de Wald; de lo contrario, se aplicó la resonancia magnética ponderada por varianza inversa (MR-IVW), seguida de evaluaciones de heterogeneidad y pleiotropía. Se utilizó la corrección de Bonferroni para tener en cuenta múltiples pruebas, con un umbral de p < 5,63 × 10⁻⁵ para priorizar proteínas.
Filtrado de Steiger y escaneo fenotípico
Para evaluar la causalidad inversa, realizamos un filtrado Steiger. Un resultado de "TRUE" con p < 0,05 indicó que no hubo causalidad inversa. El escaneo fenotípico se realizó utilizando LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 para examinar las asociaciones de pQTLs con otros rasgos. Los umbrales eran R² = 0,1 y una ventana de ±500.000 pares de bases. Se asignaron efectos pleiotrópicos a pQTLs que cumplen ambos: (i) asociación significativa a nivel genómico (p < 5 × 10⁻⁸), y (ii) asociación con factores de riesgo conocidos de la artrosis (OA).
Estudio de asociación a nivel de fenómenos
Para explicar la pleiotropía génica y los efectos fuera del objetivo, realizamos un estudio de asociación a nivel de fenómenos (PheWAS) utilizando el Portal PheWAS de AstraZeneca (https://azphewas.com/), que contiene 15.500 fenotipos binarios y 1.500 fenotipos continuos de ~450.000 participantes del Biobancodel Reino Unido 23. Los umbrales se establecieron a valores predeterminados para minimizar falsos positivos.
Red de interacción proteína-proteína (IBP)
Para visualizar las interacciones entre posibles objetivos proteicos identificados por RM, utilizamos GeneMANIA (https://genemania.org/) para el análisis de interacción proteína-proteína y la visualizaciónde resultados 24.
Análisis de enriquecimiento
Para investigar la relevancia biológica, realizamos análisis de enriquecimiento utilizando herramientas de bioinformática de https://www.bioinformatics.com.cn para el análisis y visualización de datos.
Flujo de trabajo transcriptómico
El ARN total se extrajo utilizando el kit de reactivos de extracción de ARN siguiendo las directrices del fabricante. La calidad del ARN se evaluó mediante un sistema automatizado de evaluación de la calidad del ARN; solo se usaron muestras con RIN ≥7.0. La calidad se confirmó mediante electroforesis en gel de agarosa libre de RNasa (gel al 1,5%). El ARNm eucariota se enriqueció utilizando perlas de Oligo(dT); ARNm procariota se enriqueció utilizando el Kit Magnético de eliminación de ARN. El ARNm se fragmentó (200-700 nt) y se convirtió en ADNc utilizando el Kit de Preparación de Bibliotecas de ARN. La biblioteca de ADNc se reparaba por extremos, se llevaba cola en A, se ligaba a adaptadores, se purificaba usando perlas magnéticas purificadoras de ADN (1,0×) y se amplificaba por PCR. La secuenciación se realizó en una plataforma de secuenciación de nueva generación de alto rendimiento. Los genes expresados diferencialmente se definieron por log₂FC > 1 y ajustaron p < 0,05.
Farmacología de redes
Para identificar posibles fármacos como proteínas objetivo, utilizamos BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Se utilizó un corte de puntuación de 0,74 (LR = 32,5) para seleccionar compuestos conocidos y predichos. Se extrajeron componentes herbales de TCMSP (https://old.tcmsp-e.com/index.php) y se filtraron con OB > 30% y DL > 0,1826.
Acoplamiento molecular
Se utilizó acoplamiento molecular para evaluar las interacciones de unión. Las estructuras proteicas se extrajeron del PDB (https://www.rcsb.org/). UCSF Chimera se utilizó para preprocesar estructuras eliminando ligandos y disolventes. AutoDock Tools se utilizó para calcular cargas de Gasteiger y definir centros y tamaños de cajas. Las estructuras de los fármacos se obtuvieron de PubChem (https://pubchem.ncbi.nlm.nih.gov/) y se preprocesaron de forma similar. El acoplamiento se realizó usando AutoDock Vina. Las dimensiones de la caja variaban según el objetivo. Se calcularon afinidades de unión y se visualizaron resultados en UCSF Chimera.