Aby zbadać możliwość zróżnicowania społeczności w mikrokosmosach glebowych, porównaliśmy społeczności mikrobowe w mikrokosmosach kompostowych przygotowanych poprzez wzbogacenie tej samej gleby początkowej — kompostu przemysłowego dostępnego w mieście Berkeley w Kalifornii — różnymi produktami roślinnymi: jabłkami, papryką, pomarańczami lub ziemniakami (każdy zestaw w trzech powtórzeniach). Następnie porównaliśmy społeczności mikrobowe każdego środowiska kompostowego z mikrobiomem jelitowym dzikiego typu C. elegans hodowanego w odpowiednim mikrokosmosie. Analizę przeprowadzono z próbkami DNA wyekstrahowanego z około 500 wysterylizowanych powierzchniowo osobników dorosłych na mikrokosmos oraz z próbek 250 mg kompostu z odpowiednich mikrokosmosów.
Charakterystyka mikrobiomów gleby środowiskowej oraz jelit dżdżownic opierała się na sekwencjonowaniu nowej generacji regionu V4 genu 16S rRNA bakterii. Przygotowanie bibliotek do sekwencjonowania przeprowadzono przy użyciu standardowych zestawów zgodnie z instrukcjami producentów, a samo sekwencjonowanie wykonano na komercyjnym sekwenatorze (patrz Tabela materiałów). Zdemultipleksowane sekwencje przetworzono za pomocą programu DADA2, przypisano im taksonomię na podstawie bazy referencyjnej SILVA v132 i przeanalizowano za pomocą pakietu phyloseq16,17,18 (szczegółowy opis sekwencjonowania i analizy znajduje się w Pliku uzupełniającym 1, Rysunku uzupełniającym S1, Rysunku uzupełniającym S2, Rysunku uzupełniającym S3, Tabeli uzupełniającej S1 oraz Tabeli uzupełniającej S2; pełny schemat obliczeniowy jest dostępny w serwisie GitHub [https://github.com/kennytrang/CompostMicrocosms]). Surowe dane są dostępne w bazie NCBI Sequence Read Archive (Bioproject ID PRJNA856419).
Średnio uzyskano 73 220 sekwencji na próbkę. Sekwencje te reprezentują 15 027 wariantów sekwencji amplikonów (ASV), obejmujących 27 typów i 216 rodzin, w tym rodziny uznawane za część rdzennego mikrobiomu jelitowego C. elegans13, takie jak Rhizobiaceae, Burkholderiaceae oraz Bacillaceae. Enterobacteriaceae i Pseudomonadaceae, które wcześniej uznano za dominujące elementy, tym razem stanowiły mniejszość, lecz nadal były wzbogacone (2–10-krotnie) w porównaniu do ich odpowiednich środowisk glebowych. Porównania oparte na nieważonych i ważonych dystansach UniFrac19,20 wykazały dobrą powtarzalność między triplikatami mikrokosmosów wzbogaconych tym samym produktem roślinnym, co wskazuje bliskie grupowanie się próbek. Z kolei mikrobiomy glebowe z różnych środowisk, wzbogacone różnymi produktami, grupowały się oddzielnie, co demonstruje możliwość zróżnicowania początkowej społeczności mikrobowej poprzez dodanie różnych produktów roślinnych (Rycina 2).
W porównaniach mikrobiomów jelit dżdżownic i społeczności środowiskowych, analiza głównych współrzędnych (PCoA) z zastosowaniem nieważonych lub ważonych odległości UniFrac wykazała wyraźne grupowanie się mikrobiomów jelit dżdżownic oddzielnie od ich odpowiednich środowisk dla każdego typu mikrokosmosu (Rysunek 2). Podczas gdy PCoA oparta na nieważonych odległościach UniFrac nie odróżniła mikrobiomów glebowych od mikrobiomów dżdżownic (Rysunek 2A), grupowanie oparte na odległościach ważonych ujawniło wyraźną separację mikrobiomów jelit dżdżownic i kompostu (Rysunek 2B). Wyniki te wspierają proces, w którym filtracja przez gospodarza oddziałuje na dostępność środowiskową, kształtując mikrobiom jelitowy, który nie jest całkowicie odrębny od swojego źródła środowiskowego pod względem obecności taksonów, lecz moduluje ich obfitość poprzez wzbogacanie o podzbiór dostępnych taksonów, co ostatecznie skutkuje powstaniem rdzennego mikrobiomu jelit dżdżownic, wspólnego dla osobników hodowanych w różnych środowiskach.

Rysunek 2: Mikrobiomy jelit dżdżownic grupujące się oddzielnie od odpowiadających im środowisk mikrobiologicznych zdywersyfikowanych produktami roślinnymi. Skład mikrobiomu określono za pomocą sekwencjonowania 16S, a społeczności z mikrokosmosów wzbogaconych wyznaczonymi produktami roślinnymi lub z dżdżownic w nich hodowanych zgrupowano za pomocą PCoA w oparciu o (A) nieważone lub (B) ważone odległości UniFrac. Przedstawione osie to te, które wyjaśniają największą zmienność w składzie społeczności między próbkami (N = 3 dla każdego typu mikrokosmosu). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Plik uzupełniający 1: Sekwencjonowanie nowej generacji i analiza danych. Przedstawiono tutaj etapy przygotowania bibliotek, sekwencjonowania w laboratorium oraz analizy danych. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S1: Przykład wykresu kontroli jakości dla odczytów odwrotnych z jednej próbki. Oś X (cykl) przedstawia pozycję nukleotydu wzdłuż odczytu sekwencji. Lewa oś Y przedstawia wynik jakości. Skala szarości mapy ciepła reprezentuje częstotliwość występowania wyniku jakości w każdej pozycji nukleotydu; zielona linia obrazuje medianę wyniku jakości w każdej pozycji nukleotydu; górna pomarańczowa linia obrazuje kwartyle rozkładu wyniku jakości; dolna czerwona linia obrazuje procent odczytów sekwencji, które sięgnęły danej pozycji nukleotydu (prawa oś Y, tutaj 100%). Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S2: Wskaźniki błędów dla różnych próbek. Częstotliwość błędów w różnych próbkach (czarne kropki) powinna maleć wraz ze wzrostem wyniku jakości (quality score) dla każdej przedstawionej możliwej substytucji par zasad, co odzwierciedla oczekiwaną tendencję. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S3: Przykład PCoA opartego na ważonych dystansach UniFrac. Nazwy grup przedstawione w legendzie reprezentują produkty użyte do wzbogacenia kompostu w różnych mikrokosmosach. Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S1: Sekwencyjne filtrowanie sekwencji. aLiczba odczytów sekwencji przed filtrowaniem. b-dKażda kolumna przedstawia liczbę odczytów sekwencji pozostałych po etapie filtracji: odfiltrowaniu odczytów o niskiej jakości (krok 2.5), algorytmie odszumiającym wykonanym przez dada() (krok 2.8), łączeniu odczytów forward i reverse (krok 2.9) oraz usuwaniu chimer (krok 2.11). Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S2: Tabela metadanych. Kliknij tutaj, aby pobrać ten plik.