W niniejszej pracy wykorzystano bezpośrednie wyliczanie czasowe produkcji fagów z kultury lizogenu PAO1 LESΦ2 hodowanej w warunkach nieindukujących, aby określić wpływ spontanicznej indukcji LESΦ2. Gęstość fagów była najniższa, osiągając średnią wartość ~2,61 x 106 jednostek tworzących plaki (PFU)·mL−1 2 h po przesiewie do świeżego medium w wczesnej fazie wzrostu wykładniczego, co sugeruje, że stanem dominującym była lizogenia. Miano LESΦ2 szybko wzrosło do średniej wartości ~2,4 x 108 PFU·mL−1 w ciągu 4 h i osiągnęło najwyższą gęstość po 6 h (średnio ~5,83 x 109 PFU·mL−1; Rycina 4).
Minimalna indukcja spontaniczna była obserwowana podczas wczesnej fazy logarytmicznej wzrostu lizogenu (po 2 h). Jednak mierzalna obecność fagów w pożywce była wynikiem wielu wcześniejszych zdarzeń, w tym: pakowania kwasów nukleinowych do białkowych główek, składania białek w cząsteczki fagowe oraz ekspresji późnych genów fagowych, genów fazy środkowej i wczesnych regulatorowych genów fagowych. Istotne było uchwycenie zainfekowanych komórek przed wystąpieniem zdarzeń replikacyjnych związanych z fagiem; dlatego wybrano czas 90 min na wzrost kultury przed indukcją. Aby uchwycić profil ekspresji genów w próbkach lizogenu PAO1, LESΦ2, pobrano materiał z kultury przed indukcją oraz po indukcji w okresie 90 min, zgodnie z opisem w kroku 3.4. Ten punkt czasowy (90 min) przypada znacznie wcześniej niż wykrycie wysokiego poziomu spontanicznej indukcji rezydentowego profaga za pomocą testu plackowego z kroku 2.3.2. Ponieważ gęstość komórek bakteryjnych we wczesnej fazie wzrostu wykładniczego była niska, objętości kultur zwiększono do 800 mL, aby zapewnić wystarczającą ilość materiału do badań ekspresji genów. Próbki z kultury nieindukowanej i kultur indukowanych pobierano co 10 min, a następnie wyizolowano RNA w celu zmapowania profilu ekspresji kluczowych markerów lizogenii i replikacji litycznej podczas wzrostu bakterii. Całkowite RNA oczyszczono i zweryfikowano pod kątem braku genomowego DNA za pomocą analiz qPCR ukierunkowanych na gen 16S rRNA (krok 6.1). Próbki, które osiągnęły RIN ≥ 9, przeszły kontrolę jakości i zostały przekształcone w cDNA.
Przeanalizowano adnotowany genom LESΦ2 w celu zidentyfikowania genów znanych z udziału w lizogennych i litycznych cyklach replikacji fagów umiarkowanych. Zidentyfikowane geny wykorzystano następnie do walidacji qRT-PCR w celu profilowania ekspresji genów ograniczonych do cyklu lizogennego oraz genów związanych z cyklem litycznym w kulturach indukowanych i nieindukowanych. Ilościowo określono absolutną liczbę kopii DNA i przeprowadzono test znaków Wilcoxona przy użyciu programu R36, aby porównać poziomy ekspresji w kulturach nieindukowanych i indukowanych (Rysunek 5). Zaobserwowano wyraźny wzrost ekspresji genu cro (wczesnego markera replikacji litycznej) z poziomu ~2,31 x 109 kopii w kulturach nieindukowanych do ~3,02 x 1011 kopii 30 min po indukcji (test znaków Wilcoxona: p < 0,01). Podobnie białka O i białka P, które są markerami średniego etapu replikacji litycznej (a według przewidywań uczestniczą w replikacji genomu faga), wykazały również istotny wzrost ekspresji z ~1,74 x 108 do ~1,25 x 1010 kopii (test znaków Wilcoxona: p < 0,01) oraz z ~6,05 x 102 do ~5,68 x 105 kopii (test znaków Wilcoxona: p < 0,01), odpowiednio. Na koniec, strukturalne geny związane z ogonkiem wykorzystano jako późne markery litycznego cyklu replikacji. Ponownie zaobserwowano istotny wzrost ekspresji z ~2,31 x 106 kopii w kulturach nieindukowanych do ~4,38 x 108 kopii 30 min po indukcji (test znaków Wilcoxona: p < 0,01). Zatem dane z ilościowego RT-PCR potwierdziły, że ekspresja genów dobrze ugruntowanych markerów replikacji litycznej przebiegała zgodnie z oczekiwanym trendem, przy czym markery wczesne, średnie i późne wykazywały wielokrotnie różnicową ekspresję w przewidzianej kolejności (Rysunek 5). Ponieważ ekspresja markerów replikacji litycznej uległa zwiększeniu 30 min po regeneracji, punkt ten uznano za odpowiedni reprezentatywny moment czasowy do badania krajobrazu transkryptomicznego aktywnych fagów umiarkowanych i ich gospodarzy bakteryjnych podczas cyklu litycznego.
Zaobserwowaliśmy pewną ekspresję genów litycznych w warunkach bez indukcji, co potwierdza, że zawsze dochodzi do pewnej indukcji spontanicznej, nawet w zoptymalizowanych kulturach, w których liczba lizogenów jest reprezentowana przez najwyższy stosunek CFU do uwolnionych PFU we wczesnej fazie logarytmicznej. Oznacza to, że w danych transkryptomicznych zawsze będzie występował pewien poziom „szumu”, co podkreśla znaczenie starannie przygotowanych kontroli, w tym kultur indukowanych i nieindukowanych. Właściwy wybór genów kontroli wewnętrznej w celu określenia krotności zmian w ekspresji opiera się na dokładnej analizie danych transkryptomicznych w celu zidentyfikowania genów, które są wyrażane na tym samym poziomie zarówno w próbkach nieindukowanych, jak i indukowanych. Nasze wstępne wyniki sugerują, że rpoD był najbardziej wiarygodnym przetestowanym genem kontrolnym i wykazywał najstabilniejszą ekspresję (~1,71 x 105 kopii przed indukcją i ~3,33 x 105 kopii 30 min po indukcji; test znaków rang Wilcoxona: p = 0,3594) w porównaniu do rrn16S lub genów proC (Rysunek 5). Zmienność ekspresji kontroli wewnętrznych doprowadziła do pomiaru bezwzględnej liczby transkryptów. Dalsza analiza danych transkryptomicznych pomoże w wyborze odpowiednich kontroli wewnętrznych do późniejszej walidacji.
W naszym ćwiczeniu z profilowania genów wykorzystano gen cI, ponieważ jest on uznanym markerem lizogenii. W porównaniu do markerów replikacji litycznej, ekspresja genu cI była stosunkowo stabilna (Rycina 5), jednak liczba kopii tego genu w kulturach nieindukowanych była satysfakcjonująco wysoka w porównaniu do markerów replikacji litycznej. Dane te są zgodne z niskimi wartościami PFU w tych samych próbkach, co potwierdza, że wysoka ekspresja represora wiązała się z niższym poziomem produkcji fagów. Przedstawione tutaj dane wykazują, że ekspresja transkryptu cI dla tego konkretnego faga nie ulega istotnemu obniżeniu po indukcji, co ma miejsce w przypadku fagów Stx11,17. Aktywność represora jest zazwyczaj kontrolowana zarówno na poziomie transkrypcyjnym, jak i potranslacyjnym, zatem gen represora może być transkrybowany, ale powstałe białko ulega natychmiastowemu autokatalitycznemu cięciu. Wymagane są dalsze eksperymenty w celu walidacji kontroli transkrypcyjnej i potranslacyjnej. Ponadto, na podstawie naszej krzywej wzorcowej, minimalna granica wykrywalności qPCR wydaje się wynosić ~102 kopii.
Łącznie wyniki uzyskane w testach płytkowych i za pomocą qRT-PCR potwierdzają słuszność naszej strategii hodowli i przygotowania próbek RNA w celu uzyskania dobrze kontrolowanego materiału wejściowego do eksperymentów RNA-Seq. Hodowle nieindukowane w wczesnej fazie wykładniczej wykazywały niski poziom spontanicznej indukcji i ekspresji genów litycznych, co sugeruje dominację lizogenii. Z kolei hodowle izolowane 30 min po indukcji wykazały znaczny wzrost ekspresji genów markerowych wskazujących na dominację replikacji litycznej.

