Association entre la longueur des télomères et la thyrotoxicose : perspectives issues d’une étude de randomisation mendélienne à deux échantillons

Résumé

Ce protocole décrit un pipeline de randomisation mendélienne à deux échantillons évaluant si la longueur des télomères affecte causalement le risque de thyrotoxicose à partir de données génétiques publiques résumées. Il couvre la sélection des instruments, l’harmonisation, l’estimation primaire, les analyses de sensibilité et le code R reproductible avec des sorties prêtes à la figure pour soutenir un reporting transparent.

Résumé

La thyrotoxicose est un trouble endocrinien caractérisé par un excès d’hormones thyroïdiennes, mais ses liens étiologiques avec la biologie du vieillissement systémique restent incomplètement définis. La longueur des télomères (TL) reflète la sénescence cellulaire et la stabilité du génome, et a été impliquée dans de multiples maladies complexes. Nous avons mené une étude de randomisation mendélienne (RM) sur deux échantillons afin d’évaluer l’effet causal de la TL génétiquement prédite sur le risque de thyrotoxicose. Les instruments génétiques pour la langue différente ont été dérivés d’une vaste étude d’association génomique (GWAS) d’ascendance européenne (n > 470 000). Les statistiques résumées de la thyrotoxicose ont été obtenues à partir de la dernière publication de FinnGen (≈4 000 cas et >210 000 témoins de témoins). Les analyses primaires pondérées par l’inverse de la variance ont indiqué que la TL génétiquement proxy plus longue est associée à un risque plus faible de thyrotoxicose, et la direction ainsi que l’ampleur de l’effet étaient cohérentes entre les estimateurs complémentaires (MR-Egger, médiane pondérée/maximum de vraisemblance, MR-PRESO et MR-RAPS). Les analyses de sensibilité n’ont montré aucune preuve de pléiotropie directionnelle, et le Q de Cochran a été utilisé pour évaluer l’hétérogénéité. Un test de directionnalité de Steiger a soutenu le flux causal de la langue de la langue vers la thyrotoxicose.

À notre connaissance, ce travail fait partie des premières analyses IRM à évaluer la relation causale entre le risque global de thyrotoxicose et la TL utilisant des ressources GWAS contemporaines, en étendant les preuves antérieures portant sur les phénotypes liés à l’hyperthyroïdie. Ces résultats suggèrent que les processus de vieillissement cellulaire indexés par la TL pourraient contribuer à la susceptibilité à la thyrotoxicose et motiver de futures études longitudinales et mécanistes sur la biologie des télomères dans les dysfonctionnements thyroïdiens.

Introduction

La thyrotoxicose est une affection endocrinienne courante caractérisée par un excès d’hormones thyroïdiennes circulantes. Les manifestations typiques incluent la perte de poids, la tachycardie, l’irritabilité, l’anxiété et les tremblements, les cas sévères prédisposant à des arythmies, une insuffisance cardiaque et une perte^{osseuse 1}. Bien que les caractéristiques cliniques et les stratégies thérapeutiques soient relativement bien établies, les mécanismes biologiques qui façonnent la susceptibilité et la progression interindividuelles restent incomplètement définis. Les étiologies vont de l’hyperfonction auto-immune (par exemple, la maladie de Basedowdoid) aux nodules toxiques, à la thyroïdite, à l’exposition aux hormones exogènes, aux réactions médicamenteuses et aux causes^{tumorales 2}. À travers ces entités, l’activation immunitaire, le stress oxydatif et la biogenèse mitochondriale altérée réapparaissent comme des thèmes mécanistes plausibles reliant l’excès hormonal à des lésions tissulaires et des complications systémiques.

La longueur des télomères (TL) — la structure ADN-protéine en tandem qui couvre les extrémités des chromosomes — intègre le vieillissement cellulaire, l’histoire réplicative et la stabilité du génome. La TL raccourcit avec la division cellulaire et est accélérée par le stress oxydatif et l’inflammation chronique, des processus également présents dans la dysfonction thyroïdienne ^3,4. Une TL plus courte a été associée à des troubles cardiométaboliques, certains cancers et des affections^{neurodégénératives 5}, et les données émergentes suggèrent des liens avec des phénotypes endocriniens tels que l’obésité et le diabète de type 2, dans lesquels les axes hormonaux et immunométaboliques s’entrecroisent avec le stress^{réplicatif 6}. Étant donné que les hormones thyroïdiennes régulent le taux métabolique basal, la fonction mitochondriale et la signalisation proliférative, la dynamique de la TL peut être particulièrement pertinente dans les états thyrotoxiques, où l’hypermétabolisme systémique et la dysrégulation immunitaire pourraient à la fois influencer et être influencés par la biologie des^{télomères 7,8}.

