Większość genów peroksydaz w borealnej próchnicy włączają zasłonakowate

Metatranskryptom zamiast pomiaru aktywności

13 sierpnia 2026 roku ukazała się w wersji online praca, którą podpisali Florian Barbi, Uwe Menzel, Domenico Simone, Tuula Niskanen i Björn D. Lindahl; numer drukowany „New Phytologist” datowany jest na październik 2026 roku — tom 252, zeszyt 1, strony 400–413 (rekord Crossref). Autorzy pracują w Instytucie Mikrobiologii Czeskiej Akademii Nauk, na Szwedzkim Uniwersytecie Nauk Rolniczych w Uppsali i w Fińskim Muzeum Historii Naturalnej przy Uniwersytecie Helsińskim. Pytanie postawili proste w brzmieniu, a trudne do rozstrzygnięcia: które grzyby w glebie borealnego lasu naprawdę włączają geny peroksydaz klasy II. To rodzina potężnych enzymów utleniających, które wyewoluowały wśród pieczarniaków (Agaricomycetes) i rozbijają ligninę za pośrednictwem silnie utleniających jonów Mn3+, odsłaniając ukryte za nią celulozę i hemicelulozy. Przypisywano je dotąd przede wszystkim saprotroficznym podstawczakom rozkładającym drewno i ściółkę. Pełny tekst jest otwarty na licencji CC BY 4.0.

Materiał pochodzi z rezerwatu Fiby Urskog: 64 hektary starodrzewu 16 km na zachód od Uppsali, od kilku stuleci bez większej ingerencji człowieka. Wyznaczono szesnaście poletek 8 × 8 m — osiem na żyznych, bogatych w azot fragmentach zdominowanych przez świerk pospolity (Picea abies) i osiem na ubogich, zdominowanych przez sosnę zwyczajną (Pinus sylvestris). Rdzenie o średnicy 3 cm pobierano przez całą miąższość poziomu organicznego w cztery dni: 4, 7, 17 i 27 października 2016 roku, po cztery rdzenie z poletka za każdym razem. Szesnaście rdzeni z jednego poletka łączono w jedną próbkę zbiorczą. Żywe mchy, świeżo opadłą ściółkę i glebę mineralną usuwano, a rdzeń zamrażano na suchym lodzie w ciągu sekund. Prace terenowe są zatem o dekadę starsze niż publikacja.

Z każdej z szesnastu próbek odczytano średnio około 125 milionów sekwencji, w sumie blisko 2 miliardy odczytów. Zespół nie składał całego metatranskryptomu. Własnym narzędziem RING (Reads Isolator from Nucleotidic Giga-dataset) wyławiał odczyty trzech rodzin genów — AA2 (peroksydazy klasy II), GH6 (celobiohydrolazy, czyli enzymy rozkładające celulozę) i GT48 (syntazy β-1,3-glukanu, wskaźnik żywej grzybni) — oraz β-tubuliny, która posłużyła za mianownik normalizujący. Do identyfikacji wzięto najsilniej reprezentowane kontigi: 146 dla AA2, 63 dla GH6 i 267 dla GT48, obejmujące po 80% odczytów zmapowanych w danej rodzinie. Zespół szedł przy tym za jasno postawioną hipotezą: że w glebie ubogiej ekspresja peroksydaz okaże się niższa, bo saprotroficznych pieczarniaków jest tam mniej.

Ile naprawdę znaczy 78 procent

Zdanie, z którego wyrasta tytuł pracy, brzmi: średnio 78% transkryptów AA2 pochodziło z rodzajów ektomikoryzowych, w większości z rodziny zasłonakowatych (Cortinariaceae), po czym w nawiasie podano 70% dla poletek żyznych i 73% dla ubogich. Nawias należy do zasłonakowatych, a nie do sumy grzybów ektomikoryzowych — tak samo zbudowane jest bliźniacze zdanie o GT48, w którym udział zasłonakowatych wynosi 28% i 25%. Saprotrofy dały średnio 25% zidentyfikowanych transkryptów AA2, na czele z grzybówkami (Mycena): 10% w drzewostanie żyznym i 9,5% w ubogim. Poza zasłonakowatymi transkrypty AA2 wykryto jeszcze — wszędzie poniżej 5% udziału — u gołąbków (Russula), włośnianek (Hebeloma) i u podziemnie owocnikującego rodzaju Hysterangium.

Te udziały przytaczamy dokładnie tak, jak podaje je tekst, i celowo nie wykonujemy na nich żadnych działań. Nie zamykają się w stu procentach: 78% plus 25% daje 103%, a 70 i 73 uśredniają się do 71,5, nie do 78. Dotyczą przy tym wyłącznie odczytów mapujących się na zidentyfikowane, najsilniej reprezentowane kontigi — nie całej puli peroksydaz obecnej w glebie.

Druga liczba mówi o czym innym i jest równie ważna. Ogólna ekspresja AA2, znormalizowana względem β-tubuliny, była ponad dwukrotnie wyższa w glebie ubogiej (0,67 ± 0,12) niż w żyznej (0,30 ± 0,06), przy P = 0,007 — czyli odwrotnie, niż zakładała pierwsza hipoteza autorów. Rozrzut podano w tekście bez etykiety; podpis ryciny 2 mówi o odchyleniu standardowym. Ogólna ekspresja GH6 i GT48 nie różniła się między typami lasu. Zmieniła się natomiast przynależność taksonomiczna transkryptów GH6: na poletkach żyznych 63% sekwencji mapujących się na zidentyfikowane kontigi pochodziło od saprotroficznych pieczarniaków (grzybówki 29%, monetnice z rodzaju Rhodocollybia 15%), a na ubogich podstawczaki dawały już tylko 25%, wobec 73% workowców, z samą klasą Leotiomycetes na poziomie 45%. To zarazem najsłabiej podparty wynik w pracy: GH6 opiera się na 50 150 zmapowanych odczytach, zachowanych przy najniższym z czterech użytych progów jakości mapowania (samtools Q 1), podczas gdy AA2 miało ich 237 605 przy Q 19, a GT48 430 598 przy Q 22.

