Zusammenhang zwischen Telomerlänge und Schilddrüsentoxikose: Erkenntnisse aus einer zweistichprobigen Mendelschen Randomisierungsstudie

Zusammenfassung

Dieses Protokoll beschreibt eine zweistichprobige Mendelsche Randomisierungspipeline, die anhand öffentlicher genetischer Zusammenfassungsdaten bewertet, ob die Telomerlänge kausal das Risiko der Thyrotoxikose beeinflusst. Es umfasst Instrumentenauswahl, Harmonisierung, Primärschätzung, Sensitivitätsanalysen und reproduzierbaren R-Code mit figurreifen Ausgaben zur Unterstützung transparenter Berichterstattung.

Zusammenfassung

Schilddrüsentoxikose ist eine endokrine Erkrankung, die durch einen Überschuss an Schilddrüsenhormonen gekennzeichnet ist, doch ihre ätiologischen Verbindungen zur systemischen Alterungsbiologie bleiben unvollständig definiert. Die Telomerlänge (TL) spiegelt die zelluläre Seneszenz und Genomstabilität wider und wurde mit mehreren komplexen Erkrankungen in Verbindung gebracht. Wir führten eine zweistichprobige Mendelsche Randomisierungsstudie (MR) durch, um die kausale Wirkung genetisch vorhergesagter TL auf das Risiko einer Thyreotoxikose zu bewerten. Genetische Instrumente für TL stammen aus einer großen, genomweiten Assoziationsstudie (GWAS) europäischer Abstammung (n > 470.000). Zusammenfassende Statistiken zur Schilddrüsentoxikose wurden aus der neuesten FinnGen-Veröffentlichung gewonnen (≈4.000 Fälle und >210.000 Kontrollen). Primäre invers-varianzgewichtete Analysen zeigten, dass längere genetisch proxiierte TL mit einem geringeren Risiko für Thyrotoxikose assoziiert ist und die Richtung sowie das Ausmaß des Effekts über komplementäre Schätzer hinweg konsistent waren (MR-Egger, gewichtete Median/maximale Wahrscheinlichkeit, MR-PRESSO und MR-RAPS). Sensitivitätsanalysen zeigten keine Hinweise auf gerichtete Pleiotropie, und Cochrans Q wurde zur Beurteilung der Heterogenität verwendet. Ein Steiger-Richtungstest unterstützte den kausalen Fluss von TL zur Thyreotoxose.

Soweit wir wissen, gehört diese Arbeit zu den ersten MRT-Analysen, die den kausalen Zusammenhang zwischen dem Gesamtrisiko für Thyrotoxikose und dem TL mit modernen GWAS-Ressourcen bewerten und damit die bisherigen Evidenzen auf hyperthyreose-bezogene Phänotypen erweitern. Diese Ergebnisse deuten darauf hin, dass zelluläre Alterungsprozesse, die von TL indiziert werden, zur Anfälligkeit für Schilddrüsentoxikose beitragen und zukünftige longitudinale und mechanistische Studien zur Telomerbiologie bei Schilddrüsenfunktionsstörungen anregen können.

Einleitung

Schilddrüsentoxikose ist eine häufige endokrine Erkrankung, die durch einen Überschuss an zirkulierenden Schilddrüsenhormonen gekennzeichnet ist. Typische Manifestationen sind Gewichtsverlust, Tachykardie, Reizbarkeit, Angst und Zittern, wobei schwere Fälle zu Arrhythmien, Herzinsufficienz und Knochenverlust^anfällig sind. Obwohl klinische Merkmale und therapeutische Strategien relativ gut etabliert sind, bleiben die biologischen Mechanismen, die die Inter-Individual-Anfälligkeit und -Entwicklung prägen, noch unvollständig definiert. Die Ätiologien reichen von Autoimmunhyperfunktion (z. B. Morbus Basedow) über toxische Knoten, Thyreoiditis, exogene Hormonexposition, Arzneimittelreaktionen bis hin zu tumorbedingten Ursachen². Zwischen diesen Entitäten wiederkehren Immunaktivierung, oxidativer Stress und veränderte mitochondriale Biogenese als plausible mechanistische Themen, die hormonellen Überschuss mit Gewebeschäden und systemischen Komplikationen in Verbindung bringen.

