Aby zwalidować metodę Circle-Seq, przesiewowo zbadano trzy populacje S. cerevisiae CEN.PK o liczebności 1 x 1010 komórek, po ich oddzielnej hodowli w pożywce YPD przez dziesięć generacji. Eliminację liniowego DNA chromosomowego potwierdzono brakiem sygnału qPCR dla genu ACT1 , zgodnie z wcześniejszym opisem20 (dane niepokazane). Oczyszczone i wzbogacone eccDNA poddano sekwencjonowaniu, uzyskując do 68 milionów odczytów (odczyty jednostronne o długości 141 nukleotydów), a następnie zmapowano je do genomu referencyjnego CEN.PK113-7D (wersja z 19 czerwca 2012 r.). Odczyty prawdopodobnych eccDNA z trzech próbek nazwanych C1, C2 i C4 przypisano do regionów genomowych zmapowanych przez ciągłe odczyty dłuższe niż 1 kb. Na podstawie 10 00 symulacji Monte Carlo oszacowano istotność każdego regionu zmapowanego przez ciągłe odczyty dłuższe niż 1 kb. W ten sposób 79, 159 i 56 regionów adnotowano jako prawdopodobne sekwencje eccDNA (p < 0.1, Dataset 1). Liczba zarejestrowanych ciągłych odczytów > 1 kb rosła w funkcji głębokości sekwencjonowania, co sugeruje, że przy głębszym sekwencjonowaniu próbek wykryto by jeszcze więcej elementów eccDNA (Figure 2). Zgodnie z oczekiwaniami, metoda Circle-Seq pozwoliła na wyodrębnienie licznych odczytów z wielu znanych kolistych elementów DNA, w tym z plazmidu 2µ, DNA mitochondrialnego, genów rRNA na chromosomie XII oraz trzech wewnętrznych plazmidów kontrolnych pBR32, pUC19 i pUG72, które dodano do próbek bezpośrednio przed oczyszczaniem kolumnowym (Figure 3).
Wideo przedstawia przykład sąsiednich odczytów zmapowanych do locus HXT7_ARS432_HXT6 na chromosomie IV. Wcześniej [HXT6/7circle] został wykryty za pomocą metody Circle-Seq w dziesięciu populacjach S28c (każda zawierająca 1 x 1010 komórek), a kolista struktura DNA została potwierdzona analizą odwrotnej PCR20. [HXT6/7circle] odnotowano również w każdej z trzech populacji CEN.PK (Rysunek 4A). Ponadto większość wspólnych genów eccDNA w próbkach powtórzonych CEN.PK pokrywała się z genami eccDNA z zestawów danych S28c (Rysunek 4B).
Aby przetestować specyficzność protokołu Circle-Seq w zakresie oczyszczania kolistego DNA, przygotowano dwie próbki, z których każda zawierała 30 µg przebadano genomowe DNA. Jedną z próbek uzupełniono o 10 ng plazmidowego DNA, a następnie z obu próbek oczyszczono eccDNA zgodnie z protokołem Circle-Seq. Po rozdziale kolumnowym wydajność DNA wyniosła 1,27% (380 ng) dla próbki bez plazmidu (GD) oraz 1,60% (480 ng) dla próbki z plazmidem (GD+P). Efektywność traktowania egzonukleazą w odniesieniu do zawartości liniowego DNA przetestowano po 29 h oraz 72 h za pomocą reakcji PCR przeciwko ACT1. W żadnych próbkach nie stwierdzono amplifikacji ACT1 (dane nie pokazano). Część każdej próbki traktowanej egzonukleazą poddano dalszej amplifikacji za pomocą ø29 polimerazę oraz produkty reakcji enzymatycznych przeanalizowano za pomocą barwienia jodkiem propidyny (Rysunek 5A-F) oraz elektroforezę w żelu agarozowym (Rysunek 5G). Próbki po traktowaniu egzonukleazą wykazały minimalne barwienie jodkiem propidyny (Rysunek 5A-B). ø29-próbka wzbogacona, zawierająca wyłącznie DNA genomowe, ujawniła struktury przypominające nici (Rycina 5C) podobne do próbki kontrolnej (Rycina 5E). The ø29-wz wzmocnionej próbce z dodatkiem plazmidu zaobserwowano ogniska (Rycina 5D) przypominający kontrolę plazmidową (Rycina 5F). Obrazy wskazały, że polimeraza ø29 wzbogaca frakcję kolistego DNA w stosunku do liniowego DNA. Większość liniowego DNA chromosomalnego została usunięta z próbek po 29-godzinnej inkubacji z egzonukleazą (Rysunek 5A-B, G). Jednakże, aby usunąć całe chromosomalne liniowe DNA, wymagana była intensywna obróbka egzonukleazą przez ponad 100 h z użyciem więcej niż 100 jednostek, ponieważ ø29-próbki po amplifikacji nadal wykazywały obecność w tle struktur przypominających nici po 72-godzinnej inkubacji z egzonukleazą (Rycina 5C-D).

