Genomic signatures of fifth autotrophic carbon assimilation pathway in bathypelagic Crenarchaeota

Metagenome-assembled genomes reveal unique metabolic adaptations of a basal marine Thaumarchaeota lineage

Thaumarchaeota constitute an abundant and ubiquitous phylum of Archaea that play critical roles in the global nitrogen and carbon cycles. Most well-characterized members of the phylum are chemolithoautotrophic ammonia-oxidizing archaea (AOA), which comprise up to 5 and 20% of the total single-celled life in soil and marine systems, respectively. Using two high-quality metagenome-assembled genomes (MAGs), here we describe a divergent marine thaumarchaeal clade that is devoid of the ammonia-oxidation machinery and the AOA-specific carbon-fixation pathway. Phylogenomic analyses placed these genomes within the uncultivated and largely understudied marine pSL12-like thaumarchaeal clade. The predominant mode of nutrient acquisition appears to be aerobic heterotrophy, evidenced by the presence of respiratory complexes and various organic carbon degradation pathways. Both genomes encoded several pyrroloquinoline quinone (PQQ)-dependent alcohol dehydrogenases, as well as a form III RuBisCO. Metabolic reconstructions suggest anaplerotic CO₂ assimilation mediated by RuBisCO, which may be linked to the central carbon metabolism. We conclude that these genomes represent a hitherto unrecognized evolutionary link between predominantly anaerobic basal thaumarchaeal lineages and mesophilic marine AOA, with important implications for diversification within the phylum Thaumarchaeota.

Contribution of Bicarbonate Assimilation to Carbon Pool Dynamics in the Deep Mediterranean Sea and Cultivation of Actively Nitrifying and CO2-Fixing Bathypelagic Prokaryotic Consortia

Covering two-thirds of our planet, the global deep ocean plays a central role in supporting life on Earth. Among other processes, this biggest ecosystem buffers the rise of atmospheric CO2. Despite carbon sequestration in the deep ocean has been known for a long time, microbial activity in the meso- and bathypelagic realm via the "assimilation of bicarbonate in the dark" (ABD) has only recently been described in more details. Based on recent findings, this process seems primarily the result of chemosynthetic and anaplerotic reactions driven by different groups of deep-sea prokaryoplankton. We quantified bicarbonate assimilation in relation to total prokaryotic abundance, prokaryotic heterotrophic production and respiration in the meso- and bathypelagic Mediterranean Sea. The measured ABD values, ranging from 133 to 370 μg C m-3 d-1, were among the highest ones reported worldwide for similar depths, likely due to the elevated temperature of the deep Mediterranean Sea (13-14°C also at abyssal depths). Integrated over the dark water column (≥200 m depth), bicarbonate assimilation in the deep-sea ranged from 396 to 873 mg C m-2 d-1. This quantity of produced de novo organic carbon amounts to about 85-424% of the phytoplankton primary production and covers up to 62% of deep-sea prokaryotic total carbon demand. Hence, the ABD process in the meso- and bathypelagic Mediterranean Sea might substantially contribute to the inorganic and organic pool and significantly sustain the deep-sea microbial food web. To elucidate the ABD key-players, we established three actively nitrifying and CO2-fixing prokaryotic enrichments. Consortia were characterized by the co-occurrence of chemolithoautotrophic Thaumarchaeota and chemoheterotrophic proteobacteria. One of the enrichments, originated from Ionian bathypelagic waters (3,000 m depth) and supplemented with low concentrations of ammonia, was dominated by the Thaumarchaeota "low-ammonia-concentration" deep-sea ecotype, an enigmatic and ecologically important group of organisms, uncultured until this study.

Shifts in the meso- and bathypelagic archaea communities composition during recovery and short-term handling of decompressed deep-sea samples

