W celu opracowania tej metody oceniono fenotypy ilościowe roślin kukurydzy B73 w temperaturach kontrolnych i podwyższonych, w warunkach odpowiedniego nawodnienia oraz suszy, a także w ich kombinacjach, co łącznie dało cztery warianty traktowania. B73 jest dobrze znaną linią inbredową kukurydzy z bogatymi danymi genetycznymi21.
Analiza obrazu w programie PlantCV pomyślnie skorygowała kolory zdjęć, przypisała wyniki ilościowe do metadanych roślin oraz przeprowadziła segmentację (innymi słowy, oddzielenie) rośliny od tła, co przedstawiono na Rysunku 3. Nieudana analiza dostarczyłaby obrazy wynikowe, na których nie widać niebieskiego obrysu ani różowego kształtu wokół rośliny, lecz występuje dodatkowy szum lub brakujące fragmenty rośliny (Rysunek 3, przykład niskiej jakości). Aby uwzględnić w masce jak największą część rośliny, konieczne było zastosowanie progu dwukanałowego (alternatywne strategie jednokanałowe opisano w pliku „photo-studio-SV-notebook.ipynb”). Niemniej jednak, u niektórych roślin, szczególnie tych z łodygami o fioletowym zabarwieniu, segmentacja nie była idealna, lecz okazała się wystarczająca do dalszej analizy (Rysunek 3).
Analiza w programie PlantCV dostarczyła 16 ilościowych cech o pojedynczej wartości, przedstawionych podsumowująco na Rysunku 4. Cecha o pojedynczej wartości dostarcza jednej wartości na roślinę – na przykład powierzchni liścia, wysokości, szerokości oraz średniej kołowej odcienia. Powierzchnię liścia określono przy użyciu funkcji „analyze.shape” (Rysunek 4), uwzględniając całą roślinę. Wysokość i szerokość określono natomiast za pomocą funkcji „analyze.bound_horizontal”, z linią ograniczającą wyznaczoną na poziomie linii gleby. Wynika to z faktu, że rośliny zwisające poza krawędź doniczki miałyby większą wysokość, gdyby traktować roślinę jako obiekt w całości; zamiast tego pomiar powinien być prowadzony wyłącznie od podstawy rośliny w miejscu jej styku z glebą. Oprócz cech o pojedynczej wartości, PlantCV generuje cechy wielowartościowe w formie histogramów, w szczególności dla wartości kolorów, reprezentujących każdy piksel rośliny. Średnia kołowa odcienia jest cechą o pojedynczej wartości, która stanowi średnią wartość odcienia dla całej rośliny.
Aby określić, które cechy są najistotniejsze, obliczono wariancję wyjaśnioną przez czynnik traktowania dla każdej cechy, a także korelację między poszczególnymi cechami, wykorzystując funkcję „frem” w pakiecie pcvr dla środowiska R Studio.20 (Rysunek 5A). Do dalszych analiz wybrano powierzchnię liścia, wysokość, szerokość oraz średnią kołową odcienia, ponieważ parametry te wyjaśniały ponad 50% wariancji wynikającej z zastosowanej próby i stanowią istotne miary wzrostu roślin (powierzchnia, wysokość i szerokość) oraz ich stanu zdrowotnego (odcień). Reprezentatywne zdjęcia roślin o największej i najmniejszej powierzchni liścia obrazują zróżnicowanie fenotypów zebranych w tym eksperymencie (Rycina 5B,C, odpowiednio). Roślina o największej powierzchni liści wykazywała również najwyższą średnią wartość odcienia i była rośliną linii B73, dobrze nawodnioną i poddaną stresowi termicznemu (Rycina 5B). Najmniejsza powierzchnia liścia charakteryzowała się najniższą średnią barwą (Rysunek 5C), i była to roślina odmiany B73 poddana stresowi suszy oraz stresowi cieplnemu.
Obróbka wodna miała istotny wpływ na powierzchnię liścia (F(1,13) = 226.5, p = 1.32 x 10-9), wysokość (F(1,13) = 21.1, p = 0.0005), szerokość (F(1,13) = 75.5, p = 8.92 x 10-7) oraz średnią kołową barwy (F(1,13) = 27.8, p = 0.0002), co określono za pomocą modelu regresji liniowej dla efektu obróbki oraz dwuczynnikowej analizy wariancji (ANOVA) w programie R Studio (pełne wyniki w Tabeli 1). Obróbka temperaturą miała istotny wpływ na wysokość (F(1,13) = 5.94, p = 0.03), lecz nie na powierzchnię liścia, szerokość ani średnią kołową barwy (p > 0.05, Tabela 1). Interakcja między temperaturą a stanem wodnym nie była istotnym czynnikiem dla żadnej z mierzonych cech (p > 0.05, Tabela 1).
Następnie poszczególne warianty zabiegów porównano z grupą kontrolną oraz między sobą, wykorzystując regresję liniową i test post-hoc (szacowane średnie marginalne) z korektą Sidaka dla porównań wielokrotnych (Rycina 6). Powierzchnia liści, wysokość rośliny, szerokość rośliny oraz średnia kołowa odcienia były znacząco zredukowane u roślin poddanych stresowi suszy w porównaniu do roślin dobrze podlanych w obu warunkach temperaturowych (p < 0,05, Rycina 6). Stres cieplny przy odpowiednim podlewaniu spowodował jedynie zmniejszenie wysokości rośliny (p < 0,05, Rycina 6), lecz nie wpłynął znacząco na powierzchnię liści, szerokość rośliny ani średnią kołową odcienia. Różnice w odcieniu można dokładniej zbadać poprzez analizę wielowartościowej cechy koloru – odcienia, przedstawionego w formie histogramu na Rycynie 7, zamiast opierać się wyłącznie na średniej. Obniżenie średniego odcienia spowodowane suszą wynikało z przesunięcia z pikseli zielonych na żółte, co jest znane jako chloroza (Rycina 7). Choć w przypadku wysokiej temperatury nie odnotowano znaczącej różnicy w średnim odcieniu, zaobserwowano zarówno redukcję, jak i wzrost wartości wynikające z żółknięcia (chlorozy) oraz ciemnienia koloru zielonego (Rycina 7).
Po obrazowaniu rośliny odcięto w miejscu przejścia korzenia w pęd, a następnie odnotowano masę części nadziemnej. Korzenie delikatnie przemyto i osuszono ręcznikiem, a następnie zważono, aby obliczyć całkowitą biomasę rośliny (zwaną również masą) w celu porównania z opisaną metodą fenotypowania opartą na analizie obrazu. Powierzchnia liści zmierzona poprzez analizę obrazu była silnie skorelowana z biomasą rośliny (R2 = 0,84, Rycina 8).

