The microplankton organisms at the oxic-anoxic interface in the pelagial of the Black Sea

Transport-limited reactions in microbial systems

Predicting microbial metabolic rates and emergent biogeochemical fluxes remains challenging due to the many unknown population dynamical, physiological and reaction-kinetic parameters and uncertainties in species composition. Here, we show that the need for these parameters can be eliminated when population dynamics and reaction kinetics operate at much shorter time scales than physical mixing processes. Such scenarios are widespread in poorly mixed water columns and sediments. In this 'fast-reaction-transport' (FRT) limit, all that is required for predictions are chemical boundary conditions, the physical mixing processes and reaction stoichiometries, while no knowledge of species composition, physiology or population/reaction kinetic parameters is needed. Using time-series data spanning years 2001-2014 and depths 180-900 m across the permanently anoxic Cariaco Basin, we demonstrate that the FRT approach can accurately predict the dynamics of major electron donors and acceptors (Pearson r ≥ 0.9 in all cases). Hence, many microbial processes in this system are largely transport limited and thus predictable regardless of species composition, population dynamics and kinetics. Our approach enables predictions for many systems in which microbial community dynamics and kinetics are unknown. Our findings also reveal a mechanism for the frequently observed decoupling between function and taxonomy in microbial systems.

Flow cytometric sorting of loricate choanoflagellates from the oligotrophic ocean

It is challenging to study protists with extensive, loosely-associated extracellular structures because of the problems with keeping specimens intact. Here we have tested the suitability of high-speed flow cytometric sorting as a tool for studying such protists using oceanic loricate choanoflagellates as a model. We chose choanoflagellates because their lorica-to-cell volume ratio is > 10 and the voluminous loricae, i.e., the siliceous cell baskets essential for taxonomic identification, only loosely enclose the cells. Besides, owing to low concentrations, choanoflagellates are grossly under-sampled in the oligotrophic ocean. On four research cruises the small heterotrophic protists from samples collected in the photic layer of the South Atlantic and South Pacific oligotrophic (sub)tropical gyres and adjacent mesotrophic waters were flow sorted at sea for electron microscopy ashore. Among the flow-sorted protozoa we were able to select loricate choanoflagellates to assess their species diversity and concentrations. The well-preserved loricae of flow-sorted choanoflagellates made identification of 29 species from 14 genera possible. In the oligotrophic waters, we found neither endemic species nor evident morphological adaptations other than a tendency for lighter silicification of loricae. Common sightings of specimens storing extra costae in preparation for division, indicate choanoflagellates thriving in oligotrophic waters rather than enduring them. Thus, this case study demonstrates that high-speed flow sorting can assist in studying protists with extracellular structures 16-78× bigger than the enclosed cell.

Notable predominant morphology of the smallest most abundant protozoa of the open ocean revealed by electron microscopy

In the microbe-driven ecosystems of the open ocean, the small heterotrophic flagellates (sHF) are the chief microbial predators and recyclers of essential nutrients to phototrophic microbes. Even with intensive molecular phylogenetic studies of the sHF, the origins of their feeding success remain obscure because of limited understanding of their morphological adaptations to feeding. Here, we examined the sHF morphologies in the largest, most oligotrophic South Pacific and Atlantic (sub)tropical gyres and adjacent mesotrophic waters. On four research cruises, the sHF cells were flow cytometrically sorted from bacterioplankton and phytoplankton for electron microscopy. The sorted sHF comprised chiefly heterokont (HK) biflagellates and unikont choanoflagellates numerically at around 10-to-1 ratio. Of the four differentiated morphological types of HK omnipresent in the open ocean, the short-tinsel heterokont (stHK), whose tinsel flagellum is too short to propagate a complete wave, is predominant and a likely candidate to be the most abundant predator on Earth. Modeling shows that the described stHK propulsion is effective in feeding on bacterioplankton cells at low concentrations; however, owing to general prey scarcity in the oligotrophic ocean, selective feeding is unsustainable and omnivory is equally obligatory for the seven examined sHF types irrespective of their mode of propulsion.

Circumventing kinetics in biogeochemical modeling

