Reprezentatywne wyniki przedstawione w niniejszym manuskrypcie pochodzą z obszaru badawczego ustanowionego w 2012 roku na farmie University of Nebraska-Lincoln Agriculture Research Division w pobliżu Mead, NE. Przed rozpoczęciem eksperymentu teren ten był zarządzany w systemie płodozmianu kukurydza-soja. Obszar badań obejmował trzy różne typy gleb, jednak dane analizowano przy założeniu, że wszystkie zmiany w zmierzonych właściwościach gleby wynikały z zastosowanych zabiegów.
Teren badawczy obejmował dwa czyste drzewostany rozsadnika (P. virgatum cv Liberty) i trzęślicy wielkiej (A. gerardii), a także mieszankę traw o niskiej różnorodności, zawierającą trzęślicę wielką, trawę indyjską (S. nutans) oraz trawę sideoats grama 'Butte' (B. curtipendula). Trzy działki traw ciepłosezonowych zostały rozmieszczone w zrandomizowanym schemacie bloków zupełnych z trzykrotnymi powtórzeniami. W obrębie trzech różnych działek traw zastosowano dwie warianty nawożenia azotowego (N), polegające na podaniu 56 (N1) i 112 (N2) kg N ha-1 w postaci mocznika. W momencie pobierania próbek mikrobiomu na koniec sezonu wegetacyjnego gleba zawierała 8,0 ± 1,1 (średnia ± SD) ppm azotanów w działkach nawożonych dawką 112 kg N ha-1 oraz 6,8 ± 0,7 (średnia + SD) ppm azotanów w działkach nawożonych dawką 56 kg N ha-1. Działki były nawożone raz w roku. Działki traw ciepłosezonowych wyznaczono jako działki główne (8000 m2), a warianty N jako pod-działki (4000 m2). Trzęślica wielka została wysiana jako mieszanka 'Bonanza' i 'Goldmine' w stosunku 50:50, a trawa indyjska jako mieszanka 'Scout' i 'Warrior' w stosunku 50:50. Działki zostały obsiane w 2012 roku, a pierwsza aplikacja azotu N miała miejsce wiosną 2013 roku.
Pobieranie próbek gleby i korzeni przeprowadzono 15 września 2014 roku. Prace opisane poniżej wykonano na polu, na którym zastosowano randomizowany układ dzielonych poletek z trzema powtórzeniami (Rysunek 1). Średnia głębokość sekwencjonowania wszystkich próbek dla endosfery wynosiła: 4871 ± 5711 (średnia ± SD), dla ryzosfery: 40726 ± 14684, dla gleby: 38184 ± 9043. Jednym z największych źródeł zmienności w tych eksperymentach, przy zastosowaniu opisanych metod, jest różnica w społecznościach mikrobiologicznych występujących pomiędzy typami próbek (Rysunek 2). W tym reprezentatywnym zbiorze danych ryzosfera i gleba wydają się być bardziej podobne pod względem składu do siebie nawzajem niż do endosfery (Rysunek 2A). Jednakże odnotowano również wysoce istotne (p = 0,001) różnice w składzie społeczności mikrobiologicznych pomiędzy ryzosferą a glebą (Rysunek 2B). Całkowita zmienność uwzględniona w tych eksperymentach analizowanych według typu próbki wyniosła 26%.
Analiza różnorodności alfa wykazała, że społeczności mikrobiologiczne w endosferze charakteryzowały się niższą różnorodnością próbek w porównaniu z glebą i ryzosferą (Rysunek 3). Jedyne istotne różnice w różnorodności między gatunkami traw w dowolnym przedziale wystąpiły pomiędzy próbkami endosfery rozsadnicy wielkiej i trawy preriowej, przy czym trawa preriowa wykazywała istotnie wyższą różnorodność gatunkową mikroorganizmów (Rysunek 3). Analiza względnej liczebności (Rysunek 4) wskazuje na dominację Proteobacteria, a następnie Actinobacteria we wszystkich typach próbek. Gleba i ryzosfera są również zdominowane przez Acidobacteria i Chloroflexi, podczas gdy w endosferze odnotowano większą względną liczebność Bacteriodetes.
W tym eksperymencie rośliny hodowano przy zastosowaniu dwóch różnych dawek nawozu azotowego (N), w związku z czym przeanalizowano dane w celu ustalenia, czy wystąpiły efekty zabiegowe. Efekty zabiegowe odpowiadały za 12% całkowitej zmienności, jednak nie były one istotnie różne, mimo że w ordynacji oba zabiegi wyglądają odmiennie (Rysunek 5). Podkreśla to znaczenie analiz statystycznych dla tego typu zbiorów danych w przeciwieństwie do inspekcji wizualnej lub ocen jakościowych.
Różnice w mikrobiomie tkanek roślinnych i gleby pod wpływem roślin zostały zwizualizowane przy użyciu ograniczonej metody ordynacji. Różnice statystyczne określono za pomocą analizy PERMANOVA, aby sprawdzić, czy konkretne zmienne, takie jak gatunek, prowadzą do istotnych różnic w składzie społeczności mikrobiologicznych pomiędzy próbkami. Gdy wszystkie typy próbek analizowano wspólnie, stwierdzono wysoce istotną różnicę w składzie społeczności mikrobiologicznych wynikającą z gatunku rośliny (Rycina 6). W tym eksperymencie stopień zmienności wyjaśniony przez gatunek rośliny wyniósł 6,7%. Na koniec każdy typ próbki analizowano indywidualnie, aby ustalić, który z nich może odpowiadać za istotny wpływ gatunku rośliny. Wyłącznie w endosferze odnotowano wysoce istotną różnicę (p = 0,001) między składem społeczności mikrobiologicznych różnych gatunków roślin (Rycina 7). W pozostałych typach próbek wpływ gatunku nie był istotny podczas analizy indywidualnej. W endosferze procent zmienności wynikający z gatunku wyniósł 27%, podczas gdy w ryzosferze (18%) i glebie (15%) był on niższy. Podkreśla to dodatkowo znaczenie indywidualnej analizy każdego typu tkanki.