Telomerlänge (TL) – die Tandem-Repeat-DNA-Protein-Struktur, die die Chromosomenenden bedeckt – integriert zelluläre Alterung, replikative Geschichte und Genomstabilität. TL verkürzt sich mit der Zellteilung und wird durch oxidativen Stress und chronische Entzündungen beschleunigt, Prozesse, die ebenfalls bei der Schilddrüsenfunktionsstörung^{ausgeprägt sind 3,4}. Die kürzere TL wurde mit kardiometabolischen Störungen, bestimmten Krebsarten und neurodegenerativen Erkrankungen^{in Verbindung gebracht}, und neue Daten deuten auf Zusammenhänge mit endokrinen Phänotypen wie Fettleibigkeit und Typ-2-Diabetes hin, bei denen hormonelle und immunmetabolische Achsen mit replikativem Stress⁶ überschneiden. Da Schilddrüsenhormone den basalen Stoffwechsel, die mitochondriale Funktion und die proliferative Signalübertragung regulieren, könnte die TL-Dynamik besonders relevant in thyreotoxischen Zuständen sein, wo systemischer Hypermetabolismus und Immundysregulation sowohl Einfluss haben als auch von der Telomerbiologie ^7,8 beeinflusst werden könnten.

Ob TL eine ursächliche Rolle bei der Schilddrüsentoxikose spielt, ist jedoch weiterhin ungewiss. Konventionelle Beobachtungsstudien sind anfällig für Störungen (z. B. Lebensstil, Begleiterkrankungen) und umgekehrte Kausalität (Krankheit, die TL betrifft). Um diese Einschränkungen zu beheben, verwenden wir eine zweistichprobige Mendelsche Randomisierung (MR), die Keimlinien-genetische Varianten, die robust mit TL assoziiert sind, als unverwirrte Instrumente nutzt, um deren Wirkung auf das^{Krankheitsrisiko 9} zu testen. Genomweite Assoziationsstudien liefern nun ausreichend starke TL-Instrumente und große, gut phänotypisierte Enddatensätze, um eine rigorose kausale Bewertung^zu ermöglichen. Durch die Integration der genetischen Epidemiologie mit der Telomerbiologie bewertet diese Studie, ob genetisch proxied TL das Risiko einer Thyreotoxikose beeinflusst, klärt mechanistische Signalwege und untersucht TL als Kandidat-Biomarker für Risikostratifizierung und zukünftige Interventionsentwicklung.

Soweit wir wissen, gehört diese Analyse zu den ersten MRT-Studien, die das Gesamtrisiko einer Thyreotoxikose im Zusammenhang mit TL anhand moderner großflächiger GWAS-Ressourcen untersucht und damit frühere Arbeiten erweitern, die sich hauptsächlich auf hyperthyreosebedingte Phänotypen oder Autoimmunerkrankungen der Schilddrüse¹¹ konzentrierten.

Protokoll

Diese Studie verwendet ein zweistichprobiges Mendelsche Randomisierungsdesign (MR), um zu bewerten, ob genetisch proklamierte Telomerlänge (TL) kausal das Risiko einer Thyrotoxose beeinflusst. Es wurden nur anonymisierte, öffentlich zugängliche GWAS-Zusammenfassungsstatistiken verwendet, und es wurden keine Daten auf Einzelebene abgerufen. Das Institutional Review Board des Qingdao Municipal Hospital hat entschieden, dass die Analyse von einer weiteren Überprüfung ausgenommen ist, da sie ausschließlich auf öffentlichen Zusammenfassungsdaten basiert. Alle teilnehmenden genomweiten Assoziationsstudien erhielten im Rahmen ihrer ursprünglichen Protokolle eine informierte Einwilligung und ethische Genehmigungen. Die Analyse ist so strukturiert, dass sie Störungen und umgekehrte Kausalitäten minimiert, indem Keimbahnvarianten, die mit TL assoziiert sind, als unverwirrte Instrumente verwendet und Sensitivitätsverfahren implementiert werden, um Pleiotropie, Heterogenität und Wirkungsrichtung zu untersuchen.

Datenquellen
Der Expositionsdatensatz für TL wurde von IEU OpenGWAS unter der Identifikatorin ieu-b-4879 gewonnen und umfasst etwa 472.174 Teilnehmer europäischer Abstammung. Der Ergebnisdatensatz für Thyrotoxikose stammt aus dem FinnGen-Konsortium unter Verwendung des Endpunkts von 2021 finn-b-thyROTOXICOSIS mit 4.142 Fällen und 213.693 Kontrollen. Diese Ressourcen liefern effektive Schätzungen und Standardfehler, die erforderlich sind, um SNP-Instrumente zu konstruieren, Allele zu harmonisieren und kausale Effekte mit etablierten MR-Schätzern zu schätzen.

