Procedura eksperymentalna SiC-seq obejmuje 3 mikroprzepływowe urządzenia PDMS wykonane metodą miękkiej litografii (Rysunek 1). Generator kropel współpłynnych (Rysunek 3A) wytwarza krople cyfrowych kodów kreskowych o rozmiarze 25 µm służące do znakowania DNA genomowego unikalnym identyfikatorem pojedynczej komórki. Oligonukleotydy kodów kreskowych składają się z 15 bp sekwencji zdegenerowanej, flankowanej przez uchwyty PCR do amplifikacji (Tabela 2, starter BAR). Kody kreskowe rozcieńcza się do stężenia femtomolarnego, aby uzyskać enkapsulację pojedynczych cząsteczek, tak aby każda kropla zawierała 0 lub 1 fragment kodu. Krople zawierające kod są amplifikowane, co daje wiele kopii dwuniciowych amplikonów kodów kreskowych. Barwnik kwasów nukleinowych jest wykorzystywany do weryfikacji pomyślnej amplifikacji oraz ilościowego określenia stopnia enkapsulacji fragmentów kodów kreskowych (Rysunek 4B). Mikrzele są generowane poprzez współpłynięcie zawiesiny komórek bakteryjnych i stopionego żelu agarozy przy równych prędkościach przepływu (Rysunek 2A). Agarozę przygotowuje się w dwukrotności pożądanej stężenia końcowego, ponieważ proces generowania kropel współpłynnych skutecznie rozcieńcza roztwory wodne dwukrotnie. W miarę stygnięcia agaroza krzepnie, tworząc mikrzele o średnicy 25 µm zajmujące sferyczną objętość kropli.
Seria etapów przemywania i lizy służy do oczyszczania wysokocząsteczkowego DNA genomowego w mikrogelach (Rysunek 2B). Po rozbiciu emulsji przeprowadza się przemywania wodne w dużych objętościach, aby rozcieńczyć śladowe ilości rozpuszczalników organicznych, które mogą hamować późniejsze procesy enzymatyczne. Przemyte mikrogela obserwuje się pod mikroskopem w celu zweryfikowania stopnia enkapsulacji komórek (Rysunek 4A). Do zawiesiny mikrogeli dodaje się koktajl enzymów o szerokim spektrum aktywności litycznej w celu strawienia ścian komórkowych bakterii i mikroorganizmów eukariotycznych19. Druga obróbka z użyciem Proteinazy K i detergentu degraduje białka i solubilizuje szczątki komórkowe.
Tagmentację oczyszczonego DNA przeprowadza się w kropelkach, aby uniknąć potencjalnej kontaminacji krzyżowej wynikającej z dyfuzji niewielkich fragmentów tagmentowanego DNA między mikrogelami18. Urządzenie do enkapsulacji w kropelkach (Ryc. 3B) wydziela każdy mikrogel wraz z buforem i enzymem do tagmentacji, który jednocześnie fragmentuje dwuniciowe DNA i „oznacza” je załadowanym wcześniej oligonukleotydem20. Mikrogele są ładowane do kropelek jako gęsto upakowane cząsteczki, co pozwala osiągnąć stopień enkapsulacji bliski 1 mikrogelowi na każdą kropelkę, z niewielką liczbą dubletów21 (Ryc. 4C).
W ostatnim kroku procesu mikroprzepływowego (Rysunek 2C), urządzenie wykonuje operację podwójnego łączenia, łącząc 1 kroplę z kodem kreskowym, 1 kroplę zawierającą mikrogel oraz mieszaninę do amplifikacji w kontrolowanym, dwuetapowym procesie. Najpierw kropla zawierająca odczynnik PCR zostaje sparowana i połączona z kroplą z kodem kreskowym w obszarze zaznaczonym na żółto (Rysunek 5). Elektrody solankowe w kanale mikroprzepływowym wytwarzają wysoki gradient pola elektrycznego, który inicjuje łączenie kropel. W podobny sposób pierwsza połączona kropla zostaje sparowana z kroplą mikrogelu i połączona po raz drugi w obszarze zaznaczonym na czerwono. Krople są zbierane i poddawane cyklom termicznym poza układem (off-chip) w reakcji PCR z jednokrotnym przedłużeniem nakładki (SOE PCR). Nakładające się komplementarne końce kodu kreskowego i tagmentowanego DNA genomowego umożliwiają fuzję i wykładniczą amplifikację wyłącznie prawidłowo okodowanych konstruktów.
Dane z sekwencjonowania są najpierw filtrowane pod kątem jakości odczytów, a następnie analizowane poprzez grupowanie odczytów zgodnie z ich 15-bp sekwencją kodów kreskowych dla pojedynczych komórek. Aby grupa kodów kreskowych została uznana za prawidłową, musi zawierać minimalną liczbę odczytów; takie progowanie ogranicza analizę do komórek z użyteczną ilością danych sekwencyjnych i usuwa z zestawu danych „sieroty”, czyli kody kreskowe zmutowane w wyniku PCR. W tym przykładowym przebiegu wartość minimalna została ustawiona na 7,5 kbps na grupę (50 odczytów po 150 bp każdy). Histogram liczby kodów kreskowych w stosunku do wielkości grupy pokazuje, że znaczna część prawidłowych grup kodów kreskowych znajduje się tuż powyżej wartości progowej (Rysunek 6A).
W eksperymencie kontrolnym, w którym skład społeczności mikrobiologicznej jest znany, do oceny jakości przebiegu SiC-seq wykorzystuje się wskaźniki czystości i względnej liczebności. W tym przypadku analizowana jest syntetyczna społeczność 10 komórek składająca się z 3 bakterii gram-ujemnych, 5 bakterii gram-dodatnich i 2 drożdży. Czystość danej grupy kodów kreskowych definiuje się jako liczbę odczytów mapujących do najczęstszego genomu w grupie podzieloną przez całkowitą liczbę odczytów w grupie. Zdecydowana większość grup kodów kreskowych wykazuje czystość przekraczającą 0,95 (Rysunek 6B). Względna liczebność typów komórek jest obliczana poprzez zliczanie surowych odczytów oraz zliczanie grup kodów kreskowych, gdzie grupom przypisuje się typ komórki odpowiadający konsensusowi odczytów będących ich członkami (Rysunek 6C). Liczebność odczytów i grup kodów kreskowych przebiega w zbliżonych proporcjach, co wskazuje, że populacje komórek są próbkowane w taki sposób, że niektóre gatunki nie są zawarte w nieproporcjonalnie małych lub dużych grupach kodów kreskowych. Wykreślenie zagregowanego pokrycia wszystkich grup kodów kreskowych z jednego gatunku wskazuje na wysokie pokrycie w obrębie całego genomu, z niewielką liczbą regionów z brakami (dropout) lub ich całkowitym brakiem (Rysunek 7). Jednolitość pokrycia można zweryfikować za pomocą rozkładu częstości znormalizowanych wartości pokrycia, w którym większość wartości skupia się wokół średniej (Rysunek 7, wstawka).