Rycina 1: Przykład projektu pola doświadczalnego. Projekt pola doświadczalnego ilustrujący randomizowany układ bloków kompletnych w trzech powtórzeniach na stanowisku badawczym zlokalizowanym w University of Nebraska-Lincoln Eastern Nebraska Research and Extension Center w pobliżu Mead, NE. Pełny opis stanowiska znajduje się w sekcji Wyniki. N1 oznacza niską dawkę azotu (56 kg N ha-1 mocznika), a N2 (112 kg N ha-1 mocznika) wyższą dawkę azotu, które zostały zastosowane. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Analiza beta-różnorodności porównująca skład mikrobiologiczny w różnych typach próbek, w tym endosfery, ryzosfery i gleby, pobranych z traw wieloletnich w 2014 roku. Analizę przeprowadzono za pomocą skryptu Python w QIIME1.9.1 w celu utworzenia macierzy dystansu Bray-Curtis. Analiza głównych współrzędnych (PCoA) oparta na macierzy dystansu Bray-Curtis została zwizualizowana w programie RStudio. PCoA1 i PCoA2 wskazują odpowiednio pierwszą i drugą największą wariancję wyjaśnioną przez analizę PCoA. Przeprowadzono statystyczną analizę PERMANOVA w celu określenia istotności różnic między typami próbek; wartość p jest wyświetlana w prawym górnym rogu. Każdy symbol na rysunkach reprezentuje całą społeczność mikrobiologiczną dla danej próbki. (A) Typy próbek z endosfery, ryzosfery i gleby zostały przeanalizowane wspólnie. Wszystkie 87 próbek zostało zrzadkowanych do 486 sekwencji na próbkę. (B) Próbki z ryzosfery i gleby zostały przeanalizowane wspólnie. Wszystkie 59 próbek zostało zrzadkowanych do 8231 sekwencji. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Analiza różnorodności alfa z wykorzystaniem indeksu Shannona dla każdego gatunku w endosferze, ryzosferze i glebie. Analizę przeprowadzono za pomocą skryptu Python w programie QIIME1.9.1. Rzadkość (rarefakcję) wykonano odpowiednio dla próbek z endosfery, ryzosfery i gleby przy 486, 17154 i 8231 sekwencjach na próbkę. Prostokąty reprezentują 25. i 75. percentyl (pierwszy i trzeci kwartyl). Pozioma linia wewnątrz prostokąta oznacza medianę, a czerwony plus wskazuje średnią. Wąsy pokazują zakres danych z wyłączeniem wartości odstających (przedstawionych jako czarne kropki), które wykraczały poza 1,5-krotność rozstępu międzykwartylowego (n = 6 dla każdej próbki, z wyjątkiem mieszanki sideoats grama, gdzie n = 5). Indeks Shannona dla wszystkich pięciu gatunków w endosferze był niższy niż w ryzosferze i glebie. Do określenia istotności różnic między gatunkami zastosowano nieparametryczny test suma rang Wilcoxona; nad prostokątami zaznaczono jedynie różnice istotne statystycznie między gatunkami. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 4: Względna liczebność na poziomie typu w endosferze, ryzosferze i glebie. Przebadano próbki w celu porównania liczebności typów mikroorganizmów w różnych rodzajach próbek (n = 29 dla każdego rodzaju próbki). Analizę przeprowadzono za pomocą skryptu Python w programie QIIME1.9.1 na podstawie tabeli OTU. Różne kolory na wykresie kołowym oznaczają typy. Procent wskazuje względną liczebność każdego typu w każdym rodzaju próbki. Informacje o typach zostały przypisane przy użyciu klasyfikatora Ribosomal Database Project (RDP)25. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5: Analiza z zastosowaniem traktowania jako czynnika ograniczającego pomiędzy wszystkimi typami próbek. Przeprowadzono kanoniczną analizę głównych współrzędnych (CAP), aby ustalić, czy występują różnice w składzie społeczności drobnoustrojów pomiędzy wariantami traktowania. Dla każdego traktowania azotem (N), n = 42 dla N1 (56 kg N ha-1) oraz n = 45 dla N2 (112 kg N ha-1). Macierz niezależności Braya-Curtisa wygenerowano za pomocą skryptu python w QIIME1.9.1. Analizę CAP opartą na macierzy niezależności Braya-Curtisa wykonano w RStudio, przyjmując traktowanie jako czynnik ograniczający. Przeprowadzono analizę PERMANOVA, aby określić, czy różnice między wariantami traktowania były istotne; wartość p przedstawiono w prawym górnym rogu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 6: Analiza z gatunkiem rośliny jako czynnikiem ograniczającym pomiędzy wszystkimi typami próbek. Przeprowadzono analizę w celu ustalenia, czy występują różnice w składzie społeczności drobnoustrojów pomiędzy gatunkami roślin we wszystkich typach próbek. Ordynacja głównych współrzędnych oraz analiza CAP dla wszystkich typów próbek (endosfera, ryzosfera i gleba) zostały wykonane przy użyciu macierzy niepodobieństwa Bray-Curtis. Macierz niepodobieństwa Bray-Curtis została wygenerowana za pomocą skryptu Python w programie QIIME1.9.1. Analizę CAP opartą na macierzy niepodobieństwa Bray-Curtis przeprowadzono w programie RStudio, przyjmując gatunek rośliny jako czynnik ograniczający. Przeprowadzono statystyczną analizę PERMANOVA w celu określenia istotności różnic pomiędzy gatunkami roślin; wartość P jest przedstawiona w prawym górnym rogu. Każdy symbol na rycinach reprezentuje całą społeczność drobnoustrojów dla danej próbki. n = 18 dla każdego gatunku we wszystkich typach próbek, z wyjątkiem n = 15 dla mieszanki sideoats grama. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 7: Przykład analizy CAP z wykorzystaniem gatunków jako czynnika ograniczającego dla każdego typu próbki z osobna. Ordynacja głównych współrzędnych i analiza CAP dla każdego typu próbki (endosfera, ryzosfera i gleba) z wykorzystaniem macierzy niezróżnicowania Bray-Curtis. Każdy typ próbki został poddany rzedukcji (rarefakcji) do odpowiednio 486, 17154 i 8231 odczytów na próbkę w endosferze, ryzosferze i glebie. Jako czynnik ograniczający ordynację wykorzystano gatunki. Przeprowadzono analizę statystyczną PERMANOVA w celu określenia istotności różnic między gatunkami roślin w każdym typie próbki, a wartość p przedstawiono w prawym górnym rogu. Każdy symbol na rysunku reprezentuje całą społeczność mikrobiologiczną dla danej próbki. Wielkość próby wynosi n = 29 dla każdego typu próbki, n = 6 dla każdego gatunku rośliny w każdym typie próbki, z wyjątkiem mieszanki sideoats grama (n = 5). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.