Do badania edycji genów w komórkach ssaków wykorzystaliśmy model oparty na korekcji mutacji punktowej w obrębie genu eGFP, zintegrowanego jako pojedyncza kopia w komórkach HCT116. Należy zaznaczyć, że jest to gen pojedynczej kopii, co pozwala na mniej złożoną analizę zmian w DNA lub mutagenezy. Rysunek 1A przedstawia sekwencję zmutowanego genu eGFP z zaznaczoną na czerwono docelową zasadą, którą jest trzecia zasada kodonu stop TAG. Przedstawiono również 72-zasadowy oligonukleotyd (72NT), który jest częściowo komplementarny do nici nietranskrybowanej genu eGFP i został zaprojektowany w celu wywołania zamiany zasady G na C. Ponadto na Rysunku 1A ukazano CRISPR, oznaczony jako 2C, z wskazaną sekwencją protospacera w orientacji 5' do 3'. Aby przeprowadzić tę reakcję edycji genów, zamiast wektora ekspresyjnego dla ssaków zawierającego gen Cas9 i odpowiednią specyficzną sekwencję przewodniego RNA, zastosowaliśmy rybonukleoproteinę (RNP) składającą się z CRISPR (cr) RNA oraz tracr (tr) RNA połączonych z oczyszczonym białkiem Cas9 (Rysunek 1B).
Kiedy specjalnie zaprojektowana cząstka RNP zostanie dostarczona wraz z jednoniciowym oligonukleotydem do komórek HCT116 poprzez elektroporację, edycję genów, potwierdzoną naprawą pojedynczej mutacji zasadowej w eGFP, obserwuje się po 72 h inkubacji przy użyciu cytometru przepływowego. Funkcjonalną naprawę można zaobserwować poprzez pojawienie się zielonej fluorescencji w komórkach docelowych, które dzięki tej fluorescencji mogą zostać oddzielone od całej populacji poprzez sortowanie. Jak pokazano na Rysunku 2, widoczna jest stopniowa odpowiedź na dawkę wraz ze wzrostem skoordynowanych poziomów RNP i 72NT. Stosunek molowy, pikomole RNP oraz stężenie mikromolarne 72NT przedstawione na rysunku opierają się na optymalnych dawkach stosowanych w reakcjach edycji genów zależnych od wprowadzenia Cas9 i specyficznego przewodzącego RNA z transfekcyjnych wektorów ekspresyjnych. Wstawka na Rysunku 2 przedstawia reakcję edycji genów przeprowadzoną w obecności cząstki RNP. W tym przypadku około 1% komórek docelowych zostaje skorygowanych, gdy w jednoskładnikowej reakcji edycji genów zastosowano 10-krotnie wyższe stężenie oligonukleotydu 72NT.
Aby ustalić, czy w miejscu docelowym występuje heterogeniczność genetyczna — tzw. efekt mutagenezy on-site — postanowiliśmy zbadać wyniki aktywności edycji genów w poszczególnych komórkach. Choć jest to znacznie bardziej pracochłonne niż badanie całej populacji, rzeczywistą miarę śladów genetycznych lub uszkodzeń można ustalić poprzez analizę genomu komórek poddanych ekspansji klonalnej. Powtórzyliśmy eksperyment opisany na Rysunku 2, tym razem stosując jedynie 100 pmol kompleksu RNP oraz 2,0 µmol oligonukleotydu 72NT. Podobnie jak wcześniej, komórki HCT 116 synchronizowano przez 24 h za pomocą afidicholiny, zatrzymując je na granicy faz G1/S. Po 4 h uwolniono komórki i wprowadzono narzędzia do edycji genów metodą elektroporacji. Po 72 h komórki przeanalizowano za pomocą FACS i posortowano pojedynczo do płytek 96-dołkowych (proces eksperymentalny przedstawiono na Rysunku 3). Komórki wykazujące zieloną fluorescencję posortowano za pomocą cytometrii przepływowej do poszczególnych dołków płytki 96-dołkowej w celu ekspansji klonalnej. Co ważne, wyizolowano i posortowano w podobny sposób również komórki pozbawione ekspresji eGFP, aby poddać je ekspansji w tych samych warunkach.
Po 14 dniach wzrostu większość poszczególnych klonów uległa wystarczającemu namnożeniu, aby umożliwić izolację DNA. W związku z tym wybrano 16 klonów z próbek eGFP-dodatnich i przeanalizowano integralność genetyczną w obrębie miejsca docelowego za pomocą sekwencjonowania DNA. Informacje dotyczące sekwencji DNA alleli w populacji uzyskano za pomocą sekwencjonowania metodą Sangera, a następnie zestawiono przy użyciu oprogramowania do wizualizacji sekwencji, aby porównać sekwencję z allelem typu dzikiego (Rycina 4A). Miejsce cięcia kompleksu RNP zaznaczono czarną strzałką, znajdującą się na zielonym pasku (2C crRNA). Jak również pokazano na Rycynie 4A, wszystkie 16 komórek eGFP-dodatnich zawiera przewidywaną zamianę nukleotydów w miejscu docelowym. Przekształcony reszty C wyróżniono na czerwono, a profil piku odzwierciedlający tę precyzyjną zmianę przedstawiono pod sekwencją eGFP-dodatnią. W podobny sposób 15 niezielonych izolatów klonalnych, które namnożyły się wystarczająco, aby umożliwić ekstrakcję DNA i sekwencjonowanie, przeanalizowano pod kątem heterogeniczności w miejscu docelowym. Zgodnie z przewidywaniami, w około połowie próbek nie zaobserwowano zamiany zasad DNA. Odzwierciedla to utrzymanie reszty G w miejscu docelowym, co pokazano na Rycynie 4B. Pozostałe badane w tych eksperymentach namnożone klony wykazały heterogenną populację mutacji delecyjnych, co wyjaśnia brak zielonej fluorescencji. Wielkość delecji wynosiła od jednej do 19 zasad. Należy zaznaczyć, że zbadano jedynie 15 próbek komórek eGFP-dodatnich i choć uważamy, że jest to reprezentatywne dla rodzaju uszkodzeń genetycznych pozostawianych przez aktywność CRISPR/Cas9, istnieje możliwość, że w populacji docelowej mogą występować inne typy lub formy indeli.
Wzięte razem, wyniki przedstawione na Rysunku 4 potwierdzają odczyt fenotypowy w systemie targetowania eGFP. Konwersja nukleotydu G na C umożliwia pojawienie się zielonej fluorescencji w skorygowanych komórkach HCT116. Nie zaobserwowano żadnej substytucji zasad w sąsiedztwie miejsca docelowego w klonach wyizolowanych do tego eksperymentu ani w poprzednich badaniach27,28. Dane wykazują również, że komórki, w których nie doszło do edycji genów poprzez naprawę mutacji punktowej, pozostają nieskorygowane, ale w niektórych przypadkach nie pozostają niezmienione, wykazując zakres heterogeniczności genetycznej w sąsiedztwie miejsca docelowego.