Rycina 1: Wytwarzanie urządzeń mikroprzepływowych metodą fotolitografii. (A) Formy matrycowe o pojedynczej wysokości struktur wytwarza się poprzez nanoszenie warstwy fotorezystu SU-8 na wafer krzemowy metodą spin coatingu. Następnie fotorezyst jest naświetlany przez maskę fotolitograficzną światłem UV, co powoduje sieciowanie naświetlonego SU-8. Na koniec nienasiękowany SU-8 zostaje rozpuszczony w kąpieli wywoływacza. Powstałą formę wykorzystuje się do odlewania PDMS, który jest następnie wiązany z szkiełkiem przedmiotowym w celu uzyskania kompletnego urządzenia mikroprzepływowego. (B) W przypadku urządzenia dwuwarstwowego proces wytwarzania analogicznie rozpoczyna się od etapów nanoszenia warstwy i naświetlania. Kroki te są następnie powtarzane w celu stworzenia urządzenia dwuwarstwowego. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Przegląd schematu pracy SiC-seq. (A) Zawiesina drobnoustrojów jest wprowadzana wspólnie z roztopioną agarozą do urządzenia do tworzenia kropli w celu enkapsulacji pojedynczych komórek w mikrożelach. (B) Mikrożele poddawane są serii płukań w celu oczyszczenia bakteryjnego DNA genomicznego. Enzymy lityczne trawią ściany komórkowe bakterii gram-dodatnich i drożdżaków, a detergent solubilizuje szczątki komórkowe. Mikrożele są ponownie enkapsulowane w kroplach w celu przeprowadzenia tagmentacji, aby zredukować zanieczyszczenia krzyżowe. (C) Mikrofluidyczny merger łączy kod kreskowy cyfrowego PCR, tagmentowany genom w mikrożelu oraz mieszaninę do amplifikacji z częstotliwością >1 kHz. Poza układem (off-chip) SOE-PCR łączy unikalny kod kreskowy pojedynczej komórki z tagmentowanym genomem i selektywnie amplifikuje w pełni obarkodowane konstrukty. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 3: Urządzenia mikroprzepływowe do wytwarzania kropel i ponownej enkapsulacji mikrożeli. (A) Ten panel przedstawia wytwórcę kropel w układzie współpływu (wysokość kanału 25 µm). Komórki i stopiona agaroza są wprowadzane do urządzenia z takimi samymi prędkościami przepływu, aby wytworzyć krople o średnicy 25 µm w złączu o wymiarach 25 µm x 25 µm. W przypadku wytwarzania kropel z cyfrowym kodem kreskowym wlot komórkowy zostaje zablokowany, a do wlotu agarozy wprowadzana jest mieszanina PCR. (B) Ten panel przedstawia urządzenie do ponownej enkapsulacji mikrożeli (wysokość kanału 25 µm). Mikrożele przepływają przez lejkowaty wlot, aby zachować ich gęste upakowanie, a następnie otrzymują objętość mieszaniny do tagmentacji przed ponowną enkapsulacją w złączu o wymiarach 25 µm x 30 µm. Proszę kliknąć tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 4: Mikrografie kropli i mikrożeli. (A) Ten panel przedstawia wypłukane mikrożele o średnicy 25 µm przed lizą enzymatyczną. Bakterie zostały zabarwione fluorescencyjnie w celu ilościowego oznaczenia stopnia enkapsulacji. Statystyka rozkładu Poissona wskazuje, że komórki powinny być enkapsulowane w tempie 1 na 10 kropli lub mniej, aby zminimalizować częstotliwość zdarzeń wielokrotnej enkapsulacji. (B) Ten panel przedstawia obraz z mikroskopii fluorescencyjnej kropli z cyfrowym kodem kreskowym o średnicy 25 µm, poddanych zabarwieniu kwasami nukleinowymi. Krople zawierające amplifikowane fragmenty kodów kreskowych emitują silny sygnał fluorescencyjny. (C) Ten panel przedstawia mikrożele ponownie enkapsulowane w kroplach o średnicy 50 µm. Gęste upakowanie mikrożeli pozwala na osiągnięcie stopnia enkapsulacji bliskiego 1 żelowi na kroplę przy niewielkiej liczbie dubletów. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rysunek 5: Mikroprzepływowe urządzenie do podwójnego łączenia kropel przeznaczone do kodowania paskowego genomu pojedynczych komórek. Dwuetapowa operacja łączenia pozwala na wysokoprzepustowe parowanie kropel z kodem kreskowym z tagmentowanymi genomami. Najpierw generowana jest kropla mieszaniny PCR, która zostaje połączona z kroplą kodu kreskowego w obszarze zaznaczonym na żółto za pomocą elektrod solankowych. Następnie wprowadzana jest kropla zawierająca mikrogel, która zostaje połączona po raz drugi w obszarze zaznaczonym na czerwono. Wloty oleju umożliwiają precyzyjną kontrolę odstępów między ponownie wstrzykiwanymi kroplami. Komora ponownego wstrzykiwania kodu kreskowego oraz olej separatorowy są umieszczone w cieńszej warstwie o grubości 25 µm, zaznaczonej na niebiesko. Pozostałe elementy urządzenia znajdują się w warstwie grubszej o całkowitej wysokości 45 µm. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 6: Metryki grup kodów kreskowych dla syntetycznej społeczności mikrobiologicznej składającej się z 10 komórek. (A) Ten panel przedstawia rozkład wielkości grup kodów kreskowych. Liczba grup danej wielkości maleje wykładniczo wraz ze wzrostem wielkości grupy. Minimalny próg 7,5 kbps na grupę ogranicza analizę do grup zawierających wystarczającą ilość informacji i usuwa „sieroty”, czyli sekwencje zmutowane podczas PCR. (B) Ten panel przedstawia rozkład czystości grup kodów kreskowych. Zdecydowana większość (>90%) grup wykazuje bardzo wysoką czystość (>95%). (C) Ten panel przedstawia względną liczność 10 gatunków obliczoną na poziomie odczytów oraz na poziomie grup kodów kreskowych. Obie metody zliczania dają podobne wyniki, co wskazuje na spójność wielkości grup kodów kreskowych pomiędzy gatunkami. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 7Zagregowane pokrycie genomowe dla Bacillus subtilis grupy kodów kreskowych Odczyty ze wszystkich grup kodów kreskowych mapujące do bakterii B. subtilis (N = 9 398) zostają połączone i przeanalizowane zbiorczo. Kołowa mapa pokrycia ilustruje jednorodność pokrycia odczytów SiC-seq, bez widocznych regionów luk (dropout). Przerywana linia wokół obwodu wskazuje średnie pokrycie (5,55x). Wstawiony histogram względnych częstotliwości pokrycia pokazuje, że większość zasad jest pokryta na głębokości zbliżonej do średniej dla całego genomu, reprezentowanej przez linię przerywaną. Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Urządzenie | Wysokość pierwszej warstwy (µm) | Prędkość wirowania pierwszej warstwy (obr./min) | Wysokość drugiej warstwy (µm) | Prędkość wirowania drugiej warstwy (obr./min) |
| Współpływowy generator kropel | 25 | 4000 | N/A | Nie dotyczy |
| Urządzenie do ponownej enkapsulacji w żelu | 25 | 4000 | Brak danych | Ponieważ nie dostarczono tekstu źródłowego do tłumaczenia, proszę o przesłanie treści, którą należy przełożyć na język polski. |
| Podwójna fuzja | 25 | 4000 | 20 | 5000 |
Tabela 1: Parametry wytwarzania urządzeń mikroprzepływowych. Tabela ta zawiera wykaz urządzeń mikroprzepływowych wykorzystywanych w procesie SiC-seq wraz z wymaganymi prędkościami nanoszenia fotorezystu metodą spin coatingu (na podstawie specyfikacji producenta dla SU-8 3025).
| Etykieta | Sekwencja (5' > 3') |
| BAR | GCAGCTGGCGTAATAGCGAGTACAATCTGCTCTGATGCCGCATAGNNNNNNNNNNNNNNNTAAGCCAGCCCCGACACT |
| KOD_KRESKOWY_DNA | CTGTCTCTTATACACATCTCCGAGCCCACGAGACGTGTCGGGGCTGGCTTA |
| P7_BAR | CAAGCAGAAGACGGCATACGAGATCAGCTGGCGTAATAGCG |
| P5_DNA | AATGATACGGCGACCACCGAGATCTACACTCGTCGGCAGCGTC |
| I7_CZYTAJ | GCCCACGAGACGTGTCGGGGCTGGCTTA |
Tabela 2: Sekwencje starterów.
Plik uzupełniający 1: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Aby pobrać ten plik, kliknij tutaj.