Software und Informatikumgebung
Alle Analysen wurden in R (Version 4.3.1) durchgeführt. Das wichtigste Analysepaket war TwoSampleMR (Version 0.5.7), ergänzt durch ieugwasr für programmatischen Zugriff auf GWAS-Ressourcen, MRPRESSO für Ausreißererkennung und Verzerrungstests, mr.raps für robuste Schätzungen unter schwachen Instrumenten und idiosynkratischer Pleiotropie, RadialMR für radiale Visualisierung sowie allgemeine Pakete wie tidyverse und data.table. Sitzungsinformationen und Paketversionen werden in eine Datei geschrieben, um eine strikte Reproduzierbarkeit sicherzustellen.

Kernannahmen der Mendelschen Randomisierung
Der MR-Rahmen setzt voraus, dass genetische Instrumente stark mit TL assoziiert sind, unabhängig von Faktoren sind, die die Beziehung zwischen TL und Thyrotoxose beeinträchtigen, und die Thyrotoksiose nur über TL und nicht über alternative Wege beeinflussen. Der Analyseplan operationalisiert diese Annahmen, indem er die Festigkeit pro Variant quantifiziert, eine unausgeglichene horizontale Pleiotropie mittels interzeptionbasierter Methoden und Ausreißer-Screening testet, die Heterogenität zwischen den Varianten bewertet und die Wirkungsrichtung von TL zu Thyrotoxikose mit einem formalen Richtungstest¹² bestätigt.

Instrumentenauswahl und Qualitätskontrolle
Genetische Instrumente wurden aus der TL GWAS mit einer genomweiten Signifikanzschwelle von p < 5×^10-8 ausgewählt. Um die Unabhängigkeit zu gewährleisten, wurde das Verknüpfungsungleichgewicht durch Klumpen bei r² = 0,001 innerhalb eines 10.000-kb-Fensters unter Verwendung europäischer Referenzdaten des 1000 Genomes Project angegangen; wo Varianten korreliert waren, wurde die Variante mit dem kleineren p-Wert für die Assoziation mit TL beibehalten. Die Instrumentenstärke wurde für jede Variante mit der Statistik F_i = \beta_{E,i}^{2}/\mathrm{SE}(\beta_{E,i})^{2} zusammengefasst und anhand des konventionellen F > 10-Kriteriums^13,14 bewertet. Wo angebracht, wurde die Sanderson-Windmeijer-Näherung verwendet, um die aggregierte Instrumentenstärke über mehrere Varianten in einer Stichprobe der Größe N zu beschreiben, wobei K-Instrumente die Expositionsvarianz R^{2} erklären: F_{\mathrm{SW}} = \{[R^{2}/(1-R^{2})]\,(N-K-1)\}/K ¹⁵. Mehrdeutige palindromische Varianten mit Allelfrequenzen nahe 0,5 wurden ausgeschlossen, sofern Allelfrequenzinformationen keine eindeutige Ausrichtung erlaubten.

Harmonisierung
Für jedes Instrument wurde die entsprechende Verbindung zur Thyreotoxikose aus FinnGen extrahiert und so ausgerichtet, dass die Effektgrößen das gleiche Effektallel über Exposition und Ergebnis repräsentieren. Die Harmonisierung entfernte Allel-Fehlanpassungen, korrigierte Strängprobleme und schloss palindromische SNPs mit ungelöster Mehrdeutigkeit aus, wodurch ein Datensatz für eine gültige Wald-Verhältnis-Konstruktion geeignet war.

Primäre kausale Schätzung
Die Primäranalyse verwendete den invers-varianzgewichteten (IVW) Schätzer, um SNP-spezifische Wald-Verhältnisse meta-zu analysieren und einen gepoolten kausalen Effekt¹⁶ zu erhalten. Schätzungen werden als Odds Ratios pro einer Standardabweichungserhöhung in TL angegeben, wobei 95%-Konfidenzintervalle unter Fest- und Zufallseffektmodellen abgeleitet werden. Die Modellwahl wurde durch Heterogenitätsdiagnostik geleitet, und beide Spezifikationen werden präsentiert, um eine robuste Interpretation zu ermöglichen.

