Lateral Gene Transfer Shapes the Distribution of RuBisCO among Candidate Phyla Radiation Bacteria and DPANN Archaea

Intimate microbe-water-mineral interactions mediate alkalization in the pyroxene-rich iron ore mines in Panxi area, Southwest China

In contrast to acid mine drainage, the microbial assembly and (bio)geochemical processes in alkaline mine conditions remain under-investigated. Here, microbe-water-mineral interactions were systematically investigated in two representative iron mines with alkaline conditions in the Panxi mining area, Southwest China. Compared to reference riverine samples less interfered by mining activities, the iron ore samples, composed of vanadium-titanium magnetite and pyroxene-rich bedrocks, exhibited elevated levels of Fe, HCl-extractable Fe(II), total sulfur, nitrate and sulfate, but lower total carbon (TC). Meanwhile, the mine drainage showed significantly higher sulfate, but lower TC concentrations than the riverine samples. Intriguingly, the Serpentinimonas spp., typically reported in serpentinites, prevailed in the microbial communities from the mine samples exhibiting higher pH. This suggests that the alkaline environments in Panxi mines result from serpentinization-like reactions. Enrichment of Thiobacillus spp. was observed in the mine-dwelling microbial communities, positively correlated with total sulfur, sulfate, nitrate, and Fe(II). Genome-resolved metagenomics suggested a chemoautotrophic lifestyle for the Thiobacillus species (e.g., carbon fixation, sulfur oxidation, and oxygen respiration), which may generate H+ and mitigate alkalization. This study provides valuable insights into progressive development of alkaline mine ecosystems and offers guidance for developing appropriate engineering strategies to restore the abandoned alkaline mines.

Metagenomic analysis sheds light on the mixotrophic lifestyle of bacterial phylum Zhuqueibacterota

Zhuqueibacterota is a novel bacterial phylum proposed based on hot spring metagenomes and public metagenome-assembled genomes, classified within the Fibrobacterota-Chlorobiota-Bacteroidota superphylum. This globally distributed phylum consists of one class and five orders, with the majority of its members being facultative anaerobes. Notably, the order Zhuqueibacterales utilizes hydrogen as an electron donor for carbon fixation through the Calvin Benson Bassham cycle. Phylogenetic and metabolic analyses reveal the phylum's key role in the carbon cycle, with frequent horizontal gene transfer events influencing its evolutionary trajectory.

Insights into chemoautotrophic traits of a prevalent bacterial phylum CSP1-3, herein Sysuimicrobiota

Candidate bacterial phylum CSP1-3 has not been cultivated and is poorly understood. Here, we analyzed 112 CSP1-3 metagenome-assembled genomes and showed they are likely facultative anaerobes, with 3 of 5 families encoding autotrophy through the reductive glycine pathway (RGP), Wood-Ljungdahl pathway (WLP) or Calvin-Benson-Bassham (CBB), with hydrogen or sulfide as electron donors. Chemoautotrophic enrichments from hot spring sediments and fluorescence in situ hybridization revealed enrichment of six CSP1-3 genera, and both transcribed genes and DNA-stable isotope probing were consistent with proposed chemoautotrophic metabolisms. Ancestral state reconstructions showed that the ancestors of phylum CSP1-3 may have been acetogens that were autotrophic via the RGP, whereas the WLP and CBB were acquired by horizontal gene transfer. Our results reveal that CSP1-3 is a widely distributed phylum with the potential to contribute to the cycling of carbon, sulfur and nitrogen. The name Sysuimicrobiota phy. nov. is proposed.

Ribulose 1,5-Bisphosphate Carboxylase/Oxygenase Is Required in Bradyrhizobium diazoefficiens for Efficient Soybean Root Colonization and Competition for Nodulation

The Hyphomicrobiales (Rhizobiales) order contains soil bacteria with an irregular distribution of the Calvin-Benson-Bassham cycle (CBB). Key enzymes in the CBB cycle are ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO), whose large and small subunits are encoded in cbbL and cbbS, and phosphoribulokinase (PRK), encoded by cbbP. These genes are often found in cbb operons, regulated by the LysR-type regulator CbbR. In Bradyrhizobium, pertaining to this order and bearing photosynthetic and non-photosynthetic species, the number of cbbL and cbbS copies varies, for example: zero in B. manausense, one in B. diazoefficiens, two in B. japonicum, and three in Bradyrhizobium sp. BTAi. Few studies addressed the role of CBB in Bradyrhizobium spp. symbiosis with leguminous plants. To investigate the horizontal transfer of the cbb operon among Hyphomicrobiales, we compared phylogenetic trees for concatenated cbbL-cbbP-cbbR and housekeeping genes (atpD-gyrB-recA-rpoB-rpoD). The distribution was consistent, indicating no horizontal transfer of the cbb operon in Hyphomicrobiales. We constructed a ΔcbbLS mutant in B. diazoefficiens, which lost most of the coding sequence of cbbL and has a frameshift creating a stop codon at the N-terminus of cbbS. This mutant nodulated normally but had reduced competitiveness for nodulation and long-term adhesion to soybean (Glycine max (L.) Merr.) roots, indicating a CBB requirement for colonizing soybean rhizosphere.

Analysis of nearly 3000 archaeal genomes from terrestrial geothermal springs sheds light on interconnected biogeochemical processes

Terrestrial geothermal springs are physicochemically diverse and host abundant populations of Archaea. However, the diversity, functionality, and geological influences of these Archaea are not well understood. Here we explore the genomic diversity of Archaea in 152 metagenomes from 48 geothermal springs in Tengchong, China, collected from 2016 to 2021. Our dataset is comprised of 2949 archaeal metagenome-assembled genomes spanning 12 phyla and 392 newly identified species, which increases the known species diversity of Archaea by ~48.6%. The structures and potential functions of the archaeal communities are strongly influenced by temperature and pH, with high-temperature acidic and alkaline springs favoring archaeal abundance over Bacteria. Genome-resolved metagenomics and metatranscriptomics provide insights into the potential ecological niches of these Archaea and their potential roles in carbon, sulfur, nitrogen, and hydrogen metabolism. Furthermore, our findings illustrate the interplay of competition and cooperation among Archaea in biogeochemical cycles, possibly arising from overlapping functional niches and metabolic handoffs. Taken together, our study expands the genomic diversity of Archaea inhabiting geothermal springs and provides a foundation for more incisive study of biogeochemical processes mediated by Archaea in geothermal ecosystems.

Considerations for Detecting Organic Indicators of Metabolism on Enceladus

Enceladus is of interest to astrobiology and the search for life since it is thought to host active hydrothermal activity and habitable conditions. It is also possible that the organics detected on Enceladus may indicate an active prebiotic or biotic system; in particular, the conditions on Enceladus may favor mineral-driven protometabolic reactions. When including metabolism-related biosignatures in Enceladus mission concepts, it is necessary to base these in a clearer understanding of how these signatures could also be produced prebiotically. In addition, postulating which biological metabolisms to look for on Enceladus requires a non-Earth-centric approach since the details of biological metabolic pathways are heavily shaped by adaptation to geochemical conditions over the planet's history. Creating metabolism-related organic detection objectives for Enceladus missions, therefore, requires consideration of how metabolic systems may operate differently on another world, while basing these speculations on observed Earth-specific microbial processes. In addition, advances in origin-of-life research can play a critical role in distinguishing between interpretations of any future organic detections on Enceladus, and the discovery of an extant prebiotic system would be a transformative astrobiological event in its own right.