Cependant, la question de savoir si la TL joue un rôle causal dans la thyrotoxicose reste incertaine. Les études observationnelles conventionnelles sont vulnérables à la confusion (par exemple, mode de vie, comorbidités) et à la causalité inverse (maladie affectant la langue des langues). Pour remédier à ces limitations, nous utilisons la randomisation mendélienne (RM) à deux échantillons, qui exploite robustement les variants génétiques de la lignée germinale associés à la TL comme instruments non confondus pour tester son effet sur le risque^{de maladie 9}. Les études d’association à l’échelle du génome fournissent désormais des instruments de traduction suffisamment solides et de grands ensembles de données de résultats bien phénotypés pour permettre une évaluation causale^{rigoureuse 10}. En intégrant l’épidémiologie génétique à la biologie des télomères, cette étude évalue si la TL génétiquement proxy influence le risque de thyrotoxicose, clarifiant ainsi les voies mécanistes et explorant la TL comme biomarqueur candidat pour la stratification du risque et le développement d’interventions futures.

À notre connaissance, cette analyse fait partie des premières études IRM à examiner le risque global de thyrotoxicose en relation avec la TL en utilisant des ressources GWAS contemporaines à grande échelle, poursuivant les travaux antérieurs portant principalement sur les phénotypes liés à l’hyperthyroïdie ou les maladies auto-immunes^{de la thyroïde 11}.

Protocole

Cette étude applique un plan de randomisation mendélienne (RM) en deux échantillons pour évaluer si la longueur des télomères (TL) génétiquement proxy influence causalement le risque de thyrotoxicose. Seules des statistiques résumées GWAS désidentifiées et accessibles au public ont été utilisées, et aucune donnée individuelle n’a été consultée. Le Conseil d’examen institutionnel de l’hôpital municipal de Qingdao a déterminé que l’analyse est exemptée de toute révision supplémentaire car elle repose exclusivement sur des données publiques résumées. Toutes les études d’association génomique ayant contribué ont obtenu un consentement éclairé et des approbations éthiques dans le cadre de leurs protocoles d’origine. L’analyse est structurée pour minimiser la confusion et la causalité inverse en utilisant des variantes germinales associées à la langue de langue comme instruments non confondus et en mettant en œuvre des procédures de sensibilité pour interroger la pléiotropie, l’hétérogénéité et la direction de l’effet.

Sources de données
L’ensemble d’expositions pour TL a été obtenu auprès de l’IEU OpenGWAS sous l’identifiant ieu-b-4879, regroupant environ 472 174 participants d’ascendance européenne. L’ensemble de données de résultats pour la thyrotoxicose a été extrait du consortium FinnGen, en utilisant le critère final de 2021 finn-b-thyROTOXICOSIS avec 4 142 cas et 213 693 témoins. Ces ressources fournissent des estimations efficaces et des erreurs standard nécessaires pour construire des instruments au niveau SNP, harmoniser les allèles et estimer les effets causaux avec des estimateurs MR établis.

Environnement logiciel et informatique
Toutes les analyses ont été réalisées dans R (version 4.3.1). Le principal paquet analytique était TwoSampleMR (version 0.5.7), complété par ieugwasr pour l’accès programmatique aux ressources GSA, MRPRESSO pour la détection des valeurs aberrantes et les tests de distorsion, mr.raps pour l’estimation robuste sous instruments faibles et pléiotropie idiosyncratique, RadialMR pour la visualisation radiale, et des paquets polyvalents incluant tidyverse et data.table. Les informations de session et les versions des packages sont écrites dans un fichier pour garantir une reproductibilité stricte.

