$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Obtivemos aprovação ética e consentimento informado do Comitê de Ética em Pesquisa Biomédica do Primeiro Hospital Afiliado da Universidade de Nanchang. Número de Ética: (2025)CDYFYYLK(08-007).
Análise de RM
Recuperação de dados
Os dados plasmáticos pQTL foram obtidos a partir do estudo de Zheng et al.14, que integrou cinco conjuntos de dadosGWAS 15, 16, 17, 18, 19, e do estudo de Ferkingstad et al. Os critérios de inclusão para os dados foram os seguintes: (i) associações significativas em todo o genoma (p < 5 × 10⁻⁸); e (ii) proteínas plasmáticas como potenciais alvos terapêuticos para a OA. O desenho do estudo é resumido na Figura 1. Primeiramente, identificamos alvos terapêuticos candidatos usando dados GWAS do OpenGWAS do IEU e dados plasmáticos de pQTL dos estudos de Zheng14 e Ferkingstad20 (Tabela Suplementar S1 e Tabela Suplementar S2). Foram então realizados filtros Steiger e varreduras de fenótipo para validar a robustez dos resultados. O OpenGWAS do IEU (https://gwas.mrcieu.ac.uk/) foi usado para obter estatísticas resumidas para OA de quadril ou joelho (n = 417.596), OA do joelho (n = 403.124) e OA do quadril (n = 393.873)21.
Comandos de filtragem SNP
SNPs com significância genômica (p < 5 × 10⁻⁸) foram submetidos a um processo de agrupamento (r² < 0,001, estatísticas F > 10, tamanho da janela = 10.000 kb) antes da análise por RM.
Análise de RM
Para investigar potenciais alvos de fármacos, a análise por RM foi realizada usando proteínas plasmáticas como exposições e OA como resultado, implementada via o pacote "TwoSampleMR" em R (v4.3.1). Quando apenas um único pQTL estava disponível para uma proteína, a razão de Wald era utilizada; caso contrário, foi aplicada a RM ponderada por variância inversa (RM-IVW), seguida por avaliações de heterogeneidade e pleiotropia. A correção de Bonferroni foi usada para considerar múltiplos testes, com um limiar de p < 5,63 × 10⁻⁵ para priorizar proteínas.
Filtragem de Steiger e varredura de fenótipos
Para avaliar a causalidade reversa, realizamos filtragem de Steiger. Um resultado de "TRUE" com p < 0,05 indicou que não houve causalidade reversa. A varredura de fenótipo foi realizada usando LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 para examinar associações de pQTLs com outros traços. Os limiares eram R² = 0,1 e uma janela de ±500.000 pares base. Efeitos pleiotrópicos foram atribuídos a pQTLs que atendem tanto aos seguintes: (i) associação significativa genômica (p < 5 × 10⁻⁸) quanto (ii) associação com fatores de risco conhecidos para OA.
Estudo de associação em ampla perfícia
Para considerar a pleiotropia gênica e os efeitos fora do alvo, realizamos um estudo de associação em todo o fenômico (PheWAS) usando o Portal PheWAS da AstraZeneca (https://azphewas.com/), que contém 15.500 fenótipos binários e 1.500 fenótipos contínuos de ~450.000 participantes do Biobancodo Reino Unido 23. Os limiares foram definidos para valores padrão para minimizar falsos positivos.
Rede de interação proteína - proteína (IBP)
Para visualizar interações entre potenciais alvos proteicos identificados pela RM, utilizamos o GeneMANIA (https://genemania.org/) para análise de interações proteína-proteína e visualizaçãode resultados 24.
Análise de enriquecimento
Para investigar a relevância biológica, realizamos análises de enriquecimento utilizando ferramentas de bioinformática da https://www.bioinformatics.com.cn para análise e visualização de dados.
Fluxo de trabalho transcriptômico
O RNA total foi extraído usando o kit de reagentes de extração de RNA seguindo as diretrizes do fabricante. A qualidade do RNA foi avaliada usando um sistema automatizado de avaliação da qualidade do RNA; apenas amostras com RIN ≥7.0 foram usadas. A qualidade foi confirmada por eletroforese em gel de agarose livre de RNase (gel de 1,5%). O mRNA eucariótico foi enriquecido usando esferas de Oligo(dT); o mRNA procariótico foi enriquecido usando o Kit Magnético de Eliminação de RNA. O mRNA foi fragmentado (200-700 nt) e convertido em cDNA usando o RNA Library Prep Kit. A biblioteca de cDNA era reparada por extremidades, com cauda em A, ligada a adaptadores, purificada usando esferas magnéticas purificadoras de DNA (1.0×) e amplificada por PCR. O sequenciamento era realizado em uma plataforma de sequenciamento de próxima geração de alta produtividade. Genes diferencialmente expressos foram definidos por log₂FC > 1 e ajustados p < 0,05.
Farmacologia de redes
Para identificar potenciais medicamentos para proteínas-alvo, usamos BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. Um corte de pontuação de 0,74 (LR = 32,5) foi usado para selecionar compostos conhecidos e previstos. Componentes herbais foram retirados do TCMSP (https://old.tcmsp-e.com/index.php) e filtrados com OB > 30% e DL > 0,1826.
Acoplamento molecular
O acoplamento molecular foi usado para avaliar interações de ligação. As estruturas proteicas foram retiradas do PDB (https://www.rcsb.org/). O UCSF Chimera foi usado para pré-processar estruturas removendo ligantes e solventes. O AutoDock Tools era usado para calcular cargas de Gasteiger e definir centros e tamanhos de caixas. As estruturas dos medicamentos foram obtidas do PubChem (https://pubchem.ncbi.nlm.nih.gov/) e pré-processadas de forma semelhante. O acoplamento foi realizado usando AutoDock Vina. As dimensões das caixas variavam conforme o alvo. As afinidades de ligação foram calculadas e os resultados visualizados na UCSF Chimera.