Conception de l'étude
Cette recherche a utilisé une approche de MR à deux échantillons afin d'explorer la relation causale entre les traits inflammatoires associés aux NETs et le DT2 ainsi que ses complications. Le cadre de la MR repose sur trois hypothèses fondamentales20 : 1) la variante génétique doit être fortement associée à l'exposition, 2) la variante génétique ne doit influencer le résultat qu'à travers l'exposition (exclusion de la pléiotropie), et 3) la variante génétique ne doit pas être associée à des facteurs de confusion. Pour un aperçu détaillé du processus de cette étude, veuillez vous référer à la Figure 1. Le manuscrit suit les recommandations MR-STROBE pour la présentation des études de MR, garantissant ainsi des normes rigoureuses de transparence et de reproductibilité21.
Sources de données
Les données d’études d’association génome entier (GWAS) utilisées dans cette analyse MR provenaient de bases de données publiques GWAS. En ce qui concerne les GWAS au niveau résumé portant sur le diabète de type 2 (T2DM), elles comprenaient la plus vaste méta-analyse réalisée sur une population européenne, incluant 74 124 cas et 824 006 témoins, ainsi qu’une étude prospective en cas-cohorte imbriquée menée en Europe, comportant 9 978 cas et 12 348 témoins22. Pour les complications diabétiques, notamment le T2DM avec complications rénales et le T2DM avec complications circulatoires périphériques, cette étude a extrait les statistiques sommaires GWAS d’études dans lesquelles les cas étaient des patients atteints de T2DM présentant une complication spécifique, tandis que les témoins étaient des individus sans T2DM. Des informations détaillées sur les sources de données sont disponibles dans le Tableau supplémentaire 1.
Les données GWAS relatives aux caractères inflammatoires associés aux NET ont été obtenues à partir du catalogue GWAS (voir Tableau supplémentaire 2). Afin de tenir compte de l'hétérogénéité en termes de spécificité biologique, cette étude a classé les expositions incluses en deux groupes distincts selon leurs rôles fonctionnels en biologie des neutrophiles : (1) Facteurs liés à la NETose principaux : cette catégorie regroupe les marqueurs directement impliqués dans la formation structurale des NET ou dans le processus enzymatique de décondensation de la chromatine. Plus précisément, cette étude a inclus les NET elles-mêmes, la myéloperoxydase (MPO), l'élastase des neutrophiles (NE) et le complexe MPO-ADN. La MPO et la NE sont des enzymes essentielles à la dégradation des histones et à la décondensation de la chromatine, tandis que le complexe MPO-ADN constitue un marqueur spécifique substitut des NET23. (2) Médiateurs inflammatoires associés aux NET : cette catégorie comprend des cytokines et des médiateurs agissant comme régulateurs en amont ou effecteurs en aval étroitement liés à la NETose, mais également impliqués dans des voies inflammatoires plus larges. Ce groupe inclut l'interleukine-6 (IL-6)24, le facteur de nécrose tumorale-alpha (TNF-α)25, les taux de lipocaline associée à la gélatinase des neutrophiles (NGAL)26 et le facteur 1 du réseau de communication cellulaire (CCN1)27.
Les données utilisées dans cette étude provenaient de bases de données en libre accès ou d'études GWAS publiées dans des travaux antérieurs ; par conséquent, une approbation éthique n'était pas requise pour cette recherche.
Sélection d'instruments génétiques liés aux NET
Étant donné que la formation des NETs représente un processus biologique dynamique qui n'est pas directement mesuré dans les études d'association génomique classiques (GWAS), cette étude a utilisé une approche génique afin d'identifier des instruments génétiques pour les traits inflammatoires associés aux NETs. Un ensemble complet de 257 gènes, connus pour être impliqués de façon critique dans la formation et la régulation des NETs, tels que MPO, ainsi que des gènes codant pour des histones et des protéines des granules des neutrophiles, a été recensé à partir d'études mécanistiques établies et de la littérature publiée.28.
Sélection de la variable instrumentale
Pour identifier des instruments génétiques robustes associés aux caractères inflammatoires liés aux NET, cette étude a mis en œuvre un flux de travail de filtrage séquentiel utilisant des statistiques sommaires d'études d'association pangénomique (GWAS). Les étapes opérationnelles spécifiques étaient les suivantes :
Dépistage initial : les variants génétiques significativement associés aux caractères d'exposition ont été extraits29 selon un seuil de significativité de P < 5 × 10⁻6.
Filtrage par fréquence de l'allèle mineur (MAF) : afin d'assurer une puissance statistique suffisante, les polymorphismes nucléotidiques simples (SNPs) ayant une MAF ≤ 0,01 ont été exclus30.
Élagage par déséquilibre de liaison (LD) : afin de réduire les biais dus au déséquilibre de liaison (LD), des SNP indépendants ont été sélectionnés à l’aide de la fonction d’élagage, avec des paramètres fixés à r2 < 0,001 dans une fenêtre de 10 000 kb31.
Évaluation de la force de l'instrument : La force de chaque IV restant a été quantifiée à l'aide de la statistique F, calculée selon la formule F = R2 × (N-2) / (1-R2). Seuls les SNP présentant une statistique F > 10 ont été conservés afin de minimiser le biais dû à un instrument faible32.
Identification des SNP proxy et harmonisation des données
Afin de remédier aux SNP manquants dans l'ensemble de données de l'issue, cette étude a effectué une recherche de substituts et un processus d'harmonisation des données :
Substitution par proxy : Lorsqu'un SNP cible de l'étude d'association à l'échelle du génome (GWAS) d'exposition n'était pas disponible dans la GWAS de l'issue, cette étude a utilisé la fonction LDproxy() basée sur le panel de référence européen du projet 1 000 génomes. Un SNP proxy candidat était sélectionné uniquement s'il présentait une forte liaison d'association avec le SNP d'origine (r2> 0,8). En l'absence de proxy qualifié, le SNP était exclu.
Harmonisation : les jeux de données relatifs à l'exposition et au résultat ont été alignés à l'aide de la fonction harmonise_data() du package TwoSampleMR (R version 4.0.5). Cette étude a défini le paramètre action = 2 afin d'aligner automatiquement tous les SNP sur le brin direct et de supprimer les SNP palindromiques ayant une orientation de brin ambiguë.
Vérification : après harmonisation, cette étude a examiné manuellement le jeu de données harmonisé afin de confirmer que les fréquences de l'allèle d'effet étaient cohérentes entre les données d'exposition et celles de l'issue.
Analyse MR et tests de sensibilité
L'inférence causale a été réalisée à l'aide du package TwoSampleMR (R version 4.0.5).
Analyses primaires et secondaires : Cette étude a appliqué la méthode pondérée par la variance inverse (IVW) comme approche principale33. Des analyses complémentaires ont été réalisées à l’aide des méthodes MR-Egger34, médiane pondérée et mode pondéré afin d’assurer la robustesse des résultats35.
Évaluations de la sensibilité : l'hétérogénéité entre les IV a été évaluée à l'aide du test Q de Cochran36 via la fonction mr_heterogeneity(). La pléiotropie horizontale a été évaluée à l'aide du test d'interception MR-Egger34 (mr_pleiotropy_test()).
Détection des valeurs aberrantes : cette étude a utilisé le package MR-PRESSO pour détecter d'éventuelles valeurs aberrantes37. La fonction mr_presso() a été exécutée avec 1 000 simulations. Les SNP aberrants identifiés (P < 0,05) ont été supprimés, et les estimations causales ont été recalculées afin de vérifier la stabilité des résultats. De plus, une analyse de type « laisser-un-de-hors » (leave-one-out) a été réalisée pour s'assurer que l'association causale n'était pas due à un seul SNP38.
Correction statistique
Afin de tenir compte des tests multiples, les valeurs P obtenues à partir des analyses de RM ont été corrigées à l'aide de la méthode du taux de découverte de faux positifs (FDR). Cette correction a été mise en œuvre à l'aide de la fonction P.adjust() dans R, avec le paramètre method = "fdr". Les associations présentant une valeur P corrigée (PFDR) < 0,05 ont été considérées comme statistiquement significatives.