Dark ocean microbial communities are actively involved in chemoautotrophic and anaplerotic fixation of bicarbonate. Thus, aphotic pelagic realm of the ocean might represent a significant sink of CO2 and source of primary production. However, the estimated metabolic activities in the dark ocean are fraught with uncertainties. Typically, deep-sea samples are recovered to the sea surface for downstream processing on deck. Shifts in ambient settings, associated with such treatments, can likely change the metabolic activity and community structure of deep-sea adapted autochthonous microbial populations. To estimate influence of recovery and short-term handling of deep-sea samples, we monitored the succession of bathypelagic microbial community during its 3 days long on deck incubation. We demonstrated that at the end of exposition, the deep-sea archaeal population decreased threefold, whereas the bacterial fraction doubled in size. As revealed by phylogenetic analyses of amoA gene transcripts, dominance of the active ammonium-oxidizing bathypelagic Thaumarchaeota groups shifted over time very fast. These findings demonstrated the simultaneous existence of various 'deep-sea ecotypes', differentially reacting to the sampling and downstream handling. Our study supports the hypothesis that metabolically active members of meso- and bathypelagic Thaumarchaeota possess the habitat-specific distribution, metabolic complexity and genetic divergence at subpopulation level.

Heterotrophic bicarbonate assimilation is the main process of de novo organic carbon synthesis in hadal zone of the Hellenic Trench, the deepest part of Mediterranean Sea

Ammonium-oxidizing chemoautotrophic members of Thaumarchaea are proposed to be the key players in the assimilation of bicarbonate in the dark (ABD). However, this process may also involve heterotrophic metabolic pathways, such as fixation of carbon dioxide (CO2) via various anaplerotic reactions. We collected samples from the depth of 4900 m at the Matapan-Vavilov Deep (MVD) station (Hellenic Trench, Eastern Mediterranean) and used the multiphasic approach to study the ABD mediators in this deep-sea ecosystem. At this depth, our analysis indicated the occurrence of actively CO2-fixing heterotrophic microbial assemblages dominated by Gammaproteobacteria with virtually no Thaumarchaea present. [14C]-bicarbonate incorporation experiments combined with shotgun [14C]-proteomic analysis identified a series of proteins of gammaproteobacterial origin. More than quarter of them were closely related with Alteromonas macleodii ‘deep ecotype’ AltDE, the predominant organism in the microbial community of MVD. The present study demonstrated that in the aphotic/hadal zone of the Mediterranean Sea, the assimilation of bicarbonate is associated with both chemolithoauto- and heterotrophic ABD. In some deep-sea areas, the latter may predominantly contribute to the de novo synthesis of organic carbon which points at the important and yet underestimated role heterotrophic bacterial populations can play the in global carbon cycle/sink in the ocean interior.

Comparative analysis of deep-sea bacterioplankton OMICS revealed the occurrence of habitat-specific genomic attributes

Bathyal aphotic ocean represents the largest biotope on our planet, which sustains highly diverse but low-density microbial communities, with yet untapped genomic attributes, potentially useful for discovery of new biomolecules, industrial enzymes and pathways. In the last two decades, culture-independent approaches of high-throughput sequencing have provided new insights into structure and function of marine bacterioplankton, leading to unprecedented opportunities to accurately characterize microbial communities and their interactions with the environments. In the present review we focused on the analysis of relatively few deep-sea OMICS studies, completed thus far, to find the specific genomic patterns determining the lifeway and adaptation mechanisms of prokaryotes thriving in the dark deep ocean below the depth of 1000m. Phylogenomic and omic studies provided clear evidence that the bathyal microbial communities are distinct from the epipelagic counterparts and, along with generally larger genomes, possess their own habitat-specific genomic attributes. The high abundance in the deep ocean OMICS of the systems for environmental sensing, signal transduction and metabolic versatility as compared to the epipelagic counterparts is thought to enable the deep-sea bacterioplankton to rapidly adapt to changing environmental conditions associated with resource scarcity and high diversity of energy and carbon substrates in the bathyal biotopes. Together with a versatile heterotrophy, mixotrophy and anaplerosis are thought to enable the deep-sea bacterioplankton to cope with these environmental conditions.

Partaking of Archaea to biogeochemical cycling in oxygen-deficient zones of meromictic saline Lake Faro (Messina, Italy)

