Próbka wykorzystana w tym protokole stanowiła część badania Study of Muscle, Mobility, and Aging (SOMMA)15, które zostało zatwierdzone przez komisję Institutional Review Board Western IRB-Copernicus Group (WCG) i przeprowadzone zgodnie z Deklaracją Helsińską. Uczestnicy wyrazili pisemną świadomą zgodę na udział w badaniu.
UWAGA: Niniejszy protokół został zaadaptowany z wcześniejszego protokołu wykorzystującego 100 mg ludzkiej tkanki tłuszczowej podskórnej z okolic brzucha na platformie opartej na nanodołkach16. Obecny protokół został zoptymalizowany dla 50 mg ludzkiej IMAT i przygotowania biblioteki przy użyciu platformy opartej na kropelkach. Może być wymagana dalsza optymalizacja tego protokołu w celu izolacji jąder z nie-ludzkiej IMAT lub innych depotów tłuszczowych.
1. Przygotowanie buforów i odczynników (Tabela 1 i Tabela 2)
UWAGA: Bufor należy przygotować na świeżo w dniu eksperymentu i nie używać ponownie.
- Wstępnie schłodzić wirówkę do 4 °C.
- Przygotować bufor do homogenizacji oraz medium do izolacji jąder.
- Przygotować dwa pojemniki z lodem i wstępnie schłodzić 2 x 15 mL probówki stożkowe.
- Wymieszać wszystkie odczynniki do buforu do homogenizacji (HB) w probówce stożkowej 15 mL w kolejności wymienionej w Tabeli 1. Przechowywać w lodzie. Wymieszać za pomocą wierzchołka (vortex).
- Wymieszać wszystkie odczynniki do medium do izolacji jąder (NIM) w probówce stożkowej 15 mL zgodnie z listą w Tabeli 2. Przechowywać w lodzie. Wymieszać za pomocą wierzchołka (vortex).
- Przygotować 10% Triton-X, dodając 100 µL Triton X-100 do 900 µL wody wolnej od nukleaz. Wymieszać za pomocą wierzchołka (vortex), aby zapewnić odpowiednie połączenie składników. Przechowywać w temperaturze pokojowej (RT).
| Odczynnik | Objętość (μL) | Stężenie końcowe (mM) |
| 1x | 2x |
| 1 M MgCl22 | 10 | 20 | 5 |
| Bufor Tris 1 M, pH 8,0 | 20 | 40 | 10 |
| 2 M KCl | 25 | 50 | 25 |
| 1,5 M sacharoza (-4oC) | 334 | 668 | 250 |
| 1 mM DTT | 2 | 4 | 0.001 (~1 µM) |
| inhibitor proteaz 100x | 20 | 40 | 1x |
| Superasin 20 U/µL | 40 | 80 | 0,4 U/µL |
| Woda wolna od nukleaz | 1549 | 3098 | - |
| Objętość całkowita | 2000 | 4000 | - |
Tabela 1: Bufor do homogenizacji (HB). Przechowywać w lodzie. Mieszać za pomocą wortexa.
| Odczynnik | Objętość (μL) | Stężenie końcowe (mM) |
| 1x | 2x |
| EDTA | 0.4 | 0.8 | 0.1 |
| Inhibitor RNazy Ribolock (40U/µL) | 40 | 80 | 0.8 U/µL |
| 1% BSA-PBS (-/-) | 1959.6 | 3919.2 | - |
| Objętość całkowita | 2000 | 4000 | - |
Tabela 2: Medium do izolacji jąder (NIM). Przechowywać w lodzie. Mieszać za pomocą wortexa.
2. Pulweryzacja zamrożonej tkanki ( Rysunek 2A)
- Przygotuj stanowisko do homogenizacji.
- Napełnij pojemnik ciekłym azotem (LN2).
OSTROŻNIE: podczas pracy z LN2 zawsze należy nosić gogle i rękawice kriogeniczne.
- Przygotuj 2 moździerze, 1 tłuczek, 1 szpatułkę mikro-czerpak, 1 szklany homogenizator Dounce oraz 1 stalowy tłuczek do automatycznego homogenizatora.
- Uruchom automatyczny homogenizator.
- Wypełnij zlew lodu i schłodź w nim szklany homogenizator.
- Wypełnij 2 moździerze (zawierające tłuczek i szpatułkę) LN2, aby schłodzić narzędzia. Odczekaj do odparowania LN2 i powtórz czynność.
- Podczas chłodzenia narzędzi, dodaj 1 mL buforu HB do szklanego homogenizatora.
- Wypełnij oba moździerze LN2 po raz ostatni i wsyp 50 mg próbki IMAT do jednego z moździerzy.
- Rozdrobnij IMAT za pomocą tłuczka, delikatnie naciskając na fragment tkanki, aby rozbić go na małe kawałki. Upewnij się, że wszystkie fragmenty zostały sproszkowane.
- LN2 będzie powoli odparowywać podczas rozdrabniania tkanki. Gdy tkanka zostanie odpowiednio sproszkowana, a w moździerzu pozostanie jeszcze ¼ - ½ objętości LN2, przechyl moździerz w stronę krawędzi, aby zebrać sproszkowaną tkankę przy brzegu. Pozostaw LN2 do całkowitego odparowania.
- Natychmiast po odparowaniu resztek LN2, przełóż sproszkowaną tkankę do szklanego homogenizatora zawierającego 1 mL HB.
3. Homogenizacja sproszkowanej tkanki
- Zhomogenizować sproszkowaną tkankę za pomocą automatycznego homogenizatora typu dounce. Przesuwać szklany tłuczek w górę i w dół stalowego pestle, wykonując 10 ruchów w kierunku przednim, a następnie 10 ruchów w kierunku odwrotnym.
- Upewnić się, że po homogenizacji roztwór jest mętny i nie zawiera widocznych fragmentów tkanki. Ze względu na zanieczyszczenie tkanką mięśniową często oczekiwany jest jasnoróżowy kolor.
- Przenieść homogenizat do schłodzonej wcześniej probówki low-bind o pojemności 1,7 mL, umieszczonej w lodzie.
- Użyć 400 µL buforu HB do wypłukania homogenizatora, aby upewnić się, że cały materiał został przeniesiony, a następnie dodać go do probówki.
UWAGA: Jednocześnie można przetworzyć dwie próbki. W tym celu należy podwoić ilość HB oraz NIM. Sproszkować i zhomogenizować jedną próbkę tkanki, a bezpośrednio po niej sproszkować i zhomogenizować drugą próbkę tkanki, aby móc przeprowadzić etapy izolacji i oczyszczania równolegle.
4. Izolacja i oczyszczanie jąder komórkowych (Rysunek 2B)
- Dodaj 14 µL Triton-X (10%) do homogenatu, aby uzyskać stężenie 0,1%.
- Przechowuj probówkę na lodzie i w ciemności przez 10-15 min, mieszając na wirówce vortex co 3 min.
- Namocz jedno sitko komórkowe 100 µm i jedno sitko komórkowe 40 µm (na próbkę), używając 100 µL RT DPBS dla każdego z nich w stożkowej probówce 50 mL.
- Przefiltruj homogenat przez sitko komórkowe 100 µm.
- Wypłucz probówkę 1,7 mL za pomocą 400 µL HB i przefiltruj przez sitko komórkowe 100 µm.
- Następnie przefiltruj roztwór przez sitko komórkowe 40 µm.
- Przenieś równą ilość roztworu do dwóch schłodzonych wcześniej probówek 1,7 mL typu low-bind, co odpowiada ok. 900 µL w każdej probówce.
- Wiruj probówki przez 10 min przy 2700 x g w temperaturze 4 °C. Po wirowaniu powinien być widoczny mały pellet.
- Usuń i odrzuć górną warstwę lipidową oraz pozostały supernatant, pozostawiając ok. 50 µL roztworu w pierwszej probówce.
- Powtórz czynność dla drugiej probówki.
- Dokładnie resuspenduj pellet w pierwszej probówce, delikatnie pipetując w górę i w dół 20 razy, a następnie przenieś do nowej probówki 1,7 mL typu low-bind. Unikaj tworzenia pęcherzyków powietrza.
- Powtórz tę procedurę dla drugiej probówki i przenieś zresuspendowany roztwór do tej samej probówki.
- Dodaj 500 µL NIM i wymieszaj poprzez pipetowanie.
- Wirowuj probówkę z przeciwwagą przy 1000 x g przez 10 min w temperaturze 4 °C.
- Usuń supernatant, pozostawiając ok. 50 µL, i delikatnie pipetuj w górę i w dół, aż pellet zostanie zresuspendowany. Opcjonalnie przenieś zresuspendowany pellet do nowej, czystej probówki, jeśli na ściankach probówki pozostały resztki lipidów.
- Dodaj 200 µL NIM i wymieszaj poprzez pipetowanie.
5. Barwienie i liczenie jąder komórkowych (Rycina 2C oraz Rycina 3)
UWAGA: Aby ułatwić zliczanie, należy skonfigurować protokół „zliczania jąder” w automatycznym liczniku komórek, ponieważ regulacja kanałów jasnego pola oraz DAPI może znacząco wpłynąć na wynik zliczania. Dostrój kanały w taki sposób, aby rejestrowane były tylko jądra, a nie zanieczyszczenia. Upewnij się, że kanał jasnego pola oznacza jedynie „obiekty”, które posiadają również barwienie DAPI.
- Dodać 1 kroplę roztworu do barwienia żywych komórek i pozostawić w ciemności, w lodzie, na 15 min.
- Przefiltrować roztwór przez sitko komórkowe 30 µm.
- Wymieszać roztwór jąder metodą pipetowania i nanieść 10 µL roztworu na szkiełko z komorą do liczenia komórek.
- Policzyć jądra za pomocą automatycznego licznika komórek.
UWAGA: Optymalne stężenie wynosi 1000 jąder/µL, co odpowiada 1,0 x 106/mL.
- Upewnić się, że nie występują skupiska jąder, ponieważ mogłyby one zapchać chip do generowania kropel z pojedynczymi jądrami (Rysunek 3).
- Jeśli stężenie jąder jest zbyt niskie, odwirować roztwór przy 1000 x g przez 10 min w 4 °C w celu uzyskania osadu, usunąć nadsącz i resuspender w mniejszej objętości.
- Jeśli w roztworze znajduje się dużo zanieczyszczeń (debris), resuspender roztwór jąder w większej objętości NIM (np. 1 mL) i ponownie przefiltrować przez sitko komórkowe 30 µm. Następnie odwirować przy 1000 x g przez 10 min w 4 °C i resuspender w odpowiedniej objętości w zależności od stężenia jąder.
- Po uzyskaniu odpowiedniego stężenia jąder przejść bezpośrednio do pierwszego etapu przygotowania biblioteki.