Reversed oxidative TCA (roTCA) for carbon fixation by an Acidimicrobiia strain from a saline lake

Acidimicrobiia are widely distributed in nature and suggested to be autotrophic via the Calvin-Benson-Bassham (CBB) cycle. However, direct evidence of chemolithoautotrophy in Acidimicrobiia is lacking. Here, we report a chemolithoautotrophic enrichment from a saline lake, and the subsequent isolation and characterization of a chemolithoautotroph, Salinilacustristhrix flava EGI L10123T, which belongs to a new Acidimicrobiia family. Although strain EGI L10123T is autotrophic, neither its genome nor Acidimicrobiia metagenome-assembled genomes from the enrichment culture encode genes necessary for the CBB cycle. Instead, genomic, transcriptomic, enzymatic, and stable-isotope probing data hinted at the activity of the reversed oxidative TCA (roTCA) coupled with the oxidation of sulfide as the electron donor. Phylogenetic analysis and ancestral character reconstructions of Acidimicrobiia suggested that the essential CBB gene rbcL was acquired through multiple horizontal gene transfer events from diverse microbial taxa. In contrast, genes responsible for sulfide- or hydrogen-dependent roTCA carbon fixation were already present in the last common ancestor of extant Acidimicrobiia. These findings imply the possibility of roTCA carbon fixation in Acidimicrobiia and the ecological importance of Acidimicrobiia. Further research in the future is necessary to confirm whether these characteristics are truly widespread across the clade.

Deep-branching evolutionary intermediates reveal structural origins of form I rubisco

