Alle experimentellen Verfahren in der vorliegenden Studie wurden gemäß den Richtlinien des National Institutes of Health zum Schutz und zur Verwendung von Laborversuchstieren durchgeführt und von dem Institutionellen Ausschuss für Tierversuche genehmigt (Genehmigungsnummer IACUC-2023-0156).
Tiere
Weibliche Sprague-Dawley-Ratten (n = 24, 6 Wochen alt, Körpergewicht 180–200 g) wurden vom Laboratoriumstierzentrum der angeschlossenen Einrichtung bezogen und unter kontrollierten Umweltbedingungen (Temperatur 22 ± 2°C, Luftfeuchtigkeit 50 ± 10%, 12-Stunden-Licht/Dunkel-Wechsel) gehalten, mit freiem Zugang zu Standardfutter für Nager und gereinigtem Wasser. Nach einer einwöchigen Akklimatisierungsphase wurden die Tiere zufällig drei Versuchsgruppen zugewiesen (n = 8 pro Gruppe): Kontrolle, PCOS-Modell und ZeLan-Behandlung. Die Stichprobengröße wurde basierend auf früheren Studien zur Untersuchung von Veränderungen der Darmmikrobiota in PCOS-Modellen festgelegt; eine Power-Analyse ergab, dass 8 Tiere pro Gruppe eine statistische Power von 80 % ermöglichen, einen 1,5-fachen Unterschied in den Alpha-Diversitätsindizes bei α = 0,05 nachzuweisen.
Herstellung des PCOS-Modells und Behandlungsprotokoll
Das PCOS-Modell wurde mithilfe eines gut validierten Protokolls etabliert, das die subkutane Injektion von Dehydroepiandrosteron (DHEA; Sigma-Aldrich, USA), gelöst in Sesamöl, in den dorsalen Nackenbereich in einer Dosis von 6 mg/100 g Körpergewicht täglich über 21 aufeinanderfolgende Tage umfasst. Kontrolltiere erhielten äquivalente Volumina des Sesamöl-Fahrzeugs. Der vollständige experimentelle Zeitplan umfasste 7 Tage Akklimatisierung, 21 Tage Modellinduktion und 28 Tage Intervention, wobei die Endpunktprobenahme am Tag 56 des Protokolls erfolgte, einschließlich der Akklimatisierungsphase (Abbildung 1). Der ZeLan-Extrakt wurde aus authentifiziertem Kraut von Lycopus lucidus Turcz. durch wässrige Extraktion hergestellt und standardisiert, um gesamte Flavonoide in einer Konzentration von mindestens 2,5 % gemäß Hochleistungsflüssigkeitschromatographie zu enthalten, was als chemisches Qualitätskontrollkriterium für den pflanzlichen Extrakt diente. Nach Abschluss der DHEA-induzierten Modellbildung erhielten die Ratten der ZeLan-Gruppe täglich über 28 aufeinanderfolgende Tage eine orale Gavage mit ZeLan-Extrakt in einer Dosis von 200 mg/kg Körpergewicht, während die Kontroll- und Modellgruppen äquivalente Volumina destillierten Wassers erhielten. Die Dosis von 200 mg/kg wurde als explorative pharmakologische Dosis basierend auf früheren Tierversuchserfahrungen mit pflanzlichen Extrakten und vorläufigen Überlegungen zur Verträglichkeit gewählt; eine formale Dosis-Wirkungs-Untersuchung war nicht Bestandteil des vorliegenden Versuchsdesigns. Das Körpergewicht wurde während des gesamten Experimentzeitraums wöchentlich erfasst. In diesem experimentellen Design wurde die erfolgreiche Etablierung des PCOS-Modells anhand des DHEA-Induktionsprotokolls sowie der endokrinen und metabolischen Phänotypen am Versuchsende bewertet, einschließlich Körpergewicht, Serum-Testosteron, LH/FSH-Verhältnis und nüchtern Insulin. Die histopathologische Untersuchung der Ovarien, die Follikelzählung, die Quantifizierung des Corpus luteum sowie die kontinuierliche Überwachung des Östruszyklus wurden nicht als Endpunkte zur Bestätigung des Modells einbezogen, sodass die Interpretation einer ovariellen Dysfunktion auf endokrinen und metabolischen Messgrößen und nicht auf einer direkten Gewebeanalyse beruhte.
Probensammlung und -aufbereitung
Proben (Serum und Fäzes) wurden von allen drei Gruppen (n = 8 pro Gruppe) zur anschließenden biochemischen, mikrobiotischen und metabolomischen Analyse entnommen. Die Endpunktprobenahmen erfolgten morgens nach einer nächtlichen Fastenperiode, um zirkadiane und ernährungsbedingte Variabilität zu minimieren. Blutproben wurden unter Anästhesie aus der Aorta abdominalis entnommen und bei 4 °C 15 Minuten lang bei 3.000 × g zentrifugiert, um das Serum zu gewinnen. Dieses wurde portioniert und bis zur Hormon- und Metabolomanalyse bei -80 °C gelagert. Frische Fäkalproben wurden direkt aus dem Zäkum unter aseptischen Bedingungen in sterile Kryoröhrchen gesammelt, während des Transports auf Trockeneis aufbewahrt, innerhalb von 10 Minuten nach der Entnahme in flüssigem Stickstoff schockgefroren und bis zur 16S-rRNA-Gen-Sequenzierung bei -80 °C gelagert.
Bestimmung von Hormonen im Serum
Die Serumkonzentrationen von Testosteron, luteinisierendem Hormon (LH), follikelstimulierendem Hormon (FSH) und Insulin wurden mithilfe kommerziell erhältlicher ELISA-Kits (Cusabio Biotech, Wuhan, China) gemäß den Herstelleranweisungen quantifiziert. Das LH/FSH-Verhältnis wurde als zusammengesetzter Index einer gonadotropen Dysbalance, die für das PCOS charakteristisch ist, berechnet. Alle Proben wurden im Duplikat analysiert, und Messungen mit einem Variationskoeffizienten der Duplikate von mehr als 15 % wurden wiederholt. Die intraassay- und interassay-Variationskoeffizienten mussten unter 10 % bzw. 15 % liegen. Mögliche Ausreißer wurden anhand der Originaldatensätze und der Informationen zur Assay-Performance überprüft, und kein Wert wurde ausschließlich aus statistischen Gründen ohne dokumentierten technischen Grund ausgeschlossen.
16S-rRNA-Gen-Sequenzierung und bioinformatische Analyse
Die gesamte bakterielle genomische DNA wurde aus Zäkum-Inhalten mit dem QIAamp DNA Stool Mini Kit (Qiagen, Deutschland) gemäß dem Herstellungsprotokoll isoliert. Die DNA-Qualität und -menge wurde mittels NanoDrop-Spektrophotometrie und Agarose-Gelelektrophorese überprüft, und die DNA-Extrakte wurden vor der Bibliothekskonstruktion bei -80 °C gelagert. Die hypervariablen Regionen V3-V4 des 16S-rRNA-Gens wurden mit universellen Primern 338F (5'-ACTCCTACGGGAGGCAGCA-3') und 806R (5'-GGACTACHVGGGTWTCTAAT-3') unter Hinzufügung von Illumina-Adaptersequenzen amplifiziert. Die Amplikon-Bibliotheken wurden konstruiert, gereinigt, quantifiziert, in äquimolaren Konzentrationen gepoolt und auf der Illumina NovaSeq 6000-Plattform (Novogene, Peking, China) mit paired-end 250-bp-Lesungen sequenziert.
Rohsequenzierungsdaten wurden mit der QIIME2-Pipeline (Version 2022.2) verarbeitet. Paar-endige Reads wurden demultiplexed, zusammengeführt, anhand einer Phred-Score-Schwelle von mindestens 20 qualitätsgefiltert und mit DADA2 bearbeitet, um chimäre Sequenzen zu entfernen und Amplikonsequenzvarianten (ASVs) zu generieren. Die taxonomische Klassifizierung erfolgte anhand der Referenzdatenbank SILVA 138 unter Verwendung eines Naive-Bayes-Klassifikators. Alpha-Diversitätsmetriken, darunter Shannon-Index, Simpson-Index und Chao1-Schätzer, wurden berechnet, um die Diversität und Artenvielfalt innerhalb der Proben zu bewerten. Ergebnisse zur Sequenzierungstiefe, einschließlich Lesezusammenfassungen, ASV-Häufigkeitstabellen, Good's Coverage und Dateien mit Rarefaktionskurven, wurden als Teil des im Abschnitt „Data Availability“ beschriebenen Mikrobiomdatensatzes betrachtet. Die lineare Diskriminanzanalyse der Effektgröße (LEfSe) wurde verwendet, um krankheitsassoziierte Taxa im Vergleich zwischen Kontrolle und Modellgruppe zu identifizieren, wobei eine LDA-Score-Schwelle von mehr als 2,0 und ein P < 0,05 angewandt wurden; die mit ZeLan assoziierte Umkehrung wurde durch Vergleich der Richtung der taxonomischen Häufigkeitsverschiebungen in der ZeLan-Gruppe relativ zur Modellgruppe bewertet.
Untargetierte Metabolomik-Analyse
Die serummetabolomische Profilierung wurde mittels ultra-hochleistungsflüssigkeitschromatographiegekoppelt mit Quadrupol-Time-of-Flight-Massenspektrometrie (UHPLC-Q-TOF-MS; Agilent 1290 Infinity II/6550, USA) durchgeführt. Serumproben (100 μL) wurden mit 400 μL kaltem Methanol:Acetonitril (1:1, v/v), das interne Standards enthielt, versetzt, verwirbelt und zur Proteinfällung eine Stunde bei -20 °C inkubiert. Nach Zentrifugation bei 14.000 × g für 15 Minuten wurden die Überstände unter einem Stickstoffstrom eingeengt und zur Analyse in 100 μL Acetonitril:Wasser (1:1, v/v) rekonstituiert. Die chromatographische Trennung erfolgte an einer ACQUITY UPLC HSS T3-Säule (2,1 × 100 mm, 1,8 μm) bei 40 °C unter Verwendung von Wasser mit 0,1 % Ameisensäure als mobile Phase A und Acetonitril mit 0,1 % Ameisensäure als mobile Phase B bei einer Flussrate von 0,30 mL/min. Zur Trennung der Metaboliten wurde ein Gradientenprogramm verwendet, das den Anteil der organischen Phase von 5 % auf 95 % erhöhte, gefolgt von einer Säulenspülung und Re-Äquilibration. Die massenspektrometrischen Daten wurden sowohl im positiven als auch im negativen Elektrospray-Ionisationsmodus im Bereich von m/z 50–1200 erfasst; zu den Quellparametern gehörten eine Kapillarspannung von 3,5 kV, ein Trockengasfluss von 10 L/min, eine Trockengastemperatur von 325 °C, ein Nebuliserdruck von 35 psi sowie eine Fragmentor-Spannung von 120 V.
Rohdatendateien wurden in das mzXML-Format konvertiert und mit dem XCMS-Paket in R zur Peakdetektion, Retentionszeit-Korrektur und Integration der Peakfläche verarbeitet. Qualitätskontrollproben, die durch Mischen gleicher Anteile aller Serum-Extrakte hergestellt wurden, wurden während der gesamten Analytesequenz eingefügt, um die Gerätestabilität zu überwachen. Merkmale mit übermäßigen fehlenden Werten oder instabiler Reproduzierbarkeit in der Qualitätskontrolle wurden vor der statistischen Modellierung entfernt. Die Metabolitenidentifizierung erfolgte durch Abgleich der exakten Masse, der Retentionszeit und der MS/MS-Fragmentierungsmuster mit der Human Metabolome Database (HMDB), METLIN und internen Spektralbibliotheken. Multivariate statistische Analysen wurden mit MetaboAnalyst 5.0 durchgeführt. Metabolitenmerkmale, die für die nachfolgende Berichterstattung beibehalten wurden, wurden anhand eines kombinierten Kriteriums priorisiert, das den multivariaten Diskriminierungsbeitrag, die Richtung der Fold-Change und die statistische Signifikanz einbezog; P-Werte wurden mittels der Benjamini-Hochberg-Methode angepasst, um Mehrfachvergleiche zu berücksichtigen. Die vollständige annotierte Berichtstabelle mit 52 Metabolitenmerkmalen ist als ergänzende Tabelle 1 enthalten. Die Pfad-Enrichment-Analyse wurde anhand der Kyoto Encyclopedia of Genes and Genomes (KEGG) Datenbank durchgeführt.
Statistische Analyse
Alle statistischen Analysen wurden mit SPSS Version 26.0 (IBM Corporation, Armonk, NY, USA) und R-Software (Version 4.2.1) durchgeführt. Stetige Daten wurden als Mittelwert ± Standardabweichung (SD) angegeben und vor der inferenzstatistischen Prüfung mittels des Shapiro-Wilk-Tests auf Normalverteilung überprüft. Bei normalverteilten Daten wurde eine einfaktorielle Varianzanalyse (ANOVA) mit Tukey’s post-hoc-Test für den Vergleich mehrerer Gruppen verwendet; bei nicht normalverteilten Daten kam der Kruskal-Wallis-Test mit Dunn’s post-hoc-Test zum Einsatz. Zusammenhänge zwischen Darmmikroben-Taxa und serumgebundenen Metaboliten wurden mittels des Spearman-Rangkorrelationskoeffizienten bewertet und als hierarchisch geclusterte Heatmaps dargestellt. Die Korrelations-P-Werte wurden für Mehrfachvergleiche mittels der False-Discovery-Rate-Methode nach Benjamini-Hochberg korrigiert, und nur Korrelationen, die in ihrer Richtung konsistent und statistisch signifikant waren, wurden interpretiert. Ein zweiseitiger P-Wert < 0,05 galt als statistisch signifikant, sofern nicht eine angepasste P-Wert-Schwelle angegeben wurde.