Zgodnie z protokołem opisanym w sekcji 1, konające czerwone tilapia hybrydowe wykazujące kliniczne objawy infekcji TiLV (Rycina 1A) zostały uśmiercone poprzez kąpiel w wysokim stężeniu olejku goździkowego, który działa jako środek znieczulający. Opisywane objawy kliniczne są zróżnicowane, lecz do najczęstszych należą letarg, erozja i odbarwienie skóry, egzophthalmus, odwarstwienie łusek, otwarte rany/zmiany oraz nieprawidłowe zachowanie15,16,33,35,36; niektóre z nich są wyraźnie widoczne na Rycini 1A. Ścianę brzuszną usunięto, aby pobrać narządy wewnętrzne, takie jak wątroba, śledziona lub nerka główna (Rycina 1B). Na tym etapie pobrano również próbki śluzu, delikatnie zeskrobując skórę ryby od przodu do tyłu za pomocą szkiełka podstawowego lub skalpela chirurgicznego37.

Rycina 1. Sekcja tilapia i pobieranie próbek.A. Czerwona tilapia hybrydowa zainfekowana TiLV z owrzodzeniami skóry, zaczerwienieniem wokół pyska i pokrywy skrzelowej, erozją skóry oraz zmętnieniem rogówki. B. Przecięta czerwona tilapia hybrydowa w celu pobrania tkanek z wątroby (w miejscu niebieskiej strzałki), śledziony lub głowowej nerki. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Następnie zastosowano protokół szczegółowo opisany w sekcji 2 dotyczący ekstrakcji całkowitego RNA z użyciem tiocyjanianu guanidyny, fenolu i chloroformu, a następnie przeprowadzono kwantyfikację RNA zgodnie z opisem w sekcji 3, aby ocenić czystość próbek poprzez obliczenie współczynników czystości i analizę profili widmowych (Rycina 2). Rycina 2A przedstawia reprezentatywny wynik pomyślnej procedury ekstrakcji całkowitego RNA, podczas gdy Rycina 2B obrazuje niską jakość preparatu RNA. Kwasy nukleinowe wykazują maksima absorbancji przy 260 nm, natomiast białka przy 280 nm. Stosunek pomiarów przy 260 nm i 280 nm wskazuje na czystość każdej próbki, a wartości od 1,9 do 2,1 świadczą o czystości RNA, co ma miejsce w przypadku próbki na Rycynie 2A. Niższe współczynniki A260/280 zaobserwowane na Rycynie 2B wskazują na możliwe zanieczyszczenie białkami lub fenolem pozostałe po procedurze ekstrakcji RNA. Absorbancja przy 230 nm może być wynikiem zanieczyszczenia próbki, dlatego oblicza się również stosunek A260/230 nm. W przypadku czystych preparatów RNA współczynnik ten powinien mieścić się w zakresie 2,0-2,2, co ilustruje wartość 2,03 dla próbki na Rycynie 2A, podczas gdy Rycina 2B wykazuje niski stosunek A260/230 wynoszący 1,07, a profil widmowy pokazuje przesunięcie minimum przy 230 nm w stronę 240 nm, co świadczy o obecności resztkowej guanidyny lub fenolu w próbce. W przypadku próbki przedstawionej na Rycynie 2B, ponowne wytrącanie RNA w celu usunięcia zanieczyszczeń może poprawić czystość preparatu.