Hypothèses fondamentales de la randomisation mendélienne
Le cadre de la RM suppose que les instruments génétiques sont fortement associés à la TL, indépendants des facteurs qui embrouillent la relation entre la TL et la thyrotoxicose, et n’influencent la thyrotoxicose que via la TL plutôt que par des voies alternatives. Le plan d’analyse opérationnelle ces hypothèses en quantifiant la force par variante, en testant la pléiotropie horizontale déséquilibrée à l’aide de méthodes basées sur l’interception et le criblage des valeurs aberrantes, en évaluant l’hétérogénéité entre variantes, et en confirmant la direction d’effet TL-thyrotoxicose avec un test de directionnalité^{formel 12}.

Sélection des instruments et contrôle qualité
Les instruments génétiques ont été sélectionnés à partir du GWAS TL avec un seuil de signification génomique de p < 5×^10-8. Pour garantir l’indépendance, le déséquilibre de liaison a été traité en agrégeant à r² = 0,001 dans une fenêtre de 10 000 kb à l’aide de données de référence européennes du projet 1000 génomes ; lorsque les variantes étaient corrélées, la variante avec la valeur p plus petite pour l’association à la TL était conservée. La force de l’instrument a été résumée pour chaque variante à l’aide de la statistique F_i = \beta_{E,i}^{2}/\mathrm{SE}(\beta_{E,i})^{2} et jugée selon le critère conventionnel F > 10^13,14. Lorsque cela était approprié, l’approximation de Sanderson-Windmeijer a été référencée pour décrire la force globale des instruments à travers plusieurs variantes dans un échantillon de taille N, avec des instruments K expliquant la variance d’exposition R^{2} : F_{\mathrm{SW}} = \{[R^{2}/(1-R^{2})]\,(N-K-1)\}/K ¹⁵. Les variantes palindromiques ambiguës avec des fréquences allèles proches de 0,5 étaient exclues à moins que l’information sur la fréquence des allèles ne permette un alignement non ambigu.

Harmonisation
Pour chaque instrument, l’association correspondante avec la thyrotoxicose a été extraite de FinnGen et alignée de manière à ce que les tailles d’effet représentent le même allèle d’effet à travers l’exposition et le résultat. L’harmonisation a éliminé les inadéquations des allèles, corrigé les problèmes de brins et exclu les SNP palindromes avec une ambiguïté non résolue, produisant un ensemble de données adapté à la construction valide du rapport de Wald.

Estimation causale primaire
L’analyse primaire a utilisé l’estimateur pondéré en inverse de la variance (IVW) pour méta-analyser les rapports de Wald spécifiques aux SNP en un effet causal^{regroupé 16}. Les estimations sont rapportées en termes de rapports de cotes par une augmentation d’écart-type en TL, avec des intervalles de confiance à 95 % dérivés selon des modèles à effets fixes et aléatoires. Le choix du modèle a été guidé par des diagnostics d’hétérogénéité, et les deux spécifications sont présentées pour faciliter une interprétation robuste.

Analyses de sensibilité et vérifications de robustesse
La robustesse a été évaluée à l’aide de la régression MR-Egger avec un terme d’interception pour tester la pléiotropie directionnelle déséquilibrée, l’estimation du maximum de vraisemblance pour améliorer l’efficacité sous homogénéité tout en atténuant l’erreur de mesure, MR-PRESSO pour réaliser un test global des valeurs aberrantes et pour estimer les effets corrigés par les valeurs aberrantes, et MR-RAPS pour fournir des estimations résistantes aux instruments faibles et à la pléiotropie^{idiosyncratique 17}. L’hétérogénéité entre variantes a été quantifiée à l’aide de la statistique Q de Cochran dans le cadre de la^{FIV 18}. La direction de l’effet a été examinée avec le test de Steiger, qui compare la proportion de variance expliquée dans l’exposition et le résultat afin de déterminer si les données sont plus compatibles avec la TL causant la thyrotoxicose plutôt que l’inverse. L’interception MR-Egger a été utilisée comme test formel de la pléiotropie directionnelle. Les analyses Leave-One-Out ont été inspectées afin de s’assurer que l’association globale n’était pas influencée par une seule variante.

