Dodici maschi specifici di Sprague-Dawley (SD) privi di patogeno specifico (SPF), di età compresa tra 6 e 7 settimane e con un peso di 220 g ± 10 g, sono stati acquistati dalla Beijing Vital River Laboratory Animal Technology Co., Ltd. (Numero di licenza animale: SCXK (Jing) 2025-0008). I ratti sono stati ospitati presso il Centro per Animali da Laboratorio dell'Ospedale Dongzhimen, Università di Medicina Cinese di Pechino, in conformità con lo standard nazionale "Animali da laboratorio - Requisiti per l'ambiente e le strutture abitative" (GB 14925-2010, emanato dal Comitato Tecnico Nazionale per la Standardizzazione degli Animali da Laboratorio della Cina). Le procedure sperimentali sono state approvate dal Comitato per il Benessere Animale e l'Etica dell'Ospedale Dongzhimen, Università di Medicina Cinese di Pechino (Numero di approvazione: 24-54). I reagenti e le attrezzature utilizzate sono elencati nella Tabella dei Materiali.
1. Creazione del modello e raccolta dei campioni
NOTA: Prima dell'esperimento, sottoporre tutti i ratti a un periodo di acclimatazione di 1 settimana per minimizzare gli effetti indotti dallo stress.
- Dissolvere il lipopolisaccaride di Escherichia coli (LPS) in soluzione salina allo 0,9% per ottenere una soluzione LPS da 0,4 mg/mL.
- Dividiamo casualmente 12 maschi di ratti Sprague-Dawley (SD) SPF in due gruppi: un gruppo di controllo (n = 6) e un gruppo modello ARDS (n = 6).
- Somministrare soluzione salina allo 0,9% tramite iniezione intratracheale a una dose di 0,5 mL/100g di peso corporeo al gruppo di controllo.
- Somministrare la soluzione LPS da 0,4 mg/mL tramite iniezione intratracheale allo stesso volume (0,5 mL/100g) al gruppo modello ARDS e, dopo 16 ore, anestetizzare i ratti usando isoflurano (seguendo i protocolli istituzionalmente approvati).
- Raccogli sangue tramite l'aorta addominale, poi centrifuga (4 °C, 1000 × g, 30 minuti) per ottenere il siero.
- Apri la cavità toracica ed estrai il tessuto polmonare.
- Eseguire lavaggio broncoalveolare tramite il moncone bronchiale utilizzando una siringa per ottenere il Liquido di Lavaggio Broncoalveolare (BALF)24,25.
- Preleva il lobo polmonare destro e sciacqualo con PBS per un successivo esame patologico.
- Prelevare il lobo polmonare sinistro e sciacquarlo con PBS per i successivi esperimenti di sequenziamento cellulare singola ebiologia molecolare 26,27.
2. Valutazione istopatologica del tessuto polmonare
- Immergi il tessuto polmonare destro prelevato dei ratti in paraformaldeide al 4% per fissarlo per 12 ore.
- Dopo la fissazione, disidratare il tessuto attraverso una serie gradata di etanolo (70%, 80%, 90%, 95%, 100% e 100%), permettendo 2 ore a ciascuna concentrazione.
- Tratta il tessuto disidratato con xilene finché non diventa trasparente.
- Infiltra il tessuto trasparente con paraffina fusa in un forno a 60 °C per 2 ore.
- Metti il blocco di tessuto infiltrato in uno stampo di incorporazione, versa la paraffina fresca e lascia raffreddare e solidificare in un blocco di paraffina.
- Dividere il blocco di paraffina in fette spesse 4 μm usando un microtomo, far galleggiare le sezioni su un bagno d'acqua a 45 °C per spargerle, e poi sollevarle con vetrini rivestiti di poli-L-lisina.
- Cuocere i vetri in un forno a 60 °C per garantire una solida adesione delle sezioni ai veglietti.
- Deparafinizza le sezioni tramite immersione nello xilene e reidratale attraverso un gradiente discendente di etanolo (100%, 95%, 90%, 80% e 70%).
- Risciacquare le sezioni con PBS, eseguire la colorazione di ematossilina ed eosina (HE) e montare le coperture usando resina neutra.
- Valutare la lesione polmonare su sezioni scansionate intere utilizzando un sistema di punteggio standardizzato basato sull'infiltrazione di neutrofili, infiammazione interstiziale, edema e congestione.
NOTA: Ogni parametro è stato valutato da 0 a 4 in base all'entità del coinvolgimento patologico: 0 (nessuno); 1 (≤25%); 2 (26-50%); 3 (51-75%); 4 (>75%)28,29.
3. Saggio immunoenzimatico legato agli enzimi
NOTA: I campioni possono essere mantenuti a -80 °C se l'analisi non può essere eseguita immediatamente dopo l'acquisizione.
- Sciogliere gli standard e lasciarli a temperatura ambiente per 15 minuti, poi eseguire la diluizione seriale secondo le istruzioni del produttore.
- Prepara la soluzione di lavoro degli anticorpi biotinilati, la soluzione di lavoro con coniugati enzimatici, e il tampone di lavaggio come indicato nel manuale.
- Aggiungere 100 μL di ciascun campione o 100 μL di standard diluiti serialmente ai pozzi rivestiti di anticorpi, e incubare a 37 °C per 90 minuti.
- Lava il piatto 4 volte con una rondella per piatti, aggiungi 100 μL di soluzione attiva per anticorpi biotinilati e incuba a 37 °C per 30 minuti.
- Lava la piastra 4 volte, aggiungi 100 μL di soluzione lavorativa coniugata enzimatica e incuba a 37 °C per 30 minuti.
- Aggiungere 100 μL di soluzione di substrato per pozzo, incubare protetto dalla luce a 37 °C per 10 minuti.
- Aggiungi 100 μL di soluzione di stop per pozzo, mescola delicatamente e misura la densità ottica (OD) a 450 nm usando un lettore di microplacche.
- Genera una curva standard tracciando i valori OD rispetto alle concentrazioni standard utilizzando il software progettato.
- Determinare le concentrazioni del campione basandosi sui valori di OD misurati e sulla curva standard.
4. Preparazione di sospensioni a singola per il sequenziamento
NOTA: Dopo l'ottenimento del tessuto polmonare, deve essere messo sul ghiaccio e processato rapidamente per mantenere la vitalità cellulare.
- Lava il tessuto polmonare 3 volte con PBS pre-raffreddato a 4 °C.
- Aggiungi 2 mL di PBS pre-raffreddato e tritali ulteriormente in frammenti più piccoli usando forbici sterili fino a farli dimensionare circa 0,5mm 3 .
- Usa una pipetta Pasteur per trasferire il tessuto tritato insieme al PBS in un tubo centrifuga da 15 mL.
- Aggiungere un ulteriore 1 mL di PBS pre-raffreddato a 4 °C al vaso di coltura, risciacquare e risospendere i frammenti di tessuto rimanenti e trasferire la sospensione nello stesso tubo.
- Regolare il volume a 6,5 mL con PBS pre-raffreddato a 4 °C, digerire i frammenti di tessuto in una soluzione enzimatica contenente Collagenasi I (100 U/mL) e Dispase II (1 U/mL), e incubare con lievi scuotiti a temperatura ambiente per 30 minuti.
- Quando i frammenti di tessuto diventano traslucidi e non rimangono pezzi solidi distinti, filtrare il tessuto digerito attraverso filtri cellulari da 70 μm e 40 μm, poi centrifugare (4° C, 500 × g, 5 minuti) e scartare il soprantantante.
- Sospendi immediatamente il pellet in 3 mL di tampone per la lisi dei globuli rossi (RBC), mescola delicatamente con pipetto e lisa sul ghiaccio per 10 minuti.
- Quando il soprantantante appare giallo pallido o quasi incolore, si aggiungono 7 mL di PBS per terminare la reazione, centrifuga (4 °C, 500 × g, 5 minuti) e si scarta il soprantantante.
- Aggiungere 3 mL di PBS per risospendere le cellule e determinare la concentrazione e la vitalità cellulare usando la colorazione con tripano blu.
NOTA: Se la vitalità cellulare è inferiore al 70%, esegui una procedura di rimozione delle cellule morte.
- Lavare le cellule due volte con un tampone di lavaggio contenente albumina sirica bovina (BSA) allo 0,04% preparata in DPBS.
- Sospendere le cellule per il conteggio e regolare la concentrazione a 1.000 cellule/μL per il carico e la rilevazione successiva.
5. Controllo qualità per i dati di sequenziamento a singola cella
NOTA: Le librerie di RNA a singola cellula sono state preparate da sospensioni cellulari utilizzando un metodo di costruzione di librerie a singola cellula basata su gocce. Le librerie venivano sequenziate su una piattaforma ad alto rendimento utilizzando la chimica del sequenziamento basata su nanoball di DNA. Circa 21.000 celle sono state caricate per libreria, puntando a un tasso di recupero delle celle del ~60%. Il sequenziamento generava un totale di 5,0 × 108 letture per libreria, con una profondità media di 30.000-40.000 letture per cella. I dati grezzi del sequenziamento sono stati processati utilizzando la corrispondente pipeline di analisi dell'RNA a singola cellula per produrre la matrice di espressione genica per ogni campione.
- Costruire la libreria Oligo tramite PCR e marcatura con codici a barre, garantendo concentrazione >10 ng/μL e un picco di 180 pb ±10 pb.
- Costruire la libreria di cDNA tramite frammentazione, riparazione dell'estremità, legatura degli adattatori e PCR, assicurando una concentrazione >10 ng/μL e un picco di 350–550 bp.
- Genera un oggetto Seurat per ogni campione dalla matrice di espressione genica utilizzando la funzione Read10X nel pacchetto Seurat (v4.3.0) nel software R (v4.2.0).
- Rimuovere cellule di bassa qualità filtrando quelle con proporzione genica mitocondriale >10% o proporzione gene dell'emoglobina >5%.
- Normalizzare i dati usando il metodo "LogNormalize" con un fattore di scala di 10.000.
- Identificare caratteristiche altamente variabili utilizzando il metodo "vst" e mantenere i primi 2.000 geni per l'analisi a valle.
- Eseguire la valutazione del ciclo cellulare utilizzando i geni marcatori integrati di Seurat e regredire gli effetti del ciclo cellulare (G2M. Score e S.Score) durante la scalatura dei dati.
- Rilevare e rimuovere potenziali doppietti usando DoubletFinder (v2.0.3) con pN = 0,25, pK = 0,09 e un tasso atteso di doppio di 0,06, escludendo questi dall'ulteriore elaborazione.
- Effetti batch corretti tra i campioni usando il pacchetto Harmony (v0.1.1) con theta = 2 e max.iter.harmony = 20.
- Aggiungi colonne informative per nomi di campioni individuali e dettagli di raggruppamento del gruppo modello ARDS e del gruppo di controllo ai metadati nell'oggetto Serat.
6. Riduzione della dimensionalità, clustering e annotazione dei tipi di cella
- Eseguire analisi a componente principale (PCA) su geni altamente variabili utilizzando la funzione RunPCA in Seurat.
- Determina il numero ottimale di componenti principali tramite la visualizzazione ElbowPlot.
- Identificare i principali cluster di celle applicando l'algoritmo di Leiden tramite la funzione FindClusters e la funzione FindNeighbors di Seurat con una risoluzione di 0,1.
- Effettuare la riduzione della dimensionalità UMAP utilizzando i componenti principali selezionati.
- Annota i tipi cellulari identificando geni marcatori con la funzione FindAllMarkers.
- Visualizza l'espressione genica dei marcatori tra i tipi cellulari usando la funzione DotPlot.
- Aggiungi una nuova colonna informativa dei risultati delle annotazioni delle celle ai metadati e visualizza l'annotazione del tipo di cella con Umap.
- Calcola la proporzione di ciascun gruppo sperimentale all'interno di ciascun tipo di cellula identificato usando funzioni dplyr per contare le celle e calcolare le percentuali.
- Genera grafici a barre impilate con ggplot2 per illustrare la composizione dei gruppi per tipo di cella, inclusi etichetti percentuali e confronti statistici tra gruppi.
7. Analisi delle sottopopolazioni
- Estrarre la sottopopolazione delle cellule target utilizzando la funzione sottoinsieme e creare un oggetto Seurat corrispondente per l'analisi a livello di sottopopolazione a valle.
- Rinormalizzare e ridimensionare i dati della sottopopolazione mentre si fanno regressioni ai geni mitocondriali, ai geni dell'emoglobina, ai conteggi UMI e ai punteggi del ciclo cellulare.
- Eseguire un'analisi PCA secondaria e applicare la correzione batch di Harmony specificamente per la sottopopolazione.
- Effettuare ri-clustering utilizzando funzioni FindNeighbors e FindClusters per identificare sottotipi distinti all'interno delle sottopopolazioni di cellule target.
- Generare nuovi embeddings UMAP e visualizzare i cluster di sottopopolazione con annotazioni di tipo cellulare seguendo lo stesso approccio descritto in precedenza.
8. Analisi differenziale dell'espressione genica
- Identificare geni differenzialmente espressi (DEG) tra il gruppo modello ARDS e il gruppo di controllo per ciascun sottocluster utilizzando il metodo statistico MAST.
- Incorporare la regressione delle variabili latenti per controllare gli effetti di profondità di sequenziamento durante i test di espressione differenziale.
- Filtrare i DEG basandosi su soglie di cambiamento logaritarim predefinite e valori p regolati usando la funzione FindMarkers.
- Genera grafici vulcanici con la funzione ggplot per ogni sottoammasso per visualizzare i modelli differenziali di espressione genica con l'etichettatura di significatività.
- Esportare i risultati del DEG e creare tabelle riassuntive complete per facilitare l'arricchimento funzionale e altre analisi a valle.
9. Analisi dell'arricchimento
- Mappare i simboli genici agli ID ENTREZ usando l'org. Rn.eg.db database di annotazioni per la compatibilità dell'analisi dei percorsi.
- Combinare i DEG di tutti i sottocluster e filtrare geni significativi in base alle soglie di variazione logaritarime e ai valori p aggiustati (test di Wilcoxon rank-sum, logfc.threshold = 0,25, p_val < 0,05).
- Eseguire analisi di arricchimento tramite Ontologia Genica (GO) utilizzando la funzione di arricchimento GO con categorie di ontologie di Processo Biologico, Funzione Molecolare, Componenti Cellulari e metodo di correzione di Benjamini-Hochberg.
- Condurre l'analisi di arricchimento delle vie della Kyoto Encyclopedia of Genes and Genomes (KEGG) utilizzando la funzione di arricchimento KEGG con il database degli organismi ratti e applicare soglie di significatività statistica.
- Genera dot plot per i risultati di arricchimento GO e KEGG utilizzando la funzione dotplot, mostrando i primi 15 termini significativamente arricchiti con rapporti genici e valori p corretti.
- Esporta i risultati dell'arricchimento come file CSV e salva grafici di visualizzazione per un'interpretazione completa dell'analisi dei percorsi.