L'esercizio su tapis roulant è stato associato al miglioramento di anomalie metaboliche nei topi diabetici
Rispetto al gruppo di controllo, i topi diabetici hanno mostrato valori significativamente più elevati di glicemia a digiuno (18,23 ± 0,61 contro 6,03 ± 0,21 mmol/L, P < 0,001), insulina sierica (16,30 ± 1,04 contro 8,23 ± 0,38 mIU/L, P < 0,001), HOMA-IR (13,19 ± 1,18 contro 2,21 ± 0,17, P < 0,001), colesterolo totale (4,62 ± 0,23 contro 2,01 ± 0,09 mmol/L, P < 0,001) e trigliceridi (2,12 ± 0,15 contro 0,80 ± 0,08 mmol/L, P < 0,001). Dopo l'intervento di esercizio su tapis roulant della durata di 8 settimane, questi parametri si sono ridotti rispettivamente a 12,03 ± 0,63 mmol/L, 10,61 ± 0,71 mIU/L, 5,69 ± 0,67, 3,09 ± 0,24 mmol/L e 1,20 ± 0,14 mmol/L, rispetto ai topi diabetici (P < 0,01 o P < 0,001); tuttavia, la maggior parte dei valori è rimasta più alta rispetto a quelli del gruppo di controllo (Figura 1A–E). L'analisi della varianza ad un criterio di classificazione (one-way ANOVA) ha evidenziato effetti significativi del gruppo per tutti i parametri (tutti F(2,27) > 25,0, P < 0,001).
L'allenamento fisico è stato associato a un'espressione più elevata di mRNA di Occludina nel colon e a un'espressione genica più bassa correlata all'inflammasoma cardiaco
Nei topi diabetici, l'espressione dell'mRNA di Occludin nel colon era circa il 62% inferiore rispetto al gruppo Controllo (0,38 ± 0,04 vs. 1,00 ± 0,06, P < 0,001). Contemporaneamente, l'espressione cardiaca degli mRNA di Nlrp3, Il1b e Tnfa aumentava rispettivamente di 5,6, 6,5 e 5,1 volte i valori del gruppo Controllo (tutti P < 0,001). Dopo l'intervento di esercizio fisico, l'espressione dell'mRNA di Occludin aumentava a 0,79 ± 0,05 (P < 0,001 rispetto ai topi diabetici), mentre l'espressione degli mRNA di Nlrp3, Il1b e Tnfa diminuiva rispettivamente a 2,30 ± 0,21, 2,66 ± 0,25 e 2,14 ± 0,17 (tutti P < 0,001 rispetto ai topi diabetici; Figura 2A–D).
L'esercizio fisico è stato associato a una minore espressione genica legata alla fibrosi e all'apoptosi miocardica
Il gruppo diabetico ha mostrato un'espressione cardiaca aumentata dei geni associati alla fibrosi Col1a1 (4,46 ± 0,29 volte, P < 0,001) e Tgfb1 (4,17 ± 0,24 volte, P < 0,001), insieme a un rapporto maggiore di mRNA Bax/Bcl2 (3,06 ± 0,19, P < 0,001) rispetto al gruppo Controllo. Dopo l'allenamento fisico, l'espressione di mRNA di Col1a1 e Tgfb1 è diminuita rispettivamente a 2,02 ± 0,16 e 1,92 ± 0,16 (entrambi P < 0,001 rispetto ai topi diabetici), e il rapporto di mRNA Bax/Bcl2 si è ridotto a 1,48 ± 0,10 (P < 0,001 rispetto ai topi diabetici; Figura 3A–C).
L'esercizio fisico è stato associato a cambiamenti in alcuni taxa microbici intestinali selezionati nei topi diabetici
Rispetto al gruppo di controllo, i topi diabetici hanno mostrato abbondanze relative significativamente più basse di quattro taxa batterici selezionati: 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) e Bifidobacterium (0,56 ± 0,06 vs. 1,00 ± 0,06; tutti P < 0,001). Al contrario, le abbondanze relative di E. coli (3,84 ± 0,33 vs. 0,99 ± 0,08) e Desulfovibrio (3,07 ± 0,26 vs. 0,99 ± 0,07) erano significativamente più elevate (entrambe P < 0,001). L'intervento con esercizio fisico è stato associato a un aumento significativo delle abbondanze relative di F. prausnitzii, A. muciniphila, Lactobacillus e Bifidobacterium (P < 0,01 o P < 0,001), con A. muciniphila che ha raggiunto un valore di 0,84 ± 0,07. Contestualmente, le abbondanze relative di E. coli e Desulfovibrio sono diminuite rispettivamente a 1,68 ± 0,13 e 1,52 ± 0,08 (entrambe P < 0,001 rispetto ai topi diabetici; Figura 4A–F).
Alterazioni microbiche selezionate sono state correlate all'espressione genica cardiaca e colonica
È stata eseguita un'analisi della correlazione di Spearman per valutare le associazioni tra le abbondanze relative di determinati taxa microbici e i marcatori di espressione genica nei tessuti cardiaci e colici (Figura 5A). L'abbondanza relativa di A. muciniphila è risultata negativamente correlata con l'mRNA di Nlrp3 cardiaco (R = −0,81, P = 0,0046) e con l'mRNA di Tnfa (R = −0,80, P = 0,0054). L'abbondanza relativa di E. coli è risultata positivamente correlata con l'espressione mRNA di Occludin nel colon (R = 0,66, P = 0,036; Figura 5B–D). Questi risultati hanno evidenziato associazioni tra determinati taxa microbici intestinali e marcatori trascrizionali cardiaci e colici, senza tuttavia stabilire relazioni causali.
DISPONIBILITÀ DEI DATI:
I set di dati generati e/o analizzati durante lo studio in corso sono disponibili presso l'autore corrispondente su richiesta ragionevole.

