$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Eksperymentalny przepływ pracy SiC-seq zawiera 3 urządzenia mikroprzepływowe PDMS wykonane przy użyciu procedury miękkiej litografii (Rysunek 1). Kroplarka współprzepływowa (Rysunek 3A) generuje 25 μm kropelek cyfrowego kodu kreskowego do znakowania genomowego DNA za pomocą unikalnego identyfikatora pojedynczej komórki. Oligonukleotydy kodu kreskowego składają się ze zdegenerowanej sekwencji o długości 15 pz otoczonej uchwytami PCR do amplifikacji (tabela 2, starter BAR). Kody kreskowe są rozcieńczane do stężenia femtomolowego w celu uzyskania enkapsulacji pojedynczej cząsteczki, a wszystkie kropelki otrzymują 0 lub 1 fragment(y) kodu kreskowego. Kropelki zawierające kod kreskowy są amplifikowane, co daje wiele kopii dwuniciowych amplikonów kodów kreskowych. Barwienie kwasem nukleinowym służy do weryfikacji pomyślnej amplifikacji i ilościowego określenia szybkości enkapsulacji fragmentów kodu kreskowego (Rysunek 4B). Mikrożele są wytwarzane przez współprzepływ zawiesiny komórek bakteryjnych i stopionego żelu agarozowego przy równych prędkościach przepływu (Rysunek 2A). Agaroza jest przygotowywana w dwukrotnie wyższym stężeniu końcowym niż pożądane, ponieważ proces wytwarzania kropli w przepływie skutecznie rozcieńcza roztwory wodne 2-krotnie. Gdy agaroza stygnie, krzepnie w mikrożel o średnicy 25 μm, zajmujący kulistą objętość kropli.
Seria etapów mycia i lizy oczyszcza genomowe DNA o wysokiej masie cząsteczkowej w mikrożelach (Rysunek 2B). Po rozdrobnieniu emulsji przeprowadza się przemywanie wodą w dużych ilościach w celu rozcieńczenia śladowych rozpuszczalników organicznych, które mogą hamować dalszą obróbkę enzymatyczną. Umyte mikrożele obserwuje się pod mikroskopem w celu sprawdzenia szybkości enkapsulacji komórek (Rysunek 4A). Do zawiesiny mikrożelu dodaje się koktajl enzymów o szerokiej aktywności litycznej w celu trawienia ścian komórkowych bakterii i drobnoustrojów eukariotycznych19. Druga obróbka proteinazą K i detergentem degraduje białka i rozpuszcza resztki komórek.
Znakowanie oczyszczonego DNA odbywa się w kropelkach, aby uniknąć potencjalnego zanieczyszczenia krzyżowego wynikającego z dyfuzji małych oznaczonych fragmentów DNA między mikrożelami18. Urządzenie do kapsułkowania kropelek ( Rysunek 3B) dzieli każdy mikrożel na przedziały za pomocą buforu i enzymu tagującego, który jednocześnie fragmentuje dwuniciowe DNA, jednocześnie "oznaczając" je wstępnie załadowanym oligonukleotydem20. Mikrożele są ładowane do kropelek jako gęsto upakowane cząstki, osiągając szybkość enkapsulacji zbliżającą się do 1 mikrożelu na każdą kroplę z kilkoma dubletami21 (Rysunek 4C).
W końcowym kroku przepływu pracy mikroprzepływów (Rysunek 2C), urządzenie wykonuje operację podwójnego łączenia, łącząc 1 kroplę kodu kreskowego, 1 kroplę zawierającą mikrożel i mieszankę wzmacniającą w kontrolowanym dwuetapowym procesie. Najpierw kropla zawierająca odczynnik PCR jest parowana i łączona z kroplą kodu kreskowego w obszarze pokazanym na żółto (Rysunek 5). Elektrody słonowodne w kanale mikroprzepływowym wytwarzają wysoki gradient pola elektrycznego, który powoduje łączenie się kropel. W podobny sposób pierwsza scalona kropla jest parowana z kroplą mikrożelu i łączona po raz drugi w obszarze pokazanym na czerwono. Kropelki są zbierane i poddawane cyklom termicznym poza chipem w PCR z pojedynczym nakładaniem się (SOE). Nakładające się na siebie komplementarne końce kodu kreskowego i oznakowane genomowe DNA umożliwiają fuzję i wykładniczą amplifikację tylko prawidłowo oznaczonych konstruktów oznaczonych kodem kreskowym.
Dane sekwencjonowania są najpierw filtrowane według jakości odczytu, a następnie analizowane poprzez grupowanie odczytów według sekwencji kodów kreskowych 15 bp z pojedynczą komórką. Aby grupa kodów kreskowych została uznana za ważną, powinna zawierać minimalną liczbę odczytów; ten próg ogranicza analizę do komórek z użyteczną ilością danych sekwencjonowania i usuwa "osierocone" kody kreskowe z mutacją PCR ze zbioru danych. W tym przykładowym przebiegu wartość minimalna jest ustawiona na 7,5 kb/s 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 rozmiaru progowego (Rysunek 6A).
W eksperymencie kontrolnym, w którym znany jest skład społeczności mikrobiologicznej, wskaźniki czystości i względnej obfitości są używane do oceny jakości serii SiC-seq. Tutaj analizowana jest syntetyczna 10-komórkowa społeczność 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 odwzorowujących najpowszechniejszy genom w grupie podzieloną przez całkowitą liczbę odczytów w grupie. Zdecydowana większość grup kodów kreskowych ma czystość większą niż 0,95 (Rysunek 6B). Względna obfitość typów komórek jest obliczana przez zliczanie surowych odczytów i zliczanie grup kodów kreskowych, gdzie grupom przypisuje się typ komórki odpowiadający konsensusowi odczytów ich członków (Rysunek 6C). Obfitość odczytów i grup kodów kreskowych przebiega w mniej więcej równych 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 wysoki zasięg w całym genomie, z niewielką liczbą regionów rezygnacji lub bez nich (Rysunek 7). Jednorodność pokrycia można zweryfikować za pomocą rozkładu częstości znormalizowanych wartości pokrycia, przy czym większość wartości jest wyśrodkowana wokół średniej (Rysunek 7, wstawka).

