W niniejszym badaniu wykorzystano zautomatyzowane fenotypowanie obrazowe w celu zbadania odpowiedzi morfologicznych i fizjologicznych ziemniaka (cv. Lady Rosetta) pod wpływem pojedynczego i połączonego stresu. Zastosowane podejście wykazało dynamiczne odpowiedzi roślin z wysoką rozdzielczością czasoprzestrzenną w momencie wywołania stresu na etapie inicjacji bulw. Aby ocenić wczesne i późne fazy stresu, wyniki przedstawiono w trzech okresach czasowych ([0-5 dni fenotypowania (DOP)], [6-10 DOP] oraz [11-15 DOP]) (Rysunek 1). Do 0 DOP wszystkie rośliny były hodowane w warunkach kontrolnych (C), a następnie w dniach 1-5 DOP zastosowano stres zalewowy (W) oraz stres termiczny (H). W związku z tym odpowiedzi obserwowano w następujący sposób: (i) w dniach 0-5 DOP określono początkowy stres termiczny i zalewowy; (ii) w dniach 6-10 DOP obserwowano wczesną suszę (D) oraz połączony stres termiczny i suszę (HD) oraz (iii) w dniach 11-15 DOP wykazano późny stres termiczny, suszę oraz połączony stres termiczny + susza + zalew (HDW). Regenerację po zalaniu obserwowano w okresach 6-10 DOP oraz 11-15 DOP.
Cechy morfologiczne
Obrazowanie RGB zastosowano w celu określenia wpływu różnych stresów i ich kombinacji na wzrost nadziemnych części roślin. Wyniki przedstawione na Rysunku 4 pokazują, że obróbka termiczna oraz stres spowodowany podtopieniem (0-5 DOP) już powodują zmniejszenie objętości roślin i RGR w porównaniu z grupą kontrolną. W okresie 6-10 DOP objętość roślin i RGR roślin kontrolnych stale rosły, natomiast w warunkach wysokiej temperatury oraz połączonego stresu cieplnego, suszy i podtopienia, wzrost objętości roślin był wyraźnie ograniczony (Rysunek 4A). Ze względu na wysoką podatność roślin na stres spowodowany podtopieniem, odnotowano znaczny spadek RGR (Rysunek 4B). Podczas późnego stresu suszowego (11-15 DOP), przy utrzymaniu SRWC na poziomie 20%, zaobserwowano wyraźny spadek RGR w porównaniu z grupą kontrolną. Jednak w późnej fazie połączonego stresu HDW zastosowanie zabiegu podtopienia spowodowało wzrost RGR w ostatnim dniu stresu.
Cechy fizjologiczne
W celu ujawnienia dalszych reakcji na stres zastosowano kombinację fenotypowania strukturalnego i fizjologicznego. Wykorzystanie wielu czujników obrazujących umożliwia określenie odpowiedzi fizjologicznych w wczesnej fazie stresu. Dalsza analiza danych fluorescencji chlorofilu wykazała, że zalewanie negatywnie wpływało na wydajność fotosyntetyczną, gdzie Fv'/Fm' (Fv/Fm_Lss) gwałtownie spadło w dniach 0-5 DOP i 6-10 DOP, natomiast w dniach 11-15 DOP zaobserwowano reakcję regeneracyjną, w której Fv'/Fm' nieznacznie wzrosło (Rysunek 5A). Podczas późnej fazy stresu (11-15 DOP) zaobserwowano spadek Fv'/Fm' w warunkach suszy oraz połączonego stresu cieplnego i suszy. W roślinach zalewanych wydajność operacyjna roślin (QY_Lss, tzw. φPSII) była znacznie niższa w porównaniu z innymi wariantami w dniach 0-5 DOP i 6-10 DOP, ale odnotowano niewielki wzrost w dniach 11-15 DOP, co wskazuje na regenerację roślin (Rysunek 5B). Ponadto, różne mechanizmy regulujące wydajność przyczyniające się do ochrony PSII określono poprzez obliczenie frakcji otwartych centrów reakcyjnych w PSII w stanie stacjonarnym oświetlenia (qL_Lss) (Rysunek 5C). Tylko w warunkach suszy zaobserwowano wzrost qL, co prawdopodobnie wynikało z fotoinhibicji.
Wyniki te były zgodne z danymi IR, które odzwierciedlały różne mechanizmy leżące u podstaw reakcji na stresy (Rycina 6). W warunkach zalewy zaobserwowano wzrost deltaT (ΔT), co spowodowało zmniejszenie szybkości wymiany gazowej. W przypadku późnej suszy oraz połączonego stresu wysoką temperaturą i suszą wzrost ΔT wynikał z zamknięcia aparatów szparkowych, co uznaje się za jedną z głównych odpowiedzi w celu uniknięcia nadmiernej utraty wody. Z drugiej strony, w przypadku zabiegów wysoką temperaturą zaobserwowano spadek ΔT, ponieważ aparaty szparkowe otwierają się, aby zwiększyć efektywność transpiracji i ochłodzić powierzchnię liścia.
Na podstawie analizy danych hiperspektralnych z zakresu VNIR wybrano dwa parametry do oceny indeksów odbicia liści: NDVI jako wskaźnika zawartości chlorofilu oraz PRI jako wskaźnika wydajności fotosyntezy. Wyniki wykazały spadek NDVI i PRI wyłącznie w warunkach zalewania, co korelowało z obserwowaną redukcją cech morfologicznych (Rysunek 7A,B). Ponadto w hiperspektralnych danych z zakresu SWIR, wykorzystanych do oceny zawartości wody w roślinach, zaobserwowano wzrost indeksu wody podczas zalewania w okresie 0-5 DOP (Rysunek 7C). Natomiast w przypadku traktowania temperaturą zaobserwowano odwrotną reakcję, w której indeks wody był niższy niż w grupie kontrolnej. Wyniki te były zgodne z analizą roślinności przeprowadzoną na podstawie segmentacji kolorystycznej obrazów RGB z widoku z góry. Zmiany w proporcjach odcieni wskazują na odpowiedzi stresowe w czasie (Rysunek 8). Indeks zazielenienia wykazał redukcję zawartości pigmentów w warunkach suszy oraz połączonego stresu HDW w późnej fazie stresu, a także stopniową regenerację po traktowaniu zalewaniem. Zatem zastosowanie wielu czujników obrazujących odzwierciedliło korelację cech morfofizjologicznych i umożliwiło ocenę ogólnej kondycji roślin pod wpływem stresów abiotycznych.