The enzyme rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) catalyzes the majority of biological carbon fixation on Earth. Although the vast majority of rubiscos across the tree of life assemble as homo-oligomers, the globally predominant form I enzyme-found in plants, algae, and cyanobacteria-forms a unique hetero-oligomeric complex. The recent discovery of a homo-oligomeric sister group to form I rubisco (named form I') has filled a key gap in our understanding of the enigmatic origins of the form I clade. However, to elucidate the series of molecular events leading to the evolution of form I rubisco, we must examine more distantly related sibling clades to contextualize the molecular features distinguishing form I and form I' rubiscos. Here, we present a comparative structural study retracing the evolutionary history of rubisco that reveals a complex structural trajectory leading to the ultimate hetero-oligomerization of the form I clade. We structurally characterize the oligomeric states of deep-branching form Iα and I'' rubiscos recently discovered from metagenomes, which represent key evolutionary intermediates preceding the form I clade. We further solve the structure of form I'' rubisco, revealing the molecular determinants that likely primed the enzyme core for the transition from a homo-oligomer to a hetero-oligomer. Our findings yield new insight into the evolutionary trajectory underpinning the adoption and entrenchment of the prevalent assembly of form I rubisco, providing additional context when viewing the enzyme family through the broader lens of protein evolution.

Globally distributed Myxococcota with photosynthesis gene clusters illuminate the origin and evolution of a potentially chimeric lifestyle

Photosynthesis is a fundamental biogeochemical process, thought to be restricted to a few bacterial and eukaryotic phyla. However, understanding the origin and evolution of phototrophic organisms can be impeded and biased by the difficulties of cultivation. Here, we analyzed metagenomic datasets and found potential photosynthetic abilities encoded in the genomes of uncultivated bacteria within the phylum Myxococcota. A putative photosynthesis gene cluster encoding a type-II reaction center appears in at least six Myxococcota families from three classes, suggesting vertical inheritance of these genes from an early common ancestor, with multiple independent losses in other lineages. Analysis of metatranscriptomic datasets indicate that the putative myxococcotal photosynthesis genes are actively expressed in various natural environments. Furthermore, heterologous expression of myxococcotal pigment biosynthesis genes in a purple bacterium supports that the genes can drive photosynthetic processes. Given that predatory abilities are thought to be widespread across Myxococcota, our results suggest the intriguing possibility of a chimeric lifestyle (combining predatory and photosynthetic abilities) in members of this phylum.

Metagenomic Discovery of "Candidatus Parvarchaeales"-Related Lineages Sheds Light on Adaptation and Diversification from Neutral-Thermal to Acidic-Mesothermal Environments

"Candidatus Parvarchaeales" microbes, representing a DPANN archaeal group with limited metabolic potential and reliance on hosts for their growth, were initially found in acid mine drainage (AMD). Due to the lack of representatives, however, their ecological roles and adaptation to extreme habitats such as AMD as well as how they diverge across the lineage remain largely unexplored. By applying genome-resolved metagenomics, 28 Parvarchaeales-associated metagenome-assembled genomes (MAGs) representing two orders and five genera were recovered. Among them, we identified three new genera and proposed the names "Candidatus Jingweiarchaeum," "Candidatus Haiyanarchaeum," and "Candidatus Rehaiarchaeum," with the former two belonging to a new order, "Candidatus Jingweiarchaeales." Further analyses of the metabolic potentials revealed substantial niche differentiation between Jingweiarchaeales and Parvarchaeales. Jingweiarchaeales may rely on fermentation, salvage pathways, partial glycolysis, and the pentose phosphate pathway (PPP) for energy conservation reservation, while the metabolic potentials of Parvarchaeales might be more versatile. Comparative genomic analyses suggested that Jingweiarchaeales favor habitats with higher temperatures and that Parvarchaeales are better adapted to acidic environments. We further revealed that the thermal adaptation of these lineages, especially Haiyanarchaeum, might rely on genomic features such as the usage of specific amino acids, genome streamlining, and hyperthermophile featured genes such as rgy. Notably, the adaptation of Parvarchaeales to acidic environments was possibly driven by horizontal gene transfer (HGT). The reconstruction of ancestral states demonstrated that both may have originated from thermal and neutral environments and later spread to mesothermal and acidic environments. These evolutionary processes may also be accompanied by adaptation to oxygen-rich environments via HGT. IMPORTANCE "Candidatus Parvarchaeales" microbes may represent a lineage uniquely distributed in extreme environments such as AMD and hot springs. However, little is known about the strategies and processes of how they adapted to these extreme environments. By the discovery of potential new order-level lineages, "Ca. Jingweiarchaeales," and in-depth comparative genomic analysis, we unveiled the functional differentiation of these lineages. Furthermore, we show that the adaptation of these lineages to high-temperature and acidic environments was driven by different strategies, with the former relying more on genomic characteristics such as genome streamlining and amino acid compositions and the latter relying more on the acquisition of genes associated with acid tolerance. Finally, by the reconstruction of the ancestral states of the optimal growth temperature (OGT) and isoelectric point (pI), we showed the potential evolutionary process of Parvarchaeales-related lineages with regard to the shift from the high-temperature environment of their common ancestors to low-temperature (potentially acidic) environments.

Microbial life in 25-m-deep boreholes in ancient permafrost illuminated by metagenomics

This study describes the composition and potential metabolic adaptation of microbial communities in northeastern Siberia, a repository of the oldest permafrost in the Northern Hemisphere. Samples of contrasting depth (1.75 to 25.1 m below surface), age (from ~ 10 kyr to 1.1 Myr) and salinity (from low 0.1-0.2 ppt and brackish 0.3-1.3 ppt to saline 6.1 ppt) were collected from freshwater permafrost (FP) of borehole AL1_15 on the Alazeya River, and coastal brackish permafrost (BP) overlying marine permafrost (MP) of borehole CH1_17 on the East Siberian Sea coast. To avoid the limited view provided with culturing work, we used 16S rRNA gene sequencing to show that the biodiversity decreased dramatically with permafrost age. Nonmetric multidimensional scaling (NMDS) analysis placed the samples into three groups: FP and BP together (10-100 kyr old), MP (105-120 kyr old), and FP (> 900 kyr old). Younger FP/BP deposits were distinguished by the presence of Acidobacteriota, Bacteroidota, Chloroflexota_A, and Gemmatimonadota, older FP deposits had a higher proportion of Gammaproteobacteria, and older MP deposits had much more uncultured groups within Asgardarchaeota, Crenarchaeota, Chloroflexota, Patescibacteria, and unassigned archaea. The 60 recovered metagenome-assembled genomes and un-binned metagenomic assemblies suggested that despite the large taxonomic differences between samples, they all had a wide range of taxa capable of fermentation coupled to nitrate utilization, with the exception of sulfur reduction present only in old MP deposits.

Genome-resolved metagenomics reveals depth-related patterns of microbial community structure and functions in a highly stratified, AMD overlaying mine tailings

Acid mine drainage (AMD) is a worldwide environmental problem, yet bioremediation is hampered by a limited knowledge of the reductive microbial processes in the AMD ecosystem. Here, we generate extensive metagenome and geochemical datasets to investigate how microbial populations and metabolic capacities driving major element cycles are structured in a highly stratified, AMD overlaying tailings environment. The results demonstrated an explicit depth-dependent differentiation of microbial community composition and function profiles between the surface and deeper tailings layers, paralleling the dramatic shifts in major physical and geochemical properties. Specifically, key genes involved in sulfur and iron oxidation were significantly enriched in the surface tailings, whereas those associated with reductive nitrogen, sulfur, and iron processes were enriched in the deeper layers. Genome-resolved metagenomics retrieved 406 intermediate or high-quality genomes spanning 26 phyla, including major new groups (e.g., Patescibacteria and DPANN). Metabolic models involving nitrogen, sulfur, iron, and carbon cycles were proposed based on the functional potentials of the abundant microbial genomes, emphasizing syntrophy and the importance of lesser-known taxa in the degradation of complex carbon compounds. These results have implications for in situ AMD bioremediation.

Variable impact of geochemical gradients on the functional potential of bacteria, archaea, and phages from the permanently stratified Lac Pavin

BACKGROUND

Permanently stratified lakes contain diverse microbial communities that vary with depth and so serve as useful models for studying the relationships between microbial community structure and geochemistry. Recent work has shown that these lakes can also harbor numerous bacteria and archaea from novel lineages, including those from the Candidate Phyla Radiation (CPR). However, the extent to which geochemical stratification differentially impacts carbon metabolism and overall genetic potential in CPR bacteria compared to other organisms is not well defined.

RESULTS

Here, we determine the distribution of microbial lineages along an oxygen gradient in Lac Pavin, a deep, stratified lake in central France, and examine the influence of this gradient on their metabolism. Genome-based analyses revealed an enrichment of distinct C1 and CO2 fixation pathways in the oxic lake interface and anoxic zone/sediments, suggesting that oxygen likely plays a role in structuring metabolic strategies in non-CPR bacteria and archaea. Notably, we find that the oxidation of methane and its byproducts is largely spatially separated from methane production, which is mediated by diverse communities of sediment methanogens that vary on the centimeter scale. In contrast, we detected evidence for RuBisCO throughout the water column and sediments, including form II/III and form III-related enzymes encoded by CPR bacteria in the water column and DPANN archaea in the sediments. On the whole, though, CPR bacteria and phages did not show strong signals of gene content differentiation by depth, despite the fact that distinct species groups populate different lake and sediment compartments.

CONCLUSIONS

Overall, our analyses suggest that environmental gradients in Lac Pavin select for capacities of CPR bacteria and phages to a lesser extent than for other bacteria and archaea. This may be due to the fact that selection in the former groups is indirect and depends primarily on host characteristics. Video Abstract.

Photorespiration - Rubisco's repair crew

The photorespiratory repair pathway (photorespiration in short) was set up from ancient metabolic modules about three billion years ago in cyanobacteria, the later ancestors of chloroplasts. These prokaryotes developed the capacity for oxygenic photosynthesis, i.e. the use of water as a source of electrons and protons (with O₂ as a by-product) for the sunlight-driven synthesis of ATP and NADPH for CO₂ fixation in the Calvin cycle. However, the CO₂-binding enzyme, ribulose 1,5-bisphosphate carboxylase (known under the acronym Rubisco), is not absolutely selective for CO₂ and can also use O₂ in a side reaction. It then produces 2-phosphoglycolate (2PG), the accumulation of which would inhibit and potentially stop the Calvin cycle and subsequently photosynthetic electron transport. Photorespiration removes the 2-PG and in this way prevents oxygenic photosynthesis from poisoning itself. In plants, the core of photorespiration consists of ten enzymes distributed over three different types of organelles, requiring interorganellar transport and interaction with several auxiliary enzymes. It goes together with the release and to some extent loss of freshly fixed CO₂. This disadvantageous feature can be suppressed by CO₂-concentrating mechanisms, such as those that evolved in C₄ plants thirty million years ago, which enhance CO₂ fixation and reduce 2PG synthesis. Photorespiration itself provided a pioneer variant of such mechanisms in the predecessors of C₄ plants, C₃-C₄ intermediate plants. This article is a review and update particularly on the enzyme components of plant photorespiration and their catalytic mechanisms, on the interaction of photorespiration with other metabolism and on its impact on the evolution of photosynthesis. This focus was chosen because a better knowledge of the enzymes involved and how they are embedded in overall plant metabolism can facilitate the targeted use of the now highly advanced methods of metabolic network modelling and flux analysis. Understanding photorespiration more than before as a process that enables, rather than reduces, plant photosynthesis, will help develop rational strategies for crop improvement.

Characterization of cyanobacterial isolates from freshwater and saline subtropical desert lakes

Characterization of Cyanobacteria in lakes with different physicochemical properties provides insights into the diversity of this phylum and knowledge of their features that are relevant to biotechnology applications. Six Cyanobacterial isolates were recovered from freshwater Lake Nasser and saline Lake Qarun, Egypt. The isolates were identified based on both morphology and molecular markers, 16S rRNA, and RuBisCO cbbL genes. The isolates SN1, SN2, SN3, SN4, Q1, and Q2 showed homologies with Merismopedia, Oscillatoria, Limnothrix, Persinema, and Jacksonvillea, respectively. The cbbL sequences for isolates SN1, Q1, and Q2 represented the first records for candidates relating to the genera Merismopedia and Persinema, and Jacksonvillea, respectively. Biochemical contents, carbohydrates, proteins, lipids, pigments, and ash-free dry weight were measured for each isolate. Isolate SN2 had the highest content of allophycocyanin, 71 ± 4.8 mg/g DW, and phycoerythrin, 98 ± 6.7 mg/g DW, while the isolate SN4 had the highest composition of total protein, lipid, carotenoid, and chlorophyll a, recording 364.7 ± 6.4 mg/g DW, 67.6 ± 0.2 mg/g DW, 0.261 ± 0.01 mg/g DW, and 10 ± 0.6 mg/g DW, respectively. Isolate Q1 recorded the maximum amount of phycocyanin, 114 ± 20.7 mg/g DW among isolates. The isolate Q2 was observed to have the highest carbohydrate content, 274 ± 14.5 (mg/g DW), and ash-free dry weight, 891.8 ± 2.8 mg/g DW. Thus, the study indicated that the current isolates may represent promising resources for biotechnological applications.

Functional differentiation determines the molecular basis of the symbiotic lifestyle of Ca. Nanohaloarchaeota

BACKGROUND

Candidatus Nanohaloarchaeota, an archaeal phylum within the DPANN superphylum, is characterized by limited metabolic capabilities and limited phylogenetic diversity and until recently has been considered to exclusively inhabit hypersaline environments due to an obligate association with Halobacteria. Aside from hypersaline environments, Ca. Nanohaloarchaeota can also have been discovered from deep-subsurface marine sediments.

RESULTS

Three metagenome-assembled genomes (MAGs) representing a new order within the Ca. Nanohaloarchaeota were reconstructed from a stratified salt crust and proposed to represent a novel order, Nucleotidisoterales. Genomic features reveal them to be anaerobes capable of catabolizing nucleotides by coupling nucleotide salvage pathways with lower glycolysis to yield free energy. Comparative genomics demonstrated that these and other Ca. Nanohaloarchaeota inhabiting saline habitats use a "salt-in" strategy to maintain osmotic pressure based on the high proportion of acidic amino acids. In contrast, previously described Ca. Nanohaloarchaeota MAGs from geothermal environments were enriched with basic amino acids to counter heat stress. Evolutionary history reconstruction revealed that functional differentiation of energy conservation strategies drove diversification within Ca. Nanohaloarchaeota, further leading to shifts in the catabolic strategy from nucleotide degradation within deeper lineages to polysaccharide degradation within shallow lineages.

CONCLUSIONS

This study provides deeper insight into the ecological functions and evolution of the expanded phylum Ca. Nanohaloarchaeota and further advances our understanding on the functional and genetic associations between potential symbionts and hosts. Video Abstract.

Compositional and Metabolic Responses of Autotrophic Microbial Community to Salinity in Lacustrine Environments

The compositional and physiological responses of autotrophic microbiotas to salinity in lakes remain unclear. In this study, the community composition and carbon fixation pathways of autotrophic microorganisms in lacustrine sediments with a salinity gradient (82.6 g/L to 0.54 g/L) were investigated by using metagenomic analysis. A total of 117 metagenome-assembled genomes (MAGs) with carbon fixation potentially belonging to 20 phyla were obtained. The abundance of these potential autotrophs increased significantly with decreasing salinity, and the variation of sediment autotrophic microbial communities was mainly affected by salinity, pH, and total organic carbon. Notably, along the decreasing salinity gradient, the dominant lineage shifted from Desulfobacterota to Proteobacteria. Meanwhile, the dominant carbon fixation pathway shifted from the Wood-Lungdahl pathway to the less-energy-efficient Calvin-Benson-Bassham cycle, with glycolysis shifting from the Embden-Meyerhof-Parnas pathway to the less-exergonic Entner-Doudoroff pathway. These results suggest that the physiological efficiency of autotrophic microorganisms decreased when the environmental salinity became lower. Metabolic inference of these MAGs revealed that carbon fixation may be coupled to the oxidation of reduced sulfur compounds and ferrous iron, dissimilatory nitrate reduction at low salinity, and dissimilatory sulfate reduction in hypersaline sediments. These results extend our understanding of metabolic versatility and niche diversity of autotrophic microorganisms in saline environments and shed light on the response of autotrophic microbiomes to salinity. These findings are of great significance for understanding the impact of desalination caused by climate warming on the carbon cycle of saline lake ecosystems. IMPORTANCE The Qinghai-Tibetan lakes are experiencing water increase and salinity decrease due to climate warming. However, little is known about how the salinity decrease will affect the composition of autotrophic microbial populations and their carbon fixation pathways. In this study, we used genome-resolved metagenomics to interpret the dynamic changes in the autotrophic microbial community and metabolic pathways along a salinity gradient. The results showed that desalination drove the shift of the dominant microbial lineage from Desulfobacterota to Proteobacteria, enriched autotrophs with lower physiological efficiency pathways, and enhanced coupling between the carbon cycle and other element cycles. These results can predict the future response of microbial communities to lake desalination and improve our understanding of the effect of climate warming on the carbon cycle in saline aquatic ecosystems.

Environmental Galenics: large-scale fortification of extant microbiomes with engineered bioremediation agents

Contemporary synthetic biology-based biotechnologies are generating tools and strategies for reprogramming genomes for specific purposes, including improvement and/or creation of microbial processes for tackling climate change. While such activities typically work well at a laboratory or bioreactor scale, the challenge of their extensive delivery to multiple spatio-temporal dimensions has hardly been tackled thus far. This state of affairs creates a research niche for what could be called Environmental Galenics (EG), i.e. the science and technology of releasing designed biological agents into deteriorated ecosystems for the sake of their safe and effective recovery. Such endeavour asks not just for an optimal performance of the biological activity at stake, but also the material form and formulation of the agents, their propagation and their interplay with the physico-chemical scenario where they are expected to perform. EG also encompasses adopting available physical carriers of microorganisms and channels of horizontal gene transfer as potential paths for spreading beneficial activities through environmental microbiomes. While some of these propositions may sound unsettling to anti-genetically modified organisms sensitivities, they may also fall under the tag of TINA (there is no alternative) technologies in the cases where a mere reduction of emissions will not help the revitalization of irreversibly lost ecosystems.

This article is part of the theme issue ‘Ecological complexity and the biosphere: the next 30 years’.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

Improving microalgae for biotechnology - From genetics to synthetic biology - Moving forward but not there yet

Microalgae are a diverse group of photosynthetic organisms that can be exploited for the production of different compounds, ranging from crude biomass and biofuels to high value-added biochemicals and synthetic proteins. Traditionally, algal biotechnology relies on bioprospecting to identify new highly productive strains and more recently, on forward genetics to further enhance productivity. However, it has become clear that further improvements in algal productivity for biotechnology is impossible without combining traditional tools with the arising molecular genetics toolkit. We review recent advantages in developing high throughput screening methods, preparing genome-wide mutant libraries, and establishing genome editing techniques. We discuss how algae can be improved in terms of photosynthetic efficiency, biofuel and high value-added compound production. Finally, we critically evaluate developments over recent years and explore future potential in the field.

Leave no stone unturned: individually adapted xerotolerant Thaumarchaeota sheltered below the boulders of the Atacama Desert hyperarid core

BACKGROUND

The hyperarid core of the Atacama Desert is an extremely harsh environment thought to be colonized by only a few heterotrophic bacterial species. Current concepts for understanding this extreme ecosystem are mainly based on the diversity of these few species, yet a substantial area of the Atacama Desert hyperarid topsoil is covered by expansive boulder accumulations, whose underlying microbiomes have not been investigated so far. With the hypothesis that these sheltered soils harbor uniquely adapted microbiomes, we compared metagenomes and geochemistry between soils below and beside boulders across three distantly located boulder accumulations in the Atacama Desert hyperarid core.

RESULTS

Genome-resolved metagenomics of eleven samples revealed substantially different microbial communities in soils below and beside boulders, despite the presence of shared species. Archaea were found in significantly higher relative abundance below the boulders across all samples within distances of up to 205 km. These key taxa belong to a novel genus of ammonia-oxidizing Thaumarchaeota, Candidatus Nitrosodeserticola. We resolved eight mid-to-high quality genomes of this genus and used comparative genomics to analyze its pangenome and site-specific adaptations. Ca. Nitrosodeserticola genomes contain genes for ammonia oxidation, the 3-hydroxypropionate/4-hydroxybutyrate carbon fixation pathway, and acetate utilization indicating a chemolithoautotrophic and mixotrophic lifestyle. They also possess the capacity for tolerating extreme environmental conditions as highlighted by the presence of genes against oxidative stress and DNA damage. Site-specific adaptations of the genomes included the presence of additional genes for heavy metal transporters, multiple types of ATP synthases, and divergent genes for aquaporins.

CONCLUSION

We provide the first genomic characterization of hyperarid soil microbiomes below the boulders in the Atacama Desert, and report abundant and highly adapted Thaumarchaeaota with ammonia oxidation and carbon fixation potential. Ca. Nitrosodeserticola genomes provide the first metabolic and physiological insight into a thaumarchaeal lineage found in globally distributed terrestrial habitats characterized by various environmental stresses. We consequently expand not only the known genetic repertoire of Thaumarchaeota but also the diversity and microbiome functioning in hyperarid ecosystems. Video Abstract.

Compendium of 530 metagenome-assembled bacterial and archaeal genomes from the polar Arctic Ocean

The role of the Arctic Ocean ecosystem in climate regulation may depend on the responses of marine microorganisms to environmental change. We applied genome-resolved metagenomics to 41 Arctic seawater samples, collected at various depths in different seasons during the Tara Oceans Polar Circle expedition, to evaluate the ecology, metabolic potential and activity of resident bacteria and archaea. We assembled 530 metagenome-assembled genomes (MAGs) to form the Arctic MAGs catalogue comprising 526 species. A total of 441 MAGs belonged to species that have not previously been reported and 299 genomes showed an exclusively polar distribution. Most Arctic MAGs have large genomes and the potential for fast generation times, both of which may enable adaptation to a copiotrophic lifestyle in nutrient-rich waters. We identified 38 habitat generalists and 111 specialists in the Arctic Ocean. We also found a general prevalence of 14 mixotrophs, while chemolithoautotrophs were mostly present in the mesopelagic layer during spring and autumn. We revealed 62 MAGs classified as key Arctic species, found only in the Arctic Ocean, showing the highest gene expression values and predicted to have habitat-specific traits. The Artic MAGs catalogue will inform our understanding of polar microorganisms that drive global biogeochemical cycles.

Insight into the function and evolution of the Wood-Ljungdahl pathway in Actinobacteria

Carbon fixation by chemoautotrophic microbes such as homoacetogens had a major impact on the transition from the inorganic to the organic world. Recent reports have shown the presence of genes for key enzymes associated with the Wood-Ljungdahl pathway (WLP) in the phylum Actinobacteria, which adds to the diversity of potential autotrophs. Here, we compiled 42 actinobacterial metagenome-assembled genomes (MAGs) from new and existing metagenomic datasets and propose three novel classes, Ca. Aquicultoria, Ca. Geothermincolia and Ca. Humimicrobiia. Most members of these classes contain genes coding for acetogenesis through the WLP, as well as a variety of hydrogenases (NiFe groups 1a and 3b-3d; FeFe group C; NiFe group 4-related hydrogenases). We show that the three classes acquired the hydrogenases independently, yet the carbon monoxide dehydrogenase/acetyl-CoA synthase complex (CODH/ACS) was apparently present in their last common ancestor and was inherited vertically. Furthermore, the Actinobacteria likely donated genes for CODH/ACS to multiple lineages within Nitrospirae, Deltaproteobacteria (Desulfobacterota), and Thermodesulfobacteria through multiple horizontal gene transfer events. Finally, we show the apparent growth of Ca. Geothermincolia and H2-dependent acetate production in hot spring enrichment cultures with or without the methanogenesis inhibitor 2-bromoethanesulfonate, which is consistent with the proposed homoacetogenic metabolism.

Comparative Metagenomics of the Active Layer and Permafrost from Low-Carbon Soil in the Canadian High Arctic

Approximately 87% of the Arctic consists of low-organic carbon mineral soil, but knowledge of microbial activity in low-carbon permafrost (PF) and active layer soils remains limited. This study investigated the taxonomic composition and genetic potential of microbial communities at contrasting depths of the active layer (5, 35, and 65 cm below surface, bls) and PF (80 cm bls). We showed microbial communities in PF to be taxonomically and functionally different from those in the active layer. 16S rRNA gene sequence analysis revealed higher biodiversity in the active layer than in PF, and biodiversity decreased significantly with depth. The reconstructed 91 metagenome-assembled genomes showed that PF was dominated by heterotrophic, fermenting Bacteroidota using nitrite as their main electron acceptor. Prevalent microbes identified in the active layer belonged to bacterial taxa, gaining energy via aerobic respiration. Gene abundance in metagenomes revealed enrichment of genes encoding the plant-derived polysaccharide degradation and metabolism of nitrate and sulfate in PF, whereas genes encoding methane/ammonia oxidation, cold-shock protein, and two-component systems were generally more abundant in the active layer, particularly at 5 cm bls. The results of this study deepen our understanding of the low-carbon Arctic soil microbiome and improve prediction of the impacts of thawing PF.

Comparative genomic analysis reveals metabolic flexibility of Woesearchaeota

The archaeal phylum Woesearchaeota, within the DPANN superphylum, includes phylogenetically diverse microorganisms that inhabit various environments. Their biology is poorly understood due to the lack of cultured isolates. Here, we analyze datasets of Woesearchaeota 16S rRNA gene sequences and metagenome-assembled genomes to infer global distribution patterns, ecological preferences and metabolic capabilities. Phylogenomic analyses indicate that the phylum can be classified into ten subgroups, termed A-J. While a symbiotic lifestyle is predicted for most, some members of subgroup J might be host-independent. The genomes of several Woesearchaeota, including subgroup J, encode putative [FeFe] hydrogenases (known to be important for fermentation in other organisms), suggesting that these archaea might be anaerobic fermentative heterotrophs.

Deciphering Symbiotic Interactions of "Candidatus Aenigmarchaeota" with Inferred Horizontal Gene Transfers and Co-occurrence Networks

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

The curious consistency of carbon biosignatures over billions of years of Earth-life coevolution

The oldest and most wide-ranging signal of biological activity (biosignature) on our planet is the carbon isotope composition of organic materials preserved in rocks. These biosignatures preserve the long-term evolution of the microorganism-hosted metabolic machinery responsible for producing deviations in the isotopic compositions of inorganic and organic carbon. Despite billions of years of ecosystem turnover, evolutionary innovation, organismic complexification, and geological events, the organic carbon that is a residuum of the global marine biosphere in the rock record tells an essentially static story. The ~25‰ mean deviation between inorganic and organic 13C/12C values has remained remarkably unchanged over >3.5 billion years. The bulk of this record is conventionally attributed to early-evolved, RuBisCO-mediated CO2 fixation that, in extant oxygenic phototrophs, produces comparable isotopic effects and dominates modern primary production. However, billions of years of environmental transition, for example, in the progressive oxygenation of the Earth's atmosphere, would be expected to have accompanied shifts in the predominant RuBisCO forms as well as enzyme-level adaptive responses in RuBisCO CO2-specificity. These factors would also be expected to result in preserved isotopic signatures deviating from those produced by extant RuBisCO in oxygenic phototrophs. Why does the bulk carbon isotope record not reflect these expected environmental transitions and evolutionary innovations? Here, we discuss this apparent discrepancy and highlight the need for greater quantitative understanding of carbon isotope fractionation behavior in extant metabolic pathways. We propose novel, laboratory-based approaches to reconstructing ancestral states of carbon metabolisms and associated enzymes that can constrain isotopic biosignature production in ancient biological systems. Together, these strategies are crucial for integrating the complementary toolsets of biological and geological sciences and for interpretation of the oldest record of life on Earth.

Comparative Genomics Provides Insights into the Genetic Diversity and Evolution of the DPANN Superphylum

DPANN is known as highly diverse, globally widespread, and mostly ectosymbiotic archaeal superphylum. However, this group of archaea was overlooked for a long time, and there were limited in-depth studies reported. In this investigation, 41 metagenome-assembled genomes (MAGs) belonging to the DPANN superphylum were recovered (18 MAGs had average nucleotide identity [ANI] values of <95% and a percentage of conserved proteins [POCP] of >50%, while 14 MAGs showed a POCP of <50%), which were analyzed comparatively with 515 other published DPANN genomes. Mismatches to known 16S rRNA gene primers were identified among 16S rRNA genes of DPANN archaea. Numbers of gene families lost (mostly related to energy and amino acid metabolism) were over three times greater than those gained in the evolution of DPANN archaea. Lateral gene transfer (LGT; ∼45.5% was cross-domain) had facilitated niche adaption of the DPANN archaea, ensuring a delicate equilibrium of streamlined genomes with efficient niche-adaptive strategies. For instance, LGT-derived cytochrome bd ubiquinol oxidase and arginine deiminase in the genomes of “Candidatus Micrarchaeota” could help them better adapt to aerobic acidic mine drainage habitats. In addition, most DPANN archaea acquired enzymes for biosynthesis of extracellular polymeric substances (EPS) and transketolase/transaldolase for the pentose phosphate pathway from Bacteria.

IMPORTANCE The domain Archaea is a key research model for gaining insights into the origin and evolution of life, as well as the relevant biogeochemical processes. The discovery of nanosized DPANN archaea has overthrown many aspects of microbiology. However, the DPANN superphylum still contains a vast genetic novelty and diversity that need to be explored. Comprehensively comparative genomic analysis on the DPANN superphylum was performed in this study, with an attempt to illuminate its metabolic potential, ecological distribution and evolutionary history. Many interphylum differences within the DPANN superphylum were found. For example, Altiarchaeota had the biggest genome among DPANN phyla, possessing many pathways missing in other phyla, such as formaldehyde assimilation and the Wood-Ljungdahl pathway. In addition, LGT acted as an important force to provide DPANN archaeal genetic flexibility that permitted the occupation of diverse niches. This study has advanced our understanding of the diversity and genome evolution of archaea.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Wide range of metabolic adaptations to the acquisition of the Calvin cycle revealed by comparison of microbial genomes

Knowledge of the genetic basis for autotrophic metabolism is valuable since it relates to both the emergence of life and to the metabolic engineering challenge of incorporating CO2 as a potential substrate for biorefining. The most common CO2 fixation pathway is the Calvin cycle, which utilizes Rubisco and phosphoribulokinase enzymes. We searched thousands of microbial genomes and found that 6.0% contained the Calvin cycle. We then contrasted the genomes of Calvin cycle-positive, non-cyanobacterial microbes and their closest relatives by enrichment analysis, ancestral character estimation, and random forest machine learning, to explore genetic adaptations associated with acquisition of the Calvin cycle. The Calvin cycle overlaps with the pentose phosphate pathway and glycolysis, and we could confirm positive associations with fructose-1,6-bisphosphatase, aldolase, and transketolase, constituting a conserved operon, as well as ribulose-phosphate 3-epimerase, ribose-5-phosphate isomerase, and phosphoglycerate kinase. Additionally, carbohydrate storage enzymes, carboxysome proteins (that raise CO2 concentration around Rubisco), and Rubisco activases CbbQ and CbbX accompanied the Calvin cycle. Photorespiration did not appear to be adapted specifically for the Calvin cycle in the non-cyanobacterial microbes under study. Our results suggest that chemoautotrophy in Calvin cycle-positive organisms was commonly enabled by hydrogenase, and less commonly ammonia monooxygenase (nitrification). The enrichment of specific DNA-binding domains indicated Calvin-cycle associated genetic regulation. Metabolic regulatory adaptations were illustrated by negative correlation to AraC and the enzyme arabinose-5-phosphate isomerase, which suggests a downregulation of the metabolite arabinose-5-phosphate, which may interfere with the Calvin cycle through enzyme inhibition and substrate competition. Certain domains of unknown function that were found to be important in the analysis may indicate yet unknown regulatory mechanisms in Calvin cycle-utilizing microbes. Our gene ranking provides targets for experiments seeking to improve CO2 fixation, or engineer novel CO2-fixing organisms.

Gemmatimonas groenlandica sp. nov. Is an Aerobic Anoxygenic Phototroph in the Phylum Gemmatimonadetes

The bacterial phylum Gemmatimonadetes contains members capable of performing bacteriochlorophyll-based phototrophy (chlorophototrophy). However, only one strain of chlorophototrophic Gemmatimonadetes bacteria (CGB) has been isolated to date, hampering our further understanding of their photoheterotrophic lifestyle and the evolution of phototrophy in CGB. By combining a culturomics strategy with a rapid screening technique for chlorophototrophs, we report the isolation of a new member of CGB, Gemmatimonas (G.) groenlandica sp. nov., from the surface water of a stream in the Zackenberg Valley in High Arctic Greenland. Distinct from the microaerophilic G. phototrophica strain AP64^T, G. groenlandica strain TET16^T is a strictly aerobic anoxygenic phototroph, lacking many oxygen-independent enzymes while possessing an expanded arsenal for coping with oxidative stresses. Its pigment composition and infra-red absorption properties are also different from G. phototrophica, indicating that it possesses a different photosystem apparatus. The complete genome sequence of G. groenlandica reveals unique and conserved features in the photosynthesis gene clusters of CGB. We further analyzed metagenome-assembled genomes of CGB obtained from soil and glacier metagenomes from Northeast Greenland, revealing a wide distribution pattern of CGB beyond the stream water investigated.

Marine Ultrasmall Prokaryotes Likely Affect the Cycling of Carbon, Methane, Nitrogen, and Sulfur

Recently, we uncovered the genetic components from six carbon fixation autotrophic pathways in cleaned ultrasmall size fractions from marine samples (<0.22 µm) gathered worldwide by the Tara Oceans Expedition. This first finding suggested that prokaryotic nanoorganisms, phylogenetically distantly related to the known CPR and DPANN groups, could collectively impact carbon cycling and carbon fixation across the world's ocean. To extend our mining of the functional and taxonomic microbial dark matter from the ultrasmall size fraction from the Tara Oceans Expedition, we investigated the distribution of 28 metabolic pathways associated with the cycling of carbon, methane, nitrogen, and sulfur. For all of these pathways, we report the existence not only of novel metabolic homologs in the ultrasmall size fraction of the oceanic microbiome, associated with nanoorganisms belonging to the CPR and DPANN lineages, but also of metabolic homologs exclusively found in marine host taxa belonging to other (still unassigned) microbial lineages. Therefore, we conclude that marine nanoorganisms contribute to a greater diversity of key biogeochemical cycles than currently appreciated. In particular, we suggest that oceanic nanoorganisms may be involved in a metabolic loop around Acetyl-CoA, have an underappreciated genetic potential to degrade methane, contribute to sustaining redox-reactions by producing Coenzyme F420, and affect sulfur cycling, notably as they harbor a complete suite of homologs of enzymes of the SOX system.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

Microbial dark matter filling the niche in hypersaline microbial mats

BACKGROUND

Shark Bay, Australia, harbours one of the most extensive and diverse systems of living microbial mats that are proposed to be analogs of some of the earliest ecosystems on Earth. These ecosystems have been shown to possess a substantial abundance of uncultivable microorganisms. These enigmatic microbes, jointly coined as 'microbial dark matter' (MDM), are hypothesised to play key roles in modern microbial mats.

RESULTS

We reconstructed 115 metagenome-assembled genomes (MAGs) affiliated to MDM, spanning 42 phyla. This study reports for the first time novel microorganisms (Zixibacterial order GN15) putatively taking part in dissimilatory sulfate reduction in surface hypersaline settings, as well as novel eukaryote signature proteins in the Asgard archaea. Despite possessing reduced-size genomes, the MDM MAGs are capable of fermenting and degrading organic carbon, suggesting a role in recycling organic carbon. Several forms of RuBisCo were identified, allowing putative CO₂ incorporation into nucleotide salvaging pathways, which may act as an alternative carbon and phosphorus source. High capacity of hydrogen production was found among Shark Bay MDM. Putative schizorhodopsins were also identified in Parcubacteria, Asgard archaea, DPANN archaea, and Bathyarchaeota, allowing these members to potentially capture light energy. Diversity-generating retroelements were prominent in DPANN archaea that likely facilitate the adaptation to a dynamic, host-dependent lifestyle.

CONCLUSIONS

This is the first study to reconstruct and describe in detail metagenome-assembled genomes (MAGs) affiliated with microbial dark matter in hypersaline microbial mats. Our data suggests that these microbial groups are major players in these systems. In light of our findings, we propose H₂, ribose and CO/CO₂ as the main energy currencies of the MDM community in these mat systems. Video Abstract.

Highly active rubiscos discovered by systematic interrogation of natural sequence diversity

CO2 is converted into biomass almost solely by the enzyme rubisco. The poor carboxylation properties of plant rubiscos have led to efforts that made it the most kinetically characterized enzyme, yet these studies focused on < 5% of its natural diversity. Here, we searched for fast-carboxylating variants by systematically mining genomic and metagenomic data. Approximately 33,000 unique rubisco sequences were identified and clustered into ≈ 1,000 similarity groups. We then synthesized, purified, and biochemically tested the carboxylation rates of 143 representatives, spanning all clusters of form-II and form-II/III rubiscos. Most variants (> 100) were active in vitro, with the fastest having a turnover number of 22 ± 1 s-1 -sixfold faster than the median plant rubisco and nearly twofold faster than the fastest measured rubisco to date. Unlike rubiscos from plants and cyanobacteria, the fastest variants discovered here are homodimers and exhibit a much simpler folding and activation kinetics. Our pipeline can be utilized to explore the kinetic space of other enzymes of interest, allowing us to get a better view of the biosynthetic potential of the biosphere.

Phosphoglycolate salvage in a chemolithoautotroph using the Calvin cycle

The Calvin cycle is the most important carbon fixation pathway in the biosphere. However, its carboxylating enzyme Rubisco also accepts oxygen, thus producing 2-phosphoglycolate. Phosphoglycolate salvage pathways were extensively studied in photoautotrophs but remain uncharacterized in chemolithoautotrophs using the Calvin cycle. Here, we study phosphoglycolate salvage in the chemolithoautotrophic model bacterium Cupriavidus necator H16. We demonstrate that this bacterium mainly reassimilates 2-phosphoglycolate via the glycerate pathway. Upon disruption of this pathway, a secondary route, which we term the malate cycle, supports photorespiration by completely oxidizing 2-phosphoglycolate to CO₂. While the malate cycle was not previously known to metabolize 2-phosphoglycolate in nature, a bioinformatic analysis suggests that it may support phosphoglycolate salvage in diverse chemoautotrophic bacteria.

Carbon fixation via the Calvin cycle is constrained by the side activity of Rubisco with dioxygen, generating 2-phosphoglycolate. The metabolic recycling of phosphoglycolate was extensively studied in photoautotrophic organisms, including plants, algae, and cyanobacteria, where it is referred to as photorespiration. While receiving little attention so far, aerobic chemolithoautotrophic bacteria that operate the Calvin cycle independent of light must also recycle phosphoglycolate. As the term photorespiration is inappropriate for describing phosphoglycolate recycling in these nonphotosynthetic autotrophs, we suggest the more general term “phosphoglycolate salvage.” Here, we study phosphoglycolate salvage in the model chemolithoautotroph Cupriavidus necator H16 (Ralstonia eutropha H16) by characterizing the proxy process of glycolate metabolism, performing comparative transcriptomics of autotrophic growth under low and high CO₂ concentrations, and testing autotrophic growth phenotypes of gene deletion strains at ambient CO₂. We find that the canonical plant-like C₂ cycle does not operate in this bacterium, and instead, the bacterial-like glycerate pathway is the main route for phosphoglycolate salvage. Upon disruption of the glycerate pathway, we find that an oxidative pathway, which we term the malate cycle, supports phosphoglycolate salvage. In this cycle, glyoxylate is condensed with acetyl coenzyme A (acetyl-CoA) to give malate, which undergoes two oxidative decarboxylation steps to regenerate acetyl-CoA. When both pathways are disrupted, autotrophic growth is abolished at ambient CO₂. We present bioinformatic data suggesting that the malate cycle may support phosphoglycolate salvage in diverse chemolithoautotrophic bacteria. This study thus demonstrates a so far unknown phosphoglycolate salvage pathway, highlighting important diversity in microbial carbon fixation metabolism.

Undinarchaeota illuminate DPANN phylogeny and the impact of gene transfer on archaeal evolution

The recently discovered DPANN archaea are a potentially deep-branching, monophyletic radiation of organisms with small cells and genomes. However, the monophyly and early emergence of the various DPANN clades and their role in life's evolution are debated. Here, we reconstructed and analysed genomes of an uncharacterized archaeal phylum (Candidatus Undinarchaeota), revealing that its members have small genomes and, while potentially being able to conserve energy through fermentation, likely depend on partner organisms for the acquisition of certain metabolites. Our phylogenomic analyses robustly place Undinarchaeota as an independent lineage between two highly supported 'DPANN' clans. Further, our analyses suggest that DPANN have exchanged core genes with their hosts, adding to the difficulty of placing DPANN in the tree of life. This pattern can be sufficiently dominant to allow identifying known symbiont-host clades based on routes of gene transfer. Together, our work provides insights into the origins and evolution of DPANN and their hosts.

The rise of diversity in metabolic platforms across the Candidate Phyla Radiation

BACKGROUND

A unifying feature of the bacterial Candidate Phyla Radiation (CPR) is a limited and highly variable repertoire of biosynthetic capabilities. However, the distribution of metabolic traits across the CPR and the evolutionary processes underlying them are incompletely resolved.

RESULTS

Here, we selected ~ 1000 genomes of CPR bacteria from diverse environments to construct a robust internal phylogeny that was consistent across two unlinked marker sets. Mapping of glycolysis, the pentose phosphate pathway, and pyruvate metabolism onto the tree showed that some components of these pathways are sparsely distributed and that similarity between metabolic platforms is only partially predicted by phylogenetic relationships. To evaluate the extent to which gene loss and lateral gene transfer have shaped trait distribution, we analyzed the patchiness of gene presence in a phylogenetic context, examined the phylogenetic depth of clades with shared traits, and compared the reference tree topology with those of specific metabolic proteins. While the central glycolytic pathway in CPR is widely conserved and has likely been shaped primarily by vertical transmission, there is evidence for both gene loss and transfer especially in steps that convert glucose into fructose 1,6-bisphosphate and glycerate 3P into pyruvate. Additionally, the distribution of Group 3 and Group 4-related NiFe hydrogenases is patchy and suggests multiple events of ancient gene transfer.

CONCLUSIONS

We infer that patterns of gene gain and loss in CPR, including acquisition of accessory traits in independent transfer events, could have been driven by shifts in host-derived resources and led to sparse but varied genetic inventories.

Heterotrophic Thaumarchaea with Small Genomes Are Widespread in the Dark Ocean

The Thaumarchaeota is a diverse archaeal phylum comprising numerous lineages that play key roles in global biogeochemical cycling, particularly in the ocean. To date, all genomically characterized marine thaumarchaea are reported to be chemolithoautotrophic ammonia oxidizers. In this study, we report a group of putatively heterotrophic marine thaumarchaea (HMT) with small genome sizes that is globally abundant in the mesopelagic, apparently lacking the ability to oxidize ammonia. We assembled five HMT genomes from metagenomic data and show that they form a deeply branching sister lineage to the ammonia-oxidizing archaea (AOA). We identify this group in metagenomes from mesopelagic waters in all major ocean basins, with abundances reaching up to 6% of that of AOA. Surprisingly, we predict the HMT have small genomes of ∼1 Mbp, and our ancestral state reconstruction indicates this lineage has undergone substantial genome reduction compared to other related archaea. The genomic repertoire of HMT indicates a versatile metabolism for aerobic chemoorganoheterotrophy that includes a divergent form III-a RuBisCO, a 2M respiratory complex I that has been hypothesized to increase energetic efficiency, and a three-subunit heme-copper oxidase complex IV that is absent from AOA. We also identify 21 pyrroloquinoline quinone (PQQ)-dependent dehydrogenases that are predicted to supply reducing equivalents to the electron transport chain and are among the most highly expressed HMT genes, suggesting these enzymes play an important role in the physiology of this group. Our results suggest that heterotrophic members of the Thaumarchaeota are widespread in the ocean and potentially play key roles in global chemical transformations.IMPORTANCE It has been known for many years that marine Thaumarchaeota are abundant constituents of dark ocean microbial communities, where their ability to couple ammonia oxidation and carbon fixation plays a critical role in nutrient dynamics. In this study, we describe an abundant group of putatively heterotrophic marine Thaumarchaeota (HMT) in the ocean with physiology distinct from those of their ammonia-oxidizing relatives. HMT lack the ability to oxidize ammonia and fix carbon via the 3-hydroxypropionate/4-hydroxybutyrate pathway but instead encode a form III-a RuBisCO and diverse PQQ-dependent dehydrogenases that are likely used to conserve energy in the dark ocean. Our work expands the scope of known diversity of Thaumarchaeota in the ocean and provides important insight into a widespread marine lineage.

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.

Genomic and transcriptomic evidence of light-sensing, porphyrin biosynthesis, Calvin-Benson-Bassham cycle, and urea production in Bathyarchaeota

BACKGROUND

Bathyarchaeota, a newly proposed archaeal phylum, is considered as an important driver of the global carbon cycle. However, due to the great diversity of them, there is limited genomic information that accurately encompasses the metabolic potential of the entire archaeal phylum.

RESULTS

In the current study, nine metagenome-assembled genomes of Bathyarchaeota from four subgroups were constructed from mangrove sediments, and metatranscriptomes were obtained for evaluating their in situ transcriptional activities. Comparative analyses with reference genomes and the transcripts of functional genes posit an expanded role for Bathyarchaeota in phototrophy, autotrophy, and nitrogen and sulfur cycles, respectively. Notably, the presence of genes for rhodopsins, cobalamin biosynthesis, and the oxygen-dependent metabolic pathways in some Bathyarchaeota subgroup 6 genomes suggest a light-sensing and microoxic lifestyle within this subgroup.

CONCLUSIONS

The results of this study expand our knowledge of metabolic abilities and diverse lifestyles of Bathyarchaeota, highlighting the crucial role of Bathyarchaeota in geochemical cycle. Video abstract.

Evolutionary trends in RuBisCO kinetics and their co-evolution with CO2 concentrating mechanisms

RuBisCO-catalyzed CO₂ fixation is the main source of organic carbon in the biosphere. This enzyme is present in all domains of life in different forms (III, II, and I) and its origin goes back to 3500 Mya, when the atmosphere was anoxygenic. However, the RuBisCO active site also catalyzes oxygenation of ribulose 1,5-bisphosphate, therefore, the development of oxygenic photosynthesis and the subsequent oxygen-rich atmosphere promoted the appearance of CO₂ concentrating mechanisms (CCMs) and/or the evolution of a more CO₂ -specific RuBisCO enzyme. The wide variability in RuBisCO kinetic traits of extant organisms reveals a history of adaptation to the prevailing CO₂ /O₂ concentrations and the thermal environment throughout evolution. Notable differences in the kinetic parameters are found among the different forms of RuBisCO, but the differences are also associated with the presence and type of CCMs within each form, indicative of co-evolution of RuBisCO and CCMs. Trade-offs between RuBisCO kinetic traits vary among the RuBisCO forms and also among phylogenetic groups within the same form. These results suggest that different biochemical and structural constraints have operated on each type of RuBisCO during evolution, probably reflecting different environmental selective pressures. In a similar way, variations in carbon isotopic fractionation of the enzyme point to significant differences in its relationship to the CO₂ specificity among different RuBisCO forms. A deeper knowledge of the natural variability of RuBisCO catalytic traits and the chemical mechanism of RuBisCO carboxylation and oxygenation reactions raises the possibility of finding unrevealed landscapes in RuBisCO evolution.

Computational Framework for High-Quality Production and Large-Scale Evolutionary Analysis of Metagenome Assembled Genomes

Microbial species play important roles in different environments and the production of high-quality genomes from metagenome data sets represents a major obstacle to understanding their ecological and evolutionary dynamics. Metagenome-Assembled Genomes Orchestra (MAGO) is a computational framework that integrates and simplifies metagenome assembly, binning, bin improvement, bin quality (completeness and contamination), bin annotation, and evolutionary placement of bins via detailed maximum-likelihood phylogeny based on multiple marker genes using different amino acid substitution models, next to average nucleotide identity analysis of genomes for delineation of species boundaries and operational taxonomic units. MAGO offers streamlined execution of the entire metagenomics pipeline, error checking, computational resource distribution and compatibility of data formats, governed by user-tailored pipeline processing. MAGO is an open-source-software package released in three different ways, as a singularity image and a Docker container for HPC purposes as well as for running MAGO on a commodity hardware, and a virtual machine for gaining a full access to MAGO underlying structure and source code. MAGO is open to suggestions for extensions and is amenable for use in both research and teaching of genomics and molecular evolution of genomes assembled from small single-cell projects or large-scale and complex environmental metagenomes.