Wpływ ekstraktu z ZeLan na fenotypy kliniczne PCOS
Podawanie DHEA wywołało u szczurów fenotypy endokrynne i metaboliczne przypominające PCOS, co objawiło się znacznymi wzrostami masy ciała, stężenia testosteronu w surowicy, stosunku LH/FSH oraz poziomu insuliny na czczo w porównaniu z grupą kontrolną (wszystkie P < 0,05). Ocena rozkładu z wykorzystaniem testu Shapiro-Wilka potwierdziła zasadność analizy parametrycznej dla głównych punktów końcowych dotyczących masy ciała i hormonów, zestawionych w Tabeli 1. Leczenie ekstraktem z ZeLan przez 28 dni znacząco poprawiło te zaburzenia metaboliczne i hormonalne. Masę ciała odnotowano na koniec okresu leczenia w 56. dniu protokołu, wliczając 7-dniowy okres aklimatyzacji. Jak przedstawiono w Tabeli 1, masa ciała w grupie Model (291,12 ± 15,35 g) była znacznie wyższa niż w grupie kontrolnej (227,45 ± 10,12 g, P < 0,05), natomiast leczenie ZeLan znacząco zredukowało masę ciała do 245,25 ± 10,21 g (P < 0,05 względem Model). Poziomy testosteronu w surowicy były wyraźnie podwyższone u szczurów w modelu PCOS (1,85 ± 0,21 ng/mL) w porównaniu z grupą kontrolną (0,51 ± 0,05 ng/mL), a interwencja z użyciem ZeLan obniżyła poziom testosteronu do 0,92 ± 0,12 ng/mL. Stosunek LH/FSH był znacząco zwiększony w grupie Model (3,28 ± 0,25) w stosunku do grupy kontrolnej (1,25 ± 0,11) i został częściowo znormalizowany przez leczenie ZeLan (1,85 ± 0,15). Poziomy insuliny na czczo były podwyższone u zwierząt modelowych (28,95 ± 2,35 mIU/L względem 12,55 ± 1,12 mIU/L w grupie kontrolnej) i uległy znacznemu obniżeniu po podaniu ZeLan (18,82 ± 1,55 mIU/L) (Rysunek 2). Wyniki te potwierdziły pomyślne wywołanie endokrynno-metabolicznego fenotypu przypominającego PCOS; jednak w ninianzym badaniu nie przeprowadzono histopatologii jajników, zliczania pęcherzyków, analizy ciał żółtych ani monitorowania cyklu rujowego.
Ekstrakt ZeLan przywrócił różnorodność alfa mikrobioty jelitowej
Analiza wskaźników alfa-różnorodności ujawniła istotne zmiany w strukturze mikrobiomu jelitowego związane z PCOS oraz jego leczeniem. Jak pokazano w Tabeli 2 oraz na Rysunku 3, indeks różnorodności Shannona był znacząco niższy u szczurów w modelu PCOS (4,19 ± 0,35) niż w grupie kontrolnej (6,51 ± 0,28; P < 0,05), co wskazuje na zmniejszenie alfa-różnorodności mikrobiologicznej w stanie chorobowym. Leczenie preparatem ZeLan częściowo przywróciło różnorodność mikrobiologiczną, a indeks Shannona wzrósł do 5,79 ± 0,26 (P < 0,05 w porównaniu z grupą Model). Spójne wyniki zaobserwowano dla indeksu Simpsona, który spadł z 0,98 ± 0,01 w grupie kontrolnej do 0,82 ± 0,03 u zwierząt w modelu i powrócił do poziomu 0,93 ± 0,02 po interwencji z użyciem ZeLan. Estymator Chao1, odzwierciedlający bogactwo gatunkowe, wykazał podobny wzorzec: indukcja PCOS znacząco obniżyła wartość Chao1 z 1215,3 ± 45,4 do 665,8 ± 52,1, podczas gdy leczenie preparatem ZeLan zwiększyło bogactwo do 963,5 ± 38,6. Wyniki te wspólnie wskazują, że PCOS wiąże się z dysbiozą mikrobiomu jelitowego charakteryzującą się zmniejszoną alfa-różnorodnością i bogactwem gatunkowym, a ekstrakt ZeLan może częściowo przywrócić strukturę społeczności mikrobiologicznej w kierunku profilu grupy kontrolnej.
Zmiany w składzie mikrobioty jelitowej
Profilowanie taksonomiczne na poziomie typu wykazało wyraźne różnice w składzie między grupami eksperymentalnymi. Jak przedstawiono na Rysunku 4, Firmicutes i Bacteroidetes były dominującymi typami we wszystkich grupach, stanowiąc łącznie ponad 85% całej społeczności bakteryjnej. W grupie Model odnotowano znacznie zwiększoną liczność Firmicutes (około 75%) przy jednoczesnym zmniejszeniu proporcji Bacteroidetes (około 15%) w porównaniu z grupą Kontrolną (Firmicutes około 40%, Bacteroidetes około 45%), co skutkowało wyraźnym wzrostem stosunku Firmicutes do Bacteroidetes (F/B), charakterystycznym dla dysbiozy metabolicznej. Leczenie preparatem ZeLan zmodulowało tę nierównowagę, redukując liczność Firmicutes do około 52% i zwiększając liczność Bacteroidetes do około 35%, co doprowadziło do częściowej normalizacji stosunku F/B w kierunku poziomów kontrolnych.
Do identyfikacji taksonów bakteryjnych odróżniających grupę kontrolną od grupy modelowej oraz zdefiniowania sygnatury mikrobiologicznej związanej z chorobą wykorzystano analizę LEfSe. Jak pokazano na rycinie 5, na poziomie rodzaju w grupie modelowej PCOS zaobserwowano znacznie podwyższoną liczebność Bacteroides (wynik LDA > 3,5), Prevotella (wynik LDA > 3,0) oraz Escherichia (wynik LDA > 2,5), które są zazwyczaj wiązane z dysfunkcją metaboliczną i stanami zapalnymi. Z kolei w grupie kontrolnej stwierdzono wzbogacenie w rodzaje korzystne, w tym Lactobacillus (wynik LDA > 3,0) i Ruminococcus (wynik LDA > 2,5), uznawane za istotne w utrzymaniu integralności bariery jelitowej i produkcji korzystnych metabolitów, takich jak krótkołańcuchowe kwasy tłuszczowe. Grupa ZeLan nie została przedstawiona jako oddzielna klasa LEfSe na rycinie 5, ponieważ analiza ta koncentrowała się na taksonach odróżniających stan chorobowy od zdrowego; efekty związane z ZeLan zinterpretowano na podstawie przesunięć liczebności w kierunku profilu grupy kontrolnej oraz korelacji mikrobiom-metabolom.
Profilowanie metabolomiczne i identyfikacja różnicowych metabolitów
Niekierowana analiza metabolomiczna wykazała istotne zaburzenia metaboliczne związane z PCOS oraz jego leczeniem ekstraktem z ZeLan. Wykres wulkaniczny w Rycina 6 przedstawia rozkład cech metabolitów pomiędzy grupą modelową a kontrolną, z wyróżnieniem cech o zwiększonej i zmniejszonej ekspresji na podstawie kierunku zmiany krotności (fold-change) oraz istotności statystycznej. W kompletnej tabeli raportu metabolomicznego dla porównania grupy modelowej z kontrolną zachowano łącznie 52 adnotowane cechy metabolitów. Pełna lista, obejmująca nazwę metabolitu, identyfikator HMDB, klasę biochemiczną, log2 fold change, wartość P oraz kierunek regulacji, została przedstawiona jako Tabela uzupełniająca 1. Tabela 3 przedstawia reprezentatywne adnotowane metabolity o wysokim znaczeniu biologicznym dla PCOS.
Tabela 3 podsumowuje reprezentatywne, kluczowe metabolity różnicowe o wysokich wartościach VIP i istotności biologicznej dla patofizjologii PCOS. Hormony steroidowe, w tym testosteron (VIP = 2.56, P < 0.001) i androstendion (VIP = 2.34, P < 0.001), były znacznie podwyższone u zwierząt w modelu PCOS, co jest zgodne ze stanem hiperandrogenizmu charakterystycznym dla tego schorzenia. Zaburzony był również metabolizm kwasów tłuszczowych, przy czym kwas palmitynowy (VIP = 2.12) i kwas arachidonowy (VIP = 1.89) wykazały istotne wzrosty, podczas gdy prozapalny mediator lipidowy prostaglandyna E2 był podobnie podwyższony (VIP = 2.05). Odwrotnie, metabolity związane z korzystną aktywnością mikrobioty jelitowej były zmniejszone u zwierząt z PCOS: kwas litocholowy, wtórny kwas żółciowy powstający w wyniku biotransformacji mikrobiologicznej, uległ znacznemu obniżeniu (VIP = 2.45), podobnie jak krótkołańcuchowy kwas tłuszczowy – kwas masłowy (VIP = 2.21). Dodatkowo, w grupie Model odnotowano spadek poziomu kwasu bursztynowego oraz niezbędnego aminokwasu L-tryptofanu, co sugeruje upośledzenie fermentacji mikrobiologicznej i zaburzony metabolizm aminokwasów.
Analiza wzbogacenia szlaków KEGG
Analiza wzbogacenia ścieżek z wykorzystaniem bazy danych KEGG ujawniła kilka szlaków metabolicznych w surowicy, które były istotnie zaburzone w PCOS i modulowane przez leczenie ZeLan. Jak przedstawiono na Rysunku 7, biosynteza hormonów steroidowych okazała się najbardziej wzbogaconą ścieżką (Rich Factor = 0.80, P < 0.01), co jest zgodne z fenotypem hiperandrogenowym w PCOS. Ścieżka steroidogenezy jajnikowej była podobnie silnie wzbogacona (Rich Factor = 0.75, P < 0.01), co odzwierciedla udział zaburzonej syntezy androgenów w dysfunkcji jajników. Biosynteza pierwotnych kwasów żółciowych wykazała istotne wzbogacenie (Rich Factor = 0.62, P < 0.02), co wskazuje na zaburzony krążenie jelitowo-wątrobowe kwasów żółciowych, które może być powiązane ze zmianami w mikrobiocie jelitowej. Metabolizm kwasu arachidonowego (Rich Factor = 0.56, P < 0.03) oraz ścieżka insulinooporności (Rich Factor = 0.40, P < 0.04) były również istotnie wzbogacone, co potwierdza występowanie stanów zapalnych i zaburzeń metabolicznych charakterystycznych dla PCOS. Ścieżka biosyntezy kwasów tłuszczowych wykazała umiarkowane wzbogacenie (Rich Factor = 0.30), co sugeruje udział dysregulacji metabolizmu lipidów.
Analiza korelacji mikrobiomu i metabolomu
W celu oceny zależności funkcjonalnych między składem mikrobioty jelitowej a zmianami metabolicznymi przeprowadzono analizę korelacji Spearmana z kontrolą stopy fałszywych odkryć (false-discovery rate) pomiędzy rodzajami bakterii o różnym stopniu liczebności a kluczowymi metabolitami. Hierarchicznie zgrupowana mapa ciepła korelacji przedstawiona na Rysunku 8 ujawnia wyraźne wzorce powiązań między mikrobiotą a metabolomem. Bakterie korzystne, w tym Lactobacillus i Ruminococcus, wykazały silne korelacje dodatnie z kwasem litocholowym (odpowiednio r = 0.72 i 0.68) oraz kwasem masłowym (odpowiednio r = 0.65 i 0.60), wykazując jednocześnie korelacje ujemne z testosteronem (r = -0.75 i -0.80) oraz androstendionem (r = -0.68 i -0.72). Wyniki te sugerują, że bakterie korzystne zidentyfikowane w analizie taksonomicznej były powiązane z metabolitami związanymi z androgenami oraz aktywnością metaboliczną związaną z SCFA.
Z kolei rodzaje bakterii związane z PCOS, takie jak Bacteroides, Prevotella i Escherichia, wykazały przeciwstawne wzorce korelacji. Bakterie te korelowały dodatnio z testosteronem (odpowiednio r = 0.85, 0.82 i 0.60), androstenedionem (r = 0.78, 0.75 i 0.55) oraz kwasem palmitynowym (r = 0.65, 0.60 i 0.45), wykazując jednocześnie ujemne korelacje z korzystnymi metabolitami, w tym z kwasem litocholowym i kwasem masłowym. Te wzorce korelacji wspierają istnienie powiązań funkcjonalnych między konkretnymi taksonami mikroflory jelitowej a zaburzeniami metabolicznymi charakterystycznymi dla PCOS, jednak nie dowodzą one bezpośredniego związku przyczynowo-skutkowego.
Dostępność danych:
Zbiory danych wygenerowane i przeanalizowane w ramach niniejszego badania obejmują metadane próbek, surowe dane fenotypowe, wskaźniki alfa-różnorodności 16S, zestawienia liczebności na poziomie typu, wyniki różnicowej analizy taksonów bakteryjnych, tabele różnicowej analizy metabolomicznej, reprezentatywne wyniki metabolitów, wyniki wzbogacenia ścieżek KEGG oraz macierz korelacji mikrobiota-metabolit. Te podstawowe dane źródłowe zostały udostępnione jako towarzyszący Supplementary File 1. Przetworzony zbiór danych dotyczący różnicowych metabolitów został przedstawiony w Supplementary Table 1.