Rysunek 3: Barwienie wyizolowanych jąder. Obraz z licznika komórek przedstawiający jądra zabarwione NucBlue/DAPI (obraz po lewej) oraz odpowiadający mu obraz w polu jasnym (obraz po prawej). Na obrazie w polu jasnym widoczna jest obecność niewielkiej ilości zanieczyszczeń. Użyty tutaj automatyczny licznik komórek nie posiada opcji dodawania pasków skali. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
6. Przygotowanie biblioteki i parametry sekwencjonowania
- Należy zapoznać się z dokładnym protokołem przygotowania bibliotek z wykorzystaniem metody kropelkowej dla pojedynczych jąder komórkowych, dostępnym na stronie internetowej dostawcy17.
- Należy dążyć do odzyskania docelowo 10 000 jąder. Jednak w przypadku próbek o wysokim poziomie zanieczyszczeń lub zawierających kruche jądra, należy spodziewać się uzyskania mniejszej liczby odzyskanych jąder.
- Aby umożliwić równoległe przetwarzanie większej liczby próbek, próbki można przechowywać w temperaturze 4°C do 72 h po wykonaniu kroku 2.3. protokołu przygotowania bibliotek. W tym celu należy przetwarzać dwie próbki do kroku 2.3. przez dwa kolejne dni, a w trzecim dniu przetworzyć wszystkie 4 próbki wspólnie, zaczynając od kroku 3. protokołu przygotowania bibliotek.
- Parametry sekwencjonowania: Sekwencjonowanie należy przeprowadzić na platformie sekwencyjnej, dążąc do uzyskania 50 000 odczytów par końcowych (paired-end reads) na jądro.
UWAGA: Dane przedstawione w niniejszym protokole zostały zsekwencjonowane na platformie NovaSeq 6000, przy założeniu 50 000 odczytów par końcowych na jądro.