Randomisation mendélienne inverse
Pour sonder la causalité inverse potentielle, le pipeline analytique a été répété avec la thyrotoxicose comme exposition et la TL comme résultat. Les mêmes critères de sélection des instruments, procédures d’harmonisation, estimateur primaire de la VIO et analyses de sensibilité ont été appliqués afin que des conclusions sur la directionnalité soient tirées dans un cadre causal identique.

Tests multiples, alimentation et rapports
Le test d’hypothèse principale concerne l’estimateur de la FIV pour l’effet de la TL sur la thyrotoxicose. Les estimateurs de sensibilité et les tests diagnostiques sont interprétés comme des preuves à l’appui ; Les valeurs p sont rapportées en notation scientifique pour plus de clarté, et les conclusions mettent l’accent sur la cohérence entre les méthodes plutôt que sur des seuils de signification isolés. Les résumés de la force de l’instrument et la proportion de variance expliquée indiquent une puissance approximative selon les formulations standard de non-centralité, en reconnaissant que la puissance dépend de la taille de l’échantillon, de la force de l’instrument et de l’ampleur de l’effet réel.

Reproductibilité computationnelle et commandes exactes
La reproductibilité est assurée en fournissant la séquence complète de commandes R qui recréent la sélection des instruments, l’extraction des résultats, l’harmonisation, les analyses primaires et de sensibilité, la résonance magnétique inversée, les sorties diagnostiques et l’exportation de fichiers prêts à l’analyse. Le script écrit des fichiers CSV stables et lisibles par l’humain correspondant à la liste d’instruments, au jeu de données harmonisé et au résumé des estimations et diagnostics de la RM.

# R 4.3.1 ; TwoSampleMR 0.5.7
# Installation optionnelle :
# install.packages(c(« TwoSampleMR », « ieugwasr », « MRPRESSO », « mr.raps »,
« RadialMR », « tidyverse », « data.table »))

bibliothèque(TwoSampleMR)
Bibliothèque (IEUGWASR)
bibliothèque (MRPRESSO)
bibliothèque(mr.raps)
bibliothèque (RadialMR)
bibliothèque (tidyverse)
bibliothèque(data.table)

# Exposition : longueur des télomères (IEU OpenGWAS)
exposure_id <- « IIU-B-4879 »

# Résultat : FinnGen 2021 critère thyrotoxicose utilisé pour les analyses rapportées
outcome_id <- « finn-b-thyROTOXICOSE »

# Sélection d’instruments avec seuil génomique et agglutination stricte de la LD
exp <- extract_instruments(résultats = exposure_id, p1 = 5e-8, regroupement = VRAI, r2 = 0,001, kb = 10000)
exp$F_stat <- (exp$beta.exposure^2) / (exp$se.exposure^2)
fwrite (exp, « S1_instruments_TL.csv »)

# Extraction des résultats et harmonisation
out <- extract_outcome_data(SNPS = exp$SNP, résultats = outcome_id)
dat <- harmonise_data(exp, sortie, action = 2)
fwrite(dat, « S2_harmonised_TL_vs_thyrotoxicosis.csv »)

