Dwanaście samców szczurów Sprague-Dawley (SD) wolnych od specyficznych patogenów (SPF), w wieku 6–7 tygodni i ważących 220 g ± 10 g, zostało zakupionych od Beijing Vital River Laboratory Animal Technology Co., Ltd. (numer licencji zwierzęcej: SCXK (Jing) 2025-0008). Szczury zostały umieszczone w Centrum Zwierząt Laboratoryjnych Szpitala Dongzhimen na Uniwersytecie Medycyny Chińskiej w Pekinie, zgodnie z krajowym standardem "Wymagania dotyczące środowiska i mieszkań zwierząt laboratoryjnych" (GB 14925-2010, wydany przez Narodowy Komitet Techniczny ds. Standaryzacji Zwierząt Laboratoryjnych Chin). Procedury eksperymentalne zostały zatwierdzone przez Komitet ds. Dobrostanu Zwierząt i Etyki Szpitala Dongzhimen na Uniwersytecie Medycyny Chińskiej w Pekinie (numer zatwierdzenia: 24-54). Odczynniki i używany sprzęt są wymienione w Tabeli Materiałów.
1. Ustanowienie modelu i zbieranie próbek
UWAGA: Przed eksperymentem poddaj wszystkim szczurom tygodniowy okres aklimatyzacji, aby zminimalizować efekty wywołane stresem.
- Rozpuszcz lipopolisacharyd Escherichia coli (LPS) w soli fizjologicznej 0,9%, aby uzyskać roztwór LPS o stężeniu 0,4 mg/mL.
- Losowo podziel 12 samców szczurów Sprague-Dawley (SD) z SPF na dwie grupy: grupę kontrolną (n = 6) oraz grupę modelu ARDS (n = 6).
- Podawaj 0,9% soli fizjologicznej dookołotchawicowej w dawce 0,5 mL/100g masy ciała grupie kontrolnej.
- Podawaj roztwór LPS o stężeniu 0,4 mg/mL poprzez wstrzyknięcie wewnątrztrachealnie w tej samej objętości (0,5 mL/100g) grupie modelowej ARDS, a po 16 godzinach znieczulaj szczury za pomocą izofluranu (zgodnie z zatwierdzonymi protokołami instytucjonalnymi).
- Pobieraj krew przez aortę brzuszną, a następnie wiruj (4 °C, 1000 × g, 30 min), aby uzyskać surowicę.
- Otwórz jamę piersiową i usuń tkankę płucną.
- Wykonaj płukanie oskrzelowo-pęcherzykowe przez kikut oskrzela za pomocą strzykawki, aby uzyskać płyn plawający oskrzelowo-pęcherzykowy (BALF)24,25.
- Pobierz prawy płat płuca i przepłukaj go PBS do dalszego badania patologicznego.
- Pobierz lewy płat płuca i przepłucz go PBS do kolejnych eksperymentów sekwencjonowania pojedynczych komórek oraz biologii molekularnej26,27.
2. Histopatologiczna ocena tkanki płucnej
- Zanurz pobraną tkankę prawego płuca szczurów w 4% paraformaldehydu w celu utrwalenia przez 12 godzin.
- Po utrwaleniu odwodnij tkankę przez stopniowo uporządkowaną serię etanolu (70%, 80%, 90%, 95%, 100% i 100%), z uwzględnieniem 2 godzin w każdym stężeniu.
- Traktuj odwodnioną tkankę ksylenem aż stanie się przezroczysta.
- Infiltruj przezroczystą tkankę stopioną parafiną w piekarniku o temperaturze 60 °C przez 2 godziny.
- Włóż blok tkanki do formy wsadowej, wlej świeżą parafinę i pozwól jej ostygnąć oraz stwardnieć w blok parafinowy.
- Podziel blok parafiny na plasterki o grubości 4 μm za pomocą mikrotomu, unieś je na wodzie w kąpieli wodnej 45 °C, aby je rozprowadzić, a następnie zbieraj za pomocą szkiełek pokrytych poli-L-lizyną.
- Piecz szkiełki w piekarniku w 60 °C, aby zapewnić solidne przyleganie kawałków do tych preparatów.
- Odparaffiniuj te fragmenty przez zanurzenie w ksylenie i nawilżanie je poprzez opadający gradient etanolu (100%, 95%, 90%, 80% i 70%).
- Przepłucz fragmenty PBS, wykonaj barwienie hematoksylyną i eozyną (HE) oraz zamontuj pokrywki żywicą neutralną.
- Oceń uraz płuc na całych preparatach skanowanych przez całe szkiełko, stosując standaryzowany system oceniania oparty na infiltracji neutrofili, zapaleniu śródmiąższowym, obrzęku i zatkaniu nosu.
UWAGA: Każdy parametr był oceniany od 0 do 4 na podstawie stopnia zamieszania patologicznego: 0 (brak); 1 (≤25%); 2 (26-50%); 3 (51-75%); 4 (>75%)28,29.
3. Test immunosorbentny związany z enzymami
UWAGA: Próbki mogą być przechowywane w temperaturze -80 °C, jeśli analiza nie może zostać przeprowadzona bezpośrednio po pozyskaniu.
- Rozpuszczaj standardy i pozwól im stać w temperaturze pokojowej przez 15 minut, a następnie wykonaj seryjne rozcieńczanie zgodnie z instrukcjami producenta.
- Przygotuj biotynowo-biotynowy roztwór roboczy przeciwciał, enzymatyczny sprzężony roztwór roboczy oraz bufor do płukania zgodnie z instrukcjami instrukcji.
- Dodaj 100 μL każdej próbki lub 100 μL seryjnie rozcieńczonych standardów do dołków pokrytych przeciwciałami i inkubuj w 37 °C przez 90 minut.
- Umyj płytkę 4 razy za pomocą podkładki do podkładania płyt, dodaj 100 μL biotinylowanego roztworu roboczego przeciwciał i inkubuj w temperaturze 37 °C przez 30 minut.
- Umyj płytkę 4 razy, dodaj 100 μL enzymatycznego sprzężenia roboczego i inkubuj w temperaturze 37 °C przez 30 minut.
- Dodaj 100 μL roztworu substratu na dołek, inkubuj chroniony przed światłem w temperaturze 37 °C przez 10 minut.
- Dodaj 100 μL roztworu stopowego na dołek, delikatnie wymieszaj i zmierz gęstość optyczną (OD) na 450 nm za pomocą czytnika mikropłyt.
- Wygeneruj krzywą standardową, renderując wartości OD względem standardowych stężeń za pomocą zaprojektowanego oprogramowania.
- Określ stężenia próbek na podstawie zmierzonych wartości przedawkowania i krzywej standardowej.
4. Przygotowanie zawiesin jednokomórkowych do sekwencjonowania
UWAGA: Po pozyskaniu tkanki płucnej należy ją umieścić na lodzie i szybko przetworzyć, aby utrzymać żywotność komórek.
- Przemyj tkankę płucną 3 razy za pomocą PBS schłodzonej do 4 °C.
- Dodaj 2 mL wcześniej schłodzonego PBS i dalej pokrój je na mniejsze kawałki sterylnymi nożyczkami na około 0,5 mm3 cale.
- Użyj pipety Pasteura, aby przenieść posiekaną tkankę razem z PBS do probówki wirówki o pojemności 15 mL.
- Dodaj dodatkowy 1 ml PBS schłodzonego do 4 °C do naczynia hodowniczego, przepłucz i ponownie zawiesij pozostałe fragmenty tkanki, a następnie przelej zawiesinę do tej samej probówki.
- Dostosuj objętość do 6,5 mL przy PBS schłodzonym do 4 °C, trawimy fragmenty tkanek w roztworze enzymatycznym zawierającym kolagenazę I (100 U/mL) i Dispazę II (1 U/mL), a następnie inkubuj delikatnie w temperaturze pokojowej przez 30 minut.
- Gdy fragmenty tkanek stają się przezroczyste i nie pozostaną wyraźne stałe fragmenty, przefiltruj strawioną tkankę przez sita komórkowe o pojemności 70 μm i 40 μm, następnie wiruj (4°C, 500 × g, 5 min) i wyrzucaj supernatant.
- Natychmiast ponownie zawiesić pellet w 3 mL buforu lizy czerwonych krwinek (RBC), delikatnie wymieszać pipetowaniem i lizować na lodzie przez 10 minut.
- Gdy supernatant wydaje się bladożółty lub prawie bezbarwny, dodaj 7 mL PBS, aby zakończyć reakcję, odwiąć (4 °C, 500 × g, 5 min) i wyrzucić supernatant.
- Dodaj 3 mL PBS, aby ponownie zawiesić komórki i określić stężenie oraz żywotność komórek za pomocą barwienia trypanowego na niebiesko.
UWAGA: Jeśli żywotność komórek spadnie poniżej 70%, wykonaj zabieg usunięcia martwych komórek.
- Przemyj komórki dwukrotnie buforem płuczącym zawierającym 0,04% surowicowej albuminy (BSA) przygotowanej w DPBS.
- Ponownie zawiesijcie komórki do liczenia i dostosowuję stężenie do 1 000 komórek/μL do dalszego obciążenia i wykrywania.
5. Kontrola jakości danych sekwencjonowania pojedynczych komórek
UWAGA: Biblioteki RNA pojedynczych komórek zostały przygotowane z zawiesin komórkowych przy użyciu metody budowy bibliotek pojedynczych komórek opartej na kroplach. Biblioteki były sekwencjonowane na platformie o wysokiej przepustowości, wykorzystując chemię sekwencjonowania opartą na nanokulkach DNA. Na bibliotekę załadowano około 21 000 komórek, co dążyło do tempa odzyskiwania komórek na poziomie ~60%. Sekwencjonowanie generowało łącznie 5,0 × 108 odczytów na bibliotekę, ze średnią głębokością 30 000–40 000 odczytów na komórkę. Surowe dane sekwencjonowania były przetwarzane przy użyciu odpowiadającego pociągu analizy RNA pojedynczej komórki, aby stworzyć macierz ekspresji genów dla każdej próbki.
- Zbuduj bibliotekę Oligo za pomocą PCR i oznaczania kodów kreskowych, zapewniając stężenie >10 ng/μL i szczytowy rozmiar 180 bp ±10 bp.
- Zbuduj bibliotekę cDNA poprzez fragmentację, naprawę końców, ligację adaptorów i PCR, zapewniając stężenie >10 ng/μL i szczytowy rozmiar 350–550 bp.
- Wygeneruj obiekt Seurat dla każdej próbki z macierzy ekspresji genów za pomocą funkcji Read10X w pakiecie Seurat (v4.3.0) w oprogramowaniu R (v4.2.0).
- Usuwanie komórek niskiej jakości poprzez filtrowanie tych z udziałem genów mitochondrialnych >10% lub hemoglobiną >5%.
- Normalizuj dane metodą "LogNormalize" z współczynnikiem skalowania 10 000.
- Zidentyfikuj cechy o wysokiej zmienności za pomocą metody "VST" i zachowaj 2000 najważniejszych genów do dalszej analizy.
- Przeprowadzaj punktację cyklu komórkowego z wykorzystaniem wbudowanych genów markerów cyklu komórkowego Seurat i cofaj się do efektów cyklu komórkowego (G2M). Score i S.Score) podczas skalowania danych.
- Wykrywanie i usuwanie potencjalnych dubletów za pomocą DoubletFinder (v2.0.3) z pN = 0,25, pK = 0,09 i oczekiwanym wskaźnikiem dubletów 0,06, a następnie wykluczenie ich z dalszego przetwarzania.
- Poprawne efekty wsadowe między samplami za pomocą pakietu Harmony (v0.1.1) z theta = 2 i max.iter.harmony = 20.
- Dodaj kolumny informacyjne dla nazw poszczególnych próbek oraz szczegóły grupowania grup modeli ARDS i grupy kontrolnej do metadanych w obiekcie Seurat.
6. Redukcja wymiarowości, klasteryzacja i adnotacja typów komórek
- Przeprowadz analizę głównych składowych (PCA) na genach o wysokiej zmienności z wykorzystaniem funkcji RunPCA w Seurat.
- Określ optymalną liczbę głównych komponentów za pomocą wizualizacji ElbowPlot.
- Zidentyfikuj główne klastry komórek, stosując algorytm Leiden za pomocą funkcji FindClusters i FindNeighbors firmy Seurat z rozdzielczością 0,1.
- Przeprowadz redukcję wymiarowości UMAP przy użyciu wybranych głównych komponentów.
- Anotuj typy komórek, identyfikując geny markerowe za pomocą funkcji FindAllMarkers.
- Wizualizuj ekspresję genów markerowych między typami komórek za pomocą funkcji DotPlot.
- Dodaj nową kolumnę informacji z wynikami adnotacji komórek do metadanych i wizualizuj adnotację typu komórki za pomocą Umap.
- Oblicz proporcje każdej grupy eksperymentalnej w każdym zidentyfikowanym typie komórek za pomocą funkcji dplyra do liczenia komórek i obliczania procentów.
- Wygeneruj ułożone wykresy słupkowe za pomocą ggplot2, aby zilustrować skład grup dla typu komórek, w tym etykiety procentowe i porównania statystyczne między grupami.
7. Analiza subpopulacji
- Wyodrębnij podpopulację komórek docelowych za pomocą funkcji podzbioru i stwórz odpowiadający mu obiekt Seurat do analizy na poziomie podpopulacji w dalszej fazie.
- Ponownie normalizować i skalować dane subpopulacji, jednocześnie cofając geny mitochondrialne, geny hemoglobiny, liczbę UMI oraz wyniki cyklu komórkowego.
- Wykonaj wtórną analizę PCA i zastosuj korekcję wsadową Harmony specjalnie dla subpopulacji.
- Przeprowadz ponowne klasteryzowanie za pomocą funkcji FindNeighbors i FindClusters, aby zidentyfikować odrębne podtypy w ramach podpopulacji komórek docelowych.
- Generuj nowe osadzenia UMAP i wizualizuj klastry podpopulacji z adnotacjami typów komórek, stosując tę samą wcześniej opisaną metodę.
8. Różnicowa analiza ekspresji genów
- Zidentyfikuj geny różnicowo ekspresyjne (DEG) pomiędzy grupą modelową ARDS a grupą kontrolną dla każdego podklastra za pomocą metody statystycznej MAST.
- Włącz regresję zmiennych utajonych do kontroli efektów głębokości sekwencjonowania podczas testów różnicowych ekspresji.
- Filtruj DEG na podstawie zdefiniowanych progów zmien log fold oraz skorygowanych wartości p za pomocą funkcji FindMarkers.
- Generuj wykresy wulkanów z funkcją ggplot dla każdego podklastra, aby zwizualizować różnice wzorców ekspresji genów z oznakowaniem istotności.
- Eksportuj wyniki DEG i twórz kompleksowe tabele podsumowujące, aby ułatwić wzbogacenie funkcjonalne i inne analizy dalsze.
9. Analiza wzbogacenia
- Przypisuj symbole genów do identyfikatorów ENTREZ za pomocą org. Rn.eg.db baza danych adnotacji dla kompatybilności analizy ścieżek.
- Połącz DEG ze wszystkich podklastrów i filtruj istotne geny na podstawie progów zmiennych log fałdów oraz skorygowanych wartości p (test sumy rangi Wilcoxona, logfc.threshold = 0,25, p_val < 0,05).
- Wykonaj analizę wzbogacenia ontologii genów (GO) z wykorzystaniem funkcji wzbogacenia GO z kategoriami procesów biologicznych, funkcji molekularnych, ontologii komponentów komórkowych oraz metodą korekcji Benjamini-Hochberga.
- Przeprowadz analizę wzbogacenia szlaków w Encyklopedii Genów i Genomów w Kioto (KEGG) z wykorzystaniem funkcji wzbogacania KEGG w bazę danych organizmów szczurów oraz zastosuj progi istotności statystycznej.
- Generuj wykresy krorowe dla wyników wzbogacania GO i KEGG za pomocą funkcji wykresu punktowego, pokazując 15 najwyżej wzbogaconych składników o stosunkach genów i skorygowanych wartościach p.
- Eksportuj wyniki wzbogacania w postaci plików CSV i zapisuj wykresy wizualizacyjne do kompleksowej analizy ścieżek.