Walidacja metody Worm-to-CT jako sposobu przygotowania cDNA
Aby sprawdzić, czy protokół Worm-to-CT jest prawidłową metodą ekstrakcji cDNA, porównano go ze standardowymi metodami ekstrakcji z użyciem tiocyjanianu guanidyny, fenolu i chloroformu. Wyniki przedstawiono na Rysunku 3, gdzie cDNA przygotowano ze średnio ~1 000 nicieni przy użyciu standardowych technik ekstrakcji tiocyjanianem guanidyny, fenolem i chloroformem22 oraz z 30 nicieni przy użyciu metody Worm-to-CT. Próbki poddano jednoczesnemu szokowi termicznemu (30 min w 34 °C). Ogólnie poziomy ekspresji mRNA hsp-70 na 100 ng całkowitego RNA były porównywalne przy zastosowaniu obu metod. Jednak w przypadku najwyższej ekspresji hsp-70 (tj. w szczepie N2 po szoku termicznym) poziomy ekspresji były wyższe przy zastosowaniu metody Worm-to-CT, co wskazuje na poprawioną czułość.
Aby sprawdzić, czy można odtworzyć oczekiwany spadek ekspresji hsp u hsf-1(sy441)23, mutacji głównego regulatora transkrypcyjnego chaperonów molekularnych23,24, porównano indukcję transkrypcyjną chaperonów po krótkim szoku termicznym. W obu metodach stwierdzono spadek indukcji hsp-70 u zwierząt hsf-1(sy441). Było to spodziewane, ponieważ mutanty hsf-1(sy441) wykazują obniżoną zdolność do indukcji chaperonów z powodu truncacji w domenie transaktywacyjnej HSF-1. W przypadku ekstrakcji guanidianem tiocyjanianu-fenolem-chloroformem poziom hsp70 spadł o 82,7% w porównaniu z kontrolą, a w metodzie Worm-to-CT o 92,3% w porównaniu ze zwierzętami typu dzikiego (Rycina 3). Wyniki były porównywalne dla obu metod oraz zgodne z wcześniejszymi doniesieniami23. Wyniki te wskazują, że metoda Worm-to-CT stanowi poprawną alternatywę dla standardowych technik syntezy cDNA.
Walidacja platformy nanofluidycznej PCR wykorzystywanej do amplifikacji celów mRNA
Aby przetestować spójność wyników uzyskanych przy użyciu nanofluidycznej qPCR w celu amplifikacji transkryptów, porównano wyniki PCR otrzymane metodą masową Worm-to-CT zarówno w standardowym systemie qPCR (Tabela materiałów) oraz nanofluidyczny system qPCR z wykorzystaniem chipa wielkoprzestrzennego. Krotność zmiany ekspresji trzech różnych genów, sma-3 (Rysunek 4A), sma-10 (Rycina 4B), i dnj-26, monitorowano (Rysunek 4C) u zwierząt posiadających allel null w dbl-1 (dbl-1(nk3))25 w porównaniu do odpowiedników typu dzikiego. Dbl-1 koduje jedyny ligand szlaku sygnalizacyjnego białka morfogenetycznego kości (BMP). sma-3 i sma-10 są genami kodującymi ortologi SMAD, kluczowe komponenty kaskady sygnałowej BMP. Dnj-26 koduje chaperonę molekularną, będącą celem sygnalizacji BMP. Wyniki te wykazują niewielką lub żadną różnicę w krotności zmiany (fold change) przy porównaniu rezultatów obu metod, co skutkuje nieistotnymi statystycznie wartościami P wynoszącymi odpowiednio 0,3113, 0,2635 i 0,3481 dla sma-3, sma-10i dnj-26odpowiednio. Łącznie wyniki te wykazują, że metoda Worm-to-CT zastosowana do próbek zbiorczych jest wydajnym i szybkim sposobem izolacji RNA z niewielkiej liczby nicieni i dostarcza wiarygodnych danych w połączeniu ze standardowymi systemami PCR lub wysokoprzepustowymi platformami qPCR opartymi na nanofluidyce.
Porównanie poziomów ekspresji uzyskanych z próbek zbiorczych z wartościami średnimi uzyskanymi z pojedynczych nicieni
Względne poziomy ekspresji obliczono, wykorzystując cDNA uzyskane z próbek zbiorczych (25 nicieni) lub ze średniej z 36 próbek pojedynczych nicieni (Rysunek 5). Oba rodzaje cDNA uzyskano metodą Worm-to-CT i powielono z wykorzystaniem technologii nanofluidycznej PCR. Jak zaobserwowano na Rysunku 5A–C, dla wszystkich badanych chaperonów (tj. hsp16.1, F44E5.4, hsp-70) obie metody wykryły porównywalne poziomy ekspresji. Wyniki te wskazują, że parametry uzyskane z pojedynczych nicieni są wiarygodne.
Zastosowanie technologii Worm-to-CT połączonej z nanofluidyką do szacowania parametrów ekspresji genów pojedynczych nicieni
Ponieważ chip z jedną macierzą umożliwia monitorowanie do 96 transkryptów docelowych w 96 pojedynczych próbkach, nadaje się on do badania zmienności osobniczej w ekspresji transkryptów między pojedynczymi nicieniami. Rycina 6A przedstawia reprezentatywny wynik ukazujący średnią ekspresję wielu transkryptów hsp z pojedynczych nicieni po krótkim szoku termicznym. Jak zaobserwowano na rycinie, zmienność ekspresji transkryptów różniła się drastycznie w zależności od genu (Rycina 6A). Aby uzyskać więcej informacji, obliczono współczynnik zmienności (CV), dzieląc odchylenie standardowe przez średnią poziomów ekspresji26 (Rycina 6B). Monitorowano trzy geny, których wartości CV zostały wcześniej oszacowane alternatywnymi metodami (dane niepublikowane). Dwa stabilne transkrypty (ife-1 oraz Y45F10D.4) i jeden zmienny (nlp-2927) wykazały oczekiwaną zmienność. Wykres wyraźnie przedstawia również znany odwrotny związek między wartościami zmienności a poziomami ekspresji26 (Rycina 6B).
Powtórzenia techniczne są kluczowe dla zapewnienia powtarzalności w przypadku stosowania próbek zbiorczych. Niekoniecznie ma to jednak miejsce w eksperymentach na pojedynczych komórkach.14,15,28Aby ustalić, czy w przypadku stosowania próbek pojedynczych nicieni konieczne jest wykorzystanie powtórzeń technicznych do szacowania parametrów, 28 pobrano pojedyncze nicienie, po krótkim szoku termicznym, a następnie przetworzono z wykorzystaniem triplikatów technicznych. Wartości CV obliczone z danych dla pojedynczych nicieni uzyskanych w trzech powtórzeniach (niebieskie kropki w Rycina 7, techniczny współczynnik zmienności (CV)) w porównaniu do tych dla każdego transkryptu uzyskanego z poszczególnych nicieni (czerwone kropki w Rycina 7, zmienności biologicznej) poddano porównaniu. W przypadku każdego przetestowanego transkryptu współczynniki zmienności (CV) technicznej były niższe niż współczynniki zmienności biologicznej, co wskazuje, że do oszacowania parametrów nie były wymagane powtórzenia techniczne w trzech kopiach. Brak konieczności wykonywania powtórzeń technicznych zwiększa przepustowość eksperymentu bez obniżania jego jakości.

