Protokół sekwencjonowania ultra-długiego DNA wykorzystuje HMW DNA do konstrukcji biblioteki. Dlatego w kroku zbierania komórek krytyczne jest wybranie dobrze hodowanych komórek z odsetkiem żywych komórek >85%. Liczba komórek użytych do ekstrakcji DNA wpłynie na jakość i ilość HMW DNA. Liza komórek nie przebiega prawidłowo, jeśli proces rozpocznie się od zbyt dużej liczby komórek. Zastosowanie zbyt małej liczby komórek nie pozwala na uzyskanie wystarczającej ilości DNA do konstrukcji biblioteki, ponieważ wytrącanie HMW DNA odbywa się poprzez delikatne obracanie ręczne, a nie wirowanie przy wysokich prędkościach. Przykład HMW DNA po dodaniu lodowatego 100% etanolu i obracaniu przedstawiono jako biały, przypominający watę osad na Rysunku 2.
Przed rozpoczęciem konstrukcji biblioteki ważne jest sprawdzenie jakości wejściowego DNA. Degradacja, nieprawidłowe ilościowe oznaczenie, zanieczyszczenia (np. białka, RNA, detergenty, surfaktanty oraz pozostałości fenolu lub etanolu) oraz DNA o niskiej masie cząsteczkowej mogą mieć istotny wpływ na kolejne procedury i końcową długość odczytu. Zalecamy przeprowadzenie analizy kontroli jakości (QC) z wykorzystaniem DNA pobranego z trzech różnych miejsc w probówce zawierającej HMW DNA. Z wyników pomiarów UV dla HMW DNA wynika, że wartość OD260/OD280 wynosi około 1,9, a wartość OD260/OD230 wynosi około 2,3 (Rycina 3A,B). W przypadku dobrej próbki HMW DNA wartości tych stosunków są spójne we wszystkich trzech testach. Różne metody fragmentacji wymagają różnych objętości wejściowego DNA. Stężenie HMW DNA musi być >200 ng/µL w przypadku fragmentacji mechanicznej, natomiast >1 µg/µL w przypadku fragmentacji za pomocą transpozazy. Stężenie wykryte za pomocą fluorometru jest nieco niższe niż w pomiarze UV. Niemniej jednak współczynnik zmienności stężenia tej samej próbki HMW DNA powinien być mniejszy niż 15% zarówno w analizie fluorometrycznej, jak i w pomiarze UV. Fragmentacja mechaniczna wykorzystuje strzykawkę z igłą do rozrywania HMW DNA, zatem liczba przejść przez igłę wpływa na rozmiar pofragmentowanego DNA i końcową długość odczytu. Zaleca się przeprowadzenie kontroli jakości rozmiaru po fragmentacji igłą, aby upewnić się, że większość HMW DNA ma rozmiar większy niż 50 kb, co zilustrowano na Rycini 4. W metodzie fragmentacji mechanicznej 30 przejść dało najlepsze wyniki sekwencjonowania, biorąc pod uwagę zarówno długość, jak i wydajność.
Wartość N50 biblioteki przygotowanej metodą mechanicznego ścinania wynosi 50-70 kb, natomiast biblioteki przygotowanej metodą fragmentacji za pomocą transpozazy wynosi 90-100 kb. Wyniki czterech przebiegów z wykorzystaniem linii komórkowej HG00733 przedstawiono w Tabeli 1. We wszystkich czterech przebiegach uzyskano ponad 2300 odczytów o długości przekraczającej 100 kb. Maksymalna długość jest większa w bibliotekach przygotowanych metodą fragmentacji za pomocą transpozazy (455 kb i 489 kb) w porównaniu z bibliotekami przygotowanymi metodą mechanicznego ścinania (348 kb i 387 kb), podczas gdy te ostatnie wygenerowały więcej odczytów całkowitych, co wskazuje na wyższą wydajność. Konstrukcja biblioteki metodą fragmentacji za pomocą transpozazy obejmuje mniej etapów i wymaga krótszego czasu przygotowania, dzięki czemu powstaje mniej krótkich fragmentów. Dwa przebiegi z wykorzystaniem transpozazy charakteryzowały się dłuższą średnią długością (>30 kb) i medianą długości (>10 kb). Ponadto dane wykazują spójną, wysoką jakość we wszystkich przebiegach (średni wynik jakości wynosi około 10,0, co odpowiada ok. 90% dokładności zasad). Ponad 97% wszystkich zasad zostało dopasowanych do ludzkiego genomu referencyjnego (hg19) przy użyciu programu Minimap216 z ustawieniami domyślnymi. Oczekiwane rozkłady wielkości surowych odczytów przedstawiono na Rysunku 5. We wszystkich przebiegach duża część danych przekracza 50 kb, natomiast biblioteki przygotowane metodą fragmentacji za pomocą transpozazy mają wyższy stosunek odczytów ultradługich (np. > 100 kb). Niniejszy protokół został z powodzeniem zastosowany w wielu ludzkich liniach komórkowych (Tabela uzupełniająca 1).