Rysunek 1: Wytwarzanie urządzeń mikroprzepływowych za pomocą fotolitografii. (A) Formy wzorcowe o pojedynczej wysokości cechy są wytwarzane przez wirowanie warstwy fotorezystu SU-8 na płytce krzemowej. Fotorezystor jest następnie modelowany za pomocą maski fotolitograficznej i światła UV, sieciując naświetlony SU-8. Na koniec nieusieciowany SU-8 rozpuszcza się w kąpieli wywoływaczowej. Powstała forma jest używana do odlewania PDMS, który jest połączony ze szklanym szkiełkiem w celu wytworzenia kompletnego urządzenia mikroprzepływowego. (B) W przypadku urządzenia dwuwarstwowego, produkcja podobnie rozpoczyna się od etapów powlekania wirowego i naświetlania. Te kroki są następnie powtarzane w celu stworzenia urządzenia dwuwarstwowego. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Przegląd przepływu pracy SiC-seq. (A) Zawiesina mikrobiologiczna jest współpłynna ze stopioną agarozą w urządzeniu do wytwarzania kropli, aby zamknąć pojedyncze komórki w mikrożelach. (B) Mikrożele są poddawane serii płukania w celu oczyszczenia bakteryjnego genomowego DNA. Enzymy lityczne trawią ściany komórkowe bakterii Gram-dodatnich i drożdży, a detergent rozpuszcza resztki komórkowe. Mikrożele są ponownie kapsułkowane w kropelki w celu znakowania w celu zmniejszenia zanieczyszczenia krzyżowego. (C) Fuzja mikroprzepływowa łączy cyfrowy kod kreskowy PCR, oznakowany genom mikrożelu i mieszankę amplifikacyjną z częstotliwością >1 kHz. Off-chip SOE-PCR łączy unikalny jednokomórkowy kod kreskowy z oznakowanym genomem i selektywnie amplifikuje w pełni zakodowane konstrukty. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Urządzenia mikroprzepływowe do wytwarzania kropli i re-kapsułkowania mikrożeli. (A) Ten panel przedstawia kroplówkę współprzepływową (25 μm wysokości elementu). Komórki i stopiona agaroza są wprowadzane do urządzenia przy równych natężeniach przepływu, aby wytworzyć kropelki 25 μm na złączu 25 μm x 25 μm. W przypadku cyfrowego tworzenia kodów kreskowych wlot ogniwa jest zatykany, a do wlotu agarozy wprowadzana jest mieszanina PCR. (B) Ten panel przedstawia urządzenie do rekapsułkowania mikrożelu (25 μm wysokości elementu). Mikrożele przepływają do wlotu w kształcie lejka, aby utrzymać ich ciasno upakowaną kolejność i otrzymują objętość mieszanki znakującej przed ponownym enkapsulacją w złączu 25 μm x 30 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Mikrofotografie kropelek i mikrożeli. (A) Ten panel pokazuje umyte mikrożele o wielkości 25 μm przed lizą enzymatyczną. Bakterie są barwione fluorescencyjnie w celu ilościowego określenia szybkości enkapsulacji. Statystyki ładowania Poissona dyktują, że komórki powinny być kapsułkowane z szybkością 1 na 10 kropli lub mniej, aby zminimalizować częstotliwość zdarzeń wielokrotnej enkapsulacji. (B) Ten panel pokazuje obraz z mikroskopii fluorescencyjnej kropelek cyfrowego kodu kreskowego o wielkości 25 μm potraktowanych barwieniem kwasem nukleinowym. Kropelki zawierające wzmocnione fragmenty kodu kreskowego wytwarzają silny sygnał fluorescencyjny. (C) Ten panel pokazuje mikrożele ponownie zamknięte w kroplach o wielkości 50 μm. Ścisłe upakowanie mikrożeli pozwala na szybkość kapsułkowania zbliżającą się do 1 żelu na kroplę z kilkoma dubletami. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Mikroprzepływowe urządzenie do podwójnego łączenia do kodowania kreskowego genomu pojedynczej komórki. Dwuetapowa operacja łączenia łączy kropelki kodów kreskowych z oznaczonymi genomami z dużą przepustowością. Najpierw generowana jest kropla mieszaniny PCR i łączona z kropelką kodu kreskowego w obszarze pokazanym na żółto za pomocą elektrod słonowodnych. Następnie wprowadza się kroplę zawierającą mikrożel i łączy się ją po raz drugi w obszarze pokazanym na czerwono. Wloty oleju pozwalają na precyzyjną kontrolę odstępów między ponownie wstrzykiwanymi kropelkami. Komora wtrysku kodów kreskowych i jej olej dystansowy umieszczone są na krótszej warstwie 25 μm, zacienionej na niebiesko. Wszystkie inne cechy urządzenia należą do grubszej warstwy o całkowitej wysokości 45 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Metryki grupy kodów kreskowych dla 10-komórkowej społeczności syntetycznych mikroorganizmów. (A) Ten panel pokazuje rozkład rozmiarów grup kodów kreskowych. Liczba grup o danej wielkości zmniejsza się wykładniczo wraz ze wzrostem wielkości grupy. Minimalny próg wynoszący 7,5 kb/s na grupę ogranicza analizę do grup z wystarczającą ilością informacji i usuwa "osierocone" sekwencje z mutacją PCR. (B) Ten panel pokazuje rozkład czystości grup kodów kreskowych. Zdecydowana większość (>90%) grup ma bardzo wysoką czystość (>95%). (C) Ten panel pokazuje względną liczebność 10 gatunków obliczoną na poziomie grupy odczytu i kodu kreskowego. Te 2 metody liczenia dają podobne wyniki, co wskazuje, że rozmiary grup kodów kreskowych są spójne dla wszystkich gatunków. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Zagregowane pokrycie genomowe grup kodów kreskowych Bacillus subtilis. Odczyty ze wszystkich grup kodów kreskowych odwzorowujących bakterię B. subtilis (N = 9 398) są łączone i analizowane w ujęciu zagregowanym. Okrągła mapa zasięgu ilustruje jednorodność zasięgu odczytów sekwencyjnych SiC, bez obserwowalnych obszarów zaniku. Przerywana linia na obwodzie oznacza średnie pokrycie (5,55x). Wstawka histogramu względnych częstości pokrycia pokazuje, że większość zasad jest pokryta na głębokości bliskiej średniej dla całego genomu, reprezentowanej przez linię przerywaną. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
szt.
szt.
szt.
szt.
| urządzenie | Wysokość 1. warstwy (μm) | Prędkość wirowania 1. warstwy (obr./min) | Wysokość 2. warstwy (μm) | Prędkość wirowania 2. warstwy (obr./min) |
| Kroplowarka współprzepływowa | 25 | 4000 | N/a | N/a |
| Re-enkapsulator żelowy | 25 | 4000 | N/a | N/a |
| Podwójna fuzja | 25 | 4000 | 20 | 5000 |
Tabela 1: Parametry produkcji urządzeń mikroprzepływowych. Poniższa tabela przedstawia listę urządzeń mikroprzepływowych używanych w przepływie pracy SiC-seq wraz z ich wymaganymi prędkościami do powlekania spinowego fotorezystu (w oparciu o specyfikacje producenta dla SU-8 3025).
| etykieta | >Sekwencja (5' > 3') |
| bar | GCAGCTGGCGTAATAGCGAGTACAATCTGCTCTGATGCCGCATAGNNNNNNNNNNNNNNNNTAAGCCCCCCGACACT |
| DNA_BAR | CTGTCTCTTATACACATCTCCGAGCCCCCACGAGCGAGCGGGGCTGGCTCTTA |
| P7_BAR | CAAGCAGAAGACGGCATACGAGATCAGCTGGCGTAATAGCG |
| P5_DNA | AATGATACGGCGACCGAGATCTACACTCGTCGGCAGCGTC |
| I7_READ | GCCCACGAGACGTGTCGGGGCTGGCTTA |
Tabela 2: Sekwencje podkładów.
Plik uzupełniający 1: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Kliknij tutaj, aby pobrać ten plik.