$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Cette étude a analysé des statistiques désidentifiées, au niveau résumé de l’étude d’association génomique large (GWAS), qui sont publiquement accessibles. Conformément aux politiques du dépôt et aux approbations obtenues par les chercheurs initiaux, aucune nouvelle approbation du comité d’examen institutionnel ni consentement individuel éclairé supplémentaire n’a été requis pour cette analyse secondaire. Tous les GWAS contributeurs ont rapporté les procédures d’approbation éthique et de consentement dans leurs publications sources. Toutes les analyses ont été réalisées conformément aux directives institutionnelles et à la Déclaration d’Helsinki.
Aperçu et justification
L’étude a mis en œuvre un cadre bidirectionnel de randomisation mendélienne (RM) à deux échantillons, limité aux statistiques sommaires d’ascendance européenne, afin d’évaluer les relations causales potentielles entre la sclérose en plaques (SEP) et les malignités hématologiques (HM). La conception respecte les trois hypothèses fondamentales de la MR : la pertinence de l’instrument, l’indépendance vis-à-vis des facteurs de confusion et la restriction d’exclusion. Le flux de travail comprend donc (i) l’accès et la conservation des ensembles de données, (ii) la sélection des instruments à l’échelle du génome avec agrégementation par déséquilibre de liaison (LD), (iii) le criblage des confondeurs à l’aide de PhenoScanner, (iv) l’harmonisation des allèles avec traitement explicite des variantes palindromiques, (v) l’évaluation de la directionnalité via le test de Stieger¹², (vi) l’estimation primaire de l’IRM avec des méthodes complémentaires, (vii) un ensemble complet de diagnostics de sensibilité, et (viii) la génération standardisée de figures et de tables sous Contrôle à tests multiples. Chacune de ces étapes est décrite en détail dans les sous-sections suivantes du protocole, et un aperçu du pipeline est présenté à la Figure 1.
Matériaux, logiciels et RRIDs
Les analyses ont été réalisées dans la version R 4.3.1 (RRID :SCR_001905) en utilisant RStudio/Posit 2023.12+ (RRID :SCR_000432). Le regroupement LD, lorsqu’il était effectué localement, utilisait PLINK v1.9 (build 2.3 ; RRID : SCR_001757)13. L’estimation et l’extraction des données par RM utilisaient le paquet R TwoSampleMR v0.5.7 10 ; les recherches d’instruments pour les facteurs de confusion potentiels utilisaient le phénoscanner v1.0 ; la détection et la correction des valeurs aberrantes utilisaient MRPRESSO v1.0. Les versions exactes sont rapportées pour les paquets sans RRID.
Sources de données et accès
Les statistiques résumées sur la SEP ont été obtenues à partir de la méta-analyse du Consortium international de génétique de la sclérose en plaques, comprenant 47 429 cas de SEP et 68 374 témoins, avec un contrôle qualité harmonisé sur 15 cohortes. Les statistiques résumées HM ont été obtenues à partir de FinnGen (n total = 218 792 ; >16 millions de variantes) et comprenaient le lymphome de Hodgkin (HL), le lymphome diffus à grandes cellules B (DLBCL), le lymphome folliculaire (FL), les lymphomes matures à cellules T/NK (MTNKL), d’autres lymphomes non Hodgkin (NHL) ou non spécifiés, de leucémie lymphoïde, de leucémie myéloïde, de leucémie de type cellulaire non spécifié et de néoplasies de myélome multiple(plasma cellulaire) 14. Les ensembles de données étaient consultés via le portail OpenGWAS de l’IEU en utilisant les identifiants d’accèsdocumentés 15. Toutes les analyses de cette étude étaient donc exclusivement basées sur ces ensembles de données GWAS au niveau résumé accessibles au public ; Aucune cohorte institutionnelle interne ou données individuelles de patients n’a été utilisée ou générée. Comme nous n’avons pas identifié de GWAS supplémentaires avec des définitions harmonisées de la SEP et des sous-types de malignité hématologique permettant une réplication complète du pipeline, la validation externe indépendante via un jeu de données séparé n’a pas été réalisée et est reconnue comme une limitation. Le protocole est conçu pour pouvoir être directement réappliqué aux futurs ensembles de données GWAS pour une validation indépendante.
Sélection des instruments et regroupement des LD
Pour chaque exposition, des polymorphismes mononucléotides (SNP) ont été sélectionnés selon la signification génomique (P < 5 × 10-8) en utilisant la fonction extract_instruments dans TwoSampleMR appliquée aux ensembles de données OpenGWAS. Pour garantir l’indépendance de l’instrument, le regroupement LD a ensuite été effectué contre un panneau de référence d’ascendance européenne en utilisant soit les utilitaires internes de regroupement de TwoSampleMR, soit localement avec PLINK, avec un seuil r² de 0,001 et une fenêtre physique de 10 000 kilobases. Lorsque PLINK était utilisé, les paramètres de la ligne de commande étaient fixés à un seuil de signification primaire de 5 × 10-8, r² = 0,001, et une fenêtre de 10 Mo afin que les instruments regroupés correspondent exactement à ces critères. La force de l’instrument a été évaluée à l’aide de la statistique F dérivée de l’estimation de l’effet de l’exposition et de son erreur standard (F ≈ β²/SE²) ; les variantes avec F < 10 ont été exclues des ensembles finaux d’instruments, et les SNP restants ont été transférés au dépistage du PhenoScanner.
Dépistage de confusion avec PhenoScanner
Pour minimiser la pléiotropie horizontale à travers les facteurs de risque connus, chaque instrument candidat a été interrogé dans le PhenoScanner V2 à travers le catalogue GWAS en utilisant le paquet R du phénoscanner (v1.0)16,17. Pour chaque SNP, nous avons demandé toutes les associations rapportées à P < 1 × 10⁻5 et avons inspecté manuellement les traits retournés. Les associations indiquant des liens avec des facteurs de risque hématologiques de malignité établis — tels que l’exposition liée au tabagisme ou des caractéristiques adiposées/anthropométriques (par exemple, indice de masse corporelle, circonférence de taille et mesures de graisse corporelle) — ou des associations directes avec des phénotypes hématologiques de la malignité ont conduit à l’exclusion du SNP correspondant de l’ensembled’instruments 18. Les catégories de caractères considérées comme excluantes étaient basées sur des preuves antérieures liant l’obésité et le tabagisme au risque de leucémie, de lymphome oude myélome 18,19,20. Les requêtes utilisaient des racines larges de mots-clés (par exemple, fumée, cigarette, IMC, obésité, taille, adipose, malignité hématologique, lymphome, leucémie, myélome). Toutes les suppressions étaient documentées dans un tableau de suivi accompagné du trait PhenoScanner qui déclenchait l’exclusion, et les listes d’instruments nettoyées étaient ensuite transmises à l’étape d’harmonisation.
Harmonisation et manipulation palindromique
Les allèles d’effet pour chaque SNP ont été harmonisés entre les ensembles de données d’exposition et de résultats à l’aide de la fonction harmonise_data dans le package TwoSampleMR (v0.5.7, R). Nous avons aligné tous les allèles de résultat à l’allèle d’effet d’exposition afin que des coefficients bêta positifs correspondent toujours au même allèle dans les deux ensembles de données. Les variantes palindromiques (A/T ou C/G) avec des fréquences intermédiaires d’allèle d’effet (0,42-0,58) dans le panneau de référence OpenGWAS ont été traitées comme ambiguës sur le brin et automatiquement supprimées en réglant l’action d’harmonisation pour supprimer des SNP ambiguës. Les SNP palindromiques avec des fréquences effet-allèle hors de cette plage ont été conservés et alignés en utilisant les fréquences allèles rapportées. Comme la disponibilité des allèles et le statut palindromique différaient légèrement selon les résultats de FinnGen, l’harmonisation a été effectuée séparément pour chaque phénotype HM, et le nombre final d’instruments entrant chaque analyse spécifique au résultat a été extrait des objets R harmonisés et rapporté dans les tableaux.
Évaluation de la directionnalité (filtrage Steiger)
La directionnalité a été évaluée en utilisant l’approche Steiger telle qu’implémentée dans la fonction steiger_filtering de TwoSampleMR. Pour chaque SNP, la fonction a d’abord calculé la variance expliquée (R²) dans l’exposition et le résultat à partir du coefficient bêta du GSA, de l’erreur standard et de la taille de l’échantillon. L’étude a ensuite retiré les instruments pour lesquels R² était plus élevé dans le résultat que dans l’exposition, indiquant une possible direction inverse de l’effet. Le filtrage Steiger a été appliqué séparément pour chaque jeu de données de résultats, et les instruments restants (lignes avec steiger_dir == TRUE) ont été sauvegardés et utilisés lors des analyses IRM ultérieures. Les décomptes post-Steiger des instruments ont été enregistrés pour chaque résultat et sont rapportés en même temps que les estimations de la RM.
Estimation primaire de la RM et contrôle par tests multiples
Les estimations causales primaires ont été obtenues avec une RM pondérée par l’inverse de la variance (IVW) dans un modèle à effets fixes utilisant la fonction RM dans TwoSampleMR, avec des méthodes spécifiées comme « mr_ivw », « mr_egger_regression » et « mr_weighted_median ». Pour chaque résultat de la MH, des instruments harmonisés et filtrés par Steiger ont été transmis à la RM, et des logaritariths odds ratios et des erreurs standard ont été extraits et exponentiés pour obtenir des odds ratios (OR) avec des intervalles de confiance (IC) à 95 % pour les traitsbinaires 21. Pour examiner la robustesse face à des violations modestes de l’hypothèse de non-pléiotropie, nous avons également appliqué les estimateurs de régression pondérée de la médiane pondérée etde MR-Egger 22,23, implémentés dans le même package. Lorsque le test Q de Cochran (de mr_heterogeneity) a révélé une hétérogénéité substantielle (P < 0,05), l’étude a également adapté des modèles de viveu à effets multiplicatifs aléatoires et a rapporté des résultats à effets fixes et aléatoires. L’erreur familiale sur les neuf résultats de la HM a été contrôlée par correction de Bonferroni avec α = 0,05/9 = 5,56 × 10-3 ; les associations avec des valeurs de P inférieures à ce seuil étaient considérées comme statistiquement significatives, tandis que celles avec 0,0056 ≤ P < 0,05 étaient interprétées comme suggestives et décrites avec prudence.
Diagnostic de la sensibilité : hétérogénéité, pléiotropie et valeurs aberrantes
La statistique Q de Cochran a été utilisée pour évaluer l’hétérogénéité entre instruments pour les modèles IVW et MR-Egger, implémentée via la fonction mr_heterogeneity dans TwoSampleMR. La pléiotropie horizontale directionnelle a été évaluée à l’aide du test d’interception MR-Egger (mr_pleiotropy_test) et du test global dans le paquetMR-PRESSO 24. MR-PRESSO24 a été exécuté avec les réglages recommandés en R (NbDistribution ≥ 5 000, SignifThreshold = 0,05) pour détecter les valeurs aberrantes influentes et quantifier la distorsion potentielle en comparant les estimations de la FIW avant et après la suppressiondes valeurs aberrantes 25. Des analyses sans filtre (mr_leaveoneout) ont été réalisées pour chaque paire exposition-résultat afin de déterminer si un SNP unique influençait de manière disproportionnée l’estimation globale. Pour des raisons de transparence et de reproductibilité, toutes les sorties de diagnostic ont été exportées depuis R et rapportées avec les décomptes d’instruments correspondants après harmonisation, filtrage Steiger et suppression des valeurs aberrantes MR-PRESSO.
Puissance de l’instrument et évaluation NOME
La force de l’instrument pour MR-Egger a été quantifiée à l’aide de la statistique I2GG, calculée comme 1 moins la moyenne des erreurs au carré des associations SNP-exposition divisée par leur variance entre instruments26. Des valeurs proches de 1 indiquent une meilleure conformité à l’hypothèse sans erreur de mesure (NOME) ; des valeurs plus basses suggèrent une possible dilution de régression et une interprétation prudente des résultats de la MR-Egger. I2GX a été calculé et rapporté pour chaque analyse spécifique au résultat final.
Randomisation mendélienne inverse
L’ensemble du pipeline a été répété dans le sens inverse en traitant chaque sous-type de HM comme l’exposition et le MS comme le résultat. Lorsque les instruments significatifs à l’échelle du génome étaient insuffisants pour une exposition donnée à la HM, un seuil de sélection assoupli de P < 5 × 10-6 était permis tout en maintenant les mêmes paramètres de regroupement de la lumière du diamètre, le dépistage par phénoscanner, les procédures d’harmonisation, le filtrage Steiger et les diagnostics de sensibilité. Les analyses utilisant des seuils assouplissants étaient clairement indiquées dans les tableaux et légendes de figures correspondants.
Visualisation et exportation de figures
Des graphiques dispersés, forestiers, entonnoirs et laisse-un-out ont été générés avec des légendes placées sous les panneaux et des tailles de police ajustées pour s’assurer que les labels ne masquent pas les données tracées. Les limites d’axe ont été standardisées selon des résultats comparables afin de faciliter la comparaison visuelle. Les figurines étaient exportées à un minimum de 300 dpi en formats sans perte tels que TIFF ou PNG. Toutes les valeurs numériques tracées ont été croisées avec les estimations rapportées afin d’assurer la cohérence entre le texte, les tableaux et les chiffres.
Reproductibilité et partage des données
Les graines aléatoires ont été corrigées lorsque cela était applicable, les versions logicielles ont été enregistrées, et les scripts d’analyse ainsi que les objets intermédiaires ont été archivés pour permettre la réexécution de toutes les étapes. Les identifiants d’acquisition de données et les définitions de phénotypes ont été documentés, et les listes d’instruments à chaque étape de filtrage – post-regroupement, post-harmonisation, post-filtrage Steiger et post-MR-PRESSO ont été préparées pour le téléchargement sous forme de fichiers tableur conformément aux directives de la revue.