Die Mendelian Randomization (MR)-Analysen nutzten Daten aus einem öffentlich zugänglichen GWAS (https://gwas.mrcieu.ac.uk/), das bereits während der ersten Datenerhebung die erforderlichen ethischen und informierten Einwilligungsgenehmigungen eingeholt hatte. Folglich war für die Studie keine weitere ethische oder informierte Zustimmung erforderlich. Die Software und die verwendeten Datenbanken sind in der Materialtabelle aufgeführt.
1. Expositionsdaten
GWAS-Zusammenfassungsinformationen zur relativen Makronährstoffaufnahme in der Ernährung mit Fett, Eiweiß, Kohlenhydraten und Zucker wurden von über 235.000 Personen europäischerAbstammung erhalten. Die Ernährungsaufnahme des Vortages (UKB) oder der Gewohnheitsaufnahme (alle anderen Kohorten) wurde mithilfe umfassender Fragebögen erfasst. Die relative Aufnahme jedes Makronährstoffs wurde durch die Gesamtenergieaufnahme korrigiert, was nichtlineare Effekte ermöglichte. Nach strengen Qualitätskontrollverfahren beträgt die endgültige aggregierte Stichprobe 235.391 für Zucker und 268.922 für Fett, Eiweiß undKohlenhydrate 23. Die grundlegenden Merkmale der im GWAS enthaltenen Teilnehmer der relativen Makronährstoffaufnahme in der Ernährung sind in Tabelle 1 zusammengefasst.
Genetische Assoziationsdaten für Mikronährstoffe wurden aus einer umfassenden Übersicht der vorhandenen Literatur und der Verfügbarkeit robuster GWASs in europäischen Populationen gewonnen. Kurz zusammengefasst wurden GWAS-Zusammenfassungsstatistiken für sieben essentielle Mineralien enthalten: Serumcalcium, Kupfer, Eisen, Magnesium, Phosphor, Selen (gemessen in Blut- oder Fußnagelproben) und Zink 24,25,26,27,28,29,30,31sowie sieben Vitamine, darunter 25-Hydroxyvitamin D, Vitamin A1 (Retinol), Vitamin B6, Vitamin B9 (Folsäure), Vitamin B12, Vitamin C und Vitamin E 32,33,34,35,36,37. Die grundlegenden Merkmale der im GWAS enthaltenen Teilnehmer der Mikronährstoffe sind in Tabelle 1 zusammengefasst.
2. Ergebnisdaten
Die Ergebnisdaten stammen aus dem IEU GWAS Project (https://gwas.mrcieu.ac.uk/), das an der MRC Integrative Epidemiology Unit der University of Bristol entwickelt wurde. Konkret wurde Myopie als Ergebnis mit Daten aus der Quelle: UKB-b-6353 (n = 460.536), Katarakten aus der Quelle: ebi-a-GCST90018814 (n = 491.877), Glaukom aus Quelle: ebi-a-GCST90013865 (n = 406.927) und frühaltersbedingter makuladegeneration von der Quelle: ebi-a-GCST010723 (n = 105.248) untersucht. Die grundlegenden Merkmale der im GWAS enthaltenen Teilnehmer von Myopie und anderen schwerwiegenden Augenerkrankungen sind in Tabelle 2 zusammengefasst. Die Definitionen der Phänotypen von Augenerkrankungen basierten auf den ursprünglichen GWAS-Datensätzen. Myopie (ukb-b-6353) wurde aus der UK Biobank abgeleitet und basiert auf der selbst gemeldeten Verwendung von Brillen oder Kontaktlinsen, hauptsächlich für die Fernsicht (Kurzsichtigkeit). Katarakt- und Glaukom-Ergebnisse wurden aus großen GWAS-Metaanalysen auf klinisch diagnostizierten Fällen gewonnen, die durch medizinische Unterlagen oder Krankenhausepisodenstatistiken identifiziert wurden. Die frühzeitige altersbedingte Makuladegeneration (AMD) wurde gemäß standardisierten diagnostischen Kriterien definiert, die in den ursprünglichen GWAS-Konsortiumsstudien verwendet wurden. Unterschiede in Phänotypdefinitionen zwischen den Studien können eine gewisse Heterogenität verursachen und sollten bei der Interpretation der Ergebnisse berücksichtigt werden.
3. Auswahl instrumenteller Variablen
Einzelne Nukleotid-Polymorphismen (SNPs), die mit Expositionen bei genomweiter Bedeutung assoziiert sind (P < 5 × 10⁻8), wurden als instrumentelle Variablen ausgewählt. Unabhängige SNPs wurden als instrumentelle Variablen mit einem r2-Schwellenwert von < 0,001 ausgewählt, basierend auf ihrem Abstand innerhalb von 10.000 kb in den 1000G-europäischen Daten, unter Verwendung des "TwoSampleMR"-R-Pakets. Die Stärke der ausgewählten Instrumente wurde anhand der F-Statistik38 bewertet. Eine F-Statistik größer als 10 gilt im Allgemeinen als Hinweis auf starke instrumentelle Variablen39. Alle Analysen wurden harmonisiert, um die Ausrichtung der Effektallele zwischen Exposition und Ergebnisdatensätzen sicherzustellen.
4. Mendelsche Randomisierungsschätzungen
Um die kausale Beziehung zwischen Makronährstoffen und Mikronährstoffen sowie Myopie und anderen schweren Augenerkrankungen abzuschätzen, wurde das R-Paket "TwoSampleMR" für die MR-Analyse angewendet. Genetische Varianten, die mit jedem der Makronährstoffe und Mikronährstoffe in der Ernährung assoziiert sind, wurden als instrumentelle Variablen verwendet, um deren kausale Beziehung zu Myopie und anderen schweren Augenerkrankungen zu untersuchen. MR-Analyse basiert auf drei zentralen Annahmen: (1) der Relevanzannahme, bei der die genetischen Varianten stark mit der Exposition assoziiert sind; (2) die Annahme der Unabhängigkeit, bei der die Varianten unabhängig von Störfaktoren sind; und (3) die Annahme der Ausschlussbeschränkung, bei der die Varianten das Ergebnis nur durch die Offenlegung von Interesse beeinflussen. Es wurden mehrere MR-Methoden implementiert, darunter inverse Varian-gewichtete (IVW), gewichtete Median, gewichteter Modus und MR-Egger, wobei IVW der primäre Ansatz und die anderen Methoden zur Bewertung der Robustheit der Ergebnisse verwendet wurden. Die IVW-Methode liefert die präziseste kausale Schätzung, wenn alle instrumentellen Variablen gültig sind oder wenn die horizontale Pleiotropie zwischen verschiedenen Varianten ausgewogen ist. Der gewichtete Median-Schätzer kann konsistente Schätzungen liefern, selbst wenn bis zu 50 % der genetischen Instrumente ungültig sind. Die MR-Egger-Regressionsmethode ermöglicht den Nachweis der gerichteten Pleiotropie, indem ein Interzeptterm geschätzt wird, der den durchschnittlichen pleiotropen Effekt über genetische Varianten hinweg widerspiegelt. Anfangs wurden kausale Schätzungen mit IVW-Modellen mit festen Effekten berechnet, und zufällige IVW-Modelle wurden verwendet, wenn signifikante Heterogenität festgestellt wurde (S. < 0,05). Um mehrere Vergleiche über verschiedene Expositionen und Ergebnisse hinweg zu berücksichtigen, wurde die False Discovery Rate (FDR)-Korrektur mit der Benjamini–Hochberg-Methode angewendet. Die statistische Signifikanz wurde als FDR-bereinigtes P < 0,05 definiert, während P < 0,05 als nominell signifikant galt. Alle Analysen wurden in R (Version 4.1.2) durchgeführt. Um den analytischen Rahmen zu verdeutlichen, ist in Abbildung 1 eine schematische Darstellung des MR-Designs und der Annahmen dargestellt.
5. Heterogenitäts- und Sensibilitätsanalysen
Um die Robustheit der MR-Ergebnisse zu bewerten, wurden Heterogenitätstests durchgeführt, um mögliche Inkonsistenzen zwischen den instrumentellen Variablen zu erkennen. Die Heterogenität wurde anhand der Cochran-Q-Statistik bewertet, wobei ein P-Wert von < 0,05 signifikante Heterogenität40 anzeigte. Zusätzlich wurden Sensitivitätsanalysen durchgeführt, um potenzielle Pleiotropie zu bewerten, bei der genetische Varianten das Ergebnis über biologische Wege beeinflussen könnten, die nicht über die interessierende Exposition hinausgehen. Die gerichtete Pleiotropie wurde mit dem MR-Egger-Regressionsinterceptionstest41 untersucht. Zusätzlich wurde die Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO)-Methode angewandt, um potenzielle pleiotrope Ausreißer zu identifizieren und zu entfernen (S. < 0,05) sowie um bei Bedarf korrigierte kausale Schätzungen bereitzustellen41.