Studienaufbau
Diese Studie verwendete einen Two-Sample-MR-Ansatz, um die kausale Beziehung zwischen mit NETs assoziierten entzündlichen Merkmalen und Typ-2-Diabetes (T2DM) sowie dessen Komplikationen zu untersuchen. Der MR-Ansatz basiert auf drei zentralen Annahmen20: 1) die genetische Variante muss stark mit der Exposition assoziiert sein, 2) die genetische Variante darf den Endpunkt nur über die Exposition beeinflussen (Ausschluss von Pleiotropie), und 3) die genetische Variante darf nicht mit irgendwelchen Störfaktoren assoziiert sein. Für einen detaillierten Überblick über den Ablauf dieser Studie siehe Abbildung 1. Die Arbeit folgt den MR-STROBE-Richtlinien für die Berichterstattung von MR-Studien und gewährleistet so hohe Standards an Transparenz und Reproduzierbarkeit21.
Datenquellen
Die in dieser MR-Analyse verwendeten Daten aus genomweiten Assoziationsstudien (GWAS) stammten aus öffentlichen GWAS-Datenbanken. Hinsichtlich der summarischen GWAS-Daten zu Typ-2-Diabetes (T2DM) wurden die umfangreichste Metaanalyse einer europäischen Population mit 74.124 Fällen und 824.006 Kontrollen sowie eine prospektive, nested case-cohort-Studie aus Europa mit 9.978 Fällen und 12.348 Kontrollen einbezogen22. Für diabetische Komplikationen, einschließlich T2DM mit Nierenkomplikationen und T2DM mit peripheren zirkulatorischen Komplikationen, wurden in dieser Studie die zusammenfassenden GWAS-Statistiken aus Untersuchungen bezogen, bei denen die Fälle T2DM-Patienten mit der jeweiligen spezifischen Komplikation waren und die Kontrollen Personen ohne T2DM darstellten. Detaillierte Angaben zu den Datenquellen finden sich in Ergänzungstabelle 1.
GWAS-Daten für mit NETs assoziierte entzündliche Merkmale wurden aus dem GWAS-Katalog bezogen (siehe Zusätzliche Tabelle 2). Um die Heterogenität in der biologischen Spezifität zu berücksichtigen, kategorisierte diese Studie die eingeschlossenen Expositionen in zwei unterschiedliche Gruppen basierend auf ihren funktionellen Rollen in der Neutrophilen-Biologie: (1) Kernfaktoren der NETose: Diese Kategorie umfasst Marker, die direkt an der strukturellen Bildung von NETs oder am enzymatischen Prozess der Chromatin-Entfaltung beteiligt sind. Konkret wurden in dieser Studie NETs selbst, Myeloperoxidase (MPO), Neutrophile Elastase (NE) sowie der MPO-DNA-Komplex eingeschlossen. MPO und NE sind essentielle Enzyme für den Histonabbau und die Chromatin-Entfaltung, während der MPO-DNA-Komplex ein spezifischer Surrogatmarker für NETs ist23. (2) Mit NETs assoziierte entzündliche Mediatoren: Diese Kategorie umfasst Zytokine und Mediatoren, die als upstream regulierende Faktoren oder als downstream wirksame Effektoren eng mit der NETose verknüpft sind, jedoch auch an breiteren entzündlichen Signalwegen beteiligt sind. Zu dieser Gruppe gehören Interleukin-6 (IL-6)24, Tumornekrosefaktor-alpha (TNF-α)25, die Konzentration von neutrophilen Gelatinase-assoziiertem Lipocalin (NGAL)26 und der Faktor des zellulären Kommunikationsnetzwerks 1 (CCN1)27.
Die in dieser Studie verwendeten Daten stammten aus Datenbanken mit offenen Zugang oder aus früheren Studien veröffentlichten GWAS; daher war für diese Forschung keine ethische Genehmigung erforderlich.
Auswahl genetischer Instrumente im Zusammenhang mit NETs
Da die NET-Bildung einen dynamischen biologischen Prozess darstellt, der in herkömmlichen GWAS nicht direkt gemessen wird, verwendete diese Studie einen genbasierten Ansatz, um genetische Instrumente für mit NETs assoziierte entzündliche Merkmale zu identifizieren. Ein umfassender Satz von 257 Genen, von denen bekannt ist, dass sie entscheidend an der Bildung und Regulation von NETs beteiligt sind, wie beispielsweise MPO, sowie Gene, die für Histone und neutrophile Granulaproteine kodieren, wurde aus etablierten mechanistischen Studien und der publizierten Literatur zusammengestellt.28.
Auswahl der Instrumentenvariablen
Um robuste genetische Instrumente für mit NETs assoziierte entzündliche Merkmale zu identifizieren, verwendete diese Studie einen sequenziellen Filterungsworkflow unter Nutzung von GWAS-Zusammenfassungsstatistiken. Die spezifischen Arbeitsschritte waren wie folgt:
Erste Auswahl: Genetische Varianten, die signifikant mit den Expositionseigenschaften assoziiert waren, wurden29 basierend auf einer Signifikanzschwelle von P < 5 × 10⁻6 extrahiert.
Filterung nach der Frequenz der Minor-Allele (MAF): Um die statistische Aussagekraft sicherzustellen, wurden Einzelnukleotid-Polymorphismen (SNPs) mit einer MAF ≤ 0,01 ausgeschlossen30.
LD-Klumpenbildung (LD Clumping): Um Störungen durch Kopplungsungleichgewicht (Linkage Disequilibrium, LD) zu verringern, wurden unabhängige SNPs mithilfe der Klumpenbildungsfunktion ausgewählt, wobei die Parameter auf r2 < 0,001 innerhalb eines 10.000-kb-Fensters festgelegt wurden31.
Beurteilung der Instrumentenstärke: Die Stärke jedes verbliebenen IV wurde mittels der F-Statistik quantifiziert, berechnet als F = R2 × (N-2) / (1-R2). Es wurden nur SNPs mit einer F-Statistik > 10 beibehalten, um eine Verzerrung durch schwache Instrumente zu minimieren32.
Identifizierung von Proxy-SNPs und Datenharmonisierung
Um fehlende SNPs im Ergebnisdatensatz zu beheben, wurde in dieser Studie eine Stellvertretersuche und ein Datensynchronisationsprozess durchgeführt:
Proxiesubstitution: Wenn ein Ziel-SNP aus der Expositions-GWAS in der Outcome-GWAS nicht verfügbar war, nutzte diese Studie die LDproxy()-Funktion basierend auf dem europäischen Referenzpanel des 1,000 Genomes Projects. Ein Kandidaten-Proxy-SNP wurde nur dann ausgewählt, wenn er eine hohe Kopplungsunabhängigkeit (LD) mit dem ursprünglichen SNP aufwies (r2> 0,8). Falls kein geeigneter Proxy gefunden wurde, wurde der SNP ausgeschlossen.
Harmonisierung: Expositions- und Ergebnisdatensätze wurden mithilfe der Funktion harmonise_data() aus dem TwoSampleMR-Paket (R-Version 4.0.5) angeglichen. In dieser Studie wurde der Parameter action = 2 festgelegt, um alle SNPs automatisch auf den Vorwärtstrang auszurichten und palindromische SNPs mit mehrdeutiger Strangorientierung zu entfernen.
Überprüfung: Nach der Harmonisierung hat diese Studie den harmonisierten Datensatz manuell überprüft, um sicherzustellen, dass die Effekt-Allelfrequenzen zwischen den Expositions- und Outcomedaten übereinstimmten.
MR-Analyse und Sensitivitätsprüfungen
Die kausale Inferenz wurde mit dem TwoSampleMR-Paket (R-Version 4.0.5) durchgeführt.
Primäre und sekundäre Analysen: Diese Studie verwendete die Methode der inversen Varianzgewichtung (Inverse Variance Weighted, IVW) als Hauptansatz33. Ergänzende Analysen wurden mithilfe der MR-Egger34-, gewichteten Median- und gewichteten Modus-Methode durchgeführt, um die Robustheit zu gewährleisten35.
Sensitivitätsanalysen: Die Heterogenität zwischen den Instrumentalvariablen wurde mithilfe des Cochran’schen Q-Tests36 über die Funktion mr_heterogeneity() bewertet. Horizontale Pleiotropie wurde mit dem MR-Egger34-Intercept-Test (mr_pleiotropy_test()) überprüft.
Ausreißererkennung: Diese Studie verwendete das MR-PRESSO-Paket zur Erkennung potenzieller Ausreißer37. Die Funktion mr_presso() wurde mit 1.000 Simulationen durchgeführt. Identifizierte Ausreißer-SNPs (P < 0,05) wurden entfernt, und die kausalen Schätzungen wurden erneut berechnet, um die Stabilität der Ergebnisse zu überprüfen. Zusätzlich wurde eine Leave-One-Out-Analyse durchgeführt, um sicherzustellen, dass der kausale Zusammenhang nicht durch einen einzelnen SNP verursacht wurde38.
Statistische Korrektur
Um Mehrfachprüfungen zu berücksichtigen, wurden die aus den MR-Analysen erhaltenen P-Werte mittels der False-Discovery-Rate-(FDR-)Methode korrigiert. Dies erfolgte unter Verwendung der Funktion P.adjust() in R mit dem Parameter method = „fdr“. Assoziationen mit einem korrigierten P-Wert (PFDR) < 0,05 galten als statistisch signifikant.