Rysunek 1: Oś czasu stosowania różnych zabiegów, obejmująca wiek roślin w dniach po przesadzeniu sadzonek in vitro. Dzień 0 fenotypowania (DOP) został zmierzony w warunkach kontrolnych (C), a następnie indukowano różne stresy o różnym czasie trwania. Od 1. do 5. DOP stosowano stres podtopieniowy (W) oraz początkową reakcję na zabieg wysoką temperaturą (H). W kolejnych dniach 6-10 DOP przedstawiono początkową fazę stresu suszy (D) oraz połączony stres wysoką temperaturą i suszą (HD). W dniach 11-15 DOP odzwierciedlono reakcję roślin na późną fazę zabiegów suszy i wysokiej temperatury oraz zastosowanie podtopienia w grupie HD (HDW) przez 1 dzień. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 2: Schemat podsumowujący protokół fenotypowania i analizę danych. (A) Przegląd protokołu fenotypowania. Rośliny są transportowane do systemu fenotypowania z kontrolowanych warunków w komorze wzrostu FS-WI (PSI). Przed pomiarami rośliny były aklimatyzowane do światła w komorze adaptacyjnej przez 5 min przy natężeniu 500 µmol.m-2.s-1. Do określenia cech morfologicznych i fizjologicznych wykorzystano wiele czujników obrazujących, po czym następowała stacja ważenia i podlewania. W zależności od traktowania rośliny były przenoszone z powrotem do kontrolowanych warunków, albo w temperaturze 22 °C/19 °C, albo 30 °C/28 °C. (B) Automatyczna ekstrakcja i segmentacja w ramach potoku przetwarzania obrazów z każdego czujnika obrazującego. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Przegląd krótkiego protokołu oświetleniowego dla obrazowania fluorescencji chlorofilu. Protokół pomiarowy rozpoczynał się od włączenia chłodnego białego światła aktynicznego w celu pomiaru fluorescencji stanu stacjonarnego w świetle (Ft_Lss), a następnie zastosowania impulsu nasycającego do pomiaru maksymalnej fluorescencji stanu stacjonarnego w świetle (Fm_Lss). Światło aktyniczne zostało wyłączone, a włączono światło dalekiej czerwieni w celu oznaczenia minimalnej fluorescencji stanu stacjonarnego w świetle (Fo_Lss). Czas trwania protokołu wynosił 10 s na roślinę. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4: Obrazowanie RGB wykorzystane do oceny morfologicznej. (A) Objętość rośliny obliczona na podstawie pola powierzchni widoków RGB z góry i z boku. (B) Względna szybkość wzrostu (RGR) w fazie inicjacji bulw. Dane przedstawiają wartości średnie ± odchylenie standardowe (n = 10). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Obrazowanie fluorescencji chlorofilu u roślin zaadaptowanych do światła. (A) Maksymalna wydajność fotochemii PSII próbki zaadaptowanej do światła w stanie stacjonarnym oświetlenia (Fv/Fm_Lss). (B) Wydajność kwantowa fotosystemu II lub wydajność operacyjna fotosystemu II w stanie stacjonarnym oświetlenia (QY_Lss). (C) Odsetek otwartych centrów reakcyjnych w PSII w stanie stacjonarnym oświetlenia (utleniony QA) (qL_Lss). Dane przedstawiają wartości średnie ± odchylenie standardowe (n = 10). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6: Wykorzystano obrazowanie termowizyjne IR do obliczenia różnicy między średnią temperaturą korony rośliny wyznaczoną z obrazów termowizyjnych a temperaturą powietrza (ΔT). Dane przedstawiają wartości średnie ± odchylenie standardowe (n = 10). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 7: Obrazowanie hiperspektralne w celu wyznaczenia wskaźników roślinności i zawartości wody. (A) Znormalizowany Różnicowy Wskaźnik Wegetacji (NDVI). (B) Fotochemiczny Wskaźnik Refleksji (PRI) obliczony na podstawie obrazowania VNIR. (C) Wskaźnik wody obliczony na podstawie obrazowania SWIR. Dane przedstawiają wartości średnie ± odchylenie standardowe (n = 10). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 8: Indeks zieleni dla roślin poddanych różnym zabiegom. Przetwarzanie obrazu opiera się na transformacji oryginalnego obrazu RGB w mapę kolorów składającą się z 6 zdefiniowanych odcieni. Dane reprezentują wartości średnie ± odchylenie standardowe (n = 10). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Rycina uzupełniająca 1: Natężenie światła mierzone w dniach fenotypowania (DOP). Czas trwania pomiarów od 9:00 do 12:35. LI_Buff odnosi się do danych mediany z 5 czujników światła rozmieszczonych w szklarni. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 2: Wilgotność względna (RH) mierzona w dniach fenotypowania (DOP).Czas trwania pomiarów od godziny 9:00 do 12:35. RH_Buff odnosi się do średnich danych z 5 czujników wilgotności rozmieszczonych w szklarni. RH2 odnosi się do wilgotności względnej w komorze adaptacyjnej. Proszę kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 3: Temperatura mierzona w dniach fenotypowania (DOP).Czas trwania pomiarów od godziny 9:00 do 12:35. T_Buff odnosi się do mediany danych z 5 czujników temperatury rozmieszczonych w szklarni. T2 odnosi się do temperatury w komorze adaptacyjnej. T3 odnosi się do temperatury ściany grzejnej. T4 odnosi się do temperatury w jednostce obrazowania termowizyjnego IR. Proszę kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 4: Zrzut ekranu z oprogramowania do analizy danych przedstawiający parametry dostosowane do analizy maski roślin w czujnikach obrazowania fluorescencji chlorofilu. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 5: Zrzut ekranu z oprogramowania do analizy danych przedstawiający parametry dostosowane do analizy maski roślin w czujnikach obrazowania w podczerwieni termicznej. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 6: Zrzut ekranu z oprogramowania analizatora danych przedstawiający parametry dostosowane do analizy maski rośliny w czujnikach obrazujących widok jednostronny RGB. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 7: Zrzut ekranu z oprogramowania do analizy danych przedstawiający parametry dostosowane do analizy maski roślin w czujnikach obrazujących widok z góry RGB2. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 8: Zrzut ekranu z oprogramowania analizatora danych pokazujący parametry dostosowane do analizy maski roślin w czujnikach obrazowania VNIR. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 9: Zrzut ekranu z oprogramowania do analizy danych przedstawiający parametry dostrojone do analizy maski roślin w czujnikach obrazowania SWIR. Proszę kliknąć tutaj, aby pobrać ten plik.