We used a combination of molecular and microbiological approaches to determine the activity, abundance and diversity of archaeal populations inhabiting meromictic saline Lake Faro (Messina, Italy). Analysis of archaeal 16S rRNA, amoA, accA and hbd genes and transcripts revealed that sub- and anoxic layers of Lake Faro are primarily inhabited by the organisms related to the clusters of Marine Group I.1a of Thaumarchaeota frequently recovered from oxygen-depleted marine ecosystems. These organisms dominated the metabolically active archaea down to the bottom of the lake, indicating their adaptation to recurrent changes in the levels of water column hypoxia. The upper microaerobic layer of Lake Faro redoxcline has the maximal rates of dark primary production much lower than those of other previously studied pelagic redoxclines, but comparable to the values of meso- and bathypelagic areas of Mediterranean Sea. Application of bacterial inhibitors, especially azide, significantly declined the CO2 fixation rates in the low interface and monimolimnion, whereas archaea-specific inhibitor had effect only in upper part of the redoxcline. Based on these findings, we hypothesize that dark bicarbonate fixation in suboxic zone of Lake Faro results mainly from archaeal activity which is affected by the predicted lack in oxygen in lower layers.

Metagenomic analysis of hadopelagic microbial assemblages thriving at the deepest part of Mediterranean Sea, Matapan-Vavilov Deep

The marine pelagic zone situated > 200 m below the sea level (bls) is the largest marine subsystem, comprising more than two-thirds of the oceanic volume. At the same time, it is one of the least explored ecosystems on Earth. Few large-scale environmental genomics studies have been undertaken to examine the phylogenetic diversity and functional gene repertoire of planktonic microbes present in mesopelagic and bathypelagic environments. Here, we present the description of the deep-sea microbial community thriving at > 4900 m depth in Matapan-Vavilov Deep (MVD). This canyon is the deepest site of Mediterranean Sea, with a deepest point located at approximately 5270 m, 56 km SW of city Pylos (Greece) in the Ionian Sea (36°34.00N, 21°07.44E). Comparative analysis of whole-metagenomic data revealed that unlike other deep-sea metagenomes, the prokaryotic diversity in MVD was extremely poor. The decline in the dark primary production rates, measured at 4908 m depth, was coincident with overwhelming dominance of copiotrophic Alteromonas macleodii'deep-ecotype' AltDE at the expense of other prokaryotes including those potentially involved in both autotrophic and anaplerotic CO(2) fixation. We also demonstrate the occurrence in deep-sea metagenomes of several clustered regularly interspaced short palindromic repeats systems.

Community structure and function of planktonic Crenarchaeota: changes with depth in the South China Sea

Marine Crenarchaeota represent a widespread and abundant microbial group in marine ecosystems. Here, we investigated the abundance, diversity, and distribution of planktonic Crenarchaeota in the epi-, meso-, and bathypelagic zones at three stations in the South China Sea (SCS) by analysis of crenarchaeal 16S rRNA gene, ammonia monooxygenase gene amoA involved in ammonia oxidation, and biotin carboxylase gene accA putatively involved in archaeal CO(2) fixation. Quantitative PCR analyses indicated that crenarchaeal amoA and accA gene abundances varied similarly with archaeal and crenarchaeal 16S rRNA gene abundances at all stations, except that crenarchaeal accA genes were almost absent in the epipelagic zone. Ratios of the crenarchaeal amoA gene to 16S rRNA gene abundances decreased ~2.6 times from the epi- to bathypelagic zones, whereas the ratios of crenarchaeal accA gene to marine group I crenarchaeal 16S rRNA gene or to crenarchaeal amoA gene abundances increased with depth, suggesting that the metabolism of Crenarchaeota may change from the epi- to meso- or bathypelagic zones. Denaturing gradient gel electrophoresis profiling of the 16S rRNA genes revealed depth partitioning in archaeal community structures. Clone libraries of crenarchaeal amoA and accA genes showed two clusters: the "shallow" cluster was exclusively derived from epipelagic water and the "deep" cluster was from meso- and/or bathypelagic waters, suggesting that niche partitioning may take place between the shallow and deep marine Crenarchaeota. Overall, our results show strong depth partitioning of crenarchaeal populations in the SCS and suggest a shift in their community structure and ecological function with increasing depth.

The Thaumarchaeota: an emerging view of their phylogeny and ecophysiology