Rycina 2. Spektrofotometryczna kwantyfikacja całkowitego RNA wyekstrahowanego z tkanek chorych tilapii. Awskaźniki czystości i profile widmowe z prawidłowo przeprowadzonej preparatyki RNA. BJak w przypadku A, z wyjątkiem próbki reprezentującej niską wydajność procedury ekstrakcji RNA. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.
Aby wykryć TiLV za pomocą metody RT-PCR, wykorzystuje się czyste próbki, takie jak ta przedstawiona w Rycina 2A poddano odwrotnej transkrypcji (protokół 4) do cDNA, które posłużyło jako matryca w analizie PCR opisanej w sekcji 5; reprezentatywne wyniki przedstawiono w Rysunek 3AStartery zaznaczone pogrubioną czcionką w Tabela 1 wykorzystano do amplifikacji fragmentu o długości 491 bp z segmentu 3 genomu TiLV14Produkty PCR rozdzielono metodą elektroforezy w żelu i zabarwiono EtBr w celu ich wizualizacji. Rysunek 3A przedstawia wyniki dwuetapowej reakcji RT-PCR z wykorzystaniem 4 próbek cDNA (S1-S4), wyizolowanych z wątroby chorych tilapi z Tajlandii; w każdej z próbek widoczny jest wyraźny pojedynczy prążek o wielkości około 500 bp, co oznacza, że próbki 1-4 są dodatnie w kierunku TiLV. Taki sam produkt PCR uzyskano w próbce kontroli dodatniej, zawierającej cDNA segmentu 3 TiLV sklonowane do plazmidu32 podczas gdy kontrola bez matrycy (NTC) nie dała produktów PCR. Analiza w Rycina 3B przeprowadzono przy użyciu tych samych starterów, co w Rysunek 3A lecz w innym laboratorium, stosując metodę jednoetapowej RT-PCR oraz 5 próbek RNA wyizolowanych z tkanek nerki głównej tilapi pochodzących z egipskiej akwakultury15Za pomocą tego testu detekcyjnego ustalono, że próbki 1, 3 i 5 są dodatnie w kierunku TiLV, natomiast próbki 2 i 4 są ujemne w kierunku TiLV, ponieważ nie stwierdzono produktu PCR o prawidłowej wielkości. Kontrole ujemne, w tym dwie kontrole bez odwrotnej transkryptazy oraz dwie kontrole NTC, nie wygenerowały żadnych produktów PCR. Przeprowadzono również jednokrokowy test RT-PCR ukierunkowany na gen ActinB tilapia. We wszystkich próbkach (S1-S5) zgodnie z oczekiwaniami uzyskano amplikon o wielkości 217 bp.38Analiza ta służyła jako kontrola integralności próbek RNA, a także umożliwiła półilościową ocenę próbek dodatnich w kierunku TiLV. Przy założeniu, że ilość uzyskanego produktu ActB Tilapii jest względnie stała, różnice w ilości wygenerowanego specyficznego produktu PCR TiLV można interpretować jako rzeczywiste odzwierciedlenie ilości TiLV w danej próbce tkanki.

