$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Reprezentatywne wyniki przedstawione w tym manuskrypcie pochodzą z terenu założonego w 2012 roku na farmie University of Nebraska-Lincoln Agriculture Research Division w pobliżu Mead, NE. Przed eksperymentem miejsce to było zarządzane jako płodozmian kukurydziano-soi. Stanowisko badawcze znajdowało się na trzech różnych typach gleb, ale dane analizowano tak, jakby wszystkie zmiany w mierzonych właściwościach gleby były spowodowane nałożonymi zabiegami.
Teren pola zawierał dwa, czyste, stanowiska proso rózgowatego (P. virgatum cv Liberty) i dużego bluestem (A. gerardii), a także mieszankę traw o niskiej różnorodności, zawierającą trawę dużą, trawę indyjską (S. nutans) i 'Butte' sideoats grama (B. curtipendula). Trzy poletka trawiaste w ciepłym sezonie były w losowym projekcie kompletnego bloku, który został trzykrotnie powtórzony. Na trzech różnych poletkach trawiastych zagnieżdżono dwa zabiegi nawożenia azotem (N), które wyniosły 56 (N1) i 112 (N2) kg N ha-1 zastosowanego mocznika. W momencie pobierania próbek mikrobiomu pod koniec sezonu wegetacyjnego gleba zawierała 8,0 ± 1,1 (średnio ± SD) ppm azotanów na poletkach nawożonych 112 kg N ha-1 i 6,8 ± 0,7 (średnia + SD) ppm azotanów na poletkach nawożonych 56 kg N ha-1. Poletka były nawożone raz w roku. Jako główne poletka (8000 m2) wyznaczono poletka trawiaste w ciepłym sezonie, a poletka N podzielone (4000 m2). Big bluestem został wysiew w proporcji 50:50 odmian 'Bonanza' i 'Goldmine', a Indiangrass został wysiew w proporcji 50:50 odmian 'Scout' i 'Warrior'. Działki zostały obsadzone w 2012 roku, a pierwsza aplikacja N miała miejsce wiosną 2013 roku.
Pobieranie próbek gleby i korzeni zostało przeprowadzone 15 września 2014 roku. Opisane poniżej prace zostały przeprowadzone na polu, które zostało skonfigurowane jako losowy projekt z podzielonym wykresem z trzema powtórzeniami (Rysunek 1). Średnia głębokość sekwencjonowania wszystkich próbek dla endosfery wynosiła: 4871 ± 5711 (średnia ± SD), ryzosfera: 40726 ± 14684, gleba: 38184 ± 9043. Jednym z największych źródeł zmienności w tych eksperymentach, przy użyciu opisanych metod, jest różnica w zbiorowiskach drobnoustrojów występujących między typami próbek (Rysunek 2). W tym reprezentatywnym zbiorze danych ryzosfera i gleba wydają się być bardziej podobne pod względem składu niż endosfera (Rysunek 2A). Stwierdzono jednak również bardzo istotne (p = 0,001) różnice w składzie zbiorowisk mikroorganizmów między ryzosferą a glebą (Rysunek 2B). Całkowita zmienność uwzględniona w tych eksperymentach analizowanych w zależności od rodzaju próby wyniosła 26%.
Analiza różnorodności alfa wykazała, że społeczności mikroorganizmów w endosferze były mniej zróżnicowane pod względem próbki w porównaniu z glebą i ryzosferą (Rysunek 3). Jedyne znaczące różnice w różnorodności między gatunkami traw w dowolnym przedziale dotyczyły próbek endosfery dużego bluestem i proso rózgowatego, przy czym proso rózgowate miało znacznie wyższą różnorodność gatunków mikrobiologicznych (Rysunek 3). Analiza względnej liczebności (Rysunek 4) podkreśla dominację Proteobacteria, a następnie Actinobacteria we wszystkich typach próbek. Gleba i ryzosfera są również zdominowane przez Acidobacteria i Chloroflexi, podczas gdy endosfera miała większą względną obfitość bakteriodetów.
W tym eksperymencie rośliny były uprawiane z dwiema różnymi ilościami nawozu N i dlatego przeanalizowaliśmy dane, aby ustalić, czy wystąpiły efekty leczenia. Efekty leczenia stanowiły 12% całkowitej zmienności, ale nie różniły się istotnie, chociaż w ordynacji obie terapie wyglądają inaczej (Ryc. 5). Podkreśla to znaczenie analiz statystycznych dla tych zbiorów danych, a nie oględzin lub ocen jakościowych.
Różnice w mikrobiomie tkanek roślinnych i gleby pod wpływem roślin zostały zobrazowane za pomocą metody ograniczonego ordynacji. Różnice statystyczne określono za pomocą analizy PERMANOVA w celu zbadania, czy określone zmienne, takie jak gatunki, skutkują znacząco różnym składem zbiorowisk drobnoustrojów między próbkami. Gdy wszystkie typy próbek zostały przeanalizowane razem, stwierdzono bardzo istotną różnicę w składzie zbiorowiska mikroorganizmów ze względu na gatunki roślin (Rysunek 6). W tym eksperymencie wielkość zmienności uwzględnionej w zależności od gatunku rośliny wyniosła 6,7%. Na koniec każdy typ próbki został przeanalizowany indywidualnie w celu określenia, które z typów próbek mogą napędzać znaczący wpływ na gatunki roślin. Jedynie w endosferze zaobserwowano bardzo istotną różnicę (p = 0,001) między składem zbiorowisk mikrobiologicznych różnych gatunków roślin (Ryc. 7). W innych typach próbek wpływ gatunkowy nie był istotny, gdy analizowano go indywidualnie. W endosferze procentowa zmienność ze względu na gatunki wynosiła 27%, podczas gdy była niższa w ryzosferze (18%) i glebie (15%). To dodatkowo podkreśla, jak ważna jest indywidualna analiza każdego typu tkanki.