# Estimateurs de MR primaire et de sensibilité
res <- mr(dat, method_list = c(« mr_ivw »,"mr_ivw_fe »,"mr_egger_regression »,
« mr_weighted_median », » « mr_raps" », « mr_maxlik »))

het <- mr_heterogeneity(dat) # Q de Cochran
pleio <- mr_pleiotropy_test(dat) # Egger intercept
Steiger <- directionality_test(dat) # Directionnalité Steiger

# Estimations globales et corrigées par les valeurs aberrantes MR-PRESSO
mrpresso <- mr_presso(BetaOutcome = « beta.outresult », BetaExposure = « beta.exposition »,
SdOutcome = « se.outissue », SdExposition = « se.exposition »,
OUTLIERtest = VRAI, DISTORTIONtest = VRAI,
data = dat, NbDistribution = 1000, SignifThreshold = 0,05) #

# IRM inversée : thyrotoxicose (exposition) -> TL (résultat)
rev_exp <- extract_instruments(résultats = outcome_id, p1 = 5e-8, amas = VRAI, r2 = 0,001, kb = 10000)
rev_out <- extract_outcome_data(SNPS = rev_exp$SNP, résultats = exposure_id)
rev_dat <- harmonise_data(rev_exp, rev_out, action = 2)

rev_res <- mr(rev_dat, method_list = c(« mr_ivw »,"mr_ivw_fe »,"mr_egger_regression »,

« mr_weighted_median », » « mr_raps" », « mr_maxlik »))
rev_het <- mr_heterogeneity(rev_dat)
rev_pleio <- mr_pleiotropy_test(rev_dat)
rev_steiger <- directionality_test(rev_dat)

# Exportations pour l’archivage et la génération de figures/tables
write.csv(bind_rows(res), « S3_mr_results_primary.csv », ligne.noms = FAUX)
write.csv(het, « S3_mr_heterogeneity.csv », ligne.noms = FAUX)
write.csv(pleio, « S3_mr_pleiotropy_egger.csv »,row.names = FAUX)
write.csv(steiger, « S3_mr_steiger.csv », ligne.noms = FAUX)
write.csv(bind_rows(rev_res), « S3_reverse_mr_results.csv », row.names = FAUX)
write.csv(rev_het, « S3_reverse_heterogeneity.csv », ligne.noms = FAUX)
write.csv(rev_pleio, « S3_reverse_pleiotropy_egger.csv », ligne.noms = FAUX)
write.csv(rev_steiger, « S3_reverse_steiger.csv », ligne.noms = FAUX)

Résultats

Caractéristiques de l’instrument et vérifications d’hypothèses
À partir du GWAS de longueur télomère, 129 instruments significatifs et indépendants de la LD à l’échelle du génome (p < 5×^10-8 ; r² < 0,001 ; fenêtre de 10 000 kb) ont été conservés. Les F-statistiques par SNP dépassaient le seuil conventionnel pour la grande majorité des variantes, indiquant un risque limité de biais faible pour l’instrument. L’harmonisation a supprimé les variantes palindromiques ambiguës et ...

Discussion

Cette étude de randomisation mendélienne (RM) fournit des preuves que des télomères génétiquement prédits plus longs réduisent le risque de^{thyrotoxicose 19}. Les estimations directionnellement concordantes entre IVW, MR-Egger, vraisemblance maximale, MR-PRESSO et MR-RAPS indiquent une association protectrice robuste, et les contrôles de directionnalité soutiennent une voie TL → thyrotoxicose plutôt que^{les 20} inverses. La plausibilité biologi...

Matériaux

Liste des matériaux utilisés dans cet article

Nom	Entreprise	Numéro de catalogue	Commentaires
Projet 1000 Génomes	Données de référence pour les agglomérations par déséquilibre de liaison (LD).	Projet 1000 Génomes	-
Cochran' s statistique Q	Méthode statistique utilisée pour le test d’hétérogénéité en IRM.	N/A	-
Jeu de données FinnGen	Statistiques résumées GWAS pour la thyrotoxicose (cas et témoins).	FinnGen Consortium	Sortie 2021, finn-b-thyROTOXICOSE
Jeu de données IEU OpenGWAS	Statistiques résumées GWAS pour la longueur des télomères.	IEU OpenGWAS	ID de jeu de données : ieu-b-4879
Paquet IEUGWASR	Package R utilisé pour l’accès programmatique aux données GSA.	CRAN	Version 0.1.5
Paquet Mr.Raps	R package pour l’estimation robuste en randomisation mendélienne sous instruments faibles.	CRAN	Version 0.1.2
Régression MR-Egger	Méthode d’analyse de sensibilité pour tester la pléiotropie directionnelle.	N/A	-
Test des valeurs aberrantes MR-PRESSO	Test statistique utilisé pour identifier et corriger les valeurs aberrantes dans les analyses IRM.	N/A	-
Pack MRPRESSO	R package utilisé pour la détection des valeurs aberrantes et les tests de distorsion dans l’analyse IRM.	CRAN	Version 1.0
Logiciel R	Environnement de calcul statistique et graphique utilisé pour toutes les analyses.	Fondation R pour l’informatique statistique	Version 4.3.1
Boîtier RadialMR	Package R pour la visualisation radiale des résultats de l’IRM.	CRAN	Version 1.1.2
Paquet TwoSampleMR	R package utilisé pour les analyses de randomisation mendélienne.	CRAN	Version 0.5.7