Ćwiczenia na bieżni wiązały się z poprawą nieprawidłowości metabolicznych u myszy z cukrzycą
W porównaniu z grupą kontrolną, myszy z cukrzycą wykazywały istotnie wyższy poziom glukozy w krwi na czczo (18,23 ± 0,61 vs. 6,03 ± 0,21 mmol/L, P < 0,001), insuliny w surowicy (16,30 ± 1,04 vs. 8,23 ± 0,38 mIU/L, P < 0,001), wskaźnik HOMA-IR (13,19 ± 1,18 vs. 2,21 ± 0,17, P < 0,001), cholesterol całkowity (4,62 ± 0,23 vs. 2,01 ± 0,09 mmol/L, P < 0,001) oraz trójglicerydy (2,12 ± 0,15 vs. 0,80 ± 0,08 mmol/L, P < 0,001). Po 8-tygodniowej interwencji w postaci ćwiczeń na bieżni parametry te uległy obniżeniu odpowiednio do 12,03 ± 0,63 mmol/L, 10,61 ± 0,71 mIU/L, 5,69 ± 0,67, 3,09 ± 0,24 mmol/L oraz 1,20 ± 0,14 mmol/L w porównaniu z myszami z cukrzycą (P < 0,01 lub P < 0,001); jednak większość wartości pozostała wyższa niż w grupie kontrolnej (Rycina 1A–E). Jednoczynnikowa analiza wariancji ANOVA wykazała istotne efekty grupowe dla wszystkich parametrów (wszystkie F(2,27) > 25,0, P < 0,001).
Trening fizyczny wiązał się z wyższą ekspresją mRNA okludyny w jelicie grubym oraz niższą ekspresją genów związanych z inflamasomem w sercu
U myszy z cukrzycą ekspresja mRNA Occludin w okrężnicy była o około 62% niższa niż w grupie kontrolnej (0,38 ± 0,04 vs 1,00 ± 0,06, P < 0,001). Jednocześnie ekspresja mRNA Nlrp3, Il1b oraz Tnfa w sercu wzrosła odpowiednio 5,6-krotnie, 6,5-krotnie i 5,1-krotnie w stosunku do wartości kontrolnych (wszystkie P < 0,001). Po interwencji w postaci ćwiczeń fizycznych ekspresja mRNA Occludin wzrosła do 0,79 ± 0,05 (P < 0,001 vs myszy z cukrzycą), natomiast ekspresja mRNA Nlrp3, Il1b i Tnfa spadła odpowiednio do 2,30 ± 0,21, 2,66 ± 0,25 i 2,14 ± 0,17 (wszystkie P < 0,001 vs myszy z cukrzycą; Rycina 2A–D).
Ćwiczenia fizyczne wiązały się z niższą ekspresją genów związanych z włóknieniem i apoptozą mięśnia sercowego
W grupie z cukrzycą zaobserwowano zwiększoną ekspresję w sercu genów związanych z włóknieniem Col1a1 (4,46 ± 0,29-krotny wzrost, P < 0,001) oraz Tgfb1 (4,17 ± 0,24-krotny wzrost, P < 0,001), wraz z wyższym stosunkiem mRNA Bax/Bcl2 (3,06 ± 0,19, P < 0,001) w porównaniu z grupą kontrolną. Po treningu fizycznym ekspresja mRNA Col1a1 i Tgfb1 spadła odpowiednio do 2,02 ± 0,16 i 1,92 ± 0,16 (dla obu P < 0,001 vs. myszy z cukrzycą), a stosunek mRNA Bax/Bcl2 spadł do 1,48 ± 0,10 (P < 0,001 vs. myszy z cukrzycą; Rycyna 3A–C).
Ćwiczenia fizyczne wiązały się ze zmianami w wybranych taksonach mikrobioty jelitowej u myszy z cukrzycą
W porównaniu z grupą kontrolną, u myszy z cukrzycą stwierdzono znacznie niższą względną liczebność czterech wybranych taksonów bakteryjnych: F. prausnitzii (0,25 ± 0,03 vs. 1,01 ± 0,07), A. muciniphila (0,31 ± 0,04 vs. 1,00 ± 0,06), Lactobacillus (0,48 ± 0,05 vs. 0,99 ± 0,04) oraz Bifidobacterium (0,56 ± 0,06 vs. 1,00 ± 0,06; wszystkie P < 0,001). Z kolei względna liczebność E. coli (3,84 ± 0,33 vs. 0,99 ± 0,08) i Desulfovibrio (3,07 ± 0,26 vs. 0,99 ± 0,07) była znacznie wyższa (oba P < 0,001). Interwencja w postaci ćwiczeń fizycznych wiązała się ze znacznie wyższą względną liczebnością F. prausnitzii, A. muciniphila, Lactobacillus i Bifidobacterium (P < 0,01 lub P < 0,001), przy czym wartość dla A. muciniphila osiągnęła 0,84 ± 0,07. Jednocześnie względna liczebność E. coli i Desulfovibrio spadła odpowiednio do 1,68 ± 0,13 i 1,52 ± 0,08 (oba P < 0,001 w porównaniu z myszami z cukrzycą; Rycina 4A–F).
Wybrane zmiany mikrobiologiczne skorelowano z ekspresją genów w sercu i okrężnicy
Przeprowadzono analizę korelacji Spearmana w celu oceny powiązań między względną licznością wybranych taksonów mikrobowych a markerami ekspresji genów w tkankach serca i okrężnicy (Rysunek 5A). Względna liczność A. muciniphila była ujemnie skorelowana z mRNA Nlrp3 w sercu (R = −0.81, P = 0.0046) oraz mRNA Tnfa w sercu (R = −0.80, P = 0.0054). Względna liczność E. coli była dodatnio skorelowana z ekspresją mRNA Occludin w okrężnicy (R = 0.66, P = 0.036; Rysunek 5B–D). Wyniki te wykazały powiązania między wybranymi taksonami mikrobowymi jelit a markerami transkrypcyjnymi serca i okrężnicy, lecz nie ustanowiły związków przyczynowo-skutkowych.
DOSTĘPNOŚĆ DANYCH:
Zbiory danych wygenerowane i/lub analizowane w ramach niniejszego badania są dostępne u autora korespondencyjnego na uzasadnioną prośbę.