Rycina 1: Protokół tworzenia gospodarza wskaźnikowego opornego na rifampicynę (Utworzono w BioRender.com). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Projekt eksperymentalny dla wyznaczania PFU i CFU lizogenu z tej samej próbki. (Utworzono za pomocą BioRender.com) Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 3: Schemat doświadczalny pobierania próbek z kultur indukowanych i nieindukowanych do izolacji RNA. (Utworzone w BioRender.com) Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Czasowe wyliczanie indukcji spontanicznej. Czasowe wyliczanie produkcji spontanicznego profaga LES z wykorzystaniem PFU z lizogenu PAO1 Φ2 wraz z jednoczesnym CFU, n = 8 (dwie powtórzenia biologiczne i cztery techniczne); słupki błędu reprezentują odchylenie standardowe. Ciemnoczerwone punkty wskazują CFU·mL−1 w LB; ciemnoniebieskie punkty wskazują PFU·mL−1 w LB. Spontaniczne uwalnianie zakaźnego faga φ2 przez lizogeny znajduje się na najniższym mierzalnym poziomie 2 h po inokulacji. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Bezwzględna liczba kopii docelowych genów markerowych. Bezwzględna liczba kopii genów markerowych faga potwierdza przewidywane wzorce ekspresji, uzyskane za pomocą RT-qPCR, dla genów, których oczekuje się, że odgrywają istotną rolę w cyklach lizogennych i litycznych. Kropki reprezentują zarówno trzy powtórzenia biologiczne, jak i trzy techniczne (n = 9). (A) Czerwone pudełko reprezentuje marker lizogenii, cI; (B) zielony reprezentuje wczesny marker lityczny, cro; (C,D) niebieski reprezentuje środkowe markery lityczne, geny replikacji DNA; (E) purpurowy reprezentuje późny marker lityczny, geny strukturalne ogonka; (F–H) szary reprezentuje markery gospodarza, które zostały użyte jako kontrole wewnętrzne, a (I) biały reprezentuje gyrazę DNA B, która została użyta jako kontrola indukcji. Ciągłe poziome linie wskazują medianę rozkładu. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Tabela 1: Startery zaprojektowane w niniejszym badaniu. Podano sekwencje swoistych starterów dla genów markerowych i kontroli wewnętrznych wykorzystanych w niniejszym badaniu, wraz z odpowiadającymi im numerami dostępu NCBI. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Efektywność starterów użytych w niniejszym badaniu, obliczona na podstawie krzywej standardowej qPCR. Kliknij tutaj, aby pobrać tę tabelę.