I dettagli di tutti gli animali, diete, reagenti, kit, strumenti, software e fornitori di servizi utilizzati in questo protocollo sono forniti nella Tabella dei Materiali. Tutti gli studi sugli animali sono stati condotti seguendo le Linee Guida del National Institutes of Health (NIH) per la Cura e l'Uso degli Animali da Laboratorio e hanno ricevuto l'approvazione dal Comitato per gli Animali Sperimentali dello Shandong Provincial Hospital (approvazione etica n. 2022–007, approvata il 24 gennaio 2022).
ATTENZIONE: Il sangue animale, i campioni fecali e i tessuti devono essere manipolati come materiali potenzialmente biopericolosi. Tutte le procedure che coinvolgono animali vivi, campioni biologici, oggetti taglienti, solventi organici e reagenti per analisi chimiche devono essere eseguite in conformità con le normative istituzionali sulla biosicurezza, la cura degli animali e la sicurezza chimica. Le carcasse animali, i tessuti, i materiali contaminati dal sangue e i campioni fecali devono essere raccolti in contenitori designati per rifiuti biologici e smaltiti tramite la struttura istituzionale per animali o l'ufficio biosicurezza. Aghi, aghi da gallina, tubi capillari e altri oggetti taglienti devono essere scartati immediatamente dopo l'uso nei contenitori approvati per taglienti. I rifiuti organici a base di solventi, inclusi metanolo, acetonitrile e soluzioni contenenti acido formico, dovrebbero essere raccolti in contenitori di rifiuti chimici etichettati e smaltiti tramite il programma istituzionale di rifiuti chimici pericolosi. I rifiuti provenienti da sagaggi biochimici e di stress ossidativo devono essere raccolti e smaltiti in conformità con le schede dati di sicurezza del produttore e le regole istituzionali sulla sicurezza chimica.
Animali sperimentali
I maschi di topi C57BL/6J di sei settimane (n = 4 per gruppo) sono stati mantenuti in condizioni specifiche prive di patogeni (SPF) con un ciclo chiaro/scuro di 12 ore a 22 ± 2 °C e umidità relativa del 50–60%. Gli animali avevano accesso ad libitum a cibo e acqua, con un massimo di cinque topi per gabbia. I topi di controllo sono stati nutriti con una dieta standard di chow composta dal 70% di carboidrati, 20% di proteine e 10% di grassi (kcal). Per stabilire il modello murino MASLD, i topi sono stati nutriti con una dieta ricca di grassi per 8 settimane contenente il 20% di carboidrati, il 20% di proteine e il 60% di grassi (kcal). Tutte le diete venivano conservate a 4 °C e rifornite due volte a settimana. I topi sono stati assegnati casualmente a gruppi sperimentali utilizzando un generatore di numeri casuali, e la dimensione del campione (n = 4 per gruppo) è stata determinata sulla base di studi precedenti con modelli murini simili indotti da diete ricche di grassi e su esperimenti preliminari. Tutti i trattamenti sono stati somministrati alla stessa ora ogni giorno per minimizzare la variazione circadiana.
L'olio di microalghe veniva estratto e purificato dall'Accademia di Scienze Agrarie dello Shandong, con una purezza di circa il 99%. L'olio di microalghe veniva conservato a −30 °C protetto dalla luce. Prima dell'amministrazione, l'olio veniva preparato con un emulsionante e miscelato accuratamente per garantire l'omogeneità. In particolare, l'olio di microalghe ricco di DHA è stato emulsionato in sterile carbossimetilcellulosa sodica contenente lo 0,5% di Tween-80 nella soluzione salina normale. L'olio di microalghe e il veicolo emulsionante sono stati miscelati in un rapporto di 1:9 (v/v) e sottoposti a vortici per 2 minuti immediatamente prima della gara per ottenere una sospensione uniforme. I topi del gruppo trattato con olio di microalghe ricevevano la preparazione una volta al giorno tramite gavage orale usando un ago sterile per l'alimentazione a una dose di 10 μL/g di peso corporeo. Il volume del gavage veniva regolato quotidianamente in base al peso corporeo per garantire la coerenza della dose. I topi del gruppo di controllo hanno ricevuto un volume equivalente di soluzione salina normale sterile nelle stesse condizioni. La soluzione salina utilizzata era iniezione di cloruro di sodio allo 0,9%.
Il trattamento probiotico consisteva in una formulazione composita contenente Lactobacillus plantarum (DY-1), Lactobacillus acidophilus (KDB-03), Lactobacillus casei (KDB-LC) e Bacillus coagulans (GIM 1.645). I quattro ceppi sono stati miscelati con un rapporto cellulare vitali uguale di 1:1:1:1. Ogni ceppo ha contribuito con 2,5 × 107 CFU/mL alla formulazione finale, ottenendo una concentrazione batterica totale di 1 × 108 CFU/mL. La miscela probiotica è stata sintetizzata dall'Accademia di Scienze Agrarie dello Shandong e non è disponibile commercialmente. Il trattamento probiotico consisteva in un brodo composito di fermentazione probiotico contenente Lactobacillus plantarum DY-1, Lactobacillus acidophilus KDB-03, Lactobacillus casei KDB-LC e Bacillus coagulans GIM 1.645. La concentrazione batterica della formulazione era di 1 × 108 CFU/mL. Il brodo probiotico di fermentazione veniva conservato a 4 °C prima dell'uso. I topi del gruppo trattato con probiotici ricevevano la formulazione probiotica una volta al giorno tramite gavage orale a una dose di 0,01 mL/g di peso corporeo durante tutto il periodo sperimentale.
I singoli ceppi sono stati coltivati in condizioni appropriate e raccolti nella fase di crescita a metà troncaggio (OD₆₀₀ ≈ 0,6), quindi combinati al rapporto indicato. Le colture batteriche sono state centrifugate per rimuovere il surnatante e i pellet sono stati risospesi in soluzione salina sterile tamponata con fosfato (PBS). La sospensione batterica appena preparata veniva conservata in ghiaccio prima della somministrazione per preservare la vitalità. La concentrazione batterica finale è stata regolata a 1 × 108 CFU/mL prima della somministrazione. I topi del gruppo probiotico hanno ricevuto la formulazione tramite gavage orale a una dose di 0,01 mL/g di peso corporeo una volta al giorno durante il periodo sperimentale. Non è stato applicato alcun pretrattamento antibiotico prima della somministrazione dei probiotici. Per ridurre il bias sperimentale, la valutazione degli esiti è stata effettuata da ricercatori ciechi dove applicabile. La MASLD è stata confermata tramite analisi biochimiche del sangue.
Analisi del sangue
Al termine del periodo sperimentale, i topi sono stati digiunati durante la notte (12 ore) con libero accesso all'acqua prima della misurazione della glicemia e della raccolta terminale del sangue. I topi sono stati profondamente anestetizzati con isoflurano utilizzando un vaporizzatore anestesitico fino a quando non è stata confermata la perdita del riflesso pedalale. Campioni di sangue terminale sono stati raccolti tramite sanguinamento retro-orbitale. Dopo il prelievo di sangue, i topi venivano soppressi per lussazione cervicale sotto anestesia profonda, e la morte veniva confermata con la cessazione della respirazione e del battito cardiaco. Il sangue veniva lasciato coagulare a temperatura ambiente e poi centrifugato a 3.000 × g per 10 minuti per ottenere il siero. I campioni di siero sono stati aliquotati e conservati a −80 °C fino all'analisi. Il glucosio sierico, il colesterolo totale (TC), i trigliceridi (TG), il colesterolo lipoproteico a bassa densità (LDL-c), il colesterolo lipoproteico ad alta densità (HDL-c), l'aspartato aminotransferasi (AST) e l'alanina aminotransferasi (ALT) sono stati misurati utilizzando un analizzatore biochimico automatico, seguendo le istruzioni del produttore. Questi parametri sono stati quantificati utilizzando saggi colorimetrici enzimatici con kit di reagenti standard integrati nel sistema analizzatrice. Tutti i parametri biochimici sono stati normalizzati al volume siero ed espressi nelle unità specificate dal produttore.
Saghi di stress ossidativo
I biomarcatori dello stress ossidativo sono stati misurati utilizzando kit, inclusi quelli per malondialdeide, superossido dismutasi e catalasi. I test venivano condotti seguendo le istruzioni del produttore. Per un breve tratto, i tessuti epatici furono omogeneizzati in tampone gelato e centrifugati per ottenere soprannatanti da analisi. I livelli di malondialdeide (MDA) sono stati determinati utilizzando un test di sostanze reattive all'acido tiobarbiturico (TBARS); l'attività della dismutasi del superossido (SOD) è stata misurata dalla sua capacità di inibire reazioni mediate dal superossido; e l'attività della catalasi (CAT) è stata quantificata monitorando il tasso di decomposizione del perossido di idrogeno. Le misurazioni venivano effettuate utilizzando un lettore di microplacche e i valori calcolati utilizzando curve standard o le formule fornite nei kit. I marcatori di stress ossidativo sono stati normalizzati al peso del tessuto epatico ed espressi secondo le specifiche del produttore.
Sequenziamento del rRNA 16S e analisi del microbioma
I campioni fecali venivano raccolti al termine del periodo sperimentale direttamente da singoli topi in condizioni sterili, immediatamente congelati in azoto liquido e conservati a −80 °C fino all'estrazione del DNA. Il DNA microbico è stato estratto utilizzando un kit commerciale secondo il protocollo del produttore. La regione V3–V4 del gene batterico 16S rRNA è stata amplificata utilizzando primer 341F/806R. L'amplificazione PCR è stata eseguita in un sistema di reazione da 25 μL contenente DNA modello, 2 mix master PCR e 0,2 μM di ciascun primer. Il programma di amplificazione consisteva in una denaturazione iniziale a 95 °C per 3 minuti, seguita da 25 cicli a 95 °C per 30 secondi, 55 °C per 30 secondi e 72 °C per 45 secondi, con un'estensione finale a 72 °C per 10 minuti. Le librerie Amplicon venivano costruite e sequenziate su una piattaforma Illumina MiSeq utilizzando una strategia a 2 × 250 bp a pari pari da un fornitore commerciale di servizi di sequenziamento. La qualità della libreria veniva valutata prima del sequenziamento e le librerie aggregate venivano caricate alla densità di cluster raccomandata dal produttore. Le letture di bassa qualità con basi ambigue, lunghezza insufficiente o un punteggio medio di qualità inferiore a Q20 sono state rimosse prima dell'analisi a valle.
Le letture di sequenziamento grezzo sono state processate utilizzando QIIME 2 versione 2019.420. In breve, le letture a pari sono state demultiplexate, filtrate per qualità, denoise ridotte, unite e filtrate per chimera usando DADA2 con parametri predefiniti, salvo diversa specificazione. Le varianti di sequenza amplicon non singole (ASV) sono state mantenute per l'analisi a valle. L'assegnazione tassonomica è stata eseguita utilizzando la versione 13_8 del database Greengenes implementata in QIIME221. Date le limitazioni intrinseche della risoluzione del sequenziamento genico del rRNA 16S, in particolare nell'utilizzo del database Greengenes, le annotazioni tassonomiche sono state interpretate esclusivamente a livello di genere. La classificazione a livello di specie non era considerata affidabile e quindi non veniva utilizzata in analisi a valle, confronti statistici o interpretazioni biologiche. Tutte le caratteristiche microbiche sono state compresse a livello di genere prima della visualizzazione e dell'analisi.
Sono state utilizzate metriche di alpha-diversità (Chao1, specie osservate, Shannon, Simpson, PD di Faith e uniforme di Pielou). Confronti statistici della diversità alfa tra i gruppi sono stati effettuati utilizzando test non parametrici appropriati (ad esempio, test di Kruskal–Wallis). La distanza di Jaccard e la dissimilanza di Bray-Curtis sono state utilizzate come metriche di diversità beta. Sono stati utilizzati grafici di Analisi Principale delle Coordinate (PCoA), scalatura multidimensionale non metrica (NMDS) e analisi discriminante lineare della dimensione dell'effetto (LEfSe) per illustrare i cambiamenti nelle comunità microbiome, con una soglia di punteggio LDA di 2,0. L'analisi multivariata permutazionale della varianza (PERMANOVA) è stata utilizzata per valutare la significatività statistica nella diversità beta. La visualizzazione dei dati e le analisi statistiche sono state effettuate utilizzando la piattaforma GenesCloud (https://www.genescloud.cn).
I profili funzionali microbici sono stati dedotti dai dati del sequenziamento del 16S rRNA piuttosto che misurati direttamente tramite sequenziamento metagenomico o trascrittomico a pompa. Sulla base della tabella di abbondanza dei percorsi normalizzati o dei gruppi funzionali, le unità funzionali previste sono state mappate a database comunemente utilizzati, inclusi KEGG, MetaCyc e COG. L'abbondanza dei percorsi KEGG è stata riassunta secondo le sue categorie gerarchiche funzionali, con le classificazioni di percorso di livello 2 utilizzate per il confronto a valle. L'abbondanza media di ciascuna categoria di percorso è stata calcolata utilizzando il software R.
Analisi metabolomica non mirata di campioni di fegato
I tessuti epatici venivano rapidamente estratti immediatamente dopo la conferma della morte, immediatamente congelati in azoto liquido e conservati a −80 °C prima dell'analisi. Circa 50 mg di tessuto epatico congelato sono stati omogeneizzati in metanolo:acqua gelato (4:1, v/v) integrato con standard interni. Il rapporto tessuto-solvente di estrazione era 1:10 (v/p). I campioni sono stati vortizzati per 1 minuto, sonicati sul ghiaccio per 10 minuti, e poi incubati a −20 °C per 30 minuti per migliorare la precipitazione proteica. Gli omogeneati furono accuratamente miscelati, sottoposti a sonicazione sul ghiaccio e centrifugati a 12.000 × g per 15 minuti a 4 °C. I soprannatanti sono stati accuratamente raccolti, evaporati sotto un lieve flusso di azoto e ricostituiti al 50% di metanolo per analisi successive.
I campioni di controllo qualità (QC) sono stati generati raggruppando volumi uguali da ciascun campione e sono stati analizzati a intermittenza durante la prova per valutare la stabilità analitica e la riproducibilità. L'analisi metabolomica non mirata veniva eseguita utilizzando un sistema di cromatografia liquida ad altissime prestazioni abbinato a spettrometria di massa ad alta risoluzione, gestito da un fornitore di servizi commerciale. I metaboliti sono stati separati su una colonna C18 in fase inversa utilizzando un sistema solvente binario composto da acqua con acido formico allo 0,1% e acetonitrile con acido formico allo 0,1%. Il volume di iniezione era di 2 μL, la temperatura della colonna era mantenuta a 40 °C e la portata era impostata a 0,30 mL/min. È stato utilizzato un programma di eluzione a gradiente per ottenere una vasta separazione dei metaboliti. La rilevazione spettrometrica di massa è stata condotta sia in modalità di ionizzazione elettrospray positive che negative su un intervallo m/z di 70–1.000 per ottenere una copertura metabolitica ampia. La tensione di spruzzatura ionica era impostata a 3,5 kV in modalità positiva e a −2,5 kV in modalità negativa. La temperatura capillare è stata mantenuta a 320 °C e i dati sono stati acquisiti in modalità MS a scansione completa con acquisizione MS/MS dipendente dai dati per l'annotazione dei metaboliti. I dati grezzi venivano elaborati utilizzando software dedicato alla metabolomica (MassLynx v4.1) per il rilevamento, l'allineamento e la normalizzazione delle caratteristiche. L'identificazione dei metaboliti è stata effettuata sulla base di massa, tempo di ritenzione e modelli di frammentazione MS/MS accurati, confrontandola con database pubblici, inclusi il Human Metabolome Database (HMDB) e la Kyoto Encyclopedia of Genes and Genomes (KEGG) (v2023), dove applicabile.
Analisi multivariata e di percorsi
I dati di metabolomica processati sono stati sottoposti ad analisi statistiche multivariate, inclusa l'analisi delle componenti principali (PCA) e l'analisi discriminante dei minimi quadrati parziali (PLS-DA), per visualizzare le differenze metaboliche tra i gruppi sperimentali. I metaboliti differenziali sono stati identificati in base a una combinazione di punteggi di importanza variabile nella proiezione (VIP) (>1,0) e significatività statistica (p < 0,05). I metaboliti identificati sono stati ulteriormente mappati nelle vie metaboliche utilizzando l'analisi di arricchimento delle vie KEGG per chiarire i processi biologici associati all'integrazione di olio di microalghe ricco di DHA.
Analisi statistica
I dati sono stati presentati come media ± SEM e analizzati utilizzando un'analisi unidirezionale della varianza (ANOVA) seguita dal test di confronto multiplo di Tukey. La significatività statistica è stata definita come p < 0,05. L'analisi di arricchimento dei metaboliti epatici è stata effettuata utilizzando R (versione 4.4.1). Le associazioni a coppie tra generi microbici e metaboliti epatici sono state analizzate utilizzando correlazioni di rango di Spearman. Solo correlazioni significative (tasso di false scoperte, FDR < 0,05) con |ρ| ≥ 0,6 sono stati inclusi nelle analisi integrative della rete. L'elaborazione dei dati per microbioma e metabolomica è stata effettuata come descritto sopra, inclusa l'uso della piattaforma GenesCloud e le procedure di controllo qualità.