Rysunek 1. (A) Model systemu edycji genu zmutowanego eGFP. Odpowiednie segmenty genu eGFP typu dzikiego oraz zmutowanego z docelowym kodonem, zlokalizowanym w centrum sekwencji, przedstawiono odpowiednio w kolorze zielonym i czerwonym. Nukleotyd przeznaczony do wymiany został pogrubiony i podkreślony. Bazy wyróżnione na niebiesko reprezentują sekwencję protospacera 2C CRISPR, a bazy pomarańczowe wyróżniają miejsce PAM. Oligonukleotyd użyty w tych eksperymentach ma długość 72 zasad i posiada modyfikacje fosforotiolanowe w trzech zasadach końcowych; 72-mer targetuje nietranskrybowaną nić (NT) (72NT). (B) Reakcja montażu rybonukleoproteinowego CRISPR/Cas9. crRNA zapewnia specyficzność celu (20 zasad, sekcja czerwona) odpowiadającą sekwencji protospacera 2C oraz domenę oddziaływania (niebieska) z tracrRNA (zielona). crRNA i tracrRNA są hybrydyzowane w stężeniach ekvimolarnych. Białko Cas9 (szare) jest dodawane w celu dopełnienia kompleksu RNP. Przewodniki RNA (gRNA) kierują i aktywują endonukleazę Cas9, która następnie przecina docelowe DNA. Dolna sekcja rysunku przedstawia sekwencję seed 2C oraz sekwencję tracrRNA. Rysunek zmodyfikowano na podstawie Rivera-Torres, N. et al. (2017). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2. Edycja genów jest zależna od dawki przy zastosowaniu RNP i ssODN. Zsynchronizowane i uwolnione komórki HCT 116-19 poddano elektroporacji z użyciem 24-120 pmol RNP CRISPR/Cas9 oraz 0,6-3,0 µM 72meru. Po 72-godzinnym okresie regeneracji aktywność edycji genów zmierzono za pomocą cytometru przepływowego. Edycję genów przedstawiono jako wydajność korekcji (%), wyznaczoną przez liczbę żywych komórek eGFP-dodatnich podzieloną przez całkowitą liczbę żywych komórek w populacji. Słupki błędów zostały wyznaczone z trzech zestawów punktów danych uzyskanych w trzech oddzielnych eksperymentach na podstawie obliczeń błędu standardowego. Wstawka: Edycja genów przy użyciu pojedynczego czynnika. Aktywność edycji genów kierowana przez jednoniciowy oligonukleotyd (72NT) w obecności kompleksu RNP w identycznych warunkach przedstawiono jako funkcję rosnącego stężenia. Rysunek zmodyfikowano na podstawie pracy Rivera-Torres, N. et al. (2017). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3. Strategia eksperymentalna izolacji klonów pojedynczych komórek. Komórki wykazujące ekspresję eGFP zostały ocenione jako pozytywne i posortowane za pomocą cytometru przepływowego jako pojedyncze komórki do poszczególnych studni w celu ekspansji klonalnej. Komórki pozbawione ekspresji eGFP zostały wyizolowane i posortowane w podobny sposób i poddane ekspansji w tych samych warunkach. Następnie wyizolowano DNA, a gen eGFP został powielony i poddany sekwencjonowaniu metodą Sangera w celu analizy aktywności edycji genu w obrębie miejsca docelowego. Rysunek zmodyfikowano na podstawie pracy Rivera-Torres, N. et al. (2017). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4. (A) Analiza alleliczna komórek eGFP-dodatnich powiększonych jako populacja klonalna. Próbki eGFP-dodatnie, wyizolowane klonalnie i powiększone (szesnaście klonów), przeanalizowano w obrębie miejsca wokół docelowej zasady; DNA z każdej z nich pobrano, oczyszczono, namnożono i zsekwencjonowano. Analizę alleliczną przeprowadzono za pomocą sekwencjonowania metodą Sangera, a wyniki złożono przy użyciu oprogramowania do wizualizacji sekwencji i porównano z sekwencją allelu typu dzikiego, przedstawioną w górnej części rysunku. Miejsce cięcia kompleksu RNP zaznaczono małą czarną strzałką znajdującą się na zielonym pasku (2C crRNA). (B) Analiza alleliczna komórek eGFP-ujemnych powiększonych jako populacja klonalna. Losowo wybrano piętnaście pojedynczych próbek, powiększonych z klonów pochodzących z populacji niepoprawionej, i przeanalizowano je pod kątem powstawania indeli w miejscu wokół nukleotydu docelowego. Jak wyżej, analizę alleliczną przeprowadzono za pomocą sekwencjonowania metodą Sangera i złożono przy użyciu oprogramowania do wizualizacji sekwencji. Ponownie przedstawiono sekwencję allelu typu dzikiego w górnej części rysunku wraz z miejscem cięcia RNP. Rysunek zmodyfikowano na podstawie pracy Rivera-Torres, N. et al. (2017). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.