Rycina 3. RT-PCR TiLV. A. Próbki cDNA otrzymane z tkanek wątroby chorych tilapi zebranych w Tajlandii zostały przebadane pod kątem zakażenia TiLV przy użyciu specyficznych starterów dla segmentu 3 (zaznaczonych pogrubieniem w Tabeli 1) TiLV w dwuetapowym teście RT-PCR. M = marker podany w parach zasad; S1-S4 = próbki 1-4; C1 = kontrola pozytywna z użyciem pTiLV jako matrycy PCR; C2 = kontrola bez matrycy (NTC). B. Jednoetapowy RT-PCR z użyciem tych samych starterów co w A oraz próbek z tkanek głowowej nerki chorych tilapi zebranych w Egipcie15. M = marker podany w parach zasad; S1-S5 = próbki 1-5. Kontrole C1-C2 to kontrole bez odwrotnej transkryptazy, a C3-C4 to kontrole NTC. Dolny panel przedstawia jednoetapowy RT-PCR z użyciem starterów skierowanych przeciwko ActinB tilapi38 (szczegóły w tekście), dający produkt PCR o długości 217 par zasad. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
W przeciwieństwie do reakcji PCR końcowopunktowych przedstawionych na Rysunku 3, testy qPCR omówione w protokole 6 mierzą ilość produktu PCR po każdym cyklu PCR. Amplifikację docelowego DNA wykrywa się przy użyciu cząsteczek fluorescencyjnych, które oddziałują z DNA generowanym w każdej rundzie reakcji. W tym przypadku wykorzystano barwnik SYBR Green I, który interpoluje się z dwuniciowym DNA. Sygnał fluorescencyjny jest monitorowany podczas reakcji, a jego intensywność odnosi się do ilości powstałego produktu 39,40,41,42,43. Testy qPCR dla TiLV przeprowadzono zgodnie z opisem w protokole 6 w różnych laboratoriach, stosując różne odczynniki SYBR Green, urządzenia do qPCR oraz próbki z różnych krajów. Wynikowe krzywe amplifikacji przedstawiono na Rysunku 4A i 4B. Można zauważyć, że w każdym teście przebieg eksperymentu składa się z czterech faz: liniowej fazy podstawowej, wczesnej fazy wykładniczej, późnej fazy wykładniczej oraz fazy plateau. Liniowa faza podstawowa występuje podczas wczesnych cykli, w których duplikacja DNA nie może być jeszcze zidentyfikowana ze względu na ilości DNA generujące niewystarczający stosunek sygnału do tła. W tej fazie obliczana jest fluorescencja bazowa. Następnie stężenie docelowego DNA zaczyna się podwajać w każdym cyklu, co powoduje, że sygnał staje się wykrywalny powyżej tła i wzrasta wykładniczo. Wydajność amplifikacji (E) dobrze zoptymalizowanego testu qPCR jest bardzo wysoka (bliska 100%) na początku reakcji i pozostaje stabilna podczas tej wczesnej wykładniczej fazy amplifikacji; to właśnie w tym punkcie, gdy wydajność reakcji jest nadal stała, przeprowadza się kwantyfikację. W późniejszych cyklach sygnał zaczyna osiągać plateau, a intensywność fluorescencji nie jest już skorelowana z początkową liczbą kopii matrycy, ponieważ składniki reakcji ulegają wyczerpaniu44. Nasycenie może również wystąpić w wyniku konkurencji ze strony reakcji ponownego hybrydyzowania, zmieniających się stosunków stężeń składników lub stosunku jednostek enzymu do cząsteczek substratu DNA. Prawdopodobnie takie parametry odpowiadają za różnice między krzywymi amplifikacji dla testów pokazanych na Rysunku 4A i 4B. Dołączone kontrole nie wygenerowały tych charakterystycznych krzywych amplifikacji.

