Aby zwalidować metodę Circle-Seq, przebadano trzy populacje S. cerevisiae CEN.PK zawierające 1 x 1010 komórek, po ich oddzielnej hodowli w pożywce YPD przez dziesięć pokoleń. Eliminację liniowego DNA chromosomowego potwierdzono brakiem sygnału qPCR dla ACT1 , zgodnie z wcześniejszym opisem20 (dane nie zostały przedstawione). Oczyszczone i wzbogacone eccDNA zsekwencjonowano, uzyskując do 68 milionów odczytów (odczyty typu single-end o długości 141 nukleotydów), a następnie zmapowano do genomu referencyjnego CEN.PK113-7D (wersja z 19 czerwca 2012 r.). Odczyty prawdopodobnych eccDNA z trzech próbek oznaczonych jako C1, C2 i C4 przypisano do regionów genomowych zmapowanych za pomocą ciągłych odczytów o długości przekraczającej 1 kb. Na podstawie 10 000 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 zinterpretowano jako prawdopodobne sekwencje eccDNA (p < 0,1, Dataset 1). Liczba zarejestrowanych ciągłych odczytów > 1 kb wzrastała wraz z głębokością sekwencjonowania, co sugeruje, że przy dalszym sekwencjonowaniu próbek zarejestrowano by jeszcze więcej elementów eccDNA (Figure 2). Zgodnie z oczekiwaniami, metoda Circle-Seq pozwoliła wyekstrahować liczne odczyty 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 pBR322, pUC19 i pUG72, które dodano do próbek bezpośrednio przed oczyszczaniem kolumnowym (Figure 3).
Film przedstawia przykład odczytów ciągłych zmapowanych do locus HXT7_ARS432_HXT6 na chromosomie IV. Wcześniej [HXT6/7circle] został wykryty metodą Circle-Seq w dziesięciu populacjach S288c (każda zawierająca 1 x 1010 komórek), a kolista struktura DNA została potwierdzona analizą inverse PCR20. [HXT6/7circle] odnotowano również w każdej z trzech populacji CEN.PK (Rycyna 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 S288c (Rycina 4B).
Aby przetestować swoistość protokołu Circle-Seq w zakresie oczyszczania kolistego DNA, przygotowano dwie próbki, z których każda zawierała 30 µg badano genomowe DNA. Jedną próbkę uzupełniono o 100 ng DNA plazmidowego i eccDNA z obu Próbki oczyszczono zgodnie z protokołem Circle-Seq. Po rozdziału kolumnowym wydajność uzyskanego 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ść trawienia egzonukleazą w odniesieniu do zawartości liniowego DNA sprawdzono po 29 godzinach. oraz 72 h wykorzystując PCR przeciwko AKT1. W żadnych próbkach nie stwierdzono amplifikacji ACT1 (dane nie przedstawiono). Część każdej próbki poddanej działaniu egzonukleazy została następnie poddana amplifikacji za pomocą ø29 polimerazę oraz produkty reakcji enzymatycznych przeanalizowano za pomocą barwienia jodkiem propidyny (Rycina 5A-F) oraz elektroforezę w żelu agarozowym (Rycina 5G). Próbki po traktowaniu egzonukleazą wykazały minimalne barwienie jodkiem propidyny (Rycina 5A-B). The ø29-próbka z amplifikacją zawierająca wyłącznie DNA genomowe wykazała obecność struktur nitkowatych (Rysunek 5C) podobnie jak w próbce kontrolnej (Rycina 5E). The ø29-w amplifikowanej próbce z dodatkiem plazmidu ujawniono ogniska (Rysunek 5D) przypominający kontrolę plazmidową (Rycina 5F). Obrazy wskazały, że øpolimeraza wzbogacona o DNA koliste w stosunku do DNA liniowego. Większość liniowego DNA chromosomalnego została usunięta z próbek po 29-godzinnej inkubacji z egzonukleazą (Rysunek 5A-B, G). Jednakże przedłużona inkubacja z egzonukleazą przez ponad 100 h i zastosowanie więcej niż 100 jednostek było konieczne w celu usunięcia całego liniowego DNA chromosomalnego, ponieważ ø29-wzmocniony próbki nadal wykazywały tło w postaci struktur przypominających nici po 72-godzinnej inkubacji z egzonukleazą (Rycina 5C-D).

Rycina 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 chromosomowego w frakcji eluated, 4) amplifikacji DNA za pomocą polimerazy DNA ø29 oraz 5) sekwencjonowania silnie wzbogaconego eccDNA i mapowania odczytów do genomu referencyjnego S. cerevisiae. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2. Ciągłe odczyty > 1 kb w funkcji głębokości sekwencjonowania. Ilość EccDNA z 1 x 1010 komórek rośnie w funkcji głębokości sekwencjonowania (w milionach zmapowanych odczytów). Przedstawiono: powtórzenia biologiczne z haploidalnych populacji S. cerevisiae CEN.PK (C1, C2, C4) oddzielone 1010 podziałami komórkowymi. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 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 biologicznych powtórzeniach C1, C2 i C4 szczepu CEN.PK. (A) Odczyty zmapowane do endogennych plazmidów drożdżowych to: 2µ; [rDNAcircle] (geny rRNA z chromosomu XII) oraz mtDNA (mitochondrialne DNA). (B) Unikalne odczyty zmapowane do plazmidów kontrolnych. Plazmidy kontrolne zostały dodane do próbek przed oczyszczaniem kolumnowym. Stosunki plazmidów na komórkę wynosiły: pBR322 (plusy) 1:1, pUC19 (koła) 1:50 oraz pUG72 (trójkąty) 1:2 500.

Rycina 4. Wspólne elementy eccDNA w szczepach CEN.PK i S288c. (A) Diagram Venna przedstawiający nakładanie się 476 geny na 294 elementach eccDNA w trzech próbkach CEN.PK (C1, C2, C4). Naniesiono adnotacje dla 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 zarejestrowanych genów na domniemanych eccDNA z trzech próbek CEN.PK (C1, C2, C4) w porównaniu do wszystkich zarejestrowanych genów na domniemanych eccDNA z 10 próbek S288c: 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 pokrywały minimum dwa tła szczepowe oraz co najmniej trzy konfiguracje eksperymentalne. Próbki C: CEN.PK; próbki R i Z: S288c BY4741; próbki S: S288c M3750. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5. Wizualizacja próbek DNA po zastosowaniu egzonukleazy i ø29. traktowanie. (A-F) Barwienie DNA jodkiem propidyny. Pasek skali, 10 µm. (A, C i E) Próbki z DNA genomicznym (GD); (B i D) próbki z plazmidem GD plus (GD+P). (A-B) po 29-godzinnej inkubacji z egzonukleazą (EXO 29 h); (C-D) po 72-godzinnej inkubacji z egzonukleazą, a następnie ø29 amplifikacja polimerazowa (EXO 72 h + ø29). (E) Kontrola genomowego DNA po e: ø29 amplifikacja polimerazowa; (F) kontrola plazmidowa (5,5 kb) po 29 amplifikacja polimerazowa; (Gø) elektroforeza w żelu agarozowym. Od lewej: L, 1 kb markery; P, kontrola plazmidowa (5,5 kb) po 29 h działania EXO; GD, po 29 h działania EXO (próbka jak w A); GD+P, po 29 godzinach EXO (próbka jako B); GD i GD+P, po 29 godzinach EXO + ø29; GD i GD+P, po 72 godzinach od EXO + ø29 (próbka jak w C-D). Patrz tabela S1 w celu uzyskania dodatkowych szczegółów. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

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