Rysunek 1: Zdjęcie układu studia fotograficznego z opisanymi elementami służącymi do gromadzenia obrazów roślin. (A) tkanina klasy fotograficznej, (B) karta kolorów, (C) linia z taśmy do jednorodnego rozmieszczania doniczek, (D) aparat cyfrowy, (E) komputer jednopłytkowy, (F) monitor wyświetlający proces przechwytywania obrazu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Przegląd procesu obrazowania. Przedstawiono reprezentatywny obraz o „wysokiej jakości” oraz „niskiej jakości”. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Przegląd przepływu pracy analizy obrazu w PlantCV. Przedstawiono reprezentatywne obrazy segmentacji o jakości „wysokiej”, „akceptowalnej” i „niskiej”; cechy PlantCV są reprezentowane za pomocą różowego obrysu kształtu całej rośliny oraz niebieskiego obrysu obszaru rośliny. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Podsumowanie i wizualna reprezentacja cech kształtu uzyskanych z analizy. Funkcja kształtu w PlantCV. Fioletowe pola reprezentują cechy, które nie wykorzystują innych cech w swoich obliczeniach. Pomarańczowe pola reprezentują cechy, które wykorzystują inne cechy w swoich obliczeniach; strzałki wskazują, które cechy są użyte w tych obliczeniach. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 5: Fenotypowanie oparte na obrazowaniu rejestruje wariancję fenotypową i cechy skorelowane. (A) Wariancja cech wyjaśniona przez zabieg (po prawej) oraz korelacja cech między sobą (po lewej). Każda cecha jest oznaczona, a następnie podano jej jednostkę pomiaru; „NA” oznacza pomiar bezjednostkowy. (B) Obraz rośliny, która miała największą powierzchnię liści i najwyższą kołową średnią barwy (stopnie). (C) Obraz rośliny, która miała najmniejszą powierzchnię liści i najniższą kołową średnią barwy (stopnie). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 6: Porównanie wpływu zabiegów na fenotypy roślin za pomocą wykresu pudełkowego. Porównano powierzchnię liścia (A), wysokość (B), szerokość (C) oraz średnią kołową odcienia (D) pod kątem wpływu suszy, wysokiej temperatury oraz stresu kombinowanego. n = 4. Litery oznaczają statystycznie istotne różnice między zabiegami (p < 0.05); zabiegi porównano za pomocą regresji liniowej i testu post-hoc (estymowane średnie marginalne) z korektą Sidaka dla porównań wielokrotnych. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 7: Histogramy sumy kolorów (odcienia) wszystkich pikseli dla każdej kombinacji zabiegów, znormalizowane przez całkowitą liczbę pikseli w danej kombinacji. Histogramy są zabarwione zgodnie z wartością koloru dla danego stopnia odcienia. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 8: Korelacja powierzchni liści zmierzonej z obrazów z biomasą całej rośliny. Każdy czarny punkt reprezentuje jedną roślinę. Różowa linia przedstawia liniową korelację między powierzchnią liści (X) a biomasą (Y). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Cecha | Czynnik | stopnie swobody | wartość F | wartość p |
| powierzchnia liścia (cm^2) | Temperatura | 1 | 0.5372 | 0.4766 |
| powierzchnia liścia (cm^2) | Woda | 1 | 226.5318 | 1.32E-09 |
| powierzchnia liścia (cm^2) | Temperatura:Woda | 1 | 2.1526 | 0.1661 |
| powierzchnia liścia (cm^2) | Reszty | 13 | | |
| wysokość (cm) | Temperatura | 1 | 5.9426 | 0.0298935 |
| wysokość (cm) | Woda | 1 | 21.107 | 0.0005029 |
| wysokość (cm) | Temperatura:Woda | 1 | 2.1752 | 0.1640545 |
| wysokość (cm) | Reszty | 13 | | |
| szerokość (cm) | Temperatura | 1 | 2.2495 | 0.1575 |
| szerokość (cm) | Woda | 1 | 75.5474 | 8.92E-07 |
| szerokość (cm) | Temperatura:Woda | 1 | 0.0839 | 0.7766 |
| szerokość (cm) | Reszty | 13 | | |
| średnia kołowa barwy (stopnie) | Temperatura | 1 | 0.9154 | 0.3561432 |
| średnia kołowa barwy (stopnie) | Woda | 1 | 27.7988 | 0.0001509 |
| średnia kołowa barwy (stopnie) | Temperatura:Woda | 1 | 3.4792 | 0.0848724 |
| średnia kołowa barwy (stopnie) | Reszty | 13 | | |
Tabela 1: Wyniki modelu regresji liniowej i analizy dwuczynnikowej ANOVA badających wpływ traktowania wodą, temperatury oraz ich interakcji na cechy roślin.