Sensitivitätsanalysen und Robustheitsprüfungen
Die Robustheit wurde mittels MR-Egger-Regression mit einem Interzeptterm bewertet, um unausgewogene gerichtete Pleiotropie zu testen, Maximum-Likelihood-Schätzung zur Verbesserung der Effizienz unter Homogenität bei gleichzeitiger Minderung von Messfehlern, MR-PRESSO zur Durchführung eines globalen Ausreißertests und zur Schätzung von Verzerrungen und ausreißerkorrigierten Effekten sowie MR-RAPS zur Schätzung widerstandsfähiger Schätzungen gegen schwache Instrumente und idiosynkratische Pleiotropie¹⁷. Die Zwischen-Varianten-Heterogenität wurde mithilfe der Cochranschen Q-Statistik im IVW-Rahmen¹⁸ quantifiziert. Die Wirkungsrichtung wurde mit dem Steiger-Test untersucht, der den Anteil der in der Exposition erklärten Varianz und Ergebnis vergleicht, um festzustellen, ob die Daten besser mit TL kompatibel sind, die Thyrotoxikose verursacht, als umgekehrt. Der MR-Egger-Intercept wurde als formaler Test für richtungsbezogene Pleiotropie verwendet. Leave-One-Out-Analysen wurden untersucht, um sicherzustellen, dass die Gesamtassoziation nicht durch eine einzelne Variante verursacht wurde.

Umgekehrte Mendelsche Randomisierung
Um eine mögliche Rückkausalität zu untersuchen, wurde die analytische Pipeline wiederholt, wobei Thyreotoxikose die Exposition und TL das Ergebnis war. Die gleichen Instrumentenauswahlkriterien, Harmonisierungsverfahren, IVW-Primärschätzer und Sensitivitätsanalysen wurden angewandt, sodass Folge zur Richtungsmäßigkeit innerhalb eines identischen kausalen Rahmens gezogen werden kann.

Mehrfache Tests, Stromversorgung und Berichterstattung
Der primäre Hypothesentest bezieht sich auf den IVW-Schätzer für die Wirkung von TL auf die Schilddrüsentoxikose. Sensitivitätsschätzer und diagnostische Tests werden als unterstützende Evidenz interpretiert; p-Werte werden zur Klarheit in wissenschaftlicher Notation angegeben, und Schlussfolgerungen betonen Konsistenz zwischen den Methoden statt isolierte Signifikanzschwellen. Zusammenfassungen der Instrumentenstärke und der erklärte Anteil der Varianz informieren die ungefähre Leistung unter Standardformulierungen der Nichtzentralität und erkennen an, dass die Leistung von der Stichprobengröße, der Instrumentenstärke und der tatsächlichen Effektstärke abhängt.

Berechnungsreproduzierbarkeit und exakte Befehle
Die Reproduzierbarkeit wird gewährleistet, indem die vollständige Abfolge von R-Befehlen bereitgestellt wird, die die Instrumentenauswahl, Ergebnisextraktion, Harmonisierung, Primär- und Sensitivitätsanalysen, Reverse MR, Diagnoseausgaben und den Export analysefähiger Dateien nachbilden. Das Skript schreibt stabile, menschenlesbare CSV-Dateien, die der Instrumentenliste, dem harmonisierten Datensatz und der Zusammenfassung der MR-Schätzungen und -Diagnosen entsprechen.

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

bibliothek(TwoSampleMR)
Bibliothek(IEUGWASR)
Bibliothek (MRPRESSO)
Bibliothek (Mr. Raps)
Bibliothek (RadialMR)
Bibliothek (Tidyverse)
library(data.table)

# Belichtung: Telomerlänge (IEU OpenGWAS)
exposure_id <- "IE-B-4879"

# Ergebnis: FinnGen 2021 Thyrotoxikose Endpunkt verwendet für die berichteten Analysen
outcome_id <- "finn-b-thyROTOXICOSE"

# Instrumentenauswahl mit genomweiter Schwelle und strenger LD-Klumpung
exp <- extract_instruments(Ergebnisse = exposure_id, p1 = 5e-8, klump = TRUE, r2 = 0,001, kb = 10000)
exp$F_stat <- (exp$beta.exposure^2) / (exp$se.exposure^2)
fwrite(exp, "S1_instruments_TL.csv")

# Ergebnisextraktion und Harmonisierung
Out <- extract_outcome_data(SNPS = exp$SNP, Ergebnisse = outcome_id)
dat <- harmonise_data(exp, out, action = 2)
fwrite(dat, "S2_harmonised_TL_vs_thyrotoxicosis.csv")

# Primäre MR- und Sensitivitätsschätzer
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) # Cochrans Q
Pleio <- mr_pleiotropy_test(dat) # Egger abfangen
steiger <- directionality_test(dat) # Steiger Directionality

# MR-PRESSO globale und von Ausreißern korrigierte Schätzungen
mrpresso <- mr_presso(BetaOutcome = "beta.outcome", BetaExposure = "beta.exposure",
SdOutcome = "se.outcome", SdExposure = "se.exposure",
OUTLIERtest = TRUE, DISTORTIONtest = TRUE,
data = dat, NbDistribution = 1000, SignifThreshold = 0,05) #