Microbial metabolism drives biogeochemical fluxes in virtually every ecosystem. Modeling these fluxes is challenged by the incredible diversity of microorganisms, whose kinetic parameters are largely unknown. In poorly mixed systems, such as stagnant water columns or sediments, however, long-term bulk microbial metabolism may become limited by physical transport rates of substrates across space. Here we mathematically show that under these conditions, biogeochemical fluxes are largely predictable based on the system's transport properties, chemical boundary conditions, and the stoichiometry of metabolic pathways, regardless of the precise kinetics of the resident microorganisms. We formalize these considerations into a predictive modeling framework and demonstrate its use for the Cariaco Basin subeuphotic zone, one of the largest anoxic marine basins worldwide. Using chemical concentration data solely from the upper boundary (depth 180 m) and lower boundary (depth 900 m), but without a priori knowledge of metabolite fluxes, chemical depth profiles, kinetic parameters, or microbial species composition, we predict the concentrations and vertical fluxes of biologically important substances, including oxygen, nitrate, hydrogen sulfide, and ammonium, across the entire considered depth range (180-900 m). Our predictions largely agree with concentration measurements over a period of 14 years ([Formula: see text] = 0.78-0.92) and become particularly accurate during a period where the system was near biogeochemical steady state (years 2007-2009, [Formula: see text] = 0.86-0.95). Our work enables geobiological predictions for a large class of ecosystems without knowledge of kinetic parameters or geochemical depth profiles. Conceptually, our work provides a possible explanation for the decoupling between microbial species composition and bulk metabolic function, observed in various ecosystems.

Integrating biogeochemistry with multiomic sequence information in a model oxygen minimum zone

Microorganisms are the most abundant lifeform on Earth, mediating global fluxes of matter and energy. Over the past decade, high-throughput molecular techniques generating multiomic sequence information (DNA, mRNA, and protein) have transformed our perception of this microcosmos, conceptually linking microorganisms at the individual, population, and community levels to a wide range of ecosystem functions and services. Here, we develop a biogeochemical model that describes metabolic coupling along the redox gradient in Saanich Inlet-a seasonally anoxic fjord with biogeochemistry analogous to oxygen minimum zones (OMZs). The model reproduces measured biogeochemical process rates as well as DNA, mRNA, and protein concentration profiles across the redox gradient. Simulations make predictions about the role of ubiquitous OMZ microorganisms in mediating carbon, nitrogen, and sulfur cycling. For example, nitrite "leakage" during incomplete sulfide-driven denitrification by SUP05 Gammaproteobacteria is predicted to support inorganic carbon fixation and intense nitrogen loss via anaerobic ammonium oxidation. This coupling creates a metabolic niche for nitrous oxide reduction that completes denitrification by currently unidentified community members. These results quantitatively improve previous conceptual models describing microbial metabolic networks in OMZs. Beyond OMZ-specific predictions, model results indicate that geochemical fluxes are robust indicators of microbial community structure and reciprocally, that gene abundances and geochemical conditions largely determine gene expression patterns. The integration of real observational data, including geochemical profiles and process rate measurements as well as metagenomic, metatranscriptomic and metaproteomic sequence data, into a biogeochemical model, as shown here, enables holistic insight into the microbial metabolic network driving nutrient and energy flow at ecosystem scales.

Impact of protist grazing on a key bacterial group for biogeochemical cycling in Baltic Sea pelagic oxic/anoxic interfaces

Barrier zones between oxic and anoxic water masses (redoxclines) host highly active prokaryotic communities with important roles in biogeochemical cycling. In Baltic Sea pelagic redoxclines, Epsilonproteobacteria of the genus Sulfurimonas (subgroup GD17) have been shown to dominate chemoautotrophic denitrification. However, little is known on the loss processes affecting this prokaryotic group. In the present study, the protist grazing impact on the Sulfurimonas subgroup GD17 was determined for suboxic and oxygen/hydrogen sulphide interface depths of Baltic Sea redoxclines, using predator exclusion assays and bacterial amendment with the cultured representative 'Sulfurimonas gotlandica' strain GD1. Additionally, the principal bacterivores were identified by RNA-Stable Isotope Probing (RNA-SIP). The natural Sulfurimonas subgroup GD17 population grew strongly under oxygen/hydrogen sulphide interface conditions (doubling time: 1-1.5 days), but protist grazing could consume the complete new cell production per day. In suboxic samples, little or no growth of Sulfurimonas subgroup GD17 was observed. RNA-SIP identified five active grazers, belonging to typical redoxcline ciliates (Oligohymenophorea, Prostomatea) and globally widespread marine flagellate groups (MAST-4, Chrysophyta, Cercozoa). Overall, we demonstrate for the first time that protist grazing can control the growth, and potentially the vertical distribution, of a chemolithoautotrophic key-player of oxic/anoxic interfaces.