Rysunek 1. Schemat metody Circle-Seq. Protokół składa się z 5 etapów: 1) hodowli komórek, 2) oczyszczania i wzbogacania eccDNA za pomocą chromatografii kolumnowej, 3) trawienia pozostałego liniowego DNA chromosomalnego w frakcji eluate, 4) amplifikacji DNA za pomocą polimerazy DNA ø29 oraz 5) sekwencjonowania silnie wzbogaconego eccDNA i mapowania odczytów do S. cerevisiae genom referencyjny Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2. Odczyty kontigowe > 1 kb w funkcji głębokości sekwencjonowania. eccDNA z 1 x 1010 liczba komórek wzrasta w funkcji głębokości sekwencjonowania (w milionach zmapowanych odczytów). Przedstawiono: triplikaty biologiczne z haploidalnego szczepu CEN.PK S. cerevisiae populacje (C1, C2, C4) oddzielone przez 1010 podziały komórkowe. Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 3. Detekcja znanych kolistych elementów DNA. (A-B) Wykresy rozrzutu pokrycia odczytów (gęstości odczytów) w procentach dla plazmidów w powtórzeniach biologicznych CEN.PK: C1, C2 i C4. (A) Liczba odczytów zmapowanych do endogennych plazmidów drożdżowych wynosiła: 2µ; [rDNAokrąg] (geny rRNA z chromosomu XII) oraz mtDNA (mitochondrialne DNA). (B) Unikalne odczyty zmapowane do plazmidów kontrolnych. Plazmidy kontrolne dodano do próbek przed oczyszczaniem kolumnowym. Stosunki plazmidów na komórkę wynosiły: pBR32 (plusy) 1:1, pUC19 (kółka) 1:50 oraz pUG72 (trójkąty) 1:2 50.

Rysunek 4. Wspólne elementy eccDNA w szczepach CEN.PK i S28c. (A) Diagram Venna przedstawiający nakładanie się 476 geny na 294 elementach eccDNA w trzech próbkach CEN.PK (C1, C2, C4). Adnotowano 16 wspólnych, nakładających się genów/plazmidów eccDNA (wszystkie nazwy genów znajdują się w Zbiór danych 1). (B) Diagram Venna wszystkich genów zidentyfikowanych w domniemanych ekstrachromosomalnych kołowych cząsteczkach DNA (eccDNA) z trzech próbek CEN.PK (C1, C2, C4), w porównaniu do wszystkich genów zidentyfikowanych w domniemanych eccDNA z 10 próbek S28c: S1-S2, R1-R4, Z1-Z4 (patrz odniesienie20). Przedstawiono 13 powtórzeń biologicznych (S1-S2, R1-R4, Z1-Z4, C1-C3) z genami/plazmidami oraz domniemanymi regionami eccDNA, które nakładały się na co najmniej 2 tła szczepowe oraz na 3 lub więcej układów doświadczalnych. Próbki C: CEN.PK; próbki R i Z: S28c BY4741; próbki S: S28c M3750. Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ilustracji.

Rysunek 5. Wizualizacja próbek DNA po zastosowaniu egzonukleazy i ø29 traktowanie. (A-F) Barwienie DNA jodkiem propidium. Pasek skali, 10 µm. (A, C i E) Próbki z kwasem DNA genomowym (GD); (B i D) próbki z GD oraz plazmidem (GD+P). (A-B) Po 29-godzinnej inkubacji z egzonukleazą (EXO 29 h); (C-D) po 72-godzinnej inkubacji z egzonukleazą, a następnie amplifikacji za pomocą polimerazy ø29 (EXO 72 h + ø29). (E) kontrola genomowego DNA po amplifikacji polimerazą ø29; (F) kontrola plazmidowa (5,5 kb) po 29 amplifikacja polimerazowa; (Gø) elektroforeza w żelu agarozowym. Od lewej: L, markery 1 kb; P, kontrola plazmidowa (5,5 kb) po 29 godz. działania EXO; GD, po 29 godz. działania EXO (próbka tak jak w A); GD+P, po 29 godzinach EXO (próbka jako B); GD i GD+P, po EXO 29 h + ø29; GD i GD+P, po EXO 72 h + ø29 (próbka jak w C-D). Patrz Tabela S1 po dodatkowe szczegóły. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Zbiór danych 1. Potencjalne regiony cyrkularyzacji DNA w CEN.PK. Prosimy kliknąć tutaj, aby pobrać ten plik.
Przedstawiono dane sekwencyjne i analizy dla 348 regionów. Kolumny A-D obejmują mapowanie eccDNA. A (pierwsza kolumna od lewej) to próbka, w której zidentyfikowano domniemane eccDNA; B to chromosom; C-D to współrzędne początkowe i końcowe domniemanych eccDNA. Kolumny E-H określają zawartość eccDNA. E to autonomicznie replikująca się sekwencja (ARS) w regionie; F to kompletny gen w regionie; G to fragment genu zawarty w regionie; H to gen zidentyfikowany za pomocą BLASTN. Kolumny I-O obejmują pokrycie eccDNA i wartości p. I to najdłuższy region z unikalnie adnotowaną sekwencją w bp; J to liczba wszystkich zmapowanych odczytów; K to pokrycie wszystkich zmapowanych odczytów wyrażone jako liczba fragmentów na kb z miliona zmapowanych odczytów (FPKM); L to wartość p dla domniemanych eccDNA w porównaniu do występowania przypadkowego w symulacjach Monte Carlo; M to liczba unikalnie zmapowanych odczytów; N-O to odpowiednio K i L z wykorzystaniem wyłącznie unikalnie zmapowanych odczytów (UFPKM). Parametry mapowania odczytów i symulacji Monte Carlo były zgodne z opisem20.