$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Aby zbadać możliwość zróżnicowania społeczności mikrokosmosów glebowych, porównaliśmy społeczności mikroorganizmów w mikrokosmosach kompostowych przygotowanych przez wzbogacenie tej samej początkowej gleby, kompostu klasy przemysłowej dostępnego w mieście Berkeley w Kalifornii, z różnymi produktami: jabłkami, papryką, pomarańczami lub ziemniakami (każdy zestaw w trzech egzemplarzach). Następnie porównaliśmy społeczności mikrobiologiczne w każdym środowisku kompostowym z mikrobiomem jelitowym C. elegans typu dzikiego wyhodowanego w odpowiednim mikrokosmosie. Analizę przeprowadzono na próbkach DNA wyekstrahowanych od około 500 sterylizowanych powierzchniowo dorosłych osób na mikrokosmos oraz z 250 mg próbek kompostu z odpowiednich mikrokosmosów.
Charakterystyka mikrobiomów środowiskowych, glebowych i jelitowych robaków opierała się na sekwencjonowaniu nowej generacji regionu V4 bakteryjnego genu 16S rRNA. Przygotowanie biblioteki sekwencjonowania zostało osiągnięte przy użyciu standardowych zestawów i przeprowadzone zgodnie z instrukcjami producenta, przy czym sekwencjonowanie przeprowadzono na komercyjnym sekwencerze (patrz tabela materiałów). Demultipleksowane sekwencje zostały przetworzone przy użyciu DADA2, przypisano im taksonomię na podstawie referencyjnej bazy danych SILVA v132 i przeanalizowano za pomocą phyloseq16,17,18 (patrz plik uzupełniający 1, rysunek uzupełniający S1, rysunek uzupełniający S2, rysunek uzupełniający S3, tabela uzupełniająca S1 i tabela uzupełniająca S2 szczegółowy opis sekwencjonowania i analizy; pełny potok obliczeniowy jest dostępny w GitHub [https://github.com/kennytrang/CompostMicrocosms]). Surowe dane są dostępne w NCBI Sequence Read Archive (identyfikator bioprojektu 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 uważane za część rdzenia mikrobiomu jelitowego C. elegans13, takie jak Rhizobiaceae, Burkholderiaceae i Bacillaceae. Enterobacteriaceae i Pseudomonadaceae, które wcześniej uznano za dominujących członków, tym razem stanowiły mniejszość, ale nadal były wzbogacone (2-10-krotnie) w porównaniu z odpowiednimi środowiskami glebowymi. Porównania oparte zarówno na odległościach nieważonych, jak i ważonych UniFrac19,20 wykazały dobrą odtwarzalność wśród mikrokosmicznych triplikatów wzbogaconych tym samym produktem, na co wskazuje ścisłe grupowanie. W przeciwieństwie do tego, środowiskowe mikrobiomy glebowe wzbogacone różnymi produktami grupowały się od siebie, wykazując zdolność do dywersyfikacji początkowej społeczności drobnoustrojów poprzez dodawanie różnych produktów (Rysunek 2).
W porównaniach mikrobiomów jelitowych robaków i społeczności środowiskowych, analiza głównych współrzędnych (PCoA) z nieważonymi lub ważonymi odległościami UniFrac wykazała wyraźne grupowanie mikrobiomów jelitowych robaków w stosunku do ich odpowiednich środowisk dla każdego typu mikrokosmosu (Rysunek 2). Podczas gdy PCoA oparte na nieważonych odległościach UniFrac nie rozróżniało mikrobiomów glebowych i robaków (Rysunek 2A), grupowanie oparte na odległościach ważonych ujawniło wyraźne oddzielenie mikrobiomów jelitowych i kompostowych robaków (Rysunek 2B). Wyniki te wspierają proces, w którym filtrowanie gospodarza wpływa na dostępność środowiskową w celu kształtowania mikrobiomu jelitowego, który nie jest całkowicie odrębny od swojego źródła środowiskowego pod względem obecności taksonów, ale moduluje ich liczebność poprzez wzbogacanie podzbioru dostępnych taksonów, co ostatecznie prowadzi do powstania mikrobiomu jelitowego robaka rdzeniowego współdzielonego przez robaki hodowane w różnych środowiskach.

Rysunek 2: Mikrobiomy jelitowe robaków grupujące się z dala od swoich zróżnicowanych środowisk mikrobiologicznych. Skład mikrobiomu określono za pomocą sekwencjonowania 16S, a zbiorowiska z mikrokosmosów wzbogaconych o wyznaczony produkt lub z wyhodowanych w nich robaków grupowano za pomocą PCoA w oparciu o (A) nieważone lub (B) ważone odległości UniFrac. Pokazane osie to te, które wyjaśniają największą zmienność w składzie zbiorowiska między próbkami (N = 3 dla każdego typu mikrokosmosu). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Plik uzupełniający 1: Sekwencjonowanie nowej generacji i analiza danych. Poniżej przedstawiono kroki przygotowania biblioteki, sekwencjonowania w laboratorium i analizy danych. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S1: Przykład wykresu kontroli jakości dla odwrotnych odczytów z jednej próbki. Oś X (cykl) pokazuje pozycję nukleotydu wzdłuż odczytanej sekwencji. Lewa oś Y pokazuje wynik jakości. Mapa cieplna w skali szarości reprezentuje częstotliwość oceny jakości w każdej pozycji nukleotydu; Zielona linia przedstawia medianę wyniku jakości w każdej pozycji nukleotydu; Górna pomarańczowa linia przedstawia kwartyle rozkładu Wyniku Jakości; dolna czerwona linia przedstawia procent odczytów sekwencji, które rozszerzyły tę pozycję nukleotydu (prawa oś Y, tutaj 100%). Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S2: Poziomy błędów dla różnych próbek. Częstość występowania błędów w różnych próbkach (czarne kropki) powinna zmniejszać się wraz ze wzrostem wyniku jakości dla każdego możliwego podstawienia par zasad, odzwierciedlając oczekiwany trend. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S3: Przykład PCoA oparty na ważonych odległościach UniFrac. Nazwy grup pokazane w legendzie reprezentują produkty używane do wzbogacania kompostu używanego w różnych mikrokosmosach. Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S1: Sekwencyjne filtrowanie sekwencji. ZALiczba odczytów sekwencji przed filtrowaniem. b-d Każda kolumna reprezentuje liczbę odczytów sekwencji pozostałych po kroku filtracji: odfiltrowanie odczytów o niskiej jakości (krok 2.5), algorytm odszumiania wykonywany przez dada() (krok 2.8), scalanie odczytów do przodu i do tyłu (krok 2.9) i usuwanie chimer (krok 2.11). Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S2: Tabela metadanych. Kliknij tutaj, aby pobrać ten plik.