Osobna niespodzianka dotyczy rodzaju Piloderma, zaliczanego dotąd w literaturze do grzybów „wydobywających azot”. Wśród najsilniej reprezentowanych kontigów AA2 nie pojawił się w ogóle. Autorzy sprawdzili więc cały wynik przeszukania kontigów AA2 w białkowych bazach NCBI i JGI i znaleźli kontigi pasujące do referencji Piloderma z podobieństwem powyżej 93%, ale stanowiły one 0,084% transkryptów AA2 — przy dużym udziale tego rodzaju w ekspresji GT48, którą traktują jako wskaźnik żywej grzybni. Piszą wprost, że ich wyniki podważają dotychczasową klasyfikację, „przynajmniej w zakresie naszego badania”. To nie znaczy, że Piloderma nie ma genów peroksydaz. Znaczy, że w tej próchnicy i w tym miesiącu prawie ich nie włączał.

Czego ta praca nie pokazała

Zmierzono transkrypty, nie rozkład. Nie badano aktywności enzymów, ubytku masy organicznej, uwalniania azotu ani strat węgla. Autorzy piszą ostrożnie: symbioza „może umożliwiać” grzybom ektomikoryzowym wykorzystanie fotoasymilatów drzewa do napędzania kosztownego energetycznie utleniania pod ziemią, a utlenianie napędzane przez mikoryzę „może dzięki temu pozwalać” drzewom pośrednio regulować rozkład i obieg składników. Cała niepewność tej pracy mieści się w odległości między „włącza gen” a „rozkłada próchnicę”.

Nie jest też prawdą, że dominacja grzybów ektomikoryzowych w ekspresji peroksydaz to zjawisko gleb ubogich. Sekwencje mapujące się na zidentyfikowane kontigi AA2 rozkładały się między rodzajami grzybów podobnie w obu warunkach, a w dyskusji autorzy piszą wprost, że grzyby ektomikoryzowe dominowały w ekspresji peroksydaz klasy II również na poletkach żyznych, tyle że przy niższym poziomie ekspresji. Różnił się poziom, nie obsada.

Nie wolno wreszcie przypisać całego kontrastu samemu azotowi. Wszystkie osiem poletek żyznych to świerczyny, wszystkie osiem ubogich to bory sosnowe — żyzność i gatunek drzewa są tu nierozdzielne. Autorzy poświęcają temu osobny akapit: interpretacji tych wzorców nie da się oprzeć wyłącznie na dostępności składników pokarmowych, bo gatunki drzew różnią się jakością ściółki, ilością węgla dostarczanego przez korzenie, preferencjami wobec symbiontów i dostępnością manganu. Manganu ekstrahowalnego było w glebie żyznej sześć razy więcej, pH wynosiło tam 4,29 ± 0,06 wobec 3,71 ± 0,03 (te wartości opisano jako błąd standardowy), a stosunki węgla do azotu ogółem oraz węgla do azotu amonowego były wyższe w glebie ubogiej — odpowiednio 1,2 i 1,8 raza. Sami autorzy sugerują, że przy sześciokrotnie mniejszej dostępności manganu więcej peroksydaz może być potrzebne choćby po to, by utrzymać roboczy zapas utlenionego Mn3+; to ostrożna propozycja z dyskusji, nie wynik pomiaru.

Jedno stanowisko, szesnaście próbek zbiorczych, jedna jesień i wyłącznie poziom organiczny, bez gleby mineralnej i bez świeżej ściółki: to nie jest portret borealnych lasów w ogóle. Nie jest to również historia o zasłonakowatych rozkładających drewno. Ekspresja genów celulolitycznych w glebie ubogiej wcale nie spadła — zmieniła właściciela.

Co z tego zostaje pod nogami

Warto pilnować rangi taksonu, bo łatwo ją przekręcić. Praca operuje rodziną zasłonakowatych, a nie rodzajem zasłonak: Cortinariaceae obejmowały dawniej jeden rodzaj, od 2022 roku podzielony na dziesięć. W obrębie tej rodziny do poziomu gatunku sięga jedno dopasowanie: kontig AA2_27, niemal idealnie zgodny z sekwencją referencyjną Cortinarius ominosus. Dla porządku — poza zasłonakowatymi praca wskazuje jeszcze idealne dopasowanie kontiga AA2_54 do Mycena sanguinolenta i niemal idealne AA2_104 do monetnicy maślanej (Rhodocollybia butyracea). Wszystkie trzy to identyfikacje sekwencji, a nie owocniki zebrane i oznaczone w rezerwacie.

Dla kogoś, kto chodzi po polskim borze, najciekawsza jest warstwa, o której ta praca mówi: ciemna próchnica między luźną ściółką a piaskiem, kilka albo kilkanaście centymetrów miąższości. Stamtąd pochodzi cały materiał. Różnicę między żyzną świerczyną a ubogim borem sosnowym widać w grubości i barwie tego poziomu, zanim jeszcze ktokolwiek spisze listę gatunków. Owocniki zasłonakowatych są sezonowym, drobnym fragmentem grzybni, która — jeśli wynik z Fiby Urskog przenosi się dalej — przez resztę roku odczytuje w tej samej próchnicy geny kojarzone dotąd z rozkładem drewna. Czy tak samo dzieje się w polskich borach, tego ta praca nie rozstrzyga; potrzebny byłby ten sam pomiar, w tej samej warstwie, wykonany jeszcze raz.