Unveiling microbial life in new deep-sea hypersaline Lake Thetis. Part I: Prokaryotes and environmental settings

In September 2008, an expedition of the RV Urania was devoted to exploration of the genomic richness of deep hypersaline anoxic lakes (DHALs) located in the Western part of the Mediterranean Ridge. Approximately 40 nautical miles SE from Urania Lake, the presence of anoxic hypersaline lake, which we named Thetis, was confirmed by swath bathymetry profiling and through immediate sampling casts. The brine surface of the Thetis Lake is located at a depth of 3258 m with a thickness of ≈ 157 m. Brine composition was found to be thalassohaline, saturated by NaCl with a total salinity of 348‰, which is one of highest value reported for DHALs. Similarly to other Mediterranean DHALs, seawater-brine interface of Thetis represents a steep pycno- and chemocline with gradients of salinity, electron donors and acceptors and posseses a remarkable stratification of prokaryotic communities, observed to be more metabolically active in the upper interface where redox gradient was sharper. [(14) C]-bicarbonate fixation analysis revealed that microbial communities are sustained by sulfur-oxidizing chemolithoautotrophic primary producers that thrive within upper interface. Besides microaerophilic autotrophy, heterotrophic sulfate reduction, methanogenesis and anaerobic methane oxidation are likely the predominant processes driving the ecosystem of Thetis Lake.

Phylogenetic diversity of bacterial and archaeal communities in the anoxic zone of the Cariaco Basin

Microbial community samples were collected from the anoxic zone of the Cariaco Basin at depths of 320, 500, and 1,310 m on a November 1996 cruise and were used to construct 16S ribosomal DNA libraries. Of 60 nonchimeric sequences in the 320-m library, 56 belonged to the epsilon subdivision of the Proteobacteria (epsilon-Proteobacteria) and 53 were closely related to ectosymbionts of Rimicaris exoculata and Alvinella pompejana, which are referred to here as epsilon symbiont relatives (ESR). The 500-m library contained sequences affiliated with the fibrobacteria, the Flexibacter-Cytophaga-Bacteroides division, the division Verrucomicrobia, the division Proteobacteria, and the OP3 candidate division. The Proteobacteria included members of the gamma, delta, epsilon and new candidate subdivisions, and gamma-proteobacterial sequences were dominant (25.6%) among the proteobacterial sequences. As in the 320-m library, the majority of the epsilon-proteobacteria belonged to the ESR group. The genus Fibrobacter and its relatives were the second largest group in the library (23.6%), followed by the delta-proteobacteria and the epsilon-proteobacteria. The 1,310-m library had the greatest diversity; 59 nonchimeric clones in the library contained 30 unique sequences belonging to the planctomycetes, the fibrobacteria, the Flexibacter-Cytophaga-Bacteroides division, the Proteobacteria, and the OP3 and OP8 candidate divisions. The proteobacteria included members of new candidate subdivisions and the beta, gamma, delta, and epsilon-subdivisions. ESR sequences were still present in the 1,310-m library but in a much lower proportion (8.5%). One archaeal sequence was present in the 500-m library (2% of all microorganisms in the library), and eight archaeal sequences were present in the 1,310-m library (13.6%). All archaeal sequences fell into two groups; two clones in the 1,310-m library belonged to the kingdom Crenarchaeota and the remaining sequences in both libraries belonged to the kingdom Euryarchaeota. The latter group appears to be related to the Eel-TA1f2 sequence, which belongs to an archaeon suggested to be able to oxidize methane anaerobically. Based on phylogenetic inferences and measurements of dark CO(2) fixation, we hypothesized that (i) the ESR are autotrophic anaerobic sulfide oxidizers, (ii) sulfate reduction and fermentative metabolism may be carried out by a large number of bacteria in the 500- and 1,310-m libraries, and (iii) members of the Euryarchaeota found in relatively large numbers in the 1,310-m library may be involved in anaerobic methane oxidation. Overall, the composition of microbial communities from the Cariaco Basin resembles the compositions of communities from several anaerobic sediments, supporting the hypothesis that the Cariaco Basin water column is similar to anaerobic sediments.

