$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Szacuje się, że tereny podmokłe zawierają 529 Pg węgla (C), głównie w postaci organicznego węgla zakopanego w złożach torfu1. Obecnie takie torfowiska działają jak pochłaniacz C netto, sekwestrując 29 Tg C y-1 w samej Ameryce Północnej1. Jednak zakłócenia środowiskowe, takie jak drenaż, pożary, susze i wyższe temperatury, mogą zrównoważyć ten pochłaniacz C poprzez zwiększenie rozkładu materii organicznej, co skutkuje zwiększonymi stratami węgla w wyniku produkcji gazów cieplarnianych (dwutlenku węgla [CO2] i metanu [CH4]) 1,2. Zmiana klimatu może przyczyniać się do utraty węgla, jeśli wyższe temperatury lub suche warunki stymulują szybszy rozkład węgla przez mikroorganizmy. Alternatywnie, wyższe temperatury i stężenia CO2 w powietrzu mogą stymulować produkcję pierwotną do sekwestracji większej ilości CO2 w postaci węgla organicznego (OC). To, w jakim stopniu i jak szybko to OC rozkłada się na CO2 i CH4, zależy od złożonych interakcji między substratami donorów elektronów, dostępnością akceptorów elektronów i mikroorganizmami pośredniczącymi w transformacji. W wielu przypadkach mechanizmy te nie są dobrze scharakteryzowane, a zatem ich reakcja na perturbacje środowiskowe nie jest dobrze ograniczona i nie jest jasne, jaki będzie wynik netto zmian klimatycznych dla bilansu węglowego w ekosystemach torfowiskowych.
Złożona natura naturalnej materii organicznej (NOM) sprawiła, że nawet identyfikacja związków organicznych obecnych w mieszaninach NOM była historycznie trudna. Ostatnie postępy znacznie poprawiły naszą zdolność do charakteryzowania związków, które tradycyjnie i do pewnego stopnia nadal są uważane za oporne związki humusowe lub fulwowe3,4,5. Teraz rozumiemy, że wiele z tych związków jest w rzeczywistości dostępnych mikrobiologicznie i może ulec rozkładowi, jeśli dostępny jest odpowiedni końcowy akceptor elektronów (TEA)6,7. Obliczanie nominalnego stopnia utlenienia węgla (NOSC) dla związku dostarcza metryki do przewidywania potencjału rozkładu i wymaganej wydajności energetycznej TEA. Wymaga to jednak charakterystyki materii organicznej na poziomie molekularnym7. NOSC oblicza się ze wzoru cząsteczkowego za pomocą następującego równania7: NOSC = − ((−z + 4(#C) + (#H) − 3(#N) − 2(#O) + 5(#P) − 2(#S)) / (#C)) + 4, gdzie z jest ładunkiem netto. NOSC jest skorelowany z termodynamiczną siłą napędową8, przy czym związki o wyższym NOSC są łatwiejsze do degradacji, podczas gdy związki o niższym NOSC wymagają coraz bardziej energetycznych TEA, aby mogły zostać zredukowane. Związki o NOSC mniejszym niż −2 wymagają wysokiej energii TEA, takich jakO2, azotan lub MnIV, i nie mogą być degradowane przez powszechnie występujące herbaty o niższej wartości energetycznej, takie jak FeIII lub siarczan7. Jest to ważna kwestia w podmokłych warunkach beztlenowych występujących na terenach podmokłych, gdzie O2 i inne herbaty o wysokiej wydajności energetycznej są rzadkie9, a zatem degradacja niższych związków NOSC w tych warunkach jest termodynamicznie ograniczona. Perturbacje środowiskowe mogą wpływać na stan termodynamiczny ekosystemu poprzez zmiany hydrologiczne, które wpływają naO2 (najbardziej energetyczny akceptor elektronów), zmiany w substratach organicznych i akceptorach elektronów udostępnianych przez produkcję pierwotną oraz w mniejszym stopniu przez temperaturę. Ważnym przykładem wpływu temperatury w systemach mokradeł jest kompromis zachodzący między homoacetogenezą (tj. produkcją octanu z CO2 i H2 ) a metanogenezą hydrogenotroficzną (tj. produkcjąCH 4 z CO2 i H2). Wydaje się, że w niskich temperaturach homoacetogeneza jest nieco faworyzowana, podczas gdy cieplejsze temperatury faworyzują produkcję CH410. Ten wpływ temperatury może mieć ważne implikacje dla reakcji ekosystemów na zmieniający się klimat, ponieważ CH4 jest znacznie silniejszym gazem cieplarnianym niż CO211, a zatem zwiększenie produkcji CH4 kosztem CO2 w cieplejszych temperaturach może przyczynić się do dodatniego sprzężenia zwrotnego z ociepleniem klimatu.
Torfowiska produkują globalnie znaczące ilości CO2 i CH46 poprzez mikrobiologiczne oddychanie naturalnie występującej materii organicznej. NOSC substratów węgla organicznego określa względną proporcję wytwarzanego CO2:CH4, co jest parametrem krytycznym ze względu na wyższe wymuszanie radiacyjne CH4 w porównaniu do CO211, ale także dlatego, że wysiłki modelowania zidentyfikowały ten stosunek jako krytyczny parametr do szacowania strumienia C na torfowiskach12. W przypadku braku końcowych akceptorów elektronów innych niż CO2 można wykazać za pomocą bilansu elektronowego, że organiczne substraty C o NOSC > 0 będą wytwarzać CO2:CH4 > 1, organiczne C o NOSC = 0 wytwarza CO2 i CH4 w stosunku równomolowym, a organiczne C o NOSC < 1 będą wytwarzać CO2: CH4 < 113. Rozkład OC w naturalnych ekosystemach odbywa się za pośrednictwem mikroorganizmów, tak że nawet jeśli degradacja określonego związku jest termodynamicznie wykonalna, jest ona kinetycznie ograniczona przez aktywność enzymów mikrobiologicznych oraz, w warunkach beztlenowych, przez termodynamiczną siłę napędową (tj. NOSC)7. Do tej pory pełne scharakteryzowanie materii organicznej było wyzwaniem, ponieważ różnorodność obecnych związków wymaga różnych uzupełniających się technik ich charakteryzacji. Ostatnie postępy wypełniły tę lukę; Korzystając z zestawu technik analitycznych, możemy analizować szeroki zakres związków organicznych, zapewniając charakterystykę na poziomie molekularnym, a w niektórych przypadkach kwantyfikację, od małych metabolitów pierwotnych, takich jak glukoza, do 800 Da poliheterocykli. Wcześniej takie duże złożone cząsteczki byłyby charakteryzowane po prostu jako podobne do ligniny lub podobne do taniny i zakładano, że są krnąbrne. Charakterystyka na poziomie molekularnym pozwala jednak na obliczenie NOSC nawet dla tych dużych złożonych cząsteczek. Te wartości NOSC są liniowo skorelowane z termodynamiczną siłą napędową, co pozwala na ocenę jakości materii organicznej dostępnej do rozkładu, co w wielu przypadkach ujawnia, że te złożone cząsteczki mogą w rzeczywistości ulegać degradacji mikrobiologicznej nawet w warunkach beztlenowych panujących na terenach podmokłych.
Ponieważ wprowadzenie O2 pozwala na rozkład materii organicznej o prawie wszystkich naturalnie obserwowanych wartościach NOSC, skupiamy się tutaj na zmianach w materii organicznej i proteomice mikrobiologicznej, które prawdopodobnie będą głównymi czynnikami napędzającymi systemy wodno-błotne (tj. ograniczone O2). Jednak wszystkie techniki, które omówimy, można zastosować do materii organicznej z dowolnego ekosystemu. Zazwyczaj do oceny jakości materii organicznej stosuje się pomiary masowe oparte na analizach optycznych i fluorescencyjnych3,14. Jednak podczas korzystania z pomiarów masowych, takich jak te, drobne szczegóły są tracone, ponieważ duża liczba cząsteczek jest kategoryzowana razem pod ogólnymi terminami, takimi jak humiki lub fulwy. Definicje tych kategorii nie są dobrze ograniczone i w rzeczywistości mogą się różnić w zależności od badania, co uniemożliwia porównania. Co więcej, pomiary masowe nie dostarczają szczegółów molekularnych niezbędnych do obliczenia termodynamiki rządzącej systemem, a zatem nie są w stanie naprawdę ocenić jakości materii organicznej15.
Indywidualne techniki, takie jak spektrometria mas rezonansu cyklotronowego z transformacją jonową Fouriera (FTICR-MS), spektroskopia magnetycznego rezonansu jądrowego (NMR), spektrometria mas chromatografii gazowej (GC-MS) i spektrometria mas chromatografii cieczowej (LC-MS) dostarczają takich szczegółów na poziomie molekularnym. Chociaż każda z tych technik ma swoje ograniczenia, mają one również swoje mocne strony, które można wykorzystać w zintegrowanym podejściu w celu uzyskania drobnych szczegółów molekularnych niezbędnych do ilościowego określenia jakości materii organicznej w rygorystycznym sensie termodynamicznym. GC-MS jest przydatny do identyfikacji krytycznych małych metabolitów, które mogą mieć proksymalny wpływ na produkcję CO2 i CH4 (np. glukozy, octanu itp.); jednak GC-MS wymaga weryfikacji w odniesieniu do wzorca i dlatego jest ograniczona do już znanych związków obecnych w bazie danych, co uniemożliwia identyfikację nowych związków. Co więcej, GC-MS jest techniką półjakościową pozwalającą na wnioskowanie o zmianach stężeń względnych, ale nie dostarczającą rzeczywistych informacji o stężeniach niezbędnych na przykład do obliczenia energii swobodnych Gibba. Wreszcie, GC-MS wymaga derywatyzacji cząsteczek przed analizą, co ogranicza rozdzielczość do związków mniejszych niż ~400 Da, a lotne alkohole są tracone na etapie suszenia.
Jednowymiarowy (1D) 1H płynny NMR pozwala na wysoce ilościową charakterystykę małych metabolitów (w tym pierwotnych metabolitów o małej masie cząsteczkowej i substancji lotnych, takich jak alkohole, octan, aceton, mrówczan, pirogronian, bursztynian, krótkołańcuchowe kwasy tłuszczowe, a także szereg węglowodanów notorycznie nieobecnych lub zagrożonych w metodach opartych na MS), a ich stężenia są szczególnie przydatne do obliczania parametrów termodynamicznych. Jednak, podobnie jak GC-MS, 1D NMR złożonych mieszanin wymaga standaryzacji w stosunku do bazy danych, a zatem sam w sobie nie pozwala na łatwą identyfikację nowych związków, które prawdopodobnie będą obfitować w złożonych naturalnych i zmieniających się ekosystemach. Ponadto NMR jest mniej czuły niż techniki oparte na MS, a zatem ilościowe profilowanie metabolitów uzyskuje się tylko powyżej 1 μM przy użyciu systemów NMR wyposażonych w zimne sondy chłodzone helem. Nie jest to powszechnie doceniane, niektóre zimne sondy NMR są odporne na sól i umożliwiają analizę mieszaniny środowiskowej w obecności milimolowych stężeń soli, gdy są używane w probówkach na próbki o mniejszej średnicy (< 3 mm średnicy zewnętrznej)16. Jednak kolejną komplikacją NMR jest to, że duże ilości metali i minerałów paramagnetycznych (np. Fe i Mn powyżej 1 - 3% wag.), które mogą występować w dużych ilościach w glebach wyżynnych, mogą poszerzać cechy spektralne i komplikować interpretację widm NMR. Zastosowanie ekstrakcji do fazy stałej (SPE) może pomóc w interpretacji metod metabolomicznych opartych zarówno na NMR, jak i MS poprzez redukcję soli mineralnych i poprawę jakości widma.
FTICR-MS poprzez bezpośrednie wstrzykiwanie jest bardzo czułą techniką, zdolną do wykrywania dziesiątek tysięcy metabolitów z pojedynczej próbki, ale nie wychwytuje krytycznych małych metabolitów, takich jak octan, pirogronian i bursztynian, i jest niezwykle trudna do użycia dla cukrów i innych węglowodanów17, ani nie dostarcza informacji ilościowych. Jednak w przeciwieństwie do innych technik, FTICR-MS doskonale radzi sobie z identyfikacją i przypisywaniem wzoru cząsteczkowego nowym związkom, a zatem identyfikuje największą liczbę związków dostarczających więcej informacji molekularnych niż którakolwiek z innych opisanych technik. Jest to przydatne, ponieważ informacje molekularne dostarczone przez FTICR-MS (i inne techniki) mogą być wykorzystane do obliczenia NOSC, która jest związana z termodynamiczną siłą napędową rządzącą prawdopodobieństwem określonych reakcji8 i ich szybkości w określonych warunkach7. Ponadto, poprzez sprzężenie FTICR-MS z technikami separacji, takimi jak LC wraz z tandemowym MS, można uzyskać ilościowe informacje strukturalne, kompensując niektóre wady tej techniki. LC-MS jest przydatny do identyfikacji związków lipidopodobnych i innych metabolitów, które nie są dobrze scharakteryzowane przez żadną z innych metod. Sprzężenie LC FTICR-MS lub LC-MS z kolektorem frakcji i zebranie frakcji specyficznych niewiadomych będących przedmiotem zainteresowania w celu wyjaśnienia strukturalnego za pomocą dwuwymiarowego (2D) NMR w stanie ciekłym to idealna sytuacja do identyfikacji i ilościowego określania nieznanych związków18,19. Jest to jednak bardzo czasochłonny krok, który można wykorzystać w razie potrzeby. Każda z tych technik, rozpatrywana indywidualnie, zapewnia inny obraz materii organicznej, a integrując je, możemy osiągnąć pełniejsze zrozumienie niż przy użyciu którejkolwiek z technik w izolacji.
Podczas gdy rozważania termodynamiczne ustalają ostateczne ograniczenia dotyczące możliwych przemian w systemie, rozkład materii organicznej jest spowodowany przez mikroorganizmy, których aktywność enzymatyczna kontroluje szybkość reakcji. W związku z tym pełne zrozumienie kontroli rozkładu materii organicznej i ostatecznie produkcji gazów cieplarnianych (CO2 i CH4) z terenów podmokłych wymaga zintegrowanego podejścia omicznego do scharakteryzowania aktywności enzymów mikrobiologicznych, a także metabolitów. W tym artykule opisujemy metodę uzyskiwania takiej kompleksowej analizy z pojedynczej próbki przy użyciu podejścia sekwencyjnego, które skutkuje w pełni sparowaną analizą. Podejście to rozszerza protokół ekstrakcji metabolitów, białek i lipidów (MPLex), w którym proteomika została sprzężona z GC-MS i LC-MS20 w celu identyfikacji małych metabolitów, białek i lipidów poprzez włączenie ilościowych informacji o metabolitach za pomocą NMR i identyfikację większych metabolitów wtórnych za pomocą FTICR-MS. Nieco inaczej niż w przypadku MPLex, protokół rozpoczynamy od ekstrakcji wodą, a następnie stosujemy ekstrakcję sekwencyjną z coraz bardziej niepolarnymi rozpuszczalnikami. Wszystkie ekstrakcje przeprowadza się na pojedynczej próbce, co pozwala zachować próbkę, gdy objętości są ograniczone lub trudne do uzyskania, i zmniejsza błąd doświadczalny wynikający z różnic między podwielokrotnościami z niejednorodnych matryc próbek (np. gleby i torfu) lub różnic w warunkach i czasie przechowywania.
Na koniec, łącząc analizy OM z analizami proteomicznymi społeczności mikrobiologicznej, możemy zbudować sieci metaboliczne, które opisują szlaki i przemiany rozkładu materii organicznej. Pozwala nam to przetestować konkretne hipotezy dotyczące tego, w jaki sposób zakłócenia systemu wpłyną na ostateczną produkcję CO2 i CH4 poprzez zmianę dostępnych substratów organicznych, akceptorów elektronów i społeczności mikroorganizmów pośredniczących w reakcjach poprzez aktywność katalizatorów enzymatycznych.
Ogólnym celem tej metody jest dostarczenie pojedynczego protokołu przepustowości do analizy metabolitów, lipidów i białek mikrobiologicznych z pojedynczej próbki, tworząc w ten sposób w pełni sparowany zestaw danych do budowy sieci metabolicznych przy jednoczesnym ograniczeniu błędów analitycznych.