$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Celem tego protokołu jest zapewnienie ogólnego potoku, od eksperymentów polowych po analizy bioinformatyczne, który określa ilościowo potencjalną wrażliwość organizmów na herbicyd glifosat. W eksperymencie 2 średnie stężenie glifosatu w paszy dla przepiórek wynosiło 164 mg/kg, a średnie stężenie glifosatu w próbkach odchodów (moczu i kału łącznie) wynosiło 199 mg/kg. Ściółki zebrane od przepiórek karmionych paszą zanieczyszczoną GBP zawierały średnio 158 mg/kg, a ściółki kontrolne 0,17 mg/kg glifosatu (Tabela 3). W doświadczeniach polowych gatunki roślin różnie reagowały na pozostałości glifosatu w glebie (sekcja 1). Biomasa rzepaku i owsa i rzepiku była większa w glebach kontrolnych w porównaniu z glebami, w których zastosowano GBP. Wydaje się jednak, że bób i ziemniaki skorzystały na traktowaniu GBP pod koniec sezonu wegetacyjnego15. Glifosat w oborniku drobiowym zmniejszył wzrost roślin traw (Festuca pratensis) i truskawek (Fragaria x vescana) (sekcja 1). Analizy mikrobioty uzyskane w ramach doświadczeń terenowych nie zostały jeszcze w pełni przeanalizowane i nie zostały przedstawione w niniejszym dokumencie (sekcja 2). Wyniki tego protokołu, odczytywane bezpośrednio (jak pokazano w sekcjach 3 i 4) lub pośrednio (sekcja 5), stanowią miarę proporcji potencjalnie wrażliwych i odpornych organizmów na glifosat w zestawie danych (Rysunek 9). Zastosowanie tej metody zostało przetestowane z kolekcją sekwencji białek EPSPS pochodzących od gatunków drobnoustrojów rdzenia ludzkiego mikrobiomu jelitowego, które uzyskano z publicznych repozytoriów12. W badaniu 890 szczepów ze 101 najliczniej występujących gatunków bakterii przeanalizowano metodą EPSPSClass w celu ilościowego określenia proporcji wrażliwych i opornych bakterii. Wyniki pokazały, że 54% gatunków w rdzeniu ludzkiego mikrobiomu jelitowego jest potencjalnie wrażliwych na glifosat12. Trend ten obserwuje się również w większości świata prokariotycznego; Dodatkowo, u eukariontów (głównie roślin i grzybów) odsetek potencjalnie wrażliwych gatunków jest jeszcze wyższy12. Co więcej, wykorzystaliśmy tę metodę do ilościowego określenia zmian wrażliwości w białku EPSPS na poziomie mikroewolucyjnym (Rysunek 10) 14. Zidentyfikowaliśmy zmiany w statusie wrażliwości w 12 z 32 ściśle spokrewnionych analizowanych grup prokariontów (Tabela 4)14. W związku z tym ciągłe stosowanie GBP może powodować dysbiozę mikrobiologiczną (tj. brak równowagi wrażliwych i odpornych gatunków bakterii) w mikrobiomach roślinnych, zwierzęcych i glebowych. Ponadto wysunięto hipotezę, że wzrost liczby bakterii opornych na glifosat może sprzyjać mikrobiomom wielolekoopornym14,41,42. W związku z tym protokół ten rzuca światło na interpretację wszystkich tych scenariuszy, ponieważ metoda klasyfikacji EPSPS zapewnia bezpośrednie oszacowanie wewnętrznej wrażliwości mikrobiomów na glifosat. Ze względu na wewnętrzną wrażliwość białka EPSPS na glifosat jest filogenetycznie konserwowana14, możliwa jest ekstrapolacja wyników z istniejących zestawów danych na nieznane mikrobiomy (Rysunek 8).