Rysunek 1: Schemat badania i harmonogram leczenia. Schematyczna ilustracja przebiegu eksperymentu przedstawiająca 7-dniowy okres aklimatyzacji, 21-dniowy okres iniekcji DHEA lub nośnika, 28-dniowe podawanie ekstraktu ZeLan lub nośnika oraz końcowe pobranie próbek surowicy i kału w 56. dniu protokołu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 2: Wpływ ekstraktu z ZeLan na parametry endokrynne i metaboliczne. Wykresy pudełkowe przedstawiające (A) stężenie testosteronu w surowicy, (B) stosunek LH/FSH oraz (C) masę ciała w grupie kontrolnej, modelowej i grupie leczonej ZeLan po okresie interwencji. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 3: Alfa-różnorodność mikrobioty jelitowej. Porównanie indeksu różnorodności Shannona pomiędzy grupami Control, Model i ZeLan, ilustrujące różnice w alfa-różnorodności mikrobioty jelitowej. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 4: Skład mikrobioty jelitowej na poziomie typu. Skumulowany wykres słupkowy przedstawiający względną obfitość dominujących typów bakterii w mikrobiocie jelitowej grup Control, Model i ZeLan. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5: Różnicowe taksony bakteryjne zidentyfikowane za pomocą analizy LEfSe. Analiza wielkości efektu liniowego dyskryminanta (LEfSe) identyfikująca rodzaje bakterii odróżniające grupę kontrolną od grupy modelowej. Wyniki LDA wskazują wielkość efektu taksonów o różnej liczebności. Zmiany związane z ZeLan opisano w sekcji Wyniki na podstawie przesunięć liczebności w stosunku do grupy modelowej. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 6: Różnicowe profilowanie metabolitów. Wykres wulkaniczny przedstawiający różnice w cechach metabolitów między grupą Model a grupą Kontrolną. Czerwone i niebieskie punkty reprezentują odpowiednio istotnie zwiększone i zmniejszone poziomy cech metabolitów, w oparciu o krotność zmiany (fold change) oraz istotność statystyczną. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 7: Analiza wzbogacenia ścieżek KEGG. Wykres bąbelkowy przedstawiający istotnie wzbogacone ścieżki metaboliczne zidentyfikowane na podstawie metabolitów różnicowych. Rozmiar bąbelka reprezentuje liczbę zmapowanych metabolitów, a kolor bąbelka wskazuje istotność statystyczną (P-value) wzbogacenia ścieżki. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 8: Korelacja między mikrobiotą jelitową a metabolitami w surowicy. Hierarchicznie zgrupowana mapa ciepła przedstawiająca współczynniki korelacji Spearmana pomiędzy rodzajami bakterii o zróżnicowanej liczebności a reprezentatywnymi metabolitami w surowicy. Kolor czerwony oznacza korelacje dodatnie, natomiast kolor niebieski oznacza korelacje ujemne. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Grupa | N | Masa ciała (g) | Testosteron (ng/mL) | Stosunek LH/FSH | Insulina (mIU/L) |
| Kontrola | 8 | 227.45 ± 10.12 | 0.51 ± 0.05 | 1.25 ± 0.11 | 12.55 ± 1.12 |
| Model | 8 | 291.12 ± 15.35a | 1.85 ± 0.21a | 3.28 ± 0.25a | 28.95 ± 2.35a |
| ZeLan | 8 | 245.25 ± 10.21b | 0.92 ± 0.12b | 1.85 ± 0.15b | 18.82 ± 1.55b |
| Uwaga: Dane przedstawiono jako średnia ± SD. a = p < 0,05 w porównaniu z grupą kontrolną; b = p < 0,05 w porównaniu z grupą Model. |
Tabela 1: Parametry hormonalne i metaboliczne u szczurów z PCOS. Zmierzono masę ciała, poziom testosteronu w surowicy, stosunek LH/FSH oraz poziom insuliny na czczo w grupach Kontrolnej, Modelu oraz leczonej preparatem ZeLan. Dane przedstawiono jako średnia ± SD.
| Grupa | Indeks Shannona | Indeks Simpsona | Estymator Chao1 |
| Kontrola | 6.51 ± 0.28 | 0.98 ± 0.01 | 1215.3 ± 45.4 |
| Model | 4.19 ± 0.35a | 0.82 ± 0.03a | 665.8 ± 52.1a |
| ZeLan | 5.79 ± 0.26b | 0.93 ± 0.02b | 963.5 ± 38.6b |
| Uwaga: Dane przedstawiono jako średnia ± SD. a = p < 0,05 w porównaniu z grupą kontrolną; b = p < 0,05 w porównaniu z grupą modelową. |
Tabela 2: Wskaźniki alfa-różnorodności mikrobiomu jelitowego. Porównanie indeksu Shannona, indeksu Simpsona i estymatora Chao1 pomiędzy grupą kontrolną (Control), grupą modelową (Model) a grupą leczoną ZeLan. Dane przedstawiono jako średnia ± SD.
| Nazwa metabolitu | ID HMDB | Klasa | Trend (model vs kontrola) | wartość p | Wskaźnik VIP |
| Testosteron | HMDB0000286 | Steroid | W górę ↑ | < 0.001 | 2.56 |
| Androstenedion | HMDB0000055 | Steryd | W górę ↑ | < 0.001 | 2.34 |
| Kwas palmitynowy | HMDB0000220 | Kwas tłuszczowy | W górę ↑ | < 0.001 | 2.12 |
| Kwas arachidonowy | HMDB0001043 | Kwas tłuszczowy | W górę ↑ | 0.003 | 1.89 |
| Prostaglandyna E2 | HMDB0001036 | Lipid | W górę ↑ | < 0.001 | 2.05 |
| Kwas litocholowy | HMDB0000752 | Kwas żółciowy | W dół ↓ | < 0.001 | 2.45 |
| Kwas masłowy | HMDB0000039 | krótkołańcuchowe kwasy tłuszczowe | W dół ↓ | < 0.001 | 2.21 |
| Kwas bursztynowy | HMDB0000254 | Kwas organiczny | W dół ↓ | 0.021 | 1.65 |
| L-tryptofan | HMDB0000929 | Aminokwas | W dół ↓ | 0.015 | 1.58 |
| Uwaga: VIP = Variable Importance in Projection (zmienna o znaczeniu istotnym dla projekcji) > 1 uznano za istotne). | | |
Tabela 3: Reprezentatywne różnicowe metabolity surowicy. Reprezentatywne adnotowane metabolity związane z PCOS i regulowane przez ekstrakt ZeLan, obejmujące identyfikatory HMDB, klasę metabolitu, trend regulacji, wartość P oraz wynik VIP. Pełna lista 52 adnotowanych cech metabolitów znajduje się w Tabeli uzupełniającej 1.
Tabela uzupełniająca 1: Kompletny adnotowany zestaw danych metabolicznych. Pełna lista 52 adnotowanych cech metabolitów zachowanych do raportowania metabolomicznego, zawierająca nazwę metabolitu, identyfikator HMDB, klasę biochemiczną, krotność zmiany (fold change), wartość P oraz kierunek regulacji.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 1: Dane źródłowe wspierające analizy eksperymentalne, mikrobiomowe i metabolomiczne. Ten arkusz zawiera informacje o próbkach, pomiary fenotypów klinicznych, wskaźniki różnorodności sekwencjonowania 16S rRNA, względną obfitość drobnoustrojów na poziomie typu, analizę różnicową taksonów LEfSe, różnicową analizę metabolomiczną, reprezentatywne metabolity, analizę wzbogacenia ścieżek KEGG, dane macierzy korelacji oraz podsumowanie weryfikacji danych wykorzystane do wygenerowania rycin i tabel przedstawionych w manuskrypcie.Kliknij tutaj, aby pobrać ten plik.