Rycina 4. Wykresy amplifikacji przedstawiające akumulację produktu w czasie trwania testu real-time PCR. A. Krzywe amplifikacji próbek dodatnich pod kątem TiLV pochodzących z Tajlandii, NTC oraz dodatniej kontroli plazmidowej z zastosowaniem 2-etapowego testu qPCR z SYBR-Green I. Wykres został wygenerowany poprzez naniesienie względnej fluorescencji (RFU) vs. numeru cyklu. B. Krzywe amplifikacji próbek dodatnich pod kątem TiLV pochodzących z Egiptu, analogicznie jak na Rycini 3B, oraz NTC. Krzywa amplifikacji przedstawia fluorescencję sygnału reporterowego znormalizowaną do fluorescencji pasywnego barwnika ROX zawartego w teście (Rn) versus numer cyklu. Aby zobaczyć powiększoną wersję tej ryciny, kliknij tutaj.
Po zakończeniu cykli termicznych qPCR na poszczególnych urządzeniach w każdym z laboratoriów dane zostały pozyskane i przeanalizowane. Rysunek 5A i 5B przedstaw reprezentatywne krzywe topnienia z testów przeprowadzonych w każdym laboratorium. Każdy aparat do qPCR został zaprogramowany tak, aby na koniec przeprowadzić analizę krzywej topnienia. Osiągnięto to poprzez stopniowe zwiększanie temperatury i monitorowanie fluorescencji w funkcji temperatury. Gdy temperatura jest wystarczająco wysoka, aby zdenaturować dsDNA, rejestrowany jest znaczny spadek fluorescencji, ponieważ cząsteczka fluoroforu zostaje uwolniona. Oprogramowanie każdego instrumentu qPCR obliczyło temperaturę topnienia (Tm) z danych krzywej topnienia, wykreślając ujemną pierwszą pochodną w funkcji temperatury (Rysunek 5). Można zauważyć, że w Rycina 5A i 5B że produkty powstałe w różnych zestawach próbek wykazują jednorodną temperaturę topnienia w oczekiwanej temperaturze wynoszącej około 80 °C do analizy. Nie zaobserwowano żadnych innych pików w niższych temperaturach. Ze względu na ich niewielki rozmiar, Tm temperatury topnienia dimerów starterów jest zazwyczaj niższa niż w przypadku docelowej sekwencji DNA. Dlatego ta różnica między Tmpozwala to na łatwą identyfikację potencjalnych dimerów starterów lub innych niespecyficznych produktów amplifikacji. Próby kontrolne nie wygenerowały krzywych topnienia charakterystycznych dla próbek i standardów dodatnich pod kątem TiLV i są widoczne jako niemal pozioma linia u dołu wykresów w Rycina 5A i 5B.

Rycina 5. Analiza krzywej topnienia w celu zapewnienia specyficzności testu; różne produkty PCR mogą zostać rozróżnione na podstawie ich charakterystyki topnienia. A. Analiza krzywej topnienia próbek dodatnich pod kątem TiLV pochodzących z Tajlandii, kontroli ujemnej oraz dodatniej kontroli plazmidowej. B. Analiza krzywej topnienia próbek dodatnich pod kątem TiLV pochodzących z Egiptu, standardów pTiLV oraz NTC. Wykresy w A i B przedstawiają zmianę fluorescencji podzieloną przez zmianę temperatury w funkcji temperatury, co pozwala na uzyskanie przejrzystego obrazu dynamiki topnienia. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Większość urządzeń do qPCR jest wyposażona w oprogramowanie ułatwiające dalszą analizę przebiegu qPCR, które dokonuje kwantyfikacji próbek poprzez generowanie krzywej wzorcowej poprzez automatyczne wykreślenie cyklu progowego (Ct) względem logarytmu liczby kopii matrycy standardów pTiLV, co przedstawiono na rysunku 6A oraz 6B dla dwóch niezależnych laboratoriów. W skrócie, Ct jest jednostką stosowaną do oceny wyników qPCR. Wartość Ct oznacza liczbę cykli wymaganych do osiągnięcia ustalonego poziomu sygnału fluorescencji progowej. Im większa ilość matrycy wyjściowej, tym mniej cykli potrzeba do osiągnięcia wykrywalnego poziomu fluorescencji. Rzeczywiście, próbki z wysokim mianem TiLV będą miały niższe wartości Ct niż próbki z niskim mianem TiLV, na przykład u ryb z infekcją subkliniczną. Aby wyznaczyć wartości Ct, z surowych danych najpierw odejmuje się poziomy fluorescencji tła. Następnie oprogramowanie powiązane z urządzeniem do qPCR automatycznie wybierze próg fluorescencji, analizując krzywe danych dla każdej próbki i wyznaczając Ct reprezentujące punkt, w którym próbka przekroczyła próg. Jest to wykonywane oddzielnie dla każdego testu i każdy próg powinien zostać starannie oceniony, aby upewnić się, że został on ustawiony w części logarytmicznej krzywych amplifikacji i w miejscu, w którym wszystkie krzywe są równoległe. Zatem uzyskane konkretne Ct jest wartością względną względem liczby kopii matrycy wyjściowej45, ale jest ono również specyficzne dla użytego urządzenia do qPCR i odczynników, wydajności amplifikacji PCR oraz czułości detekcji. Parametry te wpływają na różnice zaobserwowane przy użyciu tego samego testu na rysunku 6.
Na podstawie krzywych wzorcowych przedstawionych na Rysunku 6 przeprowadzono analizy regresji, w tym obliczenia nachyleń (m) i punktów przecięcia krzywych wzorcowych, wydajności amplifikacji (100 x (101/m -1))46 oraz liniowości reakcji. Analizy krzywych wzorcowych wykorzystano również do potwierdzenia czułości (granicy wykrywalności), powtarzalności i odtwarzalności testu. Teoretycznie ilość DNA podwaja się w każdym cyklu PCR, co oznacza, że wydajność (E) wynosi 100%. Jednak w praktyce tak idealna wydajność jest rzadko osiągana ze względu na suboptymalne warunki PCR, takie jak inhibicja polimerazy DNA, zanieczyszczenia, zbyt duża ilość cDNA oraz błędy pipetowania47. Typowo wydajność amplifikacji E dla dobrych testów mieści się w zakresie 90-110%; na Rysunku 6A obliczono wydajność na poziomie 94,5% przy użyciu 8 seryjnie rozcieńczonych próbek pTiLV, natomiast wydajność testu przedstawionego na Rysunku 6B przy użyciu 7 seryjnie rozcieńczonych próbek pTiLV wyniosła 101,2%. Wydajność powyżej 100% wynika zazwyczaj z obecności inhibitorów PCR w teście. Analiza regresji liniowej wykresu wzorcowego pozwala również na obliczenie liczby kopii TiLV w każdej próbce41,42,45, co można zaobserwować dla trzech próbek TiLV zaznaczonych na czerwono na Rysunku 6B, co jest zgodne z wynikami dla próbek S1, S3 i S5 przedstawionymi na Rysunku 3B.