Thaumarchaeota range among the most abundant archaea on Earth. Initially classified as ‘mesophilic Crenarchaeota’, comparative genomics has recently revealed that they form a separate and deep-branching phylum within the Archaea. This novel phylum comprises in 16S rRNA gene trees not only all known archaeal ammonia oxidizers but also several clusters of environmental sequences representing microorganisms with unknown energy metabolism. Ecophysiological studies of ammonia-oxidizing Thaumarchaeota suggest adaptation to low ammonia concentrations and an autotrophic or possibly mixotrophic lifestyle. Extrapolating from the wide substrate range of copper-containing membrane-bound monooxygenases, to which the thaumarchaeal ammonia monooxygenases belong, the use of substrates other than ammonia for generating energy by some members of the Thaumarchaeota seems likely.

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.

Contribution of crenarchaeal autotrophic ammonia oxidizers to the dark primary production in Tyrrhenian deep waters (Central Mediterranean Sea)

Mesophilic Crenarchaeota have recently been thought to be significant contributors to nitrogen (N) and carbon (C) cycling. In this study, we examined the vertical distribution of ammonia-oxidizing Crenarchaeota at offshore site in Southern Tyrrhenian Sea. The median value of the crenachaeal cell to amoA gene ratio was close to one suggesting that virtually all deep-sea Crenarchaeota possess the capacity to oxidize ammonia. Crenarchaea-specific genes, nirK and ureC, for nitrite reductase and urease were identified and their affiliation demonstrated the presence of 'deep-sea' clades distinct from 'shallow' representatives. Measured deep-sea dark CO(2) fixation estimates were comparable to the median value of photosynthetic biomass production calculated for this area of Tyrrhenian Sea, pointing to the significance of this process in the C cycle of aphotic marine ecosystems. To elucidate the pivotal organisms in this process, we targeted known marine crenarchaeal autotrophy-related genes, coding for acetyl-CoA carboxylase (accA) and 4-hydroxybutyryl-CoA dehydratase (4-hbd). As in case of nirK and ureC, these genes are grouped with deep-sea sequences being distantly related to those retrieved from the epipelagic zone. To pair the molecular data with specific functional attributes we performed [(14)C]HCO(3) incorporation experiments followed by analyses of radiolabeled proteins using shotgun proteomics approach. More than 100 oligopeptides were attributed to 40 marine crenarchaeal-specific proteins that are involved in 10 different metabolic processes, including autotrophy. Obtained results provided a clear proof of chemolithoautotrophic physiology of bathypelagic crenarchaeota and indicated that this numerically predominant group of microorganisms facilitate a hitherto unrecognized sink for inorganic C of a global importance.

Identification of missing genes and enzymes for autotrophic carbon fixation in crenarchaeota

Two autotrophic carbon fixation cycles have been identified in Crenarchaeota. The dicarboxylate/4-hydroxybutyrate cycle functions in anaerobic or microaerobic autotrophic members of the Thermoproteales and Desulfurococcales. The 3-hydroxypropionate/4-hydroxybutyrate cycle occurs in aerobic autotrophic Sulfolobales; a similar cycle may operate in autotrophic aerobic marine Crenarchaeota. Both cycles form succinyl-coenzyme A (CoA) from acetyl-CoA and two molecules of inorganic carbon, but they use different means. Both cycles have in common the (re)generation of acetyl-CoA from succinyl-CoA via identical intermediates. Here, we identified several missing enzymes/genes involved in the seven-step conversion of succinyl-CoA to two molecules of acetyl-CoA in Thermoproteus neutrophilus (Thermoproteales), Ignicoccus hospitalis (Desulfurococcales), and Metallosphaera sedula (Sulfolobales). The identified enzymes/genes include succinyl-CoA reductase, succinic semialdehyde reductase, 4-hydroxybutyrate-CoA ligase, bifunctional crotonyl-CoA hydratase/(S)-3-hydroxybutyryl-CoA dehydrogenase, and beta-ketothiolase. 4-Hydroxybutyryl-CoA dehydratase, which catalyzes a mechanistically intriguing elimination of water, is well conserved and rightly can be considered the key enzyme of these two cycles. In contrast, several of the other enzymes evolved from quite different sources, making functional predictions based solely on genome interpretation difficult, if not questionable.