Rycina 1: Schematyczny przegląd przepływu pracy sekwencjonowania długich odczytów nanoporowych (NLR-seq). Kolor pomarańczowy: kompleks transpozazy. Kolor żółto-zielony: adapter nanoporowy. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Reprezentatywne wytrącanie DNA w metodzie ekstrakcji fenol-chloroformowej. Biała strzałka wskazuje HMW DNA. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Przykład wyników kontroli jakości HMW DNA na podstawie pomiarów UV. (A) HMW DNA z etapu 1.21.1 gotowe do konstrukcji biblioteki opartej na mechanicznym ścinaniu. (B) HMW DNA z etapu 1.21.2 do konstrukcji biblioteki opartej na fragmentacji za pomocą transpozazy. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Wyniki kontroli jakości (QC) HMW DNA poddanego ścinaniu igłowym za pomocą elektroforezy w polu pulsującym. L1: drabinka DNA Quick-load 1 kb; L2: rozszerzona drabinka DNA Quick-load 1 kb. 1-8: DNA o różnych czasach przejścia przez igłę podczas ścinania. 1-3, brak ścinania; 4, 10 razy; 5, 20 razy; 6, 30 razy; 7, 40 razy; 8, 50 razy. Ten etap kontroli jakości jest opcjonalny. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Oczekiwane rozkłady rozmiarów bibliotek ultra-długiego DNA do sekwencjonowania nanoporowego. MS, biblioteki otrzymane w drodze mechanicznego ścinania. TF, biblioteki otrzymane poprzez fragmentację za pomocą transpozazy. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
| Mechaniczne ścinanie_rep1 | Mechaniczne ścinanie_rep2 | Fragmentacja transpozazą_rep1 | Fragmentacja transpozazą_rep2 |
| Linia komórkowa | HG00733 | HG00733 | HG00733 | HG00733 |
| N50 odczytów | 55,180 | 63,007 | 98,237 | 95,629 |
| Liczba odczytów dłuższych niż 100 Kb | 2,500 | 3,082 | 2,386 | 2,355 |
| Całkowita liczba odczytów | 97,859 | 80,465 | 24,166 | 21,032 |
| Maksymalna długość (bp) | 348,482 | 387,113 | 454,660 | 489,426 |
| Średnia długość (bp) | 17,861 | 20,395 | 33,528 | 38,175 |
| Mediana długości (bp) | 5,335 | 5,894 | 10,249 | 15,656 |
| Średnia jakość odczytów | 10.0 | 10.1 | 9.9 | 10.0 |
| Całkowita liczba zasad w surowych odczytach | 1,747,849,822 | 1,641,058,932 | 810,229,733 | 802,886,304 |
| Całkowita liczba zasad w dopasowanych odczytach | 1,693,300,832 | 1,607,975,925 | 791,422,077 | 778,417,627 |
| Stosunek mapowania całkowitej liczby zasad (hg19, Minimap2) | 96.9% | 98.0% | 97.7% | 97.0% |
| Liczba aktywnych porów | 1225: 480, 402, 254, 89 | 1058: 480, 356, 176, 46 | 958: 452, 328, 148, 30 | 1092: 487, 367, 195, 43 |
Tabela 1: Podsumowanie wskaźników wydajności z serii pomiarowych z zastosowaniem różnych protokołów ścinania.
| Biblioteka 1 | Biblioteka 2 |
| Linia komórkowa | K562 | GM19240 |
| Informacje dotyczące zamawiania komórek | ATCC, nr kat. CCL-243 | Coriell Institute, nr kat. GM19240 |
| Protokół | ścinanie mechaniczne | ścinanie mechaniczne |
| Wartość N50 odczytów | 60,063 | 55,295 |
| Całkowita liczba odczytów | 193,783 | 120,807 |
| Długość mediana (bp) | 1,843 | 4,688 |
| Średnia długość (bp) | 9,825 | 17,408 |
| Maksymalna długość (bp) | 548,780 | 212,338 |
| Całkowita liczba zasad surowych odczytów | 1,903,989,686 | 2,103,015,331 |
| Całkowita liczba zasad zmapowanych odczytów | 1,837,350,047 | 1,997,419,761 |
| Zmapowany stosunek całkowitej liczby zasad (hg19, Minimap2) | 96.6% | 95.0% |
| Liczba aktywnych porów | 1111: 482, 371, 203, 55 | 1032: 447, 333, 196, 56 |
Tabela uzupełniająca 1: Podsumowanie dwóch przebiegów NLR-seq z wykorzystaniem innych linii komórkowych przy zastosowaniu protokołu ścinania mechanicznego.