Rycina 6. Krzywe wzorcowe RT-qPCR. PCR w czasie rzeczywistym 10-krotnych rozcieńczeń szeregowych pTiLV, wzorca użytego w obu laboratoriach. A. Przebadano 8 próbek pTiLV w rozcieńczeniach szeregowych, wszystkie o znanym stężeniu skorelowanym z liczbą kopii TiLV/reakcję. Krzywą wzorcową wygenerowano poprzez wykreślenie logarytmu liczby kopii vs. próg cyklu (Ct). Nachylenie = -3,462, R2 = 0,9992, a wydajność wynosi 94,47%. B. Analogicznie jak w A, z wyjątkiem faktu, że przebadano 7 próbek pTiLV w rozcieńczeniach szeregowych (kolor zielony), a wykres przedstawia cykl progowy na osi y oraz liczbę kopii TiLV (Ilość) na osi x. Punkt przecięcia z osią y = 32,327, nachylenie = -3,292, R2 = 0,98, a wydajność wynosi 101,2%. W przypadku obu krzywych wzorcowych w A i B wartości nachylenia, punktu przecięcia z osią y oraz współczynnika korelacji (R2) są wykorzystywane do oceny sprawności testu. Co ważne, wartość R2 powinna być bliska 1, ponieważ jest ona miarą liniowości krzywej wzorcowej. Nachylenie służy do pomiaru wydajności PCR, gdzie wydajność 100% odpowiada nachyleniu -3,32; równanie oraz dalsze szczegóły znajdują się w tekście głównym. Dobra reakcja qPCR charakteryzuje się zazwyczaj wydajnością między 90-110%, co odpowiada nachyleniu między -3,58 a -3,10. Krzywa wzorcowa służy do bezwzględnej kwantyfikacji nieznanych próbek dodatnich pod kątem TiLV i pozwala określić dokładną liczbę kopii TiLV/reakcję, co ma miejsce w przypadku trzech próbek dodatnich TiLV zaznaczonych kolorem czerwonym w B.