Protokół ten został opracowany i przetestowany na podstawie próbek zebranych w ramach badań zatwierdzonych przez Instytucjonalną Radę Przeglądową UCSF.

Rysunek 1: Metagenomiczne sekwencjonowanie nowej generacji (mNGS) i analityczny workflow wykrywania patogenów w zakresie wykrywania patogenów. Ten diagram pokazuje podstawowe etapy w szlakach mNGS w laboratorium i bioinformatycznym mózgowo-rdzeniowym. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
1. Rozważania dotyczące obsługi próbek CSF
UWAGA: wyniki mNGS są optymalne, jeśli przygotowanie biblioteki CSF odbywa się w momencie pobierania, lub jeśli próbki CSF są zachowane przez szybkie zamrażanie w ciekłym azocie, a następnie długotrwałe przechowywanie w temperaturze -80 °C15. Alternatywnie, stosowanie roztworów konserwujących, takich jak osłona DNA/RNA, w momencie pobierania próbki może umożliwić przechowywanie w temperaturze pokojowej przez dłuższy czas16,17. Może to być korzystne w sytuacjach, gdy pojawiają się obawy dotyczące niezawodności łańcucha chłodniczego. Ponadto unikanie niepotrzebnych cykli zamrożenia i rozmrażania może również pomóc zachować integralność kwasów nukleinowych, a zaleca się wcześniejsze aliquotowanie próbek w momencie pobrania, jeśli przewiduje się wiele procesów analitycznych dla próbek CSF18,19.
2. Ekstrakcja kwasu nukleinowego z CSF
UWAGA: Wykonuj pracę w odpowiedniej szafce o poziomie bezpieczeństwa biologicznego, w rękawiczkach i fartuchu laboratoryjnym, aby zmniejszyć ryzyko przenoszenia organizmów zakaźnych. Odczynniki są szkodliwe, jeśli zostaną połknięte, wdychane lub mają kontakt ze skórą. Utylizacja zużytych materiałów jako niebezpiecznych odpadów chemicznych zgodnie z lokalnymi przepisami.
- Zestaw Quick-DNA/RNA zazwyczaj dobrze sprawdza się w próbkach o niskiej zawartości kwasów nukleinowych.
- Rozmroz próbki CSF na lodzie, następnie przelej 100 μL–1 mL każdej próbki do probówki o pojemności 1,5 mL i wiruj w 16 000 x g w temperaturze 4 °C przez 10 minut.
- Ostrożnie odpipetuj supernatant, pozostawiając za sobą granulkę (która często jest niewidoczna gołym okiem) oraz 100 μL supernatantu.
- Dodaj 100 μL tarczy DNA/RNA do pozostałych 100 μL supernatantu i pelletu, a następnie dokładnie wymieszaj pipetą.
- Następnie dodaj 600 μL buforu lizycznego DNA/RNA i mieszanki pipet, aż roztwór stanie się jednorodny.
- Postępować zgodnie z resztą protokołu ekstrakcji DNA/RNA, opisanym przez producenta w dodatkach uzupełniających A, z następującymi istotnymi modyfikacjami specyficznymi dla płynu mózgowo-rdzeniowego:
- Wszystkie etapy wirowania wykonuj przez 1 minutę z maksymalną prędkością (~21 000 x g), z wyjątkiem ostatniego etapu płukania DNA/RNA, który powinien trwać 3 minuty przy maksymalnej prędkości.
- Elautuj próbkę, dodając 22 μL (zamiast 25 μL) wody wolnej od nukleazy bezpośrednio do matrycy kolumny, a następnie inkubuj przez 3 minuty w temperaturze pokojowej. Większa objętość zazwyczaj powoduje większą całkowitą ilość kwasu nukleinowego udzielonego z kolumny.
- Po elucii pipetować przepływ z powrotem na matrycę kolumny i ponownie wirować przez minutę z maksymalną prędkością. Przełóż próbkę do nowej probówki magazynowej 1,5 mL i natychmiast przystąpi do przygotowania bibliotecznego lub przechowuj w temperaturze -80 °C.
3. Przygotowanie biblioteki RNA CSF
UWAGA: Wykonuj pracę w odpowiedniej szafce zabezpieczeń biologicznych, w rękawiczkach i fartuchu laboratoryjnym. Odczynniki są szkodliwe, jeśli zostaną połknięte, wdychane lub mają kontakt ze skórą. Utylizacja zużytych materiałów jako niebezpiecznych odpadów chemicznych zgodnie z lokalnymi przepisami.
UWAGA: Ważne kwestie przed rozpoczęciem: Redukcja rRNA jest preferowana przed wzbogaceniem poliadenylacji przy izolacji mRNA, ponieważ wyczerpanie rRNA nie wymaga dodatkowych kroków czyszczenia kulek (co może zmniejszyć wydajność użytecznego mRNA), a niektóre transkrypcje wirusa RNA neuroinwazyjnego nie mają poliadenylacji. ERCC (Konsorcjum Kontroli Zewnętrznego RNA) RNA Spike-In może być dodane jako środek kontroli jakości, choć nie jest to ściśle konieczne. Jeśli RNA było przechowywane w temperaturze -80 °C po ekstrakcji, rozmrozić na lodzie przed rozpoczęciem.
- Fragmentacja i priming RNA: Pipeta miesza następujące składniki w probówce reakcji łańcuchowej polimerazy (PCR): próbkę RNA o pojemności 3,5 μL, 0,5 μL rozcieńczonego spike-in ERCC spike-in, 4 μL buforu reakcji pierwszej nici (5x), 1 μL zapalników losowych oraz 1 μL rozcieńczonego rRNA w skali 1:100 (Tabela materiałów). W termocyklerze ustaw podgrzewaną pokrywę na 105 °C i inkubuj próbki w temperaturze 75 °C przez 2 minuty, 70 °C przez 2 minuty, 65 °C przez 2 minuty, 60 °C przez 2 minuty, 55 °C przez 2 minuty, 37 °C przez 5 minut i 25 °C przez 5 minut.
- Synteza cDNA w pierwszej nici: Pipette miesza rozdrobnione i przygotowane RNA (10 μL) z 8 μL wodą wolną od nukleazy oraz 2 μL mieszanką enzymów do syntezy pierwszej nici. W termocyklerze ustaw podgrzewaną pokrywę na 105 °C i inkubuj próbki w temperaturze 25 °C przez 10 minut, 42 °C przez 15 minut oraz 70 °C przez 15 minut.
- Synteza cDNA na drugiej nici: Pipette miesza DNA syntetyzowane z pierwszej nici (20 μL) z buforem reakcji reakcyjnej syntezy drugiej nici o pojemności 8 μL, mieszanką enzymów syntezy drugiej nici oraz 48 μL wodą wolną od nukleazy. W termocyklerze inkubuj próbki przez 1 godzinę w temperaturze 16 °C, z zamkniętą podgrzewaną pokrywką.
- Wykonaj oczyszczenie kulek magnetycznych (Supplementary Appendices B) z użyciem proporcji 1,8x kulek magnetycznych SPRI (144 μL). Elutuj do wody wolnej od nukleazy 53 μL i przelej 50 μL ostatniego supernatantu (zawierającego oczyszczone dwuniciowe cDNA) do czystej rurki PCR wolnej od nukleazy. W tym momencie próbki można bezpiecznie zamrozić w temperaturze -20 °C przez noc, jeśli zajdzie taka potrzeba.
- Końcowe przygotowanie biblioteki cDNA: Jeśli próbka była przechowywana w nocy w temperaturze -20 °C, przed ponownym uruchomieniem próbki rozmrozić na lodzie. Pipette miesza oczyszczone ds-cDNA (50 μL) z 7 μL buforu reakcji końcowej oraz 3 μL mieszanki enzymów końcowych.
- W termocyklerze inkubuj próbki w temperaturze 20 °C przez 30 minut oraz w 65 °C przez 30 minut, z zamkniętą podgrzewaną pokrywką.
- Ligacja adaptera (wykonaj ten krok na lodzie): Stwórz rozcieńczenie adaptera w stosunku 1:100 w wodzie wolnej od nukleazy. Pipetą przełóż mieszaninę reakcji końcowej (60 μL) do probówki zawierającej 30 μL mieszankę ligacyjną, 1 μL wzmacniacz ligacji oraz 2,5 μL rozcieńczony adapter 1:100. Upewnij się, że adapter jest dodawany osobno (np. nie jest wstępnie mieszany z mieszanką ligation master lub wzmacniaczem ligacji), aby uniknąć powstawania dimera adaptera.
- W termocyklerze inkubuj próbki w temperaturze 20 °C przez 15 minut, z zamkniętą podgrzewaną pokrywką.
- Natychmiast przejść do oczyszczania kulek magnetycznych (Dodatek 2) przy użyciu proporcji 0,9x kulek magnetycznych SPRI (87 μL).
- Elutuj do wody wolnej od nukleazy o pojemności 17 μL i przelej 15 μL ostatniego supernatantu do czystej, beznukleazowej rurki PCR
- PCR: Pipeta oczyszczone, adapterowo ligowane cDNA (15 μL) do probówki z enzymem USER o pojemności 3 μL, 25 μL mieszanką Q5 Master Mix oraz 10 μL unikalnymi zapalnikami z kodami kreskowymi. W termocyklerze ustaw podgrzewaną pokrywę na 105 °C, inkubuj próbki w 37 °C przez 15 minut, następnie 98 °C przez 30 s, a następnie wykonaj 19 cykli 98 °C przez 10 s i 65 °C przez 75 s. Zakończ PCR inkubacją w temperaturze 65 °C przez 5 minut.
- Wykonaj końcowe oczyszczenie kulek magnetycznych, używając proporcji 0,8x kulek magnetycznych (43 μL). Eluluj do wody wolnej od nukleazy o pojemności 23 μL, a po ostatnim etapie oddzielania kulek przelej 20 μL do czystej rurki PCR wolnej od nukleazy. W tym momencie próbki można bezpiecznie zamrozić w temperaturze -20 °C przez noc, jeśli zajdzie taka potrzeba.
4. Przygotowanie biblioteki DNA w CSF
UWAGA: Wykonuj pracę w odpowiednim szafie o poziomie bezpieczeństwa biologicznego, w rękawiczkach i fartuchu laboratoryjnym. Odczynniki są szkodliwe, jeśli zostaną połknięte, wdychane lub mają kontakt ze skórą. Utylizacja zużytych materiałów jako niebezpiecznych odpadów chemicznych zgodnie z lokalnymi przepisami.
- Jeśli wyodrębnione DNA przechowywano w nocy w temperaturze -20 °C, rozmrozić na lodzie przed ponownym uruchomieniem.
- Fragmentacja DNA i przygotowanie końcowe: Dokładnie wymieszaj bufor reakcji pierwszej nici, wykonując wiry i mieszając pipetą roztwór, aby ponownie zawiesować wytrącanie. Pipeta miesza 3,5 μL wyodrębnionego DNA, 22,5 μL wody wolnej od nukleazy, 7 μL buforu reakcji pierwszej nici oraz 2 μL mieszanki enzymów pierwszej nici.
- W termocyklerze ustaw podgrzewaną pokrywę na 105 °C i inkubuj próbki w 37 °C przez 5 minut, a następnie 65 °C przez 30 minut.
- Proces ligacji i kodowania kreskowego PCR jest identyczny dla bibliotek DNA i RNA. Powtórz sekcje 7–10 "Przygotowania biblioteki RNA CSF" z fragmentowanym i przygotowanym DNA, aby zakończyć przygotowanie biblioteki DNA.
5. Kontrola jakości bibliotek i łączenie w puli
- Ilościowo określ stężenie bibliotek próbek za pomocą zestawu do ilościowania DNA i odpowiadającego fluorometru. Kontrola wody powinna mieć znacznie niższe (lub niemierzalne) stężenia w porównaniu do pozostałych próbek.
- Oceń rozmiar biblioteki próbek, uruchamiając je na zautomatyzowanej maszynie do elektroforezy kapilarnej.
- Alternatywnie, jeśli nie jest dostępna zautomatyzowana maszyna do elektroforezy kapilarnej, można dodatkowo wzmocnić 1-2 μL próbki poprzez reakcję PCR z uniwersalnymi zapalnikami Illumina, a następnie wykonać elektroforezę 10 μL produktu końcowego na żelu agarozowym o stężeniu 2%. Utylizacja zużytych materiałów jako niebezpiecznych odpadów chemicznych zgodnie z lokalnymi przepisami.
UWAGA: Optymalna długość biblioteki to około 400–500 par zasad, dlatego przyszłe inkubacje fragmentacyjne powinny być odpowiednio dostosowane, jeśli rozmiary bibliotek będą większe lub krótsze niż pożądane. Ponadto kluczowe jest ocenenie obecności dimerów adaptera, które mają około 150 par zasad długości i występują, gdy dwie cząsteczki adaptera przypadkowo się ze sobą łączą bez sekwencji wstawek. Nawet jeśli są obecne w niewielkich ilościach, dymery adapterowe mają tendencję do efektywnego skupiania się i sekwencjonowania na komórce przepływowej, zmniejszając udział użytecznych odczytów podczas sekwencjonowania. W związku z tym, jeśli obecne są dimery adaptera, wykonaj dodatkowe oczyszczanie kulek magnetycznych z 0,8x stosunkiem kulek magnetycznych SPRI. (Warto zauważyć, że jeśli >10% próbki to ściemniacz z adaptera, może być potrzebne 2-3 rundy oczyszczania kulek magnetycznych, aby usunąć większość dimerów adaptera).
- Aby zapewnić podobną głębokość sekwencjonowania między próbkami w ramach tego samego cyklu sekwencjonowania, konieczne jest połączenie ilości równomolowych każdej próbki. Jeśli rozmiary bibliotek są porównywalne we wszystkich próbkach, wystarczy zgrupować równoważne masy każdej biblioteki RNA i powtórzyć oddzielnie dla każdej biblioteki DNA.
6. Rozważania dotyczące sekwencjonowania
- Zebrane próbki są już gotowe do sekwencjonowania na urządzeniach sekwencjonujących Illumina. Jeśli sekwencjonowanie odbywa się wewnętrznie (a nie przez specjalistyczną bazę sekwencjonową), starannie denaturuj i rozcieńczaj zgromadzone biblioteki zgodnie z instrukcjami urządzenia sekwencjonującego przed załadowaniem puli na wkład odczynnikowy (typowe stężenie wynosi 4 nM dla większości platform sekwencjonowania Illumina).
- Oblicz szacowaną głębokość sekwencjonowania, dzieląc oczekiwaną wydajność danych dla komórki przepływowej Illumina przez liczbę multipleksowanych próbek w cyklu.
UWAGA: Określenie pożądanej głębokości sekwencjonowania zależy od kilku czynników, w tym kosztu, używanego urządzenia sekwencjonującego oraz celów eksperymentu. Dwa przykłady różnych głębokości sekwencjonowania przedstawiono poniżej w sekcji Reprezentatywne Rezultaty.
7. Analiza danych w celu wykrywania patogenów
- Wykorzystaj otwartoźródłową platformę internetową Chan Zuckerberg ID (CZID, https://www.czid.org, Illumina mNGS Pipeline v8.3) do przeprowadzenia analizy metagenomowej wykrywania patogenów. CZID zawiera także opcje przesyłania próbek przez interfejs wiersza poleceń oraz bezpośredniego przesyłania z konta Illumina BaseSpace do CZID. Szczegółowy przegląd analizy danych mNGS CZID jest dostępny na https://chanzuckerberg.zendesk.com/hc/en-us/articles/13770737266196-Guide-to-mNGS-Data-Analysis
- Zaloguj się do CZID, kliknij "upload", wybierz "metagenomics" dla typu analizy i wybierz wejściowe pliki fastq z uruchomienia sekwencjonowania. Wprowadź niezbędne informacje dla każdej próbki, w tym nazwę próbki, typ próbki (CSF) oraz nukleotyd (RNA lub DNA).
UWAGA: Przesłane pliki przechodzą automatyczne przetwarzanie za pomocą pipeline'u CZID, który wykorzystuje kilka narzędzi obliczeniowych, takich jak STAR, Bowtie2, Trimmomatic, PriceSeq i GSNAP, między innymi, aby usunąć niskiej jakości odczyty i sekwencje odczytu ludzkiego. Następnie proces wyrównuje i składa filtrowane sekwencje, wykorzystując bazy danych nukleotydów (NT) i białek (NR) National Center for Biotechnology Information (NCBI) odpowiednio za pośrednictwem programów Minimap2 i DIAMOND. Ostateczny wynik CZID obejmuje liczby taksonów NT/NR oraz liczby kontigów (ciągłe, nakładające się segmenty generowane podczas montażu) dla każdej próbki.
- Po przetworzeniu kliknij folder projektu zawierający próbki i przejdź do pulpitu podsumowań.
- Przewiń w dół, aby sprawdzić liczbę odczytów na próbkę (i porównaj z oczekiwaną liczbą odczytów na podstawie użytej komórki przepływowej Illumina oraz liczby próbek multipleksowanych w sekwencjonacji), procent odczytów przechodzących przez filtrację kontroli jakości oraz duplikatowy współczynnik kompresji, aby ocenić przejęcie nadmiernego wzmocnienia PCR z użyciem stronniczości.
- Obecność skażenia fragmentami kwasów nukleinowych środowiskowych – pozyskanych podczas pobierania próbek płynu mózgowo-rdzeniowego lub podczas ekstrakcji i przygotowania biblioteki – jest częstym zjawiskiem20. Aby ustalić, czy odczyty bakterialne, grzybicze i pasożytnicze są prawdziwymi pozytywnymi wynikami, można rozważyć potencjalne metody wykorzystujące znormalizowane liczby odczytów (odczyty na milion, czyli rPM):
- Stwórz model tła (np. używając próbek wody, które zostały wyekstrahowane, przygotowane i zsekwencjonowane) bezpośrednio w CZID, aby odfiltrować taksony z rPM podobnymi lub mniejszymi pod względem liczebności względem modelu tła (a więc prawdopodobnie reprezentującymi skażenie). W folderze projektu sprawdź każdą próbkę kontroli wody, a następnie kliknij "model tła", aby utworzyć. Podczas analizy każdej próbki wybierz model tła i kliknij "filtry progowe", aby zastosować "NTZ Score" > 1 (lub wyższym), aby odfiltrować prawdopodobne zanieczyszczenia.
- Wykorzystaj wcześniej określone progi rPM, takie jakrpM próba > 10* (średnie rPM kontrolne) lub log10 (próba rPM) co najmniej o 1 log większy niż log 10 (średnia rPMcałej kohorty). Najpierw pobierz specyficzny dla taksonu wskaźnik wydajności dla każdej próbki, sprawdzając każdą próbkę na stronie projektu, klikając "Pobierz i wybierając "Combined Sample Taxon Results". Powstały plik .csv można przeanalizować za pomocą oprogramowania statystycznego takiego jak R lub Stata, w zależności od potrzeby.
UWAGA: Biorąc pod uwagę niższą względną częstość patogenów wirusowych, każdy wirus o znanym potencjale neuroinwazyjnym oraz przynajmniej jednym odczytie zgodnym z genomem wirusa powinien być uznany za wynik pozytywny (a najlepiej potwierdzony poprzez powtarzające się mNGS lub testy kliniczne specyficzne dla patogenów). Ponadto potok analityczny wymaga, aby każdy odczyt w pliku fastq był wyrównany obliczeniowo z najbardziej prawdopodobnym genomem referencyjnym w bazie danych NCBI. W rezultacie niektóre odczyty mogą być stosunkowo niespecyficzne dla konkretnego patogenu i teoretycznie mogły być powiązane z szerokim zakresem organizmów, a nie z konkretnym taksonem przypisanym do CZID. Biorąc pod uwagę to ryzyko, ręczne potwierdzenie pozytywnych wywołań powinno być przeprowadzone poprzez potwierdzenie wyrównania za pomocą internetowego narzędzia NCBI BLAST (konkretne odczyty taksonu można bezpośrednio przesłać z CZID do NCBI BLAST, klikając ikonę blast obok odpowiedniego taksonu).