Figura 1: Effetti dell'esercizio su tapis roulant sui parametri metabolici nei topi diabetici. (A) Glicemia a digiuno. (B) Insulina sierica. (C) HOMA-IR. (D) Colesterolo totale. (E) Trigliceridi nei gruppi Controllo, DM e DM+Ex. I dati sono espressi come media ± DS (n = 10 per gruppo), con i singoli punti dati sovrapposti a ciascuna barra. Il confronto tra i gruppi è stato effettuato mediante ANOVA a un fattore seguita dal test post hoc di Tukey. *P < 0,05, **P < 0,01, ***P < 0,001; ns = non significativo. Abbreviazioni: DM = modello di diabete; Ex = esercizio; HOMA-IR = valutazione del modello omeostatico della resistenza all'insulina; DS = deviazione standard. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 2: Effetti dell'esercizio su tapis roulant sull'espressione genica legata alla barriera intestinale e all'inflammasoma cardiaco nei topi diabetici. (A) Espressione relativa dell'mRNA di Occludin nel colon. (B) Espressione relativa dell'mRNA di Nlrp3 nel cuore. (C) Espressione relativa dell'mRNA di Il1b nel cuore. (D) Espressione relativa dell'mRNA di Tnfa nel cuore. L'espressione genica è stata normalizzata rispetto a Gapdh ed espressa come cambio di fold rispetto al gruppo Controllo. I dati sono riportati come media ± DS (n = 10 per gruppo). *P < 0,05, **P < 0,01, ***P < 0,001. Abbreviazioni: IL-1β = interleuchina-1β; TNF-α = fattore di necrosi tumorale-α; DS = deviazione standard. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3: Effetti dell'esercizio su tapis roulant sull'espressione genica di fibrosi miocardica e apoptosi nei topi diabetici. (A) Espressione relativa dell'mRNA cardiaco di Col1a1. (B) Espressione relativa dell'mRNA cardiaco di Tgfb1. (C) Rapporto di mRNA Bax/Bcl2. L'espressione genica è stata normalizzata rispetto a Gapdh ed espressa come variazione di fold rispetto al gruppo Controllo. Il rapporto Bax/Bcl2 è stato calcolato a partire dai valori di espressione relativa dell'mRNA. I dati sono presentati come media ± DS (n = 10 per gruppo), con i singoli punti dati sovrapposti a ciascuna barra. *P < 0,05, **P < 0,01, ***P < 0,001. Abbreviazioni: Col1a1 = collagena di tipo I alfa 1; Tgfb1 = fattore di crescita trasformante-β1; Bax = proteina associata a BCL2; Bcl2 = linfoma delle cellule B 2; DS = deviazione standard. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 4: Effetti dell'esercizio su tapis roulant sulle abbondanze relative di specifici taxa microbici intestinali nei topi diabetici. (A) Faecalibacterium prausnitzii. (B) Akkermansia muciniphila. (C) Lactobacillus. (D) Bifidobacterium. (E) Escherichia coli. (F) Desulfovibrio. Le abbondanze relative sono state determinate mediante qPCR utilizzando target specifici per il gene dell'rRNA 16S relativi all'rRNA 16S batterico totale e normalizzati al gruppo Controllo. I dati sono espressi come media ± DS (n = 10 per gruppo). *P < 0,05, **P < 0,01, ***P < 0,001. Abbreviazione: DS = deviazione standard. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 5: Analisi di correlazione tra specifici taxa microbici intestinali e parametri di espressione genica cardiaca e colonica. (A) Mappa termica che mostra i coefficienti di correlazione di Spearman tra specifici taxa batterici e i livelli di mRNA di Occludin nel colon, di Nlrp3, Il1b e Tnfa nel cuore, di Col1a1 e Tgfb1 mRNA, e del rapporto tra mRNA Bax/Bcl2. Le correlazioni sono state calcolate utilizzando i dati di tutti e tre i gruppi (n = 30) come correlazioni parziali di Spearman, correggendo per il gruppo sperimentale. La scala cromatica varia dal blu, che indica correlazioni negative, al rosso, che indica correlazioni positive. Gli asterischi all'interno della mappa termica indicano significatività dopo correzione per il tasso di falsa scoperta (FDR) (*P < 0,05, **P < 0,01, ***P < 0,001). (B) Correlazione tra l'abbondanza relativa di A. muciniphila e l'espressione di mRNA di Nlrp3 cardiaco nel gruppo DM+Ex. (C) Correlazione tra l'abbondanza relativa di E. coli e l'espressione di mRNA di Occludin colonico nel gruppo DM+Ex. (D) Correlazione tra l'abbondanza relativa di A. muciniphila e l'espressione di mRNA di Tnfa cardiaco nel gruppo DM+Ex. Ogni punto rappresenta un singolo topo. Sono riportate rette di regressione lineare e intervalli di confidenza al 95% per finalità illustrative. Abbreviazioni: DM = modello di diabete; Ex = esercizio fisico; FDR = false discovery rate. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| Gene/Microbo | Primer diretto | Primer inverso |
| Tnfa | CCCTCACACTCAGATCATCTTCT | GCTACGACGTGGGCTACAG |
| Il1b | GCAACTGTTCCTGAACTCAACT | ATCTTTTGGGGTCCGTCAACT |
| Tgfb1 | CTCCCGTGGCTTCTAGTGC | GCCTTAGTTTGGACAGGATCTG |
| Col1a1 | GCTCCTCTTAGGGGCCACT | CCACGTCTCACCATTGGGG |
| Bax | AGACAGGGGCCTTTTTGCTAC | AATTCGCCGGAGACACTCG |
| Bcl2 | GCTACCGTCGTTGACTTCGC | CCCCACCGAACTCAAAGAAGG |
| Nlrp3 | GTGAGCCCACACCACAGTTC | GACGTTCACCTCGCAGATGA |
| Gapdh | AGGTCGGTGTGAACGGATTTG | TGTAGACCATGTAGTTGAGGTCA |
| Batteri totali (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 |
Tabella 1: Sequenze dei primer utilizzate per la PCR quantitativa in tempo reale.