Molecular indicators for palaeoenvironmental change in a Messinian evaporitic sequence (Vena del Gesso, Italy). II: High-resolution variations in abundances and 13C contents of free and sulphur-bound carbon skeletons in a single marl bed

The extractable organic matter of 10 immature samples from a marl bed of one evaporitic cycle of the Vena del Gesso sediments (Gessoso-solfifera Fm., Messinian, Italy) was analyzed quantitatively for free hydrocarbons and organic sulphur compounds. Nickel boride was used as a desulphurizing agent to recover sulphur-bound lipids from the polar and asphaltene fractions. Carbon isotopic compositions (delta vs PDB) of free hydrocarbons and of S-bound hydrocarbons were also measured. Relationships between these carbon skeletons, precursor biolipids, and the organisms producing them could then be examined. Concentrations of S-bound lipids and free hydrocarbons and their delta values were plotted vs depth in the marl bed and the profiles were interpreted in terms of variations in source organisms, 13 C contents of the carbon source, and environmentally induced changes in isotopic fractionation. The overall range of delta values measured was 24.7%, from -11.6% for a component derived from green sulphur bacteria (Chlorobiaceae) to -36.3% for a lipid derived from purple sulphur bacteria (Chromatiaceae). Deconvolution of mixtures of components deriving from multiple sources (green and purple sulphur bacteria, coccolithophorids, microalgae and higher plants) was sometimes possible because both quantitative and isotopic data were available and because either the free or S-bound pool sometimes appeared to contain material from a single source. Several free n-alkanes and S-bound lipids appeared to be specific products of upper-water-column primary producers (i.e. algae and cyanobacteria). Others derived from anaerobic photoautotrophs and from heterotrophic protozoa (ciliates), which apparently fed partly on Chlorobiaceae. Four groups of n-alkanes produced by algae or cyanobacteria were also recognized based on systematic variations of abundance and isotopic composition with depth. For hydrocarbons probably derived from microalgae, isotopic variations are well correlated with those of total organic carbon. A resistant aliphatic biomacromolecule produced by microalgae is, therefore, probably an important component of the kerogen. These variations reflect changes in the depositional environment and early diagenetic transformations. Changes in the concentrations of S-bound lipids induced by variations in conditions favourable for sulphurization were discriminated from those related to variations in primary producer assemblages. The water column of the lagoonal basin was stratified and photic zone anoxia occurred during the early and middle stages of marl deposition. During the last stage of the marl deposition the stratification collapsed due to a significant shallowing of the water column. Contributions from anaerobic photoautotrophs were apparently associated with variations in depth of the chemocline.

Evidence for gammacerane as an indicator of water column stratification

A new route for the formation of gammacerane from tetrahymanol is proposed; in addition to dehydration and hydrogenation, sulphurisation and early C-S cleavage are shown to be important in the pathway of formation, especially in marine sediments. Evidence is twofold. First, relatively large amounts of the gammacerane skeleton are sequestered in S-rich macromolecular aggregates formed by natural sulphurisation of functionalised lipids. Selective cleavage of polysulphide linkages with MeLi/MeI led to formation of 3-methylthiogammacerane, indicating that the gammacerane skeleton is primarily bound via sulphur at position 3, consistent with the idea that tetrahymanol (or the corresponding ketone) is the precursor for gammacerane. Second, upon mild artificial maturation of two sediments using hydrous pyrolysis, gammacerane is released from S-rich macromolecular aggregates by cleavage of the relatively weak C-S bonds. The stable carbon isotopic compositions of gammacerane and lipids derived from primary producers and green sulphur bacteria in both the Miocene Gessoso-solfifera and Upper Jurassic Allgau Formations indicate that gammacerane is derived from bacterivorous ciliates which were partially feeding on green sulphur bacteria. This demonstrates that anaerobic ciliates living at or below the chemocline are important sources for gammacerane, consistent with the fact that ciliates only biosynthesize tetrahymanol if their diet is deprived of sterols. This leads to the conclusion that gammacerane is an indicator for water column stratification, which solves two current enigmas in gammacerane geochemistry. Firstly, it explains why gammacerane is often found in sediments deposited under hypersaline conditions but is not necessarily restricted to this type of deposits. Secondly, it explains why lacustrine deposits may contain abundant gammacerane since most lakes in the temperate climatic zones are stratified during summer.