Rysunek 2: Schemat procedury. Schematyczne przedstawienie przebiegu pracy w (A) krokach 2 i 3, (B) kroku 4 oraz (C) kroku 5 protokołu. Rysunek został przygotowany przy użyciu BioRender.com. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
7. Przetwarzanie i analiza danych
UWAGA: W niniejszym protokole krótko przedstawiono niektóre z zalecanych programów i pakietów R wykorzystywanych do przetwarzania uzyskanych danych z sekwencjonowania, koncentrując się na etapach następujących po wstępnym przetwarzaniu (Tabela 3). Niniejsze badanie dostarcza ogólnych wskaźników kontroli jakości (QC) oraz przykładu przybliżonej i rzutowanej projekcji jednostajnej rozmaitości (UMAP) w Rycina 4Jednak szczegółowy opis analizy bioinformatycznej wykracza poza zakres niniejszego protokołu. W związku z tym czytelnicy mogą zapoznać się z niedawnym przeglądem najlepszych praktyk w analizie pojedynczych komórek autorstwa Heumos i wsp.18.
- Wstępne przetwarzanie danych z sekwencjonowania
- Zmapuj odczyty z pojedynczych jąder do ludzkiego genomu referencyjnego GRCh38.
- Uwzględnij odczyty intronowe w zliczaniu.
- Przeprowadź kontrolę jakości (QC) i filtrowanie danych przy użyciu pakietu Seurat R19.
- Oblicz wskaźnik złożoności komórek, dzieląc log(10) liczby wykrytych genów przez log(10) liczby wykrytych odczytów.
- Przedstaw najważniejsze metryki QC na histogramach lub wykresach skrzypcowych, w tym liczbę genów wykrytych na jądro, procent odczytów mitochondrialnych oraz wskaźnik złożoności komórek.
- Odrzuć jądra z liczbą genów mniejszą niż 200 lub większą niż 10 000 na jądro, odsetkiem odczytów mitochondrialnych przekraczającym 10% oraz wskaźnikiem złożoności poniżej 0.8.
- Znormalizuj dane i przeprowadź redukcję wymiarowości.
- Użyj funkcji SCTransform z pakietu Seurat , aby znormalizować dane przy użyciu 2000 cech zmiennych.
- Pogrupuj dane (klasteryzacja), korzystając z następujących funkcji pakietu Seurat R: RunPCA, FindNeighbors, FindClusters oraz RunUMAP.
- Wygeneruj wykres UMAP, aby zwizualizować klastry danych.
- Odrzuć przewidywane dublety przy użyciu pakietu DoubletFinder R 20 i przeprowadź ponowną klasteryzację danych.
- Adnotuj klastry, korzystając ze znanych markerów genowych typów komórek spodziewanych w tkance (podejście nadzorowane) lub na podstawie 5 genów o największej różnicy w ekspresji pomiędzy klastrami (podejście nienadzorowane).
- Użyj narzędzia decontX21 , aby określić stopień kontaminacji tłem RNA (ambient RNA) i skorygować macierz ekspresji genów pod kątem tego zjawiska.
- Uwzględnij surową macierz genów jako tło.
- Zapisz obiekt Seurat do przyszłej analizy danych.
UWAGA: Kod do analizy QC i klasteryzacji znajduje się w Pliku uzupełniającym 1.
| Oprogramowanie/pakiety R wykorzystane w przepływie danych | Alternatywne oprogramowanie/pakiety | Etap przetwarzania |
| CellRanger | STARsolo, kallisto | Przycinanie, wyrównywanie, mapowanie |
| Seurat | SingleCellExperiment, Cellranger | Kontrola jakości, analiza i eksploracja danych |
| DoubletFinder | scds, scdblFinder, Scrublet | Wykrywanie dubletów |
| DecontX | SoupX, CellBender | Korekta RNA tła |
Tabela 3: Oprogramowanie/narzędzia do przepływu danych.