Rysunek 1: Ogólny potok To jest ogólny potok do analizy wrażliwości na GBP od eksperymentów terenowych do analizy bioinformatycznej. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Eksperyment polowy 1 mający na celu sprawdzenie wpływu pozostałości GBP na mikroorganizmy związane z roślinami uprawnymi. Pole doświadczalne składa się naprzemiennie z 10 powierzchni kontrolnych i 10 powierzchni zabiegowych GBP (23 m x 1,5 m) z 1,5 m pasami buforowymi między powierzchniami. Dwa razy w roku od 2014 r. poletka GBP były traktowane handlową GBP (stężenie glifosatu 450 g L-1, dawka 6,4 L ha-1 na 5 l wody z kranu na poletko), a poletka kontrolne taką samą ilością wody z kranu bez glifosatu. Zabiegi wykonywano za pomocą ręcznie obsługiwanego zbiornika ciśnieniowego z plastikowym kapturem w końcówce zraszacza, aby chronić GBP przed rozprzestrzenianiem się poza poletka zabiegowe. Po dwutygodniowym okresie bezpieczeństwa po zastosowaniu GBP wysiano owies (Avena sativa), bób (Vicia faba) i rzepak (Brassica rapa subsp. oleifera), a na poletkach posadzono ziemniaki (Solanum tuberosum). Od rozpoczęcia eksperymentu w 2014 r. próbki mikrobioty z badanych roślin uprawnych, liści i korzeni pobierano kilkakrotnie. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Eksperyment polowy 2 testował konsekwencje pozostałości GBP w nawozie gnojowym dla dwóch wieloletnich upraw i związanej z nimi mikroflory. Ściółki zebrane z 12-miesięcznego eksperymentu wolierowego z przepiórkami japońskimi karmionymi paszą kontrolną lub zanieczyszczoną GBP wykorzystano jako nawóz do obornika w doświadczeniu polowym. Pole doświadczalne składało się z 18 poletek kontrolnych i 18 działek GBP (1 m x 1 m) ułożonych w siatkę szachownicy 6 x 6. Ściółki zostały rozłożone na polu doświadczalnym dwukrotnie, w sierpniu 2018 r. i maju 2019 r. (25 L/działkę). Poletka kontrolne nawożono ściółkami zebranymi z przepiórek karmionych paszą kontrolną, a poletka GBP ściółkami z przepiórek karmionych paszą zanieczyszczoną GBP. Pozostałości glifosatu w ściółkach kontrolnych wynosiły 0,17 mg/kg glifosatu, a w ściółce GBP ilość glifosatu wynosiła 158 mg/kg. We wrześniu 2018 r., około miesiąc po rozłożeniu pierwszych rabat, posadzono dwie endofity symbiotyczne (E+), dwie wolne od endofitów (E-) Festuca pratensis i dwie Fragaria x vescana. Pomiary wydajności i kondycji roślin, a także pobieranie próbek mikrobioty związanej z korzeniami i liśćmi przeprowadzono w dwóch kolejnych sezonach wegetacyjnych (2019 i 2020). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Analiza taksonów drobnoustrojów przy użyciu genu 16S rRNA/regionu ITS oraz wrażliwość mikrobiomów na glifosat przy użyciu genu EPSPS. (A) Analiza sekwencji 16S rRNA lub ITS w celu identyfikacji taksonów drobnoustrojów. (B) Analiza sekwencji EPSPS w celu określenia wrażliwości drobnoustrojów na glifosat (GS-wrażliwy na glifosat/odporny na GR-glifosat) Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 5: Algorytm do identyfikacji klasy sekwencji białek EPSPS. Dane wejściowe to sekwencja białka EPSPS w formacie FAPSA. Algorytm przeprowadza porównania ze znanymi markerami aminokwasowymi w referencyjnych sekwencjach białek, które określają potencjalną wrażliwość na glifosat. Algorytm został zaimplementowany na ogólnodostępnym serwerze WWW EPSPSClass29. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Podstawowe dane wejściowe i wyjściowe serwera WWW EPSPSClass. (A) Dane wejściowe: sekwencja białka EPSPS w formacie FASTA. (B) Wyjście 1 - tożsamość: frakcja markerów aminokwasowych obecnych w sekwencjach zapytania (klasy I-IV) i motywach (klasa III). (C) Dane wyjściowe 2 — tożsamość: wyrównania sekwencji zapytania i odwołań. (D) Dane wyjściowe 3 — wyrównanie parami sekwencji zapytania i odwołania. (E) Referencyjne sekwencje EPSPS: Vibrio cholerae (vcEPSPS, klasa I), Coxiella burnetii (cbEPSPS, klasa II), Brevundimonas vesicularis (bvEPSPS, klasa III), Streptomyces davawensis (sdEPSPS, klasa IV). (F) Linki do wykonywania dodatkowych wyszukiwań blastp i identyfikacji zachowanych domen Kliknij tutaj aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Dostęp do wstępnie obliczonych zestawów danych sekwencji EPSPS. Postępuj zgodnie ze wskazówkami na rysunku, aby uzyskać dostęp do wstępnie obliczonego zestawu danych sekwencji EPSPS. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 8: Przykład oszacowania potencjalnej wrażliwości w projektach mikrobiomu bez sekwencji EPSPS. W przykładzie użyto wartości z bazy danych Alignable Tight Genomic Clusters30, która zawiera sekwencje gatunków prokariotycznych. Hipotetyczne gatunki z projektu dotyczącego mikrobiomu to Staphylococcus aureus, Corynebacterium diphtheriae, Campylobacter jejuni, Chlamydia psittaci i Sulfolobus islandicus. Wskaźnik wrażliwości na glifosat oblicza się według wzoru Number_Sensitive_Sequences/Total_Number_Of_Sequences. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 9: Schemat interpretacji wyników z tego protokołu i hipotetyczne scenariusze ewolucyjne. (A) W mikrobiomie proporcja potencjalnej wrażliwości (na zielono) i odporności (na czerwono) bakterii wynosi około 50:50. Czarne kropki oznaczają gatunki drobnoustrojów niesklasyfikowane; W związku z tym ich wrażliwość na glifosat jest nieznana. W niektórych mikrobiomach odsetek wrażliwych bakterii jest nieco wyższy, jak w ludzkim mikrobiomie jelitowym12. (B) Z biegiem czasu stosowanie glifosatu może prowadzić do dysbiozy mikrobiologicznej (tj. braku równowagi w proporcji wrażliwych i opornych bakterii), co prowadzi do różnych hipotetycznych scenariuszy. (C) Hipotetyczny przypadek 1 (brak wyboru): Stosowanie glifosatu nie wpływa na mikrobiom; Dzięki temu proporcja wrażliwych i opornych bakterii pozostaje stała. (D) Hipotetyczny przypadek 2: Zastosowanie glifosatu usuwa z populacji bakterie wrażliwe na glifosat. Spekulujemy, że ten scenariusz może być zależny od dawki. (E) Hipotetyczny przypadek 3: Presja selekcyjna spowodowana stosowaniem glifosatu nasila mutacje w genie EPSPS, które zmieniają stan wrażliwości bakterii. W ten sposób cała populacja drobnoustrojów staje się odporna na glifosat. Co więcej, w tym scenariuszu może dojść do wzrostu liczby bakterii wielolekoopornych. (F) Hipotetyczny przypadek 4: stosowanie glifosatu zmienia skład niektórych gatunków bakterii, powodując brak równowagi w stosunku do opornych bakterii, podczas gdy niektóre gatunki bakterii pozostają niezmienione, prawdopodobnie z powodu dodatkowych mechanizmów oporności, takich jak pompy wypływowe lub nadekspresji genu EPSPS 13. Scenariusz ten może również prowadzić do wzrostu liczby bakterii opornych na glifosat, a także wzrostu oporności bakterii na dodatkowe antybiotyki. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 10: Rozkład przewidywanej wrażliwości na glifosat w całym drzewie gatunków. Wykresy kołowe wskazują odsetek gatunków, które są przypuszczalnie wrażliwe (zielony) lub odporny (czerwony) na glifosat i niesklasyfikowane (). Ten rysunek został zaadaptowany za zgodą Rainio et al.14. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 11: Wejścia i wyjścia serwera WWW EPSPSClass do testowania własnej sekwencji referencyjnej użytkownika. (A) Wejście 1: sekwencja zapytania. (B) Dane wejściowe 2: sekwencja odniesienia. (C) Dane wejściowe 3: markery aminokwasowe w sekwencjach referencyjnych. (D) Dane wyjściowe: tożsamość: frakcja markerów aminokwasowych w sekwencjach zapytania (klasa I-IV i własne sekwencje referencyjne użytkownika). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Tabela 1: Lista starterów do amplifikacji PCR genu 16S rRNA i regionu ITS w analizie mikrobiomu Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Kody enzymu syntazy 5-enolopiruzyszykimo-3-fosforanowej (EPSPS) w różnych bazach danych Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Średnie stężenie glifosatu Kliknij tutaj, aby pobrać tę tabelę.
Tabela 4: Tabela podsumowująca procent gatunków wrażliwych/odpornych na glifosat. Ta tabela została zaadaptowana za zgodą Rainio et al.14. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 5: Pozycje markerów aminokwasów w sekwencjach referencyjnych Kliknij tutaj, aby pobrać tę tabelę.