Rysunek 1: Przykład projektu pola doświadczalnego. Eksperymentalny projekt terenowy ilustrujący losowy kompletny projekt blokowy w trzech egzemplarzach terenu znajdującego się w University of Nebraska-Lincoln Eastern Nebraska Research and Extension Center w pobliżu Mead, NE. Pełny opis witryny znajduje się w sekcji Wyniki. N1 to niska (56 kg N ha-1 mocznika), a N2 (112 kg N ha-1 mocznika) to wyższa dawka azotu, która została zastosowana. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Analiza różnorodności beta porównująca skład mikrobiologiczny w różnych typach próbek, w tym endosferze, ryzosferze i glebie, z próbek wieloletnich traw w 2014 r. Analizę przeprowadzono przy użyciu skryptu Pythona w QIIME1.9.1 w celu wytworzenia macierzy odmienności Braya-Curtisa. Analiza głównych współrzędnych (PCoA) w oparciu o macierz odmienności Braya-Curtisa została zwizualizowana w RStudio. PCoA1 i PCoA2 wskazują pierwszą i drugą co do wielkości wariancję wyjaśnioną przez analizę PCoA. Przeprowadzono analizę statystyczną PERMANOVA w celu określenia istotności między typami próbek, a wartość p jest pokazana w prawym górnym rogu. Każdy symbol na rysunkach reprezentuje całą społeczność drobnoustrojów dla każdej próbki. (A) Typy próbek endosfery, ryzosfery i gleby analizowano razem. Wszystkie 87 próbek rozrzedzono do 486 sekwencji na próbkę. (B) Próbki ryzosfery i gleby analizowano razem. Wszystkie 59 próbek rozrzedzono za pomocą 8231 sekwencji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Analiza różnorodności alfa z wykorzystaniem indeksu Shannona dla każdego gatunku w endosferze, ryzosferze i glebie. Analiza została przeprowadzona przy użyciu skryptu Python w QIIME1.9.1. Rozrzedzenie przeprowadzono odpowiednio dla typów próbek endosfery, ryzosfery i gleby, z 486, 17154 i 8231 sekwencjami na próbkę. Pola wskazują 25. i 75. percentyl (pierwszy i trzeci kwartyl). Pozioma linia w polu oznacza medianę, a czerwony plus pokazuje średnią. Wąsy pokazują zakres danych z wyłączeniem wartości odstających (które są pokazane jako czarne kropki), które spadły ponad 1,5 raza w stosunku do przedziału międzykwartylowego (n = 6 dla każdej próby, z wyjątkiem mieszanki Sideoats grama, gdzie n = 5). Indeks Shannona wszystkich pięciu gatunków w endosferze był niższy zarówno niż ryzosfera, jak i gleba. Do określenia istotności między gatunkami użyto nieparametrycznego testu sumy rang Wilcoxona, a na wierzchu pól pokazano tylko znaczące różnice między gatunkami. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Względna liczebność na poziomie gromady w endosferze, ryzosferze i glebie. Próbki analizowano w celu porównania liczebności gromad drobnoustrojów między różnymi typami próbek (n = 29 dla każdego typu próbki). Analizę przeprowadzono przy użyciu skryptu Python w QIIME1.9.1 z tabeli OTU. Różne kolory wewnątrz wykresu kołowego oznaczają typ. Wartość procentowa wskazuje względną liczebność każdej gromady w każdym typie próbki. Informacje o typie zostały oznaczone przy użyciu klasyfikatora projektu bazy danych rybosomalnych (RDP)25. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Analiza z zastosowaniem traktowania jako czynnika ograniczającego między wszystkimi typami próbek. Przeprowadzono analizę kanoniczną współrzędnych głównych (CAP) w celu ustalenia, czy występują różnice w składzie społeczności drobnoustrojów między traktowanymi lekami. Dla każdego zabiegu N n = 42 dla N1 (56 kg N ha-1) i n = 45 dla N2 (112 kg N ha-1). Macierz odmienności Braya-Curtisa została wygenerowana przy użyciu skryptu Pythona w QIIME1.9.1. Analizę CAP opartą na macierzy odmienności Braya-Curtisa przeprowadzono poprzez ograniczenie leczenia jako czynnika w RStudio. Analizę PERMANOVA przeprowadzono w celu określenia, czy różnice w leczeniu były istotne, a wartość p jest pokazana w prawym górnym rogu. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Analiza z wykorzystaniem gatunków roślin jako czynnika ograniczającego wszystkie typy próbek. Przeprowadzono analizę w celu ustalenia, czy występują różnice w składzie zbiorowisk mikrobiologicznych między gatunkami roślin we wszystkich typach próbek. Uporządkowanie głównych współrzędnych i analiza CAP wszystkich typów próbek (endosfera, ryzosfera i gleba) przeprowadzono przy użyciu macierzy odmienności Braya-Curtisa. Macierz odmienności Braya-Curtisa została wygenerowana przy użyciu skryptu Pythona w QIIME1.9.1. Analizę CAP w oparciu o macierz odmienności Braya-Curtisa przeprowadzono poprzez ograniczenie gatunków roślin jako czynnika w RStudio. Przeprowadzono analizę statystyczną PERMANOVA w celu określenia istotności między gatunkami roślin, a wartość P jest pokazana w prawym górnym rogu. Każdy symbol na rysunkach reprezentuje całą społeczność drobnoustrojów dla tej 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 zobaczyć większą wersję tego rysunku.

Rysunek 7: Przykład analizy CAP z wykorzystaniem gatunków jako czynnika ograniczającego dla każdego rodzaju próbki z osobna. Uporządkowanie głównych współrzędnych i analiza CAP każdego typu próbki (endosfera, ryzosfera i gleba) przy użyciu macierzy odmienności Braya-Curtisa. Każdy typ próbki został rozrzedzony do 486, 17154 i 8231 odczytów na próbkę odpowiednio w endosferze, ryzosferze i glebie. Gatunek został użyty jako czynnik ograniczający święcenia. Przeprowadzono analizę statystyczną PERMANOVA w celu określenia istotności między gatunkami roślin w każdym typie próby, a wartość p jest pokazana w prawym górnym rogu. Każdy symbol na rysunku reprezentuje całą społeczność drobnoustrojów dla każdej próbki. Wielkość próby wynosi n = 29 dla każdego rodzaju 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 zobaczyć większą wersję tego rysunku.