Rysunek 1: Przegląd protokołu Worm-to-CT.
Na tym rysunku przedstawiono krótki przegląd poszczególnych etapów wymaganych do przeprowadzenia analizy nicieni zgodnie z protokołem Worm-to-CT. Przedstawiono dwie opcjonalne metody dla etapu odwrotnej transkrypcji; są to metody zamienne, możliwe do zastosowania dla każdego typu chipa. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Przegląd przygotowania i przeprowadzania nanofluidycznego qPCR.
Na tym rysunku przedstawiono przygotowania do uruchomienia nanofluidycznego systemu qPCR z wykorzystaniem chipa wielomaćcianowego oraz chipa jednocząstkowego. Aby wyświetlić większą wersję tego rysunku, kliknij tutaj.

Rysunek 3: Protokół Worm-to-CT zastosowany do próbek zbiorczych dostarczył wiarygodnych wyników.
Porównanie protokołu Worm-to-CT z regularną ekstrakcją tiocyjanianem guanidyny, fenolem i chloroformem22 w próbkach zbiorczych. Zgodnie z wcześniejszymi ustaleniami, u mutantów hsf-1(sy441)23 poziomy transkryptów hsp w odpowiedzi na szok termiczny uległy obniżeniu. Powyższe histogramy przedstawiają indukcję hsp-70 w przypadku braku (-) lub po (+) krótkim szoku termicznym trwającym 30 min w 34 °C. cDNA uzyskano za pomocą ekstrakcji tiocyjanianem guanidyny, fenolem i chloroformem zastosowanej do 1 000 nicieni (lewo) lub za pomocą metody Worm-to-CT zastosowanej do 30 połączonych nicieni (prawo). Porównano poziomy ekspresji hsp-70 na 100 ng całkowitego RNA uzyskanego każdą z metod. Zgodnie z oczekiwaniami, u hsf-1(sy441) indukcja transkrypcyjna hsp-70 w odpowiedzi na szok termiczny znacząco spadła o 82,7% przy użyciu tiocyjanianu guanidyny, fenolu i chloroformu oraz o 92,3% przy użyciu metody Worm-to-CT. Poziomy mRNA genów docelowych znormalizowano względem średniej z trzech genów housekeepingowych: cdc-42, pmp-3 i ire-1. Każda kropka reprezentuje powtórzenie biologiczne. Dane zostały poddane transformacji logarytmicznej do analizy statystycznej, ponieważ nie spełniały założeń wymaganych do analizy parametrycznej. Analizę statystyczną przeprowadzono za pomocą RM-One-way ANOVA z zastosowaniem testu wielokrotnych porównań Sidaka. Typ dziki = N2, hsf-1 = hsf-1(sy441). Słupki oznaczają błąd standardowy średniej. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 4: Wzorce ekspresji były spójne pomiędzy standardowymi systemami qPCR a nanofluidycznymi systemami qPCR.
(A) Poziom ekspresji mRNA dla sma-3 (A), sma-10 (B) lub dnj-26 (C) został określony za pomocą standardowego qPCR oraz nanofluidycznego qPCR (chip wielkopolowy) z trzech powtórzeń biologicznych cDNA otrzymanego metodą Worm-to CT z dzikiego szczepu (N2) oraz szczepu z nokautem dbl-1(nk3)25. Względne poziomy ekspresji mRNA dla każdego szczepu określono metodą Delta-Ct21. Następnie obliczono krotność zmiany (fold change), dzieląc poziomy ekspresji uzyskane u nicieni dbl-1(nk3) przez odpowiadające im poziomy mRNA w szczepie N2. Jak pokazano na panelu A, wzorce były spójne dla obu metod w każdym poszczególnym powtórzeniu biologicznym. (B) i (C) przedstawiają analogiczne wyniki jak w przypadku (A) dla poziomów mRNA sma-10 i dnj-26 odpowiednio. Poziomy docelowego mRNA znormalizowano względem genów housekeepingowych cdc-42 i pmp-3. Analizę statystyczną dla każdego genu przeprowadzono za pomocą sparowanego testu t, porównując wyniki trzech powtórzeń biologicznych uzyskanych za pomocą standardowego qPCR z wynikami uzyskanymi za pomocą nanofluidycznego qPCR. Wartości P dla tych porównań wynosiły odpowiednio 0,3113, 0,2635 i 0,3481 dla sma-3, sma-10 i dnj-26. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Zastosowanie metody Worm-to-CT w próbkach zbiorczych lub w pojedynczych nicieniach pozwoliło uzyskać podobne poziomy ekspresji po normalizacji na jednego nicienia.
Poziomy ekspresji (A) hsp-16.1/11, (B) F44E5.4 oraz (C) hsp-70 (C12C8.1) analizowano u młodych osobników dorosłych w warunkach braku szoku termicznego, przeprowadzając procedurę Worm-to-CT albo na zbiorczej próbce 25 osobników, albo na 36 pojedynczych osobnikach. Po normalizacji danych na jednego nicienia nie stwierdzono istotnych różnic między poziomami uzyskanymi na jednego nicienia dla każdego transkryptu przy zastosowaniu obu metod. Poziomy mRNA genów docelowych znormalizowano względem średniej z trzech genów pomocniczych (housekeeping genes): cdc-42, pmp-3 oraz ire-1. Słupki reprezentują błąd standardowy średniej. Analiza statystyczna = test t dla prób zależnych. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 6: Wysokoprzepustowa analiza RT-qPCR pojedynczych nicieni z wykorzystaniem metody Worm-to-CT pozwala na monitorowanie zmienności międzyosobniczej w ekspresji genów.
(A) średnie poziomy ekspresji 53 transkryptów uzyskane po krótkotrwałym szoku termicznym (30 min w temperaturze 34 °C). Wykresy pudełkowe przedstawiają rozkład średniej ekspresji mRNA z poszczególnych nicieni (dla każdego osobnika wykorzystano średnią z trzech powtórzeń technicznych). Kropki reprezentują poziomy ekspresji u 28 poszczególnych nicieni. Poziomy mRNA genów docelowych znormalizowano względem średniej z trzech genów pomocniczych. cdc-42, pmp-3i ire-1. (B) Współczynnik zmienności26 (CV) jako funkcja średniej ekspresji mRNA dla 53 transkryptów po ekspozycji na krótki szok termiczny została obliczona na podstawie danych z 28 osobników (surowe dane przedstawiono w panelu B). Zbiór transkryptytów obejmuje zmienną nlp-29 transkrypcja27 oraz dwa stabilne transkrypty (ife-1 i Y45F10D.4; niepublikowane dane). CV (współczynnik zmienności) jest stosunkiem odchylenia standardowego do średniej. Współczynnik ten wykorzystano do oszacowania zmienności międzyosobniczej w ekspresji transkryptów pomiędzy poszczególnymi osobnikami nicieni. Zgodnie z oczekiwaniami, zmienność międzyosobnicza rosła wraz ze spadkiem średnich poziomów ekspresji. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 7: Powtórzenia techniczne nie były konieczne podczas analizy zmienności międzyosobniczej w ekspresji genów z wykorzystaniem chipa nanofluidycznego.
Dane przedstawione na tym wykresie uzyskano u 28 poszczególnych nicieni po krótkim szoku termicznym (30 min w 34 °C). Każda czerwona kropka reprezentuje współczynnik zmienności (CV) średnich poziomów ekspresji transkryptu dla jednego badanego transkryptu pomiędzy 28 poszczególnymi nicierniami (bio CV). Każda niebieska kropka reprezentuje CV poziomów ekspresji pomiędzy trzema powtórzeniami technicznymi uzyskanymi z jednego nicienia, na każdy badany transkrypt (technical CV). Wykres ten pokazuje, że zmienność techniczna (między powtórzeniami technicznymi) była znacznie niższa niż zmienność biologiczna (między poszczególnymi nicierniami), co sugeruje, że przy badaniu ekspresji genów w pojedynczych nicierniami, podobnie jak w badaniach pojedynczych komórek14,15,28, nie ma potrzeby wykonywania powtórzeń technicznych na nanofluidycznym układzie do ekspresji genów. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Tabela 1: Schemat rozmieszczenia dla chipa wielomatrycowego. Powyższa tabela przedstawia prosty układ, który można wykorzystać podczas planowania analizy na chipie wielomatrycowym. Po lewej stronie znajdują się miejsca, które należy wypełnić docelowymi starterami, a po prawej miejsca przeznaczone na analizowane próbki. Każda matryca testów i próbek jest sparowana numerycznie w obrębie chipa. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Plan układu dla chipa z pojedynczą macierzą. Powyższa tabela przedstawia prosty układ, który można wykorzystać podczas planowania pracy z chipem z pojedynczą macierzą. Po lewej stronie znajdują się miejsca, które należy wypełnić docelowymi starterami, a po prawej miejsca przeznaczone na badane próbki. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Lista starterów RT-qPCR użytych w tym badaniu. Kliknij tutaj, aby pobrać tę tabelę.
Tabela uzupełniająca 1: Primery z bazy danych primerów do RT-qPCR. Kliknij tutaj, aby pobrać tę tabelę.