Ética de aprobación y consentimiento para participar
Todos los diseños experimentales con animales se llevaron a cabo estrictamente conforme a las Medidas de Manejo de Animales de Laboratorio de la Provincia de Ningxia y fueron aprobados por el Comité de Ética del Hospital Popular de la Región Autónoma Hui de Ningxia, Universidad Médica de Ningxia.
Diseño del estudio
Este estudio investigó las relaciones causales entre 1.400 metabolitos plasmáticos y la ATI empleando una metodología de resonancia magnética de dos muestras. Para garantizar la validez de los IVs en la inferencia causal, deben cumplirse tres supuestos críticos: (1) la variación genética está directamente asociada con la exposición de interés; (2) la variación genética no se confunde por posibles factores de confusión entre la exposición y el resultado; y (3) la variación genética influye exclusivamente en el resultado a través de la exposición y no por vías alternativas.
Fuentes de datos GWAS para ATI
Los datos para el Estudio de Asociación Generalizada del Genoma (GWAS) sobre ATI se obtuvieron del proyecto IEU-OpenGWAS (https://gwas.mrcieu.ac.uk/), utilizando específicamente el conjunto de datos ebi-a-GCST90018895. El análisis GWAS abarcó a 357.658 individuos de ascendencia europea, comprendiendo 7.746 casos y 349.912 controles. Tras procedimientos de control de calidad, se identificaron un total de 19.072.056 polimorfismos de nucleótido único (SNPs). Los procedimientos de control de calidad incluyeron la eliminación de SNPs con frecuencia alelada menor<0,01, equilibrio Hardy-Weinberg P<1×10-⁶ y tasa de llamada genotípica <0,95. La estratificación poblacional se evaluó mediante análisis de componentes principales y se calculó el factor de inflación genómica (λ) para evaluar el sesgo sistemático.
Fuentes de datos GWAS para metabolitos plasmáticos
El conjunto de datos de Factores de Exposición incluye datos GWAS para 1.400 metabolitos plasmáticos, que consisten en 1.091 niveles de metabolitos y 309 proporciones de metabolitos, según se cataloga en la base de datos GWAS (https://www.ebi.ac.uk/gwas/studies/GCST90199621-90201020). Los análisis de RM se restringieron a personas de ascendencia europea. La información detallada sobre el conjunto de datos se proporciona en Materiales Suplementarios, Tabla S1.
Selección de IVs
En la identificación de IVs, el umbral de significación estadística de 19.072.056 SNPs asociados a metabolitos plasmáticos se estableció en P < 1 × 10-5. Para asegurar la independencia de los análisis, los SNPs se agregaron usando umbrales de r2 = 0,001 y kilobase (kb) = 10.000, excluyendo 19.037.213 SNPs con r2 > 0,001 y kb < 10.000 para eliminar posibles efectos del Desequilibrio de Enlace (LD). Los SNPs con valores de estadístico F < 10 fueron excluidos de análisis posteriores debido a su alta heterogeneidad, mientras que 34.843 SNPs que cumplían los criterios se mantuvieron como IVs válidos.
Animales y diseño experimental
Todos los experimentos con animales se llevaron a cabo de estricta conformidad con las Medidas de Manejo de Animales de Laboratorio de la Provincia de Ningxia y fueron aprobados por el Comité de Ética del Hospital Popular de la Región Autónoma Hui de Ningxia, Universidad Médica de Ningxia. Se adquirieron ratones BALB/c machos (8 semanas, grado SPF) a Beijing Vital River Laboratory Animal Technology Co., Ltd. (Pekín, China). Los ratones se aclimataron durante una semana en un entorno controlado (temperatura 23 ± 1°C, humedad 60 ± 10%), con libre acceso a comida y agua estándar, y un ciclo de 12 horas de luz y 12 horas de oscuridad. Tras la aclimatación, se estableció un modelo de lesión del tendón murino. Tras una incisión longitudinal cutánea, se expuso el tendón de Aquiles izquierdo y se realizó una transsección transversal estandarizada de grosor completo en el punto medio del tendón utilizando tijeras microquirúrgicas. El tendón fue reparado inmediatamente con suturas 6-0 no absorbibles bajo un estereomicroscopio. Inmediatamente después de la cirugía, se realizaron radiografías para verificar el establecimiento exitoso del modelo murino ATI. Simultáneamente, se adquirieron imágenes de rayos X para el grupo de control bajo condiciones anestésicas idénticas para servir como un controlen blanco 28. Todos los procedimientos se realizaron bajo anestesia general inducida por inyección intraperitoneal de pentobarbital sódico (50 mg·kg-1), y el análisis postoperatorio se administró buprenorfina (0,05 mg·kg-1, subcutánea) cada 12 horas durante 48 h29. Tras la inducción exitosa del modelo, los ratones fueron asignados aleatoriamente a los siguientes grupos: grupos de control, modelo y tratamiento que recibieron Dih en dosis bajas (10 mg·kg-1, LD-Hyd), medias (30 mg·kg-1, MD-Hyd) y altas (50 mg·kg-1, HD-Hyd). Cálculo del tamaño de la muestra: Basándose en datos preliminares que muestran una mejora del 40% en las puntuaciones de cicatrización con el tratamiento de metabolitos (DE = 15%), se determinó que n = 6 ratones por grupo proporcionarían un 80% de potencia para detectar esta diferencia en α = 0,05 (de dos colas). La aleatorización se realizó utilizando números aleatorios generados por ordenador por un investigador cegado ante la asignación de grupos. Los procedimientos experimentales y las evaluaciones de resultados fueron realizados por investigadores que no se veían ante las tareasgrupales 30.
Tinción por hematoxilina y eosina (HE)
Al día 7 tras el tratamiento, los ratones fueron eutanasiados bajo anestesia profunda inducida por inyección intraperitoneal de pentobarbital sódico (100 mg·kg-1), asegurando la pérdida total de reflejos antes de la recogida de tejido, en consecuencia. El tejido tendinoso lesionado, incluyendo una pequeña porción del tejido blando circundante, fue cuidadosamente extraído y procesado. Se extirpó un trozo de pielde 15 mm 2 centrado en el tendón de Aquiles. Las muestras de piel se fijaron en una solución de paraformaldehído al 4% a temperatura ambiente durante 48 horas, seguidas de deshidratación, incrustación de parafina y seccionamiento de 5 μm. Se evaluaron las características estructurales de la piel de la herida, la formación de tejido de granulación y la infiltración celular inflamatoria mediante tinción H&E. Para la evaluación histológica, se desparafinan secciones en xileno (2 × 10 minutos), se rehidrataron mediante etanol graduado (100%, 95%, 85% y 75%, 5 minutos cada una) y se enjuagaron en agua destilada. Las láminas se teñieron con hematoxilina Harris durante 5 minutos, se enjuagaron con agua del grifo corriente durante 5 minutos, se diferenciaron en alcohol ácido al 1% durante 10 segundos y se enazulearon en agua al 0,2% de amoníaco durante 30 segundos. Después del lavado, las secciones se contra-teñieron con eosina al 0,5% durante 2 minutos, seguidas de deshidratación, limpieza y montaje. Fiabilidad de la evaluación histológica: Para garantizar la precisión de la medición, el 20% de las muestras seleccionadas al azar fueron puntuadas de forma independiente por dos observadores cegados. La fiabilidad entre evaluadores se evaluó utilizando el coeficiente de correlación intraclase (ICC). La ICC para la organización del colágeno fue de 0,89 (IC 95%: 0,81-0,94), para la degeneración tisular de 0,92 (IC 95%: 0,86-0,96) y para la infiltración celular inflamatoria de 0,87 (IC 95%: 0,78-0,93), indicando una excelente fiabilidad. La fiabilidad intra-evaluadora se evaluó haciendo que un observador re-puntuara el 15% de las muestras tras un intervalo de dos semanas, dando valores ICC > 0,85 para todos los parámetros31.
Análisis de Western blot
La expresión proteica fue analizada por Western blot. La cuantificación de proteínas se realizó utilizando el kit de ensayo de concentración de proteínas de BCA. Cantidades iguales de proteína se separaron mediante electroforesis en gel de dodecil sulfato de sodio y poliacrilamida al 10% (SDS-PAGE) y se transfirieron a una membrana de nitrocelulosa. La membrana estuvo bloqueada durante 1 hora con leche desnatada al 5% en solución fisiológica con Tris-buffer que contenía 0,1% Tween-20 a temperatura ambiente. Las membranas se incubaron durante la noche a 4°C con anticuerpos primarios para VEGFA (1:1000, ab46154), TGF-β1 (1:1000, 21898-1-AP), IL-6 (1:1000, 26404-1-AP) y TNF-α (1:1000, 17590-1-AP). Tras tres lavados con PBS, las membranas fueron incubadas con anticuerpos secundarios (1:5000) durante 90 minutos a temperatura ambiente. Las bandas proteicas se visualizaron utilizando el kit de sustrato quimioluminiscencia ECL (véase la Tabla de Materiales para más detalles) siguiendo el protocolo del fabricante y se analizaron usando el software ImageJ (versión 1.8; Institutos Nacionales de Salud). El GAPDH (10494-1-AP) se utilizó como control de carga.
Análisis estadístico
Los análisis estadísticos se realizaron utilizando el software R Studio. Los paquetes TwoSampleMR e ieugwasr R se emplearon para extraer datos de exposición y filtrar IVs válidos, incluyendo agrupamiento (r² < 0,001, tamaño de ventana = 10.000 kb), armonización de alelos y cálculo de estadística F (F = β²_exposure/SE²_exposure). Para análisis de RM e IRM inversa, utilizamos los paquetes VariantAnnotation, gwasglue y TwoSampleMR, que facilitaron pruebas robustas de inferencia causal y sensibilidad. Se realizó un estudio de resonancia magnética de dos muestras para investigar rasgos relacionados con metabolitos plasmáticos. Se emplearon cinco enfoques analíticos distintos: Ponderación por Inversa de Varianza (IVW), MR-Egger, Mediana Ponderada (WM), Modo Simple y Modo Ponderado. Los resultados primarios se determinaron mediante el método de la FIV, con significación estadística fijada en P < 0,01. El método de la VV asigna pesos para minimizar el impacto de la variación, proporcionando estimaciones de efecto más precisas. La regresión de MR-Egger se utilizó para evaluar la pleiotropía (el fenómeno en el que los genes afectan a múltiples rasgos), evaluando la pleiotropía horizontal mediante el valor P del intercepto (P > 0,05), y la heterogeneidad evaluada usando el método de IVW combinado con la estadística Q de MR-Egger. No se observó heterogeneidad significativa cuando P > 0,05. Se aplicó corrección por pruebas múltiples utilizando el método de la tasa de descubrimiento de falsos errores (FDR) de Benjamini-Hochberg, considerando que FDR < 0,05 era significativo. La heterogeneidad se cuantificó utilizando la estadística Q de Cochran y la estadística I2 , donde la I2 > 50% indicó una heterogeneidad sustancial. Para los valores de I2 , se calcularon intervalos de confianza del 95% siguiendo el método descrito por Ioannidis et al. Los análisis de sensibilidad se realizaron eliminando secuencialmente cada SNP para evaluar la influencia de variantes individuales en los resultados generales.
Para evaluar la fiabilidad de las mediciones, se realizaron análisis de prueba y reprueba sobre medidas clave de resultado (cuantificación de la burbuja blanca y puntuación histológica) utilizando los coeficientes de correlación intraclase (ICC). El ICC se empleó para evaluar la fiabilidad interevaluador e intra-evaluador de la cuantificación y la puntuación histológica de Western Blot (WB). Los criterios se interpretaron de la siguiente manera: por debajo de 0,50 (pobre), entre 0,50 y 0,75 (moderado), entre 0,75 y 0,90 (bueno) y por encima de 0,90 (excelente)32. Tres investigadores independientes evaluaron ciegamente el mismo conjunto de muestras en dos momentos separados por una semana. Los ICC para la cuantificación Western blot de IL-6, TNF-α, VEGFA y TGF-β oscilaron entre 0,87 y 0,93, lo que indica una excelente fiabilidad. Las ICC para la puntuación histológica fueron 0,85 para la organización del colágeno y 0,82 para la evaluación de neovascularización.