Rycina 1: Wpływ ćwiczeń na bieżni na parametry metaboliczne u myszy z cukrzycą. (A) Glikemia na czczo. (B) Insulina w surowicy. (C) HOMA-IR. (D) Cholesterol całkowity. (E) Triglicerydy w grupach Control, DM oraz DM+Ex. Dane przedstawiono jako średnia ± SD (n = 10 na grupę), z naniesionymi na każdy słupek indywidualnymi punktami danych. Porównania między grupami przeprowadzono za pomocą jednoczynnikowej analizy wariancji ANOVA, a następnie testu post hoc Tukeya. *P < 0.05, **P < 0.01, ***P < 0.001; ns = brak istotności statystycznej. Skróty: DM = model cukrzycy; Ex = ćwiczenia; HOMA-IR = homeostatyczna ocena insulinooporności; SD = odchylenie standardowe. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 2: Wpływ ćwiczeń na bieżni na ekspresję genów związanych z barierą jelitową oraz inflamasomem sercowym u myszy z cukrzycą. (A) Względna ekspresja mRNA Occludin w okrężnicy. (B) Względna ekspresja mRNA Nlrp3 w sercu. (C) Względna ekspresja mRNA Il1b w sercu. (D) Względna ekspresja mRNA Tnfa w sercu. Ekspresję genów znormalizowano względem Gapdh i przedstawiono jako zmianę krotności (fold change) w stosunku do grupy kontrolnej. Dane przedstawiono jako średnia ± SD (n = 10 na grupę). *P < 0,05, **P < 0,01, ***P < 0,001. Skróty: IL-1β = interleukina-1β; TNF-α = czynnik martwicy nowotworów-α; SD = odchylenie standardowe. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 3: Wpływ ćwiczeń na bieżni na ekspresję genów związanych z włóknieniem mięśnia sercowego i apoptozą u myszy z cukrzycą. (A) Względna ekspresja mRNA Col1a1 w sercu. (B) Względna ekspresja mRNA Tgfb1 w sercu. (C) Stosunek mRNA Bax/Bcl2. Ekspresję genów znormalizowano względem Gapdh i przedstawiono jako krotność zmiany w stosunku do grupy kontrolnej. Stosunek Bax/Bcl2 obliczono na podstawie wartości względnej ekspresji mRNA. Dane przedstawiono jako średnia ± SD (n = 10 na grupę), z naniesionymi na każdy słupek indywidualnymi punktami danych. *P < 0.05, **P < 0.01, ***P < 0.001. Skróty: Col1a1 = kolagen typu I alfa 1; Tgfb1 = transformujący czynnik wzrostu-β1; Bax = białko X związane z BCL2; Bcl2 = chłoniak komórek B 2; SD = odchylenie standardowe. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 4: Wpływ ćwiczeń na bieżni na względną liczebność wybranych taksonów mikrobioty jelitowej u myszy z cukrzycą. (A) Faecalibacterium prausnitzii. (B) Akkermansia muciniphila. (C) Lactobacillus. (D) Bifidobacterium. (E) Escherichia coli. (F) Desulfovibrio. Względną liczebność określono metodą qPCR, wykorzystując specyficzne dla taksonów sekwencje docelowe 16S rRNA w stosunku do całkowitego bakteryjnego 16S rRNA i normalizując wyniki względem grupy kontrolnej. Dane przedstawiono jako średnia ± SD (n = 10 na grupę). *P < 0.05, **P < 0.01, ***P < 0.001. Skrót: SD = odchylenie standardowe. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Analiza korelacji między wybranymi taksonami mikrobioty jelitowej a parametrami ekspresji genów w sercu i okrężnicy. (A) Mapa ciepła przedstawiająca współczynniki korelacji Spearmana między wybranymi taksonami bakteryjnymi a mRNA Occludin w okrężnicy, mRNA Nlrp3, Il1b i Tnfa w sercu, mRNA Col1a1 i Tgfb1 oraz stosunkiem mRNA Bax/Bcl2. Korelacje obliczono na podstawie danych z trzech grup (n = 30) jako cząstkowe korelacje Spearmana z kontrolą grupy eksperymentalnej. Skala kolorystyczna obejmuje zakres od niebieskiego dla korelacji ujemnych do czerwonego dla korelacji dodatnich. Gwiazdki na mapie ciepła oznaczają istotność po korekcie FDR (*P < 0.05, **P < 0.01, ***P < 0.001). (B) Korelacja między względną licznością A. muciniphila a ekspresją mRNA Nlrp3 w sercu w grupie DM+Ex. (C) Korelacja między względną licznością E. coli a ekspresją mRNA Occludin w okrężnicy w grupie DM+Ex. (D) Korelacja między względną licznością A. muciniphila a ekspresją mRNA Tnfa w sercu w grupie DM+Ex. Każdy punkt reprezentuje pojedynczą mysz. Linie regresji liniowej i 95% przedziały ufności zostały przedstawione w celu wizualizacji. Skróty: DM = model cukrzycy; Ex = ćwiczenia; FDR = stopień fałszywych odkryć. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Gen/mikrob | Starter forward | Starter reverse |
| Tnfa | CCCTCACACTCAGATCATCTTCT | GCTACGACGTGGGCTACAG |
| Il1b | GCAACTGTTCCTGAACTCAACT | ATCTTTTGGGGTCCGTCAACT |
| Tgfb1 | CTCCCGTGGCTTCTAGTGC | GCCTTAGTTTGGACAGGATCTG |
| Col1a1 | GCTCCTCTTAGGGGCCACT | CCACGTCTCACCATTGGGG |
| Bax | AGACAGGGGCCTTTTTGCTAC | AATTCGCCGGAGACACTCG |
| Bcl2 | GCTACCGTCGTGACTTCGC | CCCCACCGAACTCAAAGAAGG |
| Nlrp3 | GTGAGCCCACACCACAGTTC | GACGTTCACCTCGCAGATGA |
| Gapdh | AGGTCGGTGTGAACGGATTTG | TGTAGACCATGTAGTTGAGGTCA |
| Całkowita liczba bakterii (16S rRNA) | ACTCCTACGGGAGGCAGCAG | ATTACCGCGGCTGCTGG |
| F. prausnitzii | GGAGGAAGAAGGTCTTCGG | AATTCCGCCTACCTCTGCACT |
| A. muciniphila | CAGCACGTGAAGGTGGGGAC | CCTTGCGGTTGGCTTCAGAT |
| Lactobacillus | AGCAGTAGGGAATCTTCCA | ATTYCACCGCTACACATG |
| Bifidobacterium | TCGCGTCCGGTGTGAAAG | CCACATCCAGCATCCAC |
| E. coli | AGGCCTTCGGGTTGTAAAGT | GTTAGCCGGTGCTTCTTCTG |
| Desulfovibrio | CCGTAGATATCTGGAGGAACATCAG | ACATCTAGCATCCATCGTTTACAGC |
Tabela 1: Sekwencje starterów użyte do ilościowej PCR w czasie rzeczywistym.