# Umgekehrte MRT: Thyreotoksiose (Exposition) -> TL (Ergebnis)
rev_exp <- extract_instruments(Ergebnisse = outcome_id, p1 = 5e-8, Klump = TRUE, r2 = 0,001, kb = 10000)
rev_out <- extract_outcome_data(SNPS = rev_exp$SNP, Ergebnisse = exposure_id)
rev_dat <- harmonise_data(rev_exp, rev_out, Aktion = 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)

# Exporte für Archivierung und Figuren-/Tabellengenerierung
write.csv(bind_rows(res), "S3_mr_results_primary.csv", reihe.namen = FALSCH)
write.csv(het, "S3_mr_heterogeneity.csv", Zeilennamen = FALSCH)
write.csv(pleio, "S3_mr_pleiotropy_egger.csv",reihe.namen = FALSCH)
write.csv(steiger, "S3_mr_steiger.csv", Zeilennamen = FALSCH)
write.csv(bind_rows(rev_res), "S3_reverse_mr_results.csv", reihe.namen = FALSCH)
write.csv(rev_het, "S3_reverse_heterogeneity.csv", zeilen.namen = FALSCH)
write.csv(rev_pleio, "S3_reverse_pleiotropy_egger.csv", zeilenannamen = FALSCH)
write.csv(rev_steiger, "S3_reverse_steiger.csv", reihe.namen = FALSCH)

Ergebnisse

Instrumenteigenschaften und Annahmeprüfungen
Aus dem telomerlangen GWAS wurden 129 genomweite, signifikante, LD-unabhängige Instrumente (p < 5×^10-8; r² < 0,001; 10.000 kb Fenster) beibehalten. Die Per-SNP-F-Statistiken überschritten die konventionelle Schwelle für die große Mehrheit der Varianten, was auf ein begrenztes Risiko für schwache Instrumentenverzerrung hinweist. Die Harmonisierung entfernte mehrdeutige palindromische Varianten und ausgerichtete Allele, sodass Expo...

Diskussion

Diese Mendelsche Randomisierungsstudie (MR) liefert Hinweise darauf, dass länger genetisch vorhergesagte Telomere das Risiko für Thyreotoxikose¹⁹ senken. Richtungsübereinstimmende Schätzungen über IVW, MR-Egger, Maximum Likelihood, MR-PRESSO und MR-RAPS deuten auf eine robuste protektive Assoziation hin, und Richtungsprüfungen unterstützen einen TL → Thyrotoxikose statt der umgekehrten²⁰. Die biologische Plausibilität entspricht der Rolle d...

Materialien

Liste der in diesem Artikel verwendeten Materialien

Name	Unternehmen	Katalognummer	Kommentare
1000 Genomes Projekt	Referenzdaten für das Klumpen von Verknüpfungsungleichgewicht (LD).	1000 Genomes Projekt	-
Cochran' s Q-Statistik	Statistische Methode zur Heterogenitätsprüfung im MR.	N/A	-
FinnGen-Datensatz	GWAS-Zusammenfassungsstatistiken für Schilddrüsentoxikose (Fälle und Kontrollen).	FinnGen-Konsortium	Veröffentlichung 2021, finn-b-thyROTOXICOSIS
IEU OpenGWAS-Datensatz	GWAS-Zusammenfassungsstatistiken für die Telomerlänge.	IEU OpenGWAS	Datensatz-ID: ieu-b-4879
IEUGWASR-Paket	R-Paket, das für programmatischen Zugriff auf GWAS-Daten verwendet wird.	CRAN	Version 0.1.5
Mr. Raps-Paket	R-Paket für robuste Schätzung in der Mendelschen Randomisierung unter schwachen Instrumenten.	CRAN	Version 0.1.2
MR-Egger-Regression	Empfindlichkeitsanalysemethode zur Untersuchung der gerichteten Pleiotropie.	N/A	-
MR-PRESSO Ausreißertest	Statistischer Test, der verwendet wird, um Ausreißer in MR-Analysen zu identifizieren und zu korrigieren.	N/A	-
MRPRESSO-Paket	R-Paket, das für Ausreißererkennung und Verzerrungstests in der MR-Analyse verwendet wird.	CRAN	Version 1.0
R-Software	Statistische Rechen- und Grafikumgebung wird für alle Analysen verwendet.	R-Stiftung für statistische Berechnung	Version 4.3.1
RadialMR-Paket	R-Paket zur radialen Visualisierung der MR-Ergebnisse.	CRAN	Version 1.1.2
TwoSampleMR-Paket	R-Paket, das für Mendelsche Randomisierungsanalysen